MA48: Sparse unsymmetric system: driver for conventional direct method

To solve a sparse unsymmetric system of \(m\) linear equations in \(n\) unknowns using Gaussian elimination. There are facilities for block triangularization, choosing a good pivot order, factorizing a matrix with a given pivot order, and solving a system of equations using the factorized matrix. Error estimates are available.