The CPUSolver class has a few routines which overlap with the mesh class routines, such as
- printFissionRates
- printNegativeSources
These should be removed from the solver class and moved to the mesh class.
Then the Solver should call the mesh class routines after every sweep as appropriate to output the desired reaction rates / source / flux / negative sources/fluxes distributions.