* Group common linear system solver code across `SparseSolverVisitor`, `SparseNonlinearSolverVisitor`, and `LinearSolverVisitor`.
* Normalize intermediate values in Gaussian elimination when the matrix has more than 5 rows.

Fixes #890