Rheolef  7.2
an efficient C++ finite element environment
mosolov_yield_surface.cc

The Mossolov problem – yield surface.

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" + to_string(4*(Th.degree()+1)) + "d");
dout << lazy_interpolate (Th1, norm(sigma_h)-Bi);
}
see the Float page for the full documentation
see the field page for the full documentation
see the catchmark page for the full documentation
Definition: catchmark.h:67
see the environment page for the full documentation
Definition: environment.h:121
see the space page for the full documentation
int main(int argc, char **argv)
This file is part of Rheolef.
field_basic< T, M > lazy_interpolate(const space_basic< T, M > &X2h, const field_basic< T, M > &u1h)
see the interpolate page for the full documentation
Definition: field.h:871
T norm(const vec< T, M > &x)
norm(x): see the expression page for the full documentation
Definition: vec.h:387
STL namespace.
rheolef - reference manual