2 from scipy.sparse
import csc_matrix
7 rng = np.random.default_rng()
9 A = rng.random((dim, dim))
10 A = (A + A.T) * 0.5 + np.diag(10.0 + rng.random(dim))
14 ldlt = eigenpy.SimplicialLDLT(A)
16 assert ldlt.info() == eigenpy.ComputationInfo.Success
20 D = csc_matrix(np.diag(ldlt.vectorD()))
25 X = rng.random((dim, 20))
31 ldlt.analyzePattern(A)
33 permutation = ldlt.permutationP()