pub fn lower_triangular(a: &CsrMatrix) -> CsrMatrix
Extract the lower triangular part of a sparse matrix (including diagonal).