52 static constexpr void interpolate (
const F& f, std::vector<C>& out)
54 out.resize(numVertices+1);
60 for (
int i = 0; i < dim; ++i) {
71 std::vector<RangeType> sfValues;
72 LB::evaluateFunction(x, sfValues);
75 for (
int i = 0; i < numVertices; ++i)
76 out[numVertices] -= out[i]*sfValues[i];