spchol - Sparse cholesky factorization (original) (raw)

X=[ 3., 0., 0., 2., 0., 0., 2., 0., 2., 0., 0. ; 0., 5., 4., 0., 0., 0., 0., 0., 0., 0., 0. ; 0., 4., 5., 0., 0., 0., 0., 0., 0., 0., 0. ; 2., 0., 0., 3., 0., 0., 2., 0., 2., 0., 0. ; 0., 0., 0., 0. , 5., 0., 0., 0., 0., 0., 4. ; 0., 0., 0., 0., 0., 4., 0., 3., 0., 3., 0. ; 2., 0., 0., 2., 0., 0., 3., 0., 2., 0., 0. ; 0., 0., 0., 0., 0., 3., 0., 4., 0., 3., 0. ; 2., 0., 0., 2., 0., 0., 2., 0., 3., 0., 0. ; 0., 0., 0., 0., 0., 3., 0., 3., 0., 4., 0. ; 0., 0., 0., 0., 4., 0., 0., 0., 0., 0., 5.]; X=sparse(X);[R,P] = spchol(X); max(PRR'*P'-X)