Function sprs_ldl::ldl_numeric

pub fn ldl_numeric<N, PStorage>(mat: CsMatView<N>, l_colptr: &[usize], parents: ParentsView, perm: &Permutation<PStorage>, l_nz: &mut [usize], l_indices: &mut [usize], l_data: &mut [N], diag: &mut [N], y_workspace: &mut [N], pattern_workspace: &mut DStack<usize>, flag_workspace: &mut [usize]) where N: Clone + Copy + PartialEq + Num + PartialOrd, PStorage: Deref<Target=[usize]>

Perform numeric LDLT decomposition

pattern_workspace is a DStack of capacity n