pub fn sparse_lower_triangular_solve(l: &CsrMatrix, b: &[f64]) -> Vec<f64>Expand description
Forward-substitution solve L x = b where L is lower-triangular (CSR).
Assumes L has non-zero diagonal entries. The diagonal entry for row i
is taken as L.get(i, i).