pub trait SparseLinearSolver { // Required method fn solve( &mut self, residuals: &Mat<f64>, jacobians: &SparseColMat<usize, f64>, ) -> Option<Mat<f64>>; }