Rheolef  7.1
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" + itos(Xh.degree()-1) + "d";
string valued = (uh.get_geo().dimension() == 3) ? "vector" : "scalar";
space Lh (uh.get_geo(), grad_approx, valued);
field curl_uh = interpolate (Lh, curl(uh));
dout << catchmark("w") << curl_uh;
}
see the field page for the full documentation
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