Re: [eigen] Re: LU precision tuning |
[ Thread Index |
Date Index
| More lists.tuxfamily.org/eigen Archives
]
- To: eigen@xxxxxxxxxxxxxxxxxxx
- Subject: Re: [eigen] Re: LU precision tuning
- From: Hauke Heibel <hauke.heibel@xxxxxxxxxxxxxx>
- Date: Mon, 11 May 2009 11:40:09 +0200
- Dkim-signature: v=1; a=rsa-sha256; c=relaxed/relaxed; d=googlemail.com; s=gamma; h=domainkey-signature:mime-version:received:in-reply-to:references :date:message-id:subject:from:to:content-type :content-transfer-encoding; bh=rJNBJCl3kDxX+DrAZm79O15kzDN9svbcOuCG80+K1Wc=; b=uB0eKd2xA6n3n10V1mPocJWmJjtdHQqBraep9DePsMs7Ulg1uW5KIMnH7whyYzz4Ru jQ0gI2HZwPtgU/6MzOntMdZXdl3tsw7VN09lkYJg0ogVkks167kaRVktBSAGeaFupgIQ HNsHhmWYINRzHdoTHTNr1Hhc/ml7nFdrSnfBc=
- Domainkey-signature: a=rsa-sha1; c=nofws; d=googlemail.com; s=gamma; h=mime-version:in-reply-to:references:date:message-id:subject:from:to :content-type:content-transfer-encoding; b=hbpgGHTckydH26dsdoo52XXsUc2PIgs6BLGjOoeNtBQ5Uq/d+ppdxanyjRYqCCbNdO DWOEU99JCcvLkkrGGuUaG+r051X32QFsgzSjytFMchA1cQ248N9Nq2AWoQ7x7LVctqSr uZLSOfLCFssa5aQbv8/zWloBhl8s3EMlx+SYQ=
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
http://rafb.net/p/K58Ni284.html
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