3#ifndef DUNE_HIERARCHICAL_PRISM_P2_LOCALINTERPOLATION_HH
4#define DUNE_HIERARCHICAL_PRISM_P2_LOCALINTERPOLATION_HH
7#include <dune/localfunctions/common/localinterpolation.hh>
19 template<
typename F,
typename C>
20 void interpolate (
const F& ff, std::vector<C>& out)
const
22 typename LB::Traits::DomainType x;
23 typename LB::Traits::RangeType y;
26 auto&& f = Impl::makeFunctionWithCallOperator<decltype(x)>(ff);
29 x[0] = 0.0; x[1] = 0.0; x[2] = 0.0; out[0] = f(x);
30 x[0] = 1.0; x[1] = 0.0; x[2] = 0.0; out[1] = f(x);
31 x[0] = 0.0; x[1] = 1.0; x[2] = 0.0; out[2] = f(x);
32 x[0] = 0.0; x[1] = 0.0; x[2] = 1.0; out[3] = f(x);
33 x[0] = 1.0; x[1] = 0.0; x[2] = 1.0; out[4] = f(x);
34 x[0] = 0.0; x[1] = 1.0; x[2] = 1.0; out[5] = f(x);
38 x[0] = 0.0; x[1] = 0.0; x[2] = 0.5; y = f(x);
39 out[6] = y - 0.5*(out[0] + out[3]);
41 x[0] = 1.0; x[1] = 0.0; x[2] = 0.5; y = f(x);
42 out[7] = y - 0.5*(out[1] + out[4]);
44 x[0] = 0.0; x[1] = 1.0; x[2] = 0.5; y = f(x);
45 out[8] = y - 0.5*(out[2] + out[5]);
47 x[0] = 0.5; x[1] = 0.0; x[2] = 0.0; y = f(x);
48 out[9] = y - 0.5*(out[0] + out[1]);
50 x[0] = 0.0; x[1] = 0.5; x[2] = 0.0; y = f(x);
51 out[10] = y - 0.5*(out[2] + out[0]);
53 x[0] = 0.5; x[1] = 0.5; x[2] = 0.0; y = f(x);
54 out[11] = y - 0.5*(out[2] + out[1]);
56 x[0] = 0.5; x[1] = 0.0; x[2] = 1.0; y = f(x);
57 out[12] = y - 0.5*(out[3] + out[4]);
59 x[0] = 0.0; x[1] = 0.5; x[2] = 1.0; y = f(x);
60 out[13] = y - 0.5*(out[3] + out[5]);
62 x[0] = 0.5; x[1] = 0.5; x[2] = 1.0; y = f(x);
63 out[14] = y - 0.5*(out[4] + out[5]);
67 x[0] = 0.5; x[1] = 0.0; x[2] = 0.5; y = f(x);
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; y = f(x);
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; y = f(x);
74 out[17] = y - 0.25*(out[2] + out[1] + out[4] + out[5] );
Definition: hierarchicalprismp2localinterpolation.hh:16
Dune namespace.
Definition: alignedallocator.hh:11