added new cmfd solver that doesnt require Petsc
added SOR
solver now works with coremap feature, loop is around rows now instead of i,j,k
default spectral radius is 0 (G-S no SOR)
added 1g solver and procedure pointer
can now set G-S tolerance, updated tests
added Wielandt shifting
added matrix copy and option to read in shift
added relative tolerancing and now can be set on input
removed copy of loss matrix for shift because we can overwrite loss matrix
updated cmfd test cases for tighter tolerances
cmfd shift cant be infinity due to FPE, making it very large and fixed matrix subtraction in solver
added res scat variables to template files
added type for 0K nuclides
added 0K xs global variables
modified template files to allow for multiple scatterers
added logical indicating resonant nuclides
first crack at weight correction method (WCM) target sampling
allow for diff E ranges and methods for each scatterer
added resonance scattering header