Function sprs::linalg::diag_solve

pub fn diag_solve<'a, N, V1, V2>(diag: V1, x: V2) where
    N: 'a + Clone + Num + DivAssign,
    N: DivAssign<&'r N>,
    V1: DenseVector<Scalar = N>,
    V2: DenseVectorMut + DenseVector<Scalar = N>, 

Diagonal solve