Rheolef  7.1
an efficient C++ finite element environment
mosolov_yield_surface.cc
Go to the documentation of this file.
1 #include "rheolef.h"
26 using namespace rheolef;
27 using namespace std;
28 int main(int argc, char**argv) {
29  environment rheolef (argc,argv);
30  Float tol = (argc > 1) ? atof(argv[1]) : 1e-15;
31  Float Bi;
32  field sigma_h;
33  din >> catchmark("Bi") >> Bi
34  >> catchmark("sigma") >> sigma_h;
35  space Th = sigma_h.get_space();
36  space Th1 (Th.get_geo(), "P" + itos(4*(Th.degree()+1)) + "d");
37  dout << interpolate (Th1, norm(sigma_h)-Bi);
38 }
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:115
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