Rheolef  7.2
an efficient C++ finite element environment
 
Loading...
Searching...
No Matches
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
see the space page for the full documentation
int main()
Definition field2bb.cc:58
This file is part of Rheolef.
STL namespace.
rheolef - reference manual