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

Extraction as a surface mesh of the zero level set – spherical geometry

#include "rheolef.h"
using namespace rheolef;
using namespace std;
#include "torus.icc"
int main (int argc, char**argv) {
environment rheolef (argc,argv);
geo lambda (argv[1]);
level_set_option opts;
opts.split_to_triangle
= (argc > 2 && argv[2] == std::string("-tq")) ? false : true;
space Xh (lambda, "P1");
field phi_h = lazy_interpolate(Xh, phi);
geo gamma = level_set (phi_h, opts);
dout << gamma;
}
see the field page for the full documentation
see the geo page for the full documentation
odiststream dout(cout)
see the diststream page for the full documentation
Definition: diststream.h:467
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
geo_basic< T, M > level_set(const field_basic< T, M > &fh, const level_set_option &opt)
Definition: level_set.cc:925
Float gamma[][pmax+1]
rheolef - reference manual
Definition: phi.h:25
The torus benchmark – level set, right-hand-side and exact solution.