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);
Definition: pk3dlocalinterpolation.hh:11
Definition: brezzidouglasmarini1cube2dlocalbasis.hh:14
void interpolate(const F &f, std::vector< C > &out) const
Definition: pk3dlocalinterpolation.hh:21