Improvement in Teuchos::SerialDenseSolver:::solve ()
Created by: bartgol
There is a check towards the beginning of the method to make sure that A and b are either both equilibrated or both not equilibrated. However, the check is done after equilibrating b and before factoring A (unless A is already factored/inverted)
Clearly, I could call the factor() method before calling solve(), and that would fix the problem. However, since solve() will automatically call factor() if needed, it would be nice if Teuchos did these checks only once it knows that the matrix is for sure factored (or inverted), since the problem of A not being factored may be solved in that step.