diff options
author | Gael Guennebaud <g.gael@free.fr> | 2013-06-17 00:14:42 +0200 |
---|---|---|
committer | Gael Guennebaud <g.gael@free.fr> | 2013-06-17 00:14:42 +0200 |
commit | 55365566b2f81656a58e3d3dad1a82ec3038ca64 (patch) | |
tree | 458e647ee49df1e3d0d48504713f4767a9b25235 /Eigen/src/SVD/UpperBidiagonalization.h | |
parent | 9f11f80db1b53d0f9b9173a777c0a644e763619a (diff) |
Fix HouseholderSequence::conjugate() and ::adjoint() and add respective unit tests.
Diffstat (limited to 'Eigen/src/SVD/UpperBidiagonalization.h')
-rw-r--r-- | Eigen/src/SVD/UpperBidiagonalization.h | 4 |
1 files changed, 2 insertions, 2 deletions
diff --git a/Eigen/src/SVD/UpperBidiagonalization.h b/Eigen/src/SVD/UpperBidiagonalization.h index 213b3100d..587de37a5 100644 --- a/Eigen/src/SVD/UpperBidiagonalization.h +++ b/Eigen/src/SVD/UpperBidiagonalization.h @@ -39,7 +39,7 @@ template<typename _MatrixType> class UpperBidiagonalization CwiseUnaryOp<internal::scalar_conjugate_op<Scalar>, const Diagonal<const MatrixType,0> > > HouseholderUSequenceType; typedef HouseholderSequence< - const MatrixType, + const typename internal::remove_all<typename MatrixType::ConjugateReturnType>::type, Diagonal<const MatrixType,1>, OnTheRight > HouseholderVSequenceType; @@ -74,7 +74,7 @@ template<typename _MatrixType> class UpperBidiagonalization const HouseholderVSequenceType householderV() // const here gives nasty errors and i'm lazy { eigen_assert(m_isInitialized && "UpperBidiagonalization is not initialized."); - return HouseholderVSequenceType(m_householder, m_householder.const_derived().template diagonal<1>()) + return HouseholderVSequenceType(m_householder.conjugate(), m_householder.const_derived().template diagonal<1>()) .setLength(m_householder.cols()-1) .setShift(1); } |