|[eigen] Re: Student contribution|
[ Thread Index |
| More lists.tuxfamily.org/eigen Archives
We have found something that could explain the bug with complexes.
When we do
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.
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?
Jean Ceccato, Pierre Zoppitelli, Gauthier Brun, Nicolas Carre