[eigen] Re: Student contribution

[ Thread Index | Date Index | More lists.tuxfamily.org/eigen Archives ]


We have found something that could explain the bug with complexes.
When we do
MatrixType result(upperbid.householderU()* upperbid.bidiagonal().toDenseMatrix() * upperbid.householderV().adjoint()); MatrixType result2(upperbid.householderV()* upperbid.bidiagonal().toDenseMatrix().adjoint() * upperbid.householderU().adjoint())

We would expect result.adjoint() == result2 but we don't. When we print them we have different results. We do not know how to solve it.

We also have difficulties to affect an UpperBidiagonalization householder to a matrix.
The lines
std::cout << UpperBid.houseHolderU() * Identity << std::endl;
MatrixType U = UpperBid.houseHolderU();
std::cout << U * Identity << std::endl;
give two different results (the first one is the one we want).

Since UpperBid.houseHolderU() is an houseHolderU type and not a simple matrix we suppose there might be some particular way to manipulate it although it gives the expected result while we don't try to affect it. What should we do for that?

Team GL27

Jean Ceccato, Pierre Zoppitelli, Gauthier Brun, Nicolas Carre

Mail converted by MHonArc 2.6.19+ http://listengine.tuxfamily.org/