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);
}
see the Float page for the full documentation
see the field page for the full documentation
idiststream din
see the diststream page for the full documentation
Definition: diststream.h:427
odiststream dout(cout)
see the diststream page for the full documentation
Definition: diststream.h:430
see the space page for the full documentation
int main(int argc, char **argv)
This file is part of Rheolef.
T norm(const vec< T, M > &x)
norm(x): see the expression page for the full documentation
Definition: vec.h:387
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
std::string itos(std::string::size_type i)
itos: see the rheostream page for the full documentation
rheolef - reference manual