29 …me MatrixType,template <typename,int> class CholType> void test_chol_update(const MatrixType& symm)  in test_chol_update()  argument
35   MatrixType symmLo = symm.template triangularView<Lower>();  in test_chol_update()
36   MatrixType symmUp = symm.template triangularView<Upper>();  in test_chol_update()
37   MatrixType symmCpy = symm;  in test_chol_update()
44     VectorType vec = VectorType::Random(symm.rows());  in test_chol_update()
78   SquareMatrixType symm =  a0 * a0.adjoint();  in cholesky()  local
83     symm += a1 * a1.adjoint();  in cholesky()
91     SquareMatrixType symmUp = symm.template triangularView<Upper>();  in cholesky()
92     SquareMatrixType symmLo = symm.template triangularView<Lower>();  in cholesky()
95     VERIFY_IS_APPROX(symm, chollo.reconstructedMatrix());  in cholesky()
97     VERIFY_IS_APPROX(symm * vecX, vecB);  in cholesky()
99     VERIFY_IS_APPROX(symm * matX, matB);  in cholesky()
103     VERIFY_IS_APPROX(symm, cholup.reconstructedMatrix());  in cholesky()
105     VERIFY_IS_APPROX(symm * vecX, vecB);  in cholesky()
107     VERIFY_IS_APPROX(symm * matX, matB);  in cholesky()
140       symm = -symm; // test a negative matrix  in cholesky()
143     SquareMatrixType symmUp = symm.template triangularView<Upper>();  in cholesky()
144     SquareMatrixType symmLo = symm.template triangularView<Lower>();  in cholesky()
147     VERIFY_IS_APPROX(symm, ldltlo.reconstructedMatrix());  in cholesky()
149     VERIFY_IS_APPROX(symm * vecX, vecB);  in cholesky()
151     VERIFY_IS_APPROX(symm * matX, matB);  in cholesky()
154     VERIFY_IS_APPROX(symm, ldltup.reconstructedMatrix());  in cholesky()
156     VERIFY_IS_APPROX(symm * vecX, vecB);  in cholesky()
158     VERIFY_IS_APPROX(symm * matX, matB);  in cholesky()
182       symm = -symm;  in cholesky()
187       SquareMatrixType A = symm;  in cholesky()
237   CALL_SUBTEST(( test_chol_update<SquareMatrixType,LLT>(symm)  ));  in cholesky()
238   CALL_SUBTEST(( test_chol_update<SquareMatrixType,LDLT>(symm) ));  in cholesky()
261   RealMatrixType symm =  a0 * a0.adjoint();  in cholesky_cplx()  local
266     symm += a1 * a1.adjoint();  in cholesky_cplx()
270     RealMatrixType symmLo = symm.template triangularView<Lower>();  in cholesky_cplx()
273     VERIFY_IS_APPROX(symm, chollo.reconstructedMatrix());  in cholesky_cplx()
275     VERIFY_IS_APPROX(symm * vecX, vecB);  in cholesky_cplx()
286       symm = -symm; // test a negative matrix  in cholesky_cplx()
289     RealMatrixType symmLo = symm.template triangularView<Lower>();  in cholesky_cplx()
292     VERIFY_IS_APPROX(symm, ldltlo.reconstructedMatrix());  in cholesky_cplx()
294     VERIFY_IS_APPROX(symm * vecX, vecB);  in cholesky_cplx()