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

Schur decomposition of 1-by-1 always converges.

Browse Source
This commit is contained in:
Jitse Niesen 2010-06-04 09:40:35 +01:00
parent 9178e2bd54
commit 1ff1bd69ac
2 changed files with 14 additions and 8 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

11
test/schur_complex.cpp
View File

@ -68,10 +68,13 @@ template<typename MatrixType> void schur(int size = MatrixType::ColsAtCompileTim
VERIFY_IS_EQUAL(cs1.matrixT(), csOnlyT.matrixT()); VERIFY_IS_EQUAL(cs1.matrixT(), csOnlyT.matrixT());
VERIFY_RAISES_ASSERT(csOnlyT.matrixU()); VERIFY_RAISES_ASSERT(csOnlyT.matrixU());
// Test matrix with NaN if (size > 1)
A(0,0) = std::numeric_limits<typename MatrixType::RealScalar>::quiet_NaN(); {
ComplexSchur<MatrixType> csNaN(A); // Test matrix with NaN
VERIFY_IS_EQUAL(csNaN.info(), NoConvergence); A(0,0) = std::numeric_limits<typename MatrixType::RealScalar>::quiet_NaN();
ComplexSchur<MatrixType> csNaN(A);
VERIFY_IS_EQUAL(csNaN.info(), NoConvergence);
}
} }
void test_schur_complex() void test_schur_complex()

11
test/schur_real.cpp
View File

@ -85,10 +85,13 @@ template<typename MatrixType> void schur(int size = MatrixType::ColsAtCompileTim
VERIFY_IS_EQUAL(rs1.matrixT(), rsOnlyT.matrixT()); VERIFY_IS_EQUAL(rs1.matrixT(), rsOnlyT.matrixT());
VERIFY_RAISES_ASSERT(rsOnlyT.matrixU()); VERIFY_RAISES_ASSERT(rsOnlyT.matrixU());
// Test matrix with NaN if (size > 1)
A(0,0) = std::numeric_limits<typename MatrixType::Scalar>::quiet_NaN(); {
RealSchur<MatrixType> rsNaN(A); // Test matrix with NaN
VERIFY_IS_EQUAL(rsNaN.info(), NoConvergence); A(0,0) = std::numeric_limits<typename MatrixType::Scalar>::quiet_NaN();
RealSchur<MatrixType> rsNaN(A);
VERIFY_IS_EQUAL(rsNaN.info(), NoConvergence);
}
} }
void test_schur_real() void test_schur_real()