aboutsummaryrefslogtreecommitdiffhomepage
diff options
context:
space:
mode:
authorGravatar Jitse Niesen <jitse@maths.leeds.ac.uk>2010-05-25 12:40:42 +0100
committerGravatar Jitse Niesen <jitse@maths.leeds.ac.uk>2010-05-25 12:40:42 +0100
commite7dc7725549a4a52b234d7c84ead58f27aa7cf76 (patch)
tree583b9c441e44795581b95447b70e5749f017868c
parente7d809d4349fd4048777be71f1c803d0b13f8fe8 (diff)
Respect MaxRowsAtCompileTime in HouseholderSequence::evalTo().
This fixes the failing test nomalloc_4. Also remove a print inserted for debugging in schur_real test.
-rw-r--r--Eigen/src/Householder/HouseholderSequence.h3
-rw-r--r--test/schur_real.cpp1
2 files changed, 2 insertions, 2 deletions
diff --git a/Eigen/src/Householder/HouseholderSequence.h b/Eigen/src/Householder/HouseholderSequence.h
index 4bac484eb..ab7702b14 100644
--- a/Eigen/src/Householder/HouseholderSequence.h
+++ b/Eigen/src/Householder/HouseholderSequence.h
@@ -156,7 +156,8 @@ template<typename VectorsType, typename CoeffsType, int Side> class HouseholderS
{
int vecs = m_actualVectors;
dst.setIdentity(rows(), rows());
- Matrix<Scalar,1,DestType::RowsAtCompileTime> temp(rows());
+ Matrix<Scalar, DestType::RowsAtCompileTime, 1,
+ AutoAlign|ColMajor, DestType::MaxRowsAtCompileTime, 1> temp(rows());
for(int k = vecs-1; k >= 0; --k)
{
int cornerSize = rows() - k - m_shift;
diff --git a/test/schur_real.cpp b/test/schur_real.cpp
index d0aca4308..bcb19c936 100644
--- a/test/schur_real.cpp
+++ b/test/schur_real.cpp
@@ -57,7 +57,6 @@ template<typename MatrixType> void schur(int size = MatrixType::ColsAtCompileTim
RealSchur<MatrixType> schurOfA(A);
MatrixType U = schurOfA.matrixU();
MatrixType T = schurOfA.matrixT();
- std::cout << "T = \n" << T << "\n\n";
verifyIsQuasiTriangular(T);
VERIFY_IS_APPROX(A, U * T * U.transpose());
}