2 eigenpy.switchToNumpyArray()
5 import numpy.linalg
as la
8 A = np.random.rand(dim,dim)
10 A = (A + A.T)*0.5 + np.diag(10. + np.random.rand(dim))
12 ldlt = eigenpy.LDLT(A)
16 P = ldlt.transpositionsP()
18 assert eigenpy.is_approx(np.transpose(P).dot(L.dot(np.diag(D).dot(np.transpose(L).dot(P)))),A)
EIGEN_DONT_INLINE bool is_approx(const Eigen::MatrixBase< MatrixType1 > &mat1, const Eigen::MatrixBase< MatrixType2 > &mat2)