5#ifndef DUNE_HIERARCHICAL_PRISM_P2_LOCALINTERPOLATION_HH
6#define DUNE_HIERARCHICAL_PRISM_P2_LOCALINTERPOLATION_HH
20 template<
typename F,
typename C>
21 void interpolate (
const F& f, std::vector<C>& out)
const
23 typename LB::Traits::DomainType x;
24 typename LB::Traits::RangeType y;
28 x[0] = 0.0; x[1] = 0.0; x[2] = 0.0; out[0] = f(x);
29 x[0] = 1.0; x[1] = 0.0; x[2] = 0.0; out[1] = f(x);
30 x[0] = 0.0; x[1] = 1.0; x[2] = 0.0; out[2] = f(x);
31 x[0] = 0.0; x[1] = 0.0; x[2] = 1.0; out[3] = f(x);
32 x[0] = 1.0; x[1] = 0.0; x[2] = 1.0; out[4] = f(x);
33 x[0] = 0.0; x[1] = 1.0; x[2] = 1.0; out[5] = f(x);
37 x[0] = 0.0; x[1] = 0.0; x[2] = 0.5; y = f(x);
38 out[6] = y - 0.5*(out[0] + out[3]);
40 x[0] = 1.0; x[1] = 0.0; x[2] = 0.5; y = f(x);
41 out[7] = y - 0.5*(out[1] + out[4]);
43 x[0] = 0.0; x[1] = 1.0; x[2] = 0.5; y = f(x);
44 out[8] = y - 0.5*(out[2] + out[5]);
46 x[0] = 0.5; x[1] = 0.0; x[2] = 0.0; y = f(x);
47 out[9] = y - 0.5*(out[0] + out[1]);
49 x[0] = 0.0; x[1] = 0.5; x[2] = 0.0; y = f(x);
50 out[10] = y - 0.5*(out[2] + out[0]);
52 x[0] = 0.5; x[1] = 0.5; x[2] = 0.0; y = f(x);
53 out[11] = y - 0.5*(out[2] + out[1]);
55 x[0] = 0.5; x[1] = 0.0; x[2] = 1.0; y = f(x);
56 out[12] = y - 0.5*(out[3] + out[4]);
58 x[0] = 0.0; x[1] = 0.5; x[2] = 1.0; y = f(x);
59 out[13] = y - 0.5*(out[3] + out[5]);
61 x[0] = 0.5; x[1] = 0.5; x[2] = 1.0; y = f(x);
62 out[14] = y - 0.5*(out[4] + out[5]);
66 x[0] = 0.5; x[1] = 0.0; x[2] = 0.5; y = f(x);
67 out[15] = y - 0.25*(out[4] + out[1] + out[0] + out[3] );
69 x[0] = 0.0; x[1] = 0.5; x[2] = 0.5; y = f(x);
70 out[16] = y - 0.25*(out[2] + out[0] + out[3] + out[5] );
72 x[0] = 0.5; x[1] = 0.5; x[2] = 0.5; y = f(x);
73 out[17] = y - 0.25*(out[2] + out[1] + out[4] + out[5] );
Definition: hierarchicalprismp2localinterpolation.hh:17
Dune namespace.
Definition: alignedallocator.hh:13