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