public static void SolveSparseLU(this SparseMatrix matrix, Vector <double> result, Vector <double> rhs, double tol = 1.0) { LU lu = matrix.SparseLU(tol: tol); lu.Solve(rhs, result); }