|Re: [eigen] Specialized QR|
[ Thread Index | Date Index | More lists.tuxfamily.org/eigen Archives ]
|> For what sizes?|
3x3 up to 150x150, inclusive rectangular matrices (e.g.: 100x50, 20x10,...). All results are similars.
> Can you send us code reproducing the problem?
> Also, 50 and 100 are HUGE as compile-time fixed sizes. For matrices of
> this size, you should definitely use dynamic size.
Here is how I measured the spent time for QR:
th = walltime(&t0);
Where MyMatrix is Matrix<double, 50, 100> (and A is a MyMatrix).
The program crashes at the initialization of qr_h with the message I posted.
> So: R rectangular and Q "really unitary".
> Notice that i say unitary instead of orthogonal because i want support
> for complex matrices.
With complex I haven't done much tests, but the code should work. I will test and eventually correct the issues.
|Mail converted by MHonArc 2.6.19+||http://listengine.tuxfamily.org/|