ShapeFunction::ShapeFunction(int nsf, const MasterElement &master) : sftable(master.ngauss_sets()), sfcache(master.ngauss_sets()), nfunctions(nsf) { for(int i=0; i<master.ngauss_sets(); i++) { sftable[i] = new ShapeFunctionTable(master.ngauss(i), nsf); sfcache[i] = new ShapeFunctionCache(master.ngauss(i), nsf); } }