Page Menu
Home
c4science
Search
Configure Global Search
Log In
Files
F68529669
bem_kato.hh
No One
Temporary
Actions
Download File
Edit File
Delete File
View Transforms
Subscribe
Mute Notifications
Award Token
Subscribers
None
File Metadata
Details
File Info
Storage
Attached
Created
Thu, Jun 27, 20:07
Size
2 KB
Mime Type
text/x-c++
Expires
Sat, Jun 29, 20:07 (2 d)
Engine
blob
Format
Raw Data
Handle
18596721
Attached To
rTAMAAS tamaas
bem_kato.hh
View Options
#ifndef BEM_KATO_H
#define BEM_KATO_H
/* -------------------------------------------------------------------------- */
#include "bem_fft_base.hh"
/* -------------------------------------------------------------------------- */
class
BemKato
:
public
BemFFTBase
{
/* ------------------------------------------------------------------------ */
/* Constructors/Destructors */
/* ------------------------------------------------------------------------ */
public
:
BemKato
(
Surface
<
Real
>
&
p
)
:
BemFFTBase
(
p
),
surface_tractions_backup
(
p
.
size
(),
p
.
getL
()){
e_star
=
1.
;
max_pressure
=
1e300
;
};
virtual
~
BemKato
(){};
/* ------------------------------------------------------------------------ */
/* Methods */
/* ------------------------------------------------------------------------ */
public
:
//! compute the equilibrium situation
Real
computeEquilibrium
(
Real
epsilon
,
Real
pressure
,
int
nthreads
=
1
);
public
:
//! perform a line scan in a given direction
Real
linescan
(
Real
scale
,
Real
pressure
);
//! perform a line search
Real
linesearch
(
Real
&
hmax
,
Real
fold
,
Real
pressure
,
int
search_flag
);
//! make a step in the gradient direction(copy in pressure2)
void
updatePressure
(
Real
scale
);
//! backup tractions before doing the linescan
void
backupTractions
();
//! find the constant pressure shift making balance with load
void
shiftPressure
(
Real
required_pressure
);
//! truncate the pressure to positive values
Real
positivePressure
(
Real
step
);
//! compute the reaction (sum positive pressures)
void
truncatePressure
();
/* ------------------------------------------------------------------------ */
/* Accessors */
/* ------------------------------------------------------------------------ */
public
:
//! get the number of iterations
UInt
getNbIterations
()
const
{
return
this
->
nb_iterations
;};
//! get the convergence steps of iterations
const
std
::
vector
<
Real
>
&
getConvergenceIterations
()
const
{
return
this
->
convergence_iterations
;};
//! set the maximal pressure admissible
void
setMaxPressure
(
Real
max_pressure
){
this
->
max_pressure
=
max_pressure
;};
/* ------------------------------------------------------------------------ */
/* Class Members */
/* ------------------------------------------------------------------------ */
protected
:
//! surface of tractions used during steepest descent
Surface
<
Real
>
surface_tractions_backup
;
//! maximal pressure admissible (for perfectly plastic purpose)
Real
max_pressure
;
};
#endif
/* BEM_KATO_H */
Event Timeline
Log In to Comment