Function sparse_diag 

pub fn sparse_diag<T: Clone + Zero>(
    diagonal: &[T],
) -> UtilsResult<(Vec<T>, Vec<usize>, Vec<usize>)>

Create sparse diagonal matrix