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

The vorticity associated to a vector field.

The vorticity associated to a vector field

#include "rheolef.h"
using namespace rheolef;
using namespace std;
int main(int argc, char** argv) {
environment rheolef (argc, argv);
field uh;
din >> uh;
const space& Xh = uh.get_space();
string grad_approx = "P" + to_string(Xh.degree()-1) + "d";
string valued = (uh.get_geo().dimension() == 3) ? "vector" : "scalar";
space Lh (uh.get_geo(), grad_approx, valued);
field curl_uh = lazy_interpolate (Lh, curl(uh));
dout << catchmark("w") << curl_uh;
}
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:121
see the space page for the full documentation
This file is part of Rheolef.
std::enable_if< details::has_field_rdof_interface< Expr >::value, details::field_expr_v2_nonlinear_terminal_field< typenameExpr::scalar_type, typenameExpr::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 > 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
STL namespace.
rheolef - reference manual
int main(int argc, char **argv)
Definition: vorticity.cc:28