Function sprs_ldl::ldl_numeric

pub fn ldl_numeric<N, I, PStorage>(
    mat: CsMatViewI<'_, N, I>, 
    l_colptr: &[I], 
    parents: ParentsView<'_>, 
    perm: &Permutation<I, PStorage>, 
    l_nz: &mut [I], 
    l_indices: &mut [I], 
    l_data: &mut [N], 
    diag: &mut [N], 
    y_workspace: &mut [N], 
    pattern_workspace: &mut DStack<I>, 
    flag_workspace: &mut [I]
) -> Result<(), LinalgError> where
    N: Clone + Copy + PartialEq + Num + PartialOrd,
    I: SpIndex,
    PStorage: Deref<Target = [I]>, 

Perform numeric LDLT decomposition

pattern_workspace is a DStack of capacity n