3 #ifndef DUNE_Pk3DLOCALINTERPOLATION_HH
4 #define DUNE_Pk3DLOCALINTERPOLATION_HH
17 static const int kdiv = (k == 0 ? 1 : k);
20 template<
typename F,
typename C>
23 typename LB::Traits::DomainType x;
24 typename LB::Traits::RangeType y;
25 typedef typename LB::Traits::DomainFieldType D;
28 for (
int i2 = 0; i2 <= k; i2++)
29 for (
int i1 = 0; i1 <= k-i2; i1++)
30 for (
int i0 = 0; i0 <= k-i1-i2; i0++)
32 x[0] = ((D)i0)/((D)kdiv);
33 x[1] = ((D)i1)/((D)kdiv);
34 x[2] = ((D)i2)/((D)kdiv);