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

Add tests for sparseQR results (value and size) covering bugs #1522 and #1544

Browse Source
This commit is contained in:
Jeff Trull 2018-04-21 10:26:30 -07:00
parent 9f0c5c3669
commit e7147f69ae
1 changed files with 22 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

22
test/sparseqr.cpp
View File

@ -54,6 +54,28 @@ template<typename Scalar> void test_sparseqr_scalar()
b = dA * DenseVector::Random(A.cols()); b = dA * DenseVector::Random(A.cols());
solver.compute(A); solver.compute(A);
// Q should be MxM
VERIFY_IS_EQUAL(solver.matrixQ().rows(), A.rows());
VERIFY_IS_EQUAL(solver.matrixQ().cols(), A.rows());
// R should be MxN
VERIFY_IS_EQUAL(solver.matrixR().rows(), A.rows());
VERIFY_IS_EQUAL(solver.matrixR().cols(), A.cols());
// Q and R can be multiplied
DenseMat recoveredA = solver.matrixQ()
* DenseMat(solver.matrixR().template triangularView<Upper>())
* solver.colsPermutation().transpose();
VERIFY_IS_EQUAL(recoveredA.rows(), A.rows());
VERIFY_IS_EQUAL(recoveredA.cols(), A.cols());
// and in the full rank case the original matrix is recovered
if (solver.rank() == A.cols())
{
VERIFY_IS_APPROX(A, recoveredA);
}
if(internal::random<float>(0,1)>0.5f) if(internal::random<float>(0,1)>0.5f)
solver.factorize(A); // this checks that calling analyzePattern is not needed if the pattern do not change. solver.factorize(A); // this checks that calling analyzePattern is not needed if the pattern do not change.
if (solver.info() != Success) if (solver.info() != Success)