CFD/eigen
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

Respect MaxRowsAtCompileTime in HouseholderSequence::evalTo().

Browse Source 
This fixes the failing test nomalloc_4.
Also remove a print inserted for debugging in schur_real test.
This commit is contained in:
Jitse Niesen 2010-05-25 12:40:42 +01:00
parent e7d809d434
commit e7dc772554
2 changed files with 2 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

3
Eigen/src/Householder/HouseholderSequence.h
View File

@ -156,7 +156,8 @@ template<typename VectorsType, typename CoeffsType, int Side> class HouseholderS
{ {
int vecs = m_actualVectors; int vecs = m_actualVectors;
dst.setIdentity(rows(), rows()); 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) for(int k = vecs-1; k >= 0; --k)
{ {
int cornerSize = rows() - k - m_shift; int cornerSize = rows() - k - m_shift;

1
test/schur_real.cpp
View File

@ -57,7 +57,6 @@ template<typename MatrixType> void schur(int size = MatrixType::ColsAtCompileTim
RealSchur<MatrixType> schurOfA(A); RealSchur<MatrixType> schurOfA(A);
MatrixType U = schurOfA.matrixU(); MatrixType U = schurOfA.matrixU();
MatrixType T = schurOfA.matrixT(); MatrixType T = schurOfA.matrixT();
std::cout << "T = \n" << T << "\n\n";
verifyIsQuasiTriangular(T); verifyIsQuasiTriangular(T);
VERIFY_IS_APPROX(A, U * T * U.transpose()); VERIFY_IS_APPROX(A, U * T * U.transpose());
} }