Function diag_matrix

pub fn diag_matrix<F>(
    diag: &[F],
    n: Option<usize>,
) -> SparseResult<CsrMatrix<F>>
where
    F: Float + NumAssign + Sum + 'static + Debug,

Create a diagonal matrix