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

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
idiststream din(cin)
see the diststream page for the full documentation
Definition: diststream.h:464
odiststream dout(cout)
see the diststream page for the full documentation
Definition: diststream.h:467
see the space page for the full documentation
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
std::enable_if< details::has_field_rdof_interface< 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
rheolef - reference manual
int main(int argc, char **argv)
Definition: vorticity.cc:28