Re: [eigen] Re: LU precision tuning

[ Thread Index | Date Index | More Archives ]

I could not help but thought about it again. The way we are currently
creating orthogonal matrices involves the QR decomposition and the
eps2 in QR was always computed as

RealScalar eps2 = precision<RealScalar>()*precision<RealScalar>();

Maybe this caused problems with the orthogonal matrices computed via
QR. I am wondering whether this code

maybe more appropriate for the creation of random orthogonal matrices.
The code uses Gram-Schmidt orthogonalization to do so.

I will give it a try when I come back from lunch.

- Hauke

Mail converted by MHonArc 2.6.19+