Defects, in some sense, are not actually part of the Q-tensor system -- they just arise from particular configurations. Hence, we should not make that system query itself for defect locations.
Instead, we should follow the lead of the deal.II tutorial and create a generic EvaluateSolution class which can be passed into a NematicSystemMPI in order to evaluate something about the solution. In the future this should be true of all the output functions just to streamline things, but for now it's good enough to just have this keep track of the defect locations and times.
Defects, in some sense, are not actually part of the Q-tensor system -- they just arise from particular configurations. Hence, we should not make that system query itself for defect locations.
Instead, we should follow the lead of the deal.II tutorial and create a generic
EvaluateSolutionclass which can be passed into aNematicSystemMPIin order to evaluate something about the solution. In the future this should be true of all the output functions just to streamline things, but for now it's good enough to just have this keep track of the defect locations and times.