3 #ifndef DUNE_HIERARCHICAL_PRISM_P2_LOCALINTERPOLATION_HH
4 #define DUNE_HIERARCHICAL_PRISM_P2_LOCALINTERPOLATION_HH
18 template<
typename F,
typename C>
21 typename LB::Traits::DomainType x;
22 typename LB::Traits::RangeType y;
29 x[0] = 0.0; x[1] = 0.0; x[2] = 0.0; f.evaluate(x, y); out[0] = y;
30 x[0] = 1.0; x[1] = 0.0; x[2] = 0.0; f.evaluate(x, y); out[1] = y;
31 x[0] = 0.0; x[1] = 1.0; x[2] = 0.0; f.evaluate(x, y); out[2] = y;
32 x[0] = 0.0; x[1] = 0.0; x[2] = 1.0; f.evaluate(x, y); out[3] = y;
33 x[0] = 1.0; x[1] = 0.0; x[2] = 1.0; f.evaluate(x, y); out[4] = y;
34 x[0] = 0.0; x[1] = 1.0; x[2] = 1.0; f.evaluate(x, y); out[5] = y;
38 x[0] = 0.0; x[1] = 0.0; x[2] = 0.5; f.evaluate(x, y);
39 out[6] = y - 0.5*(out[0] + out[3]);
41 x[0] = 1.0; x[1] = 0.0; x[2] = 0.5; f.evaluate(x, y);
42 out[7] = y - 0.5*(out[1] + out[4]);
44 x[0] = 0.0; x[1] = 1.0; x[2] = 0.5; f.evaluate(x, y);
45 out[8] = y - 0.5*(out[2] + out[5]);
47 x[0] = 0.5; x[1] = 0.0; x[2] = 0.0; f.evaluate(x, y);
48 out[9] = y - 0.5*(out[0] + out[1]);
50 x[0] = 0.0; x[1] = 0.5; x[2] = 0.0; f.evaluate(x, y);
51 out[10] = y - 0.5*(out[2] + out[0]);
53 x[0] = 0.5; x[1] = 0.5; x[2] = 0.0; f.evaluate(x, y);
54 out[11] = y - 0.5*(out[2] + out[1]);
56 x[0] = 0.5; x[1] = 0.0; x[2] = 1.0; f.evaluate(x, y);
57 out[12] = y - 0.5*(out[3] + out[4]);
59 x[0] = 0.0; x[1] = 0.5; x[2] = 1.0; f.evaluate(x, y);
60 out[13] = y - 0.5*(out[3] + out[5]);
62 x[0] = 0.5; x[1] = 0.5; x[2] = 1.0; f.evaluate(x, y);
63 out[14] = y - 0.5*(out[4] + out[5]);
67 x[0] = 0.5; x[1] = 0.0; x[2] = 0.5; f.evaluate(x, y);
68 out[15] = y - 0.25*(out[4] + out[1] + out[0] + out[3] );
70 x[0] = 0.0; x[1] = 0.5; x[2] = 0.5; f.evaluate(x, y);
71 out[16] = y - 0.25*(out[2] + out[0] + out[3] + out[5] );
73 x[0] = 0.5; x[1] = 0.5; x[2] = 0.5; f.evaluate(x, y);
74 out[17] = y - 0.25*(out[2] + out[1] + out[4] + out[5] );