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

remove dummy code in ColPivHouseholderQR::solve

Browse Source
This commit is contained in:
Gael Guennebaud 2013-01-23 18:34:29 +01:00
parent 19c78cf510
commit ee36eaefc6
1 changed files with 0 additions and 8 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

8
Eigen/src/QR/ColPivHouseholderQR.h
View File

@ -487,14 +487,6 @@ struct solve_retval<ColPivHouseholderQR<_MatrixType>, Rhs>
.template triangularView<Upper>() .template triangularView<Upper>()
.solveInPlace(c.topRows(nonzero_pivots)); .solveInPlace(c.topRows(nonzero_pivots));
typename Rhs::PlainObject d(c);
d.topRows(nonzero_pivots)
= dec().matrixQR()
.topLeftCorner(nonzero_pivots, nonzero_pivots)
.template triangularView<Upper>()
* c.topRows(nonzero_pivots);
for(Index i = 0; i < nonzero_pivots; ++i) dst.row(dec().colsPermutation().indices().coeff(i)) = c.row(i); for(Index i = 0; i < nonzero_pivots; ++i) dst.row(dec().colsPermutation().indices().coeff(i)) = c.row(i);
for(Index i = nonzero_pivots; i < cols; ++i) dst.row(dec().colsPermutation().indices().coeff(i)).setZero(); for(Index i = nonzero_pivots; i < cols; ++i) dst.row(dec().colsPermutation().indices().coeff(i)).setZero();
} }