cmfd shift cant be infinity due to FPE, making it very large and fixed matrix subtraction in solver