GitHub-Proxy/eigen
1
0
Fork 0
mirror of https://gitlab.com/libeigen/eigen.git synced 2025-08-12 19:59:05 +08:00
Code Issues Packages Projects Releases Wiki Activity

add test for absDeterminant()

Browse Source
This commit is contained in:
Benoit Jacob 2009-08-24 00:35:42 -04:00
parent c9a307f330
commit f31b5a7114
1 changed files with 11 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

11
test/qr_fullpivoting.cpp
View File

@ -68,6 +68,8 @@ template<typename MatrixType> void qr()
template<typename MatrixType> void qr_invertible() template<typename MatrixType> void qr_invertible()
{ {
typedef typename NumTraits<typename MatrixType::Scalar>::Real RealScalar; typedef typename NumTraits<typename MatrixType::Scalar>::Real RealScalar;
typedef typename MatrixType::Scalar Scalar;
int size = ei_random<int>(10,50); int size = ei_random<int>(10,50);
MatrixType m1(size, size), m2(size, size), m3(size, size); MatrixType m1(size, size), m2(size, size), m3(size, size);
@ -88,6 +90,15 @@ template<typename MatrixType> void qr_invertible()
m3 = MatrixType::Random(size,size); m3 = MatrixType::Random(size,size);
VERIFY(qr.solve(m3, &m2)); VERIFY(qr.solve(m3, &m2));
VERIFY_IS_APPROX(m3, m1*m2); VERIFY_IS_APPROX(m3, m1*m2);
// now construct a matrix with prescribed determinant
m1.setZero();
for(int i = 0; i < size; i++) m1(i,i) = ei_random<Scalar>();
RealScalar absdet = ei_abs(m1.diagonal().prod());
m3 = qr.matrixQ(); // get a unitary
m1 = m3 * m1 * m3;
qr.compute(m1);
VERIFY_IS_APPROX(absdet, qr.absDeterminant());
} }
template<typename MatrixType> void qr_verify_assert() template<typename MatrixType> void qr_verify_assert()