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

Add missing calls to tests of COD.

Browse Source 
Fix a few mistakes in 3.2 -> 3.3 port.
This commit is contained in:
Rasmus Munk Larsen 2016-02-08 08:50:34 -08:00
parent 16ec450ca1
commit 414efa47d3
2 changed files with 11 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

4
Eigen/src/QR/CompleteOrthogonalDecomposition.h
View File

@ -480,7 +480,7 @@ template <typename _MatrixType>
template <typename RhsType, typename DstType> template <typename RhsType, typename DstType>
void CompleteOrthogonalDecomposition<_MatrixType>::_solve_impl( void CompleteOrthogonalDecomposition<_MatrixType>::_solve_impl(
const RhsType& rhs, DstType& dst) const { const RhsType& rhs, DstType& dst) const {
eigen_assert(rhs().rows() == this->rows()); eigen_assert(rhs.rows() == this->rows());
const Index rank = this->rank(); const Index rank = this->rank();
if (rank == 0) { if (rank == 0) {
@ -491,7 +491,7 @@ void CompleteOrthogonalDecomposition<_MatrixType>::_solve_impl(
// Compute c = Q^* * rhs // Compute c = Q^* * rhs
// Note that the matrix Q = H_0^* H_1^*... so its inverse is // Note that the matrix Q = H_0^* H_1^*... so its inverse is
// Q^* = (H_0 H_1 ...)^T // Q^* = (H_0 H_1 ...)^T
typename RhsType::PlainObject c(rhs()); typename RhsType::PlainObject c(rhs);
c.applyOnTheLeft( c.applyOnTheLeft(
householderSequence(matrixQTZ(), hCoeffs()).setLength(rank).transpose()); householderSequence(matrixQTZ(), hCoeffs()).setLength(rank).transpose());

9
test/qr_colpivoting.cpp
View File

@ -289,6 +289,15 @@ void test_qr_colpivoting()
CALL_SUBTEST_5(( qr_fixedsize<Matrix<double,1,1>, 1 >() )); CALL_SUBTEST_5(( qr_fixedsize<Matrix<double,1,1>, 1 >() ));
} }
for(int i = 0; i < g_repeat; i++) {
CALL_SUBTEST_1( cod<MatrixXf>() );
CALL_SUBTEST_2( cod<MatrixXd>() );
CALL_SUBTEST_3( cod<MatrixXcd>() );
CALL_SUBTEST_4(( cod_fixedsize<Matrix<float,3,5>, 4 >() ));
CALL_SUBTEST_5(( cod_fixedsize<Matrix<double,6,2>, 3 >() ));
CALL_SUBTEST_5(( cod_fixedsize<Matrix<double,1,1>, 1 >() ));
}
for(int i = 0; i < g_repeat; i++) { for(int i = 0; i < g_repeat; i++) {
CALL_SUBTEST_1( qr_invertible<MatrixXf>() ); CALL_SUBTEST_1( qr_invertible<MatrixXf>() );
CALL_SUBTEST_2( qr_invertible<MatrixXd>() ); CALL_SUBTEST_2( qr_invertible<MatrixXd>() );