Rheolef
7.1
an efficient C++ finite element environment
mosolov_yield_surface.cc
The Mossolov problem – yield surface
#include "
rheolef.h
"
using namespace
rheolef
;
using namespace
std
;
int
main
(
int
argc,
char
**argv) {
environment
rheolef
(argc,argv);
Float
tol = (argc > 1) ? atof(argv[1]) : 1e-15;
Float
Bi;
field
sigma_h;
din
>>
catchmark
(
"Bi"
) >> Bi
>>
catchmark
(
"sigma"
) >> sigma_h;
space
Th = sigma_h.get_space();
space
Th1 (Th.get_geo(),
"P"
+
itos
(4*(Th.degree()+1)) +
"d"
);
dout
<<
interpolate
(Th1,
norm
(sigma_h)-Bi);
}
rheolef::catchmark
see the catchmark page for the full documentation
Definition:
catchmark.h:67
field
see the field page for the full documentation
space
see the space page for the full documentation
rheolef::norm
T norm(const vec< T, M > &x)
norm(x): see the expression page for the full documentation
Definition:
vec.h:387
rheolef.h
rheolef - reference manual
rheolef::interpolate
field_basic< T, M > interpolate(const space_basic< T, M > &V2h, const field_basic< T, M > &u1h)
see the interpolate page for the full documentation
Definition:
interpolate.cc:233
rheolef::din
idiststream din
see the diststream page for the full documentation
Definition:
diststream.h:427
rheolef::environment
see the environment page for the full documentation
Definition:
environment.h:115
rheolef
This file is part of Rheolef.
Definition:
compiler_eigen.h:37
Float
see the Float page for the full documentation
main
int main(int argc, char **argv)
Definition:
mosolov_yield_surface.cc:28
rheolef::dout
odiststream dout(cout)
see the diststream page for the full documentation
Definition:
diststream.h:430
rheolef::itos
std::string itos(std::string::size_type i)
itos: see the rheostream page for the full documentation
rheolef::std
Definition:
vec_expr_v2.h:402