3
using namespace rheolef;
4
int main (int argc, char**argv) {
5
environment rheolef (argc, argv);
7
din >> catchmark("u") >> uh
8
>> catchmark("phi") >> phi_h;
9
const space& Xh = phi_h.get_space();
11
space Wh (gh.level_set(), "P1");
12
gh.level_set().save();
13
dout << interpolate (Wh, uh);