Function sparse_matrix_trace 

pub fn sparse_matrix_trace<F>(a: &CsrMatrix<F>) -> F
where
    F: Float + NumAssign + Sum + SparseElement + Debug + 'static,

Compute the trace of a sparse matrix (sum of diagonal elements).