At this point, we refine around defects and around the center of the grid. For the former, this only includes a cylindrical refinement around the line along the z-axis through a defect point. For the former, this includes a rectangular prism refinement around the center at distances of L/(2^n).
Presumably what should happen is that each refinement scheme should correspond to further refinement around some geometric structure (point or a line). Additionally, there should be a class enum which contains the particular norm which is being computed. Then there can be a vector of pointers to the parent class GeometricGridRefinement or something. This should just take in a mutable reference to a triangulation and it should refine according to distance from the geometric object.
At this point, we refine around defects and around the center of the grid. For the former, this only includes a cylindrical refinement around the line along the z-axis through a defect point. For the former, this includes a rectangular prism refinement around the center at distances of L/(2^n).
Presumably what should happen is that each refinement scheme should correspond to further refinement around some geometric structure (point or a line). Additionally, there should be a
class enum
which contains the particular norm which is being computed. Then there can be a vector of pointers to the parent classGeometricGridRefinement
or something. This should just take in a mutable reference to a triangulation and it should refine according to distance from the geometric object.