void ShapeFunction::precompute(const MasterElement &master) {
// Use the virtual functions that evaluate the shapefunction at
// arbitrary points to store its values at the Gauss points.
// loop over integration orders (sets of gauss points)
for(int ord=0; ord<master.ngauss_sets(); ord++) {
const GaussPtTable &gptable = master.gptable(ord);
sftable[ord] = new ShapeFunctionTable(gptable.size(), nfunctions);
std::vector<DoubleVec> &f_table = sftable[ord]->f_table;
std::vector<std::vector<DoubleVec> > &df_table =
sftable[ord]->df_table;
// loop over gausspoints
for(std::vector<GaussPtData>::size_type g=0; g<gptable.size(); g++) {
MasterCoord mpos = gptable[g].position;
for(ShapeFunctionIndex n=0; n<nfunctions; ++n) { // loop over sf's
f_table[g][n] = value(n, mpos);
DoubleVec &dftemp = df_table[g][n];
for(SpaceIndex j=0; j<DIM; ++j) // loop over spatial dimensions
dftemp[j] = masterderiv(n, j, mpos);
}
}
}
}
