Design for the solver - maroonbear/MultilevelCMFD GitHub Wiki

We are going to design the solver here.

1. High order solution provides Dhat, xsecs, keff and initial flux guess as input? (Should Dhat be an input?)
2. Gauss Seidel iteration for sweeping.
3. Coarse Mesh size goes down to 1 mesh by 1 mesh.
4. "Eigenvalue" calculation is at the bottom level. 2G 1 mesh by 1 mesh. Solved by hand.