Rheolef  7.1
an efficient C++ finite element environment
vorticity.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  field uh;
31  din >> uh;
32  const space& Xh = uh.get_space();
33  string grad_approx = "P" + itos(Xh.degree()-1) + "d";
34  string valued = (uh.get_geo().dimension() == 3) ? "vector" : "scalar";
35  space Lh (uh.get_geo(), grad_approx, valued);
36  field curl_uh = interpolate (Lh, curl(uh));
37  dout << catchmark("w") << curl_uh;
38 }
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
This file is part of Rheolef.
std::enable_if< details::is_field_convertible< Expr >::value,details::field_expr_v2_nonlinear_terminal_field< typename Expr::scalar_type,typename Expr::memory_type,details::differentiate_option::curl >>::type curl(const Expr &expr)
curl(uh): see the expression page for the full documentation
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
int main(int argc, char **argv)
Definition: vorticity.cc:28