3 #ifndef DUNE_Pk2DLOCALINTERPOLATION_HH
4 #define DUNE_Pk2DLOCALINTERPOLATION_HH
20 static const int kdiv = (k == 0? 1 : k);
24 template<
typename F,
typename C>
27 typename LB::Traits::DomainType x;
28 typename LB::Traits::RangeType y;
29 typedef typename LB::Traits::DomainFieldType D;
32 for (
int j=0; j<=k; j++)
33 for (
int i=0; i<=k-j; i++)
35 x[0] = ((D)i)/((D)kdiv); x[1] = ((D)j)/((D)kdiv);