1 2 3 4 5 6 7 8 9 10 11 12
for(int f=0; f<(Element::nf)*2; f++) { if (f%2==1) p = p-1; { this->pInterpolation[i]=p; this->cInterpolation[i]=0.; this->dofInterpolation[i]=e+f; this->coefInterpolation[i]=1.; i++; p++; } }