AgglomerationHandler should not be responsible for the initialization and setup of FEValues-like objects, in the same way as DoFHandler does not evaluate any basis function. A dedicated class, with name something along the lines of FEAggloValues or FEPolyValues, should be the one that calls these methods:
Its constructor should be identical to the one of FEValuesBase with the twist that instead of using a DoFHandler, and AgglomerationHandler is working under the hood.
AgglomerationHandlershould not be responsible for the initialization and setup of FEValues-like objects, in the same way asDoFHandlerdoes not evaluate any basis function. A dedicated class, with name something along the lines ofFEAggloValuesorFEPolyValues,should be the one that calls these methods:AgglomerationHandler<dim,spacedim>::reinit(polytope),AgglomerationHandler<dim,spacedim>::reinit(polytope,face)Its constructor should be identical to the one of FEValuesBase with the twist that instead of using a DoFHandler, and
AgglomerationHandleris working under the hood.