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