For a sparse square nonsingular matrix, $A$, of order N, and an
optional vector or matrix B, this routine will factor A, solve A X =
B, \ A^T X = B, compute the reciprocal of the condition number, and
obtain the determinant and the inverse. This code has similar
functionality to the dense codes described in Chapter 4-1. The
program in Chapter 19-6 can be used together with the code here to
solve problems with the matrix given in Matrix Market format.
This software is part of the MATH77/mathc90 collection here.