diff options
author | Benoit Jacob <jacob.benoit.1@gmail.com> | 2010-01-15 00:35:26 -0500 |
---|---|---|
committer | Benoit Jacob <jacob.benoit.1@gmail.com> | 2010-01-15 00:35:26 -0500 |
commit | bfe6fdde2488c67ac8de01603d20167c61a788f2 (patch) | |
tree | 86b3d8250ff5bc707dc3761404d4ea4c575cbe00 /test/upperbidiagonalization.cpp | |
parent | ddc32adb0e9b8683a5070bc78c1a1250d9313e85 (diff) |
allow to multiply a householder sequence and a matrix when one is real and one is complex.
This is especially important as in bidiagonalization, the band matrix is real.
Diffstat (limited to 'test/upperbidiagonalization.cpp')
-rw-r--r-- | test/upperbidiagonalization.cpp | 3 |
1 files changed, 1 insertions, 2 deletions
diff --git a/test/upperbidiagonalization.cpp b/test/upperbidiagonalization.cpp index e6f66a2d5..f64c1aed0 100644 --- a/test/upperbidiagonalization.cpp +++ b/test/upperbidiagonalization.cpp @@ -24,7 +24,6 @@ #include "main.h" #include <Eigen/SVD> -#include <Eigen/LU> template<typename MatrixType> void upperbidiag(const MatrixType& m) { @@ -39,7 +38,7 @@ template<typename MatrixType> void upperbidiag(const MatrixType& m) RealMatrixType b(rows, cols); b.setZero(); b.block(0,0,cols,cols) = ubd.bidiagonal(); - MatrixType c = ubd.householderU() * b.template cast<Scalar>() * ubd.householderV().adjoint(); + MatrixType c = ubd.householderU() * b * ubd.householderV().adjoint(); VERIFY_IS_APPROX(a,c); } |