Lines Matching refs:MatrixType
29 template<typename MatrixType,template <typename,int> class CholType> void test_chol_update(const MatrixType& symm)
31 typedef typename MatrixType::Scalar Scalar;
32 typedef typename MatrixType::RealScalar RealScalar;
33 typedef Matrix<Scalar, MatrixType::RowsAtCompileTime, 1> VectorType;
35 MatrixType symmLo = symm.template triangularView<Lower>();
36 MatrixType symmUp = symm.template triangularView<Upper>();
37 MatrixType symmCpy = symm;
39 CholType<MatrixType,Lower> chollo(symmLo);
40 CholType<MatrixType,Upper> cholup(symmUp);
49 CholType<MatrixType,Lower> chol(symmCpy);
61 template<typename MatrixType> void cholesky(const MatrixType& m)
63 typedef typename MatrixType::Index Index;
70 typedef typename MatrixType::Scalar Scalar;
72 typedef Matrix<Scalar, MatrixType::RowsAtCompileTime, MatrixType::RowsAtCompileTime> SquareMatrixType;
73 typedef Matrix<Scalar, MatrixType::RowsAtCompileTime, 1> VectorType;
75 MatrixType a0 = MatrixType::Random(rows,cols);
77 MatrixType matB = MatrixType::Random(rows,cols), matX(rows,cols);
82 MatrixType a1 = MatrixType::Random(rows,cols);
109 MatrixType neg = -symmLo;
113 VERIFY_IS_APPROX(MatrixType(chollo.matrixL().transpose().conjugate()), MatrixType(chollo.matrixU()));
114 VERIFY_IS_APPROX(MatrixType(chollo.matrixU().transpose().conjugate()), MatrixType(chollo.matrixL()));
115 VERIFY_IS_APPROX(MatrixType(cholup.matrixL().transpose().conjugate()), MatrixType(cholup.matrixU()));
116 VERIFY_IS_APPROX(MatrixType(cholup.matrixU().transpose().conjugate()), MatrixType(cholup.matrixL()));
145 VERIFY_IS_APPROX(MatrixType(ldltlo.matrixL().transpose().conjugate()), MatrixType(ldltlo.matrixU()));
146 VERIFY_IS_APPROX(MatrixType(ldltlo.matrixU().transpose().conjugate()), MatrixType(ldltlo.matrixL()));
147 VERIFY_IS_APPROX(MatrixType(ldltup.matrixL().transpose().conjugate()), MatrixType(ldltup.matrixU()));
148 VERIFY_IS_APPROX(MatrixType(ldltup.matrixU().transpose().conjugate()), MatrixType(ldltup.matrixL()));
150 if(MatrixType::RowsAtCompileTime==Dynamic)
171 MatrixType m1 = MatrixType::Random(rows,cols), m2(rows,cols);
190 template<typename MatrixType> void cholesky_cplx(const MatrixType& m)
197 typedef typename MatrixType::Index Index;
202 typedef typename MatrixType::Scalar Scalar;
204 typedef Matrix<RealScalar, MatrixType::RowsAtCompileTime, MatrixType::RowsAtCompileTime> RealMatrixType;
205 typedef Matrix<Scalar, MatrixType::RowsAtCompileTime, 1> VectorType;
209 MatrixType matB = MatrixType::Random(rows,cols), matX(rows,cols);
250 template<typename MatrixType> void cholesky_bug241(const MatrixType& m)
254 typedef typename MatrixType::Scalar Scalar;
255 typedef Matrix<Scalar, MatrixType::RowsAtCompileTime, 1> VectorType;
257 MatrixType matA;
265 template<typename MatrixType> void cholesky_verify_assert()
267 MatrixType tmp;
269 LLT<MatrixType> llt;
275 LDLT<MatrixType> ldlt;