2
using namespace rheolef;
4
int main(int argc, char** argv) {
5
environment rheolef (argc, argv);
8
string grad_approx = "P" + itos(uh.get_space().degree()-1) + "d";
9
string valued = (uh.size() == 3) ? "vector" : "scalar";
10
space Lh (uh.get_geo(), grad_approx, valued);
11
form curl (uh.get_space(), Lh, "curl");
12
form inv_m (Lh, Lh, "inv_mass");
13
dout << catchmark("w") << inv_m*(curl*uh);