Re: [eigen] A not so simple product

["Benoit Jacob" <jacob.benoit.1@xxxxxxxxx>]

Yes, I confirm that the bug in my example (with the eigenvectors) is fixed too!

Thanks,
Benoit

2008/12/11 Gael Guennebaud <gael.guennebaud@xxxxxxxxx>:
> Benjamin: the segfault you found should be fixed in trunk (no vectorization
> yet).
>
> Benoit: yesterday I did not take the time to check your example, but it
> seems that my last commit also fixed that one.
>
> cheers,
> Gael.
>
> On Wed, Dec 10, 2008 at 3:46 PM, Benoit Jacob <jacob.benoit.1@xxxxxxxxx>
> wrote:
>>
>> I just discovered yet another bug in diagonal product. This part of
>> Eigen is really very buggy!!
>>
>> To grade my exam I made this program:
>>
>> #define EIGEN_DEFAULT_IO_FORMAT EIGEN_DOCS_IO_FORMAT
>> #include <Eigen/QR>
>>
>> using namespace Eigen;
>> using namespace std;
>>
>> int main()
>> {
>>  Matrix3d A;
>>  A << 0,1,1,
>>       1,0,1,
>>       1,1,0;
>>  SelfAdjointEigenSolver<Matrix3d> s(A);
>>  cout << "Eigenvalues:" << endl << s.eigenvalues() << endl;
>>  cout << "Eigenvectors:" << endl << s.eigenvectors() *
>> (s.eigenvectors().row(2).cwise().inverse().asDiagonal()).eval() <<
>> endl;
>> }
>>
>> This program produces the correct result.
>>
>> But if I remove the .eval() it produces the wrong result.
>>
>> Fixing diagonal product is really my top priority in eigen.... will
>> take me a while though as a first have to grade 400 exams.
>>
>> Cheers,
>> Benoit
>>
>> ---
>>
>
>

---