3 #ifndef DUNE_P1_LOCALINTERPOLATION_HH
4 #define DUNE_P1_LOCALINTERPOLATION_HH
10 template<
int dim,
class LB>
15 template<
typename F,
typename C>
18 typename LB::Traits::RangeType y;
19 typename LB::Traits::DomainType x;
24 for (
int i=0; i<dim; i++)
26 f.evaluate(x,y); out[0] = y;
29 for (
int i=0; i<dim; i++) {
30 for (
int j=0; j<dim; j++)
33 f.evaluate(x,y); out[i+1] = y;