Subroutines to solve a sparse linear system A*X=B - View it on GitHub