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:
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, the AgglomerationHandler class 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, the
AgglomerationHandlerclass is working under the hood.