const_iterator end(const MeshEntity& e) const
{
assert(e.getDimension() <= dimension);
assert(e.getDimension() < dimension || e.getIndex() == 0);
const typename std::map<MeshEntity, qpoints_t>::const_iterator qIter = quadratureMap.find(e);
assert(qIter != quadratureMap.end());
return qIter->second.end();
}
void setQuadrature(const MeshEntity& e, const std::map<vertex<dimension>, double> weights)
{
assert(e.getDimension() <= dimension);
// We only have one cell per cell, so if we have quadrature for MeshEntity(dimension, k)
// for k>0, something is wrong
assert(e.getDimension() < dimension || e.getIndex() == 0);
quadratureMap[e] = qpoints_t(weights.begin(), weights.end());
}
