Sparse LU factorization with partial pivoting
Rust binding to lufact FORTRAN routines
Solve sparse systems of linear equations