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*: Benoit Jacob <jacob.benoit.1@xxxxxxxxx>*Date*: Mon, 11 May 2009 19:10:29 +0200*Dkim-signature*: v=1; a=rsa-sha256; c=relaxed/relaxed; d=gmail.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=Yt8eM72+HJFUZ6L37hEub8secjAX7vAUbKRhoeHIScU=; b=Gkl7hDNt2Pa5lNXCUMPt0RgY5ChEwKMNZ8kb8in5bu3BDCh9+lQRDiWPNALyWIJHf0 2Ap7K8WUU9L9e3IGw6+qtl9Qc1os//xIa3p7LHpFIKh0j3up2ehaZUYFLT/MPyzusbHL lHsOd3bYLrpEH1HKCe14usrPLZwdEtZGF00JM=*Domainkey-signature*: a=rsa-sha1; c=nofws; d=gmail.com; s=gamma; h=mime-version:in-reply-to:references:date:message-id:subject:from:to :content-type:content-transfer-encoding; b=N9W19DInQ0GjI8Ll8A7t+FqHq95i0I1ANgbjFM3+H7Do4ANQSMvypz+njtVhTnhDjM HfdGDj35GmSBVJf4vJRvScn6YSRTJgiU0xuFrEGB6aKFt6KHieZzyhCfDieKfln7e4S0 eVbDVVLBv0hVS33wF5i6HLlUtxzJ60BGr2Lzs=

2009/5/11, Hauke Heibel <hauke.heibel@xxxxxxxxxxxxxx>: >> ... but I rather think that what is invalidated is the QR precision >> testing! >> Indeed your patch introduced m_precision in QR but it was never used >> in the computation of the QR decomposition, only in the rank() and >> solve() methods. So what's really needed is to make eps2 honor the >> m_precision. >> >> I can't tell whether the creation of random rank r matrices was >> affected by that or not. > > It seems as if it was not really affected. After fixing the precision > usage, the results did not improve. Ah, ok. > I took a short look at the code an briefly touched the literature > about RRQR (rank revealing QR) decompositions. Since the current QR > does (please correct me if I am wrong) no pivoting I assume it is > pretty safe to assume that it is not rank revealing. There are a bunch > of RRQR algorithms reviewed/described in Bischof & Quintana-Qrti, > "Computing Rank-Revealing QR Factorizations of Dense Matrices" and all > of them rely on pivoting. Indeed, the current QR doesn't do any pivoting: that completely escaped me (disclaimer -- i didn't write any of this code). So yes, it's clearly not rank-revealing and i'm amazed then how well it performed in your tests -- interesting. I have no experience with QR so I don't realize how much pivoting is essential or not for QR. So I can't answer right away the following important questions: - what is the domain of application of a no-pivoting QR? Is it restricted to only invertible matrices or something like that? (the fact that our unittest passes doesn't say much on precision, as our tests use a very low precision) - how would the domain of application be expanded if we added partial (column) pivoting? Just asking in case anyone here knows. I'm not saying we must have sorted this completely for 2.1 but at least if we don't change the QR, the API and documentation should be adapted to its domain of applicability. Thanks for the preliminary research on RRQR, it'll be useful if someone decides to implement it. Cheers, Benoit

**Follow-Ups**:**Re: [eigen] Re: LU precision tuning***From:*Jitse Niesen

**References**:**[eigen] LU precision tuning***From:*Benoit Jacob

**Re: [eigen] Re: LU precision tuning***From:*Benoit Jacob

**Re: [eigen] Re: LU precision tuning***From:*Hauke Heibel

**Re: [eigen] Re: LU precision tuning***From:*Benoit Jacob

**Re: [eigen] Re: LU precision tuning***From:*Benoit Jacob

**Re: [eigen] Re: LU precision tuning***From:*Benoit Jacob

**Re: [eigen] Re: LU precision tuning***From:*Hauke Heibel

**Re: [eigen] Re: LU precision tuning***From:*Hauke Heibel

**Re: [eigen] Re: LU precision tuning***From:*Benoit Jacob

**Re: [eigen] Re: LU precision tuning***From:*Hauke Heibel

**Messages sorted by:**[ date | thread ]- Prev by Date:
**Re: [eigen] Re: LU precision tuning** - Next by Date:
**Re: [eigen] Special Matrices: first get diagonal matrices done** - Previous by thread:
**Re: [eigen] Re: LU precision tuning** - Next by thread:
**Re: [eigen] Re: LU precision tuning**

Mail converted by MHonArc 2.6.19+ | http://listengine.tuxfamily.org/ |