11 #include <Eigen/LU>
26      LU.h  in lu_non_invertible()
71   FullPivLU<MatrixType> lu;  in lu_non_invertible()  local
76   lu.setThreshold(RealScalar(0.01));  in lu_non_invertible()
77   lu.compute(m1);  in lu_non_invertible()
80   u = lu.matrixLU().template triangularView<Upper>();  in lu_non_invertible()
83     = lu.matrixLU().block(0,0,rows,(std::min)(rows,cols));  in lu_non_invertible()
85   VERIFY_IS_APPROX(lu.permutationP() * m1 * lu.permutationQ(), l*u);  in lu_non_invertible()
87   KernelMatrixType m1kernel = lu.kernel();  in lu_non_invertible()
88   ImageMatrixType m1image = lu.image(m1);  in lu_non_invertible()
90   VERIFY_IS_APPROX(m1, lu.reconstructedMatrix());  in lu_non_invertible()
91   VERIFY(rank == lu.rank());  in lu_non_invertible()
92   VERIFY(cols - lu.rank() == lu.dimensionOfKernel());  in lu_non_invertible()
93   VERIFY(!lu.isInjective());  in lu_non_invertible()
94   VERIFY(!lu.isInvertible());  in lu_non_invertible()
95   VERIFY(!lu.isSurjective());  in lu_non_invertible()
100   check_solverbase<CMatrixType, MatrixType>(m1, lu, rows, cols, cols2);  in lu_non_invertible()
106   m2.block(0,0,m2.rows(),m2.cols()) = lu.solve(m3);  in lu_non_invertible()
121   FullPivLU<MatrixType> lu;  in lu_invertible()  local
122   lu.setThreshold(RealScalar(0.01));  in lu_invertible()
125     lu.compute(m1);  in lu_invertible()
126   } while(!lu.isInvertible());  in lu_invertible()
128   VERIFY_IS_APPROX(m1, lu.reconstructedMatrix());  in lu_invertible()
129   VERIFY(0 == lu.dimensionOfKernel());  in lu_invertible()
130   VERIFY(lu.kernel().cols() == 1); // the kernel() should consist of a single (zero) column vector  in lu_invertible()
131   VERIFY(size == lu.rank());  in lu_invertible()
132   VERIFY(lu.isInjective());  in lu_invertible()
133   VERIFY(lu.isSurjective());  in lu_invertible()
134   VERIFY(lu.isInvertible());  in lu_invertible()
135   VERIFY(lu.image(m1).fullPivLu().isInvertible());  in lu_invertible()
137   check_solverbase<MatrixType, MatrixType>(m1, lu, size, size, size);  in lu_invertible()
139   MatrixType m1_inverse = lu.inverse();  in lu_invertible()
141   m2 = lu.solve(m3);  in lu_invertible()
145   const RealScalar rcond_est = lu.rcond();  in lu_invertible()
152   VERIFY_IS_APPROX(lu.solve(m3*m4), lu.solve(m3)*m4);  in lu_invertible()
187   FullPivLU<MatrixType> lu;  in lu_verify_assert()  local
188   VERIFY_RAISES_ASSERT(lu.matrixLU())  in lu_verify_assert()
189   VERIFY_RAISES_ASSERT(lu.permutationP())  in lu_verify_assert()
190   VERIFY_RAISES_ASSERT(lu.permutationQ())  in lu_verify_assert()
191   VERIFY_RAISES_ASSERT(lu.kernel())  in lu_verify_assert()
192   VERIFY_RAISES_ASSERT(lu.image(tmp))  in lu_verify_assert()
193   VERIFY_RAISES_ASSERT(lu.solve(tmp))  in lu_verify_assert()
194   VERIFY_RAISES_ASSERT(lu.transpose().solve(tmp))  in lu_verify_assert()
195   VERIFY_RAISES_ASSERT(lu.adjoint().solve(tmp))  in lu_verify_assert()
196   VERIFY_RAISES_ASSERT(lu.determinant())  in lu_verify_assert()
197   VERIFY_RAISES_ASSERT(lu.rank())  in lu_verify_assert()
198   VERIFY_RAISES_ASSERT(lu.dimensionOfKernel())  in lu_verify_assert()
199   VERIFY_RAISES_ASSERT(lu.isInjective())  in lu_verify_assert()
200   VERIFY_RAISES_ASSERT(lu.isSurjective())  in lu_verify_assert()
201   VERIFY_RAISES_ASSERT(lu.isInvertible())  in lu_verify_assert()
202   VERIFY_RAISES_ASSERT(lu.inverse())  in lu_verify_assert()
214 EIGEN_DECLARE_TEST(lu)  in EIGEN_DECLARE_TEST()  argument