2 #ifndef DUNE_PRISM_P2_LOCALINTERPOLATION_HH
3 #define DUNE_PRISM_P2_LOCALINTERPOLATION_HH
15 template<
typename F,
typename C>
18 typename LB::Traits::RangeType y;
21 typename LB::Traits::DomainType x;
23 x[0] = 0.0; x[1] = 0.0; x[2] = 0.0;
24 f.evaluate(x,y); out[0] = y;
26 x[0] = 1.0; x[1] = 0.0; x[2] = 0.0;
27 f.evaluate(x,y); out[1] = y;
29 x[0] = 0.0; x[1] = 1.0; x[2] = 0.0;
30 f.evaluate(x,y); out[2] = y;
32 x[0] = 0.0; x[1] = 0.0; x[2] = 1.0;
33 f.evaluate(x,y); out[3] = y;
35 x[0] = 1.0; x[1] = 0.0; x[2] = 1.0;
36 f.evaluate(x,y); out[4] = y;
38 x[0] = 0.0; x[1] = 1.0; x[2] = 1.0;
39 f.evaluate(x,y); out[5] = y;
41 x[0] = 0.0; x[1] = 0.0; x[2] = 0.5;
42 f.evaluate(x,y); out[6] = y;
44 x[0] = 1.0; x[1] = 0.0; x[2] = 0.5;
45 f.evaluate(x,y); out[7] = y;
47 x[0] = 0; x[1] = 1.0; x[2] = 0.5;
48 f.evaluate(x,y); out[8] = y;
50 x[0] = 0.5; x[1] = 0.0; x[2] = 0.0;
51 f.evaluate(x,y); out[9] = y;
53 x[0] = 0.0; x[1] = 0.5; x[2] = 0.0;
54 f.evaluate(x,y); out[10] = y;
56 x[0] = 0.5; x[1] = 0.5; x[2] = 0.0;
57 f.evaluate(x,y); out[11] = y;
59 x[0] = 0.5; x[1] = 0.0; x[2] = 1.0;
60 f.evaluate(x,y); out[12] = y;
62 x[0] = 0.0; x[1] = 0.5; x[2] = 1.0;
63 f.evaluate(x,y); out[13] = y;
65 x[0] = 0.5; x[1] = 0.5; x[2] = 1.0;
66 f.evaluate(x,y); out[14] = y;
68 x[0] = 0.5; x[1] = 0.0; x[2] = 0.5;
69 f.evaluate(x,y); out[15] = y;
71 x[0] = 0.0; x[1] = 0.5; x[2] = 0.5;
72 f.evaluate(x,y); out[16] = y;
74 x[0] = 0.5; x[1] = 0.5; x[2] = 0.5;
75 f.evaluate(x,y); out[17] = y;