[eigen] Specialized QR |

[ Thread Index | Date Index | More lists.tuxfamily.org/eigen Archives ]

*To*: Eigen <eigen@xxxxxxxxxxxxxxxxxxx>*Subject*: [eigen] Specialized QR*From*: Andrea Arteaga <yo.eres@xxxxxxxxx>*Date*: Tue, 28 Apr 2009 15:40:37 +0200*Dkim-signature*: v=1; a=rsa-sha256; c=relaxed/relaxed; d=gmail.com; s=gamma; h=domainkey-signature:received:received:from:to:subject:date :user-agent:mime-version:content-type:content-transfer-encoding :content-disposition:message-id; bh=5dXisoUUk8Nh6/OsWIX7Co6ET42VePlmnqhfUE7RByw=; b=xeRui9DNBU4+Ts1xiNMDJwU7ZJN8ZXPQrWAbYn9ilWbUj8k/YQRHOEDBZmwuoE9DAs cArUg9uXpZ4OjPvRvpyvcHvxlYd5AepoT6Y2IIfhH+kweCFH0gHgsyyTfo/UBwKaEsPF nH7r6rhIrorRJTrknzJBrqgmoh+K1bytwd5fA=*Domainkey-signature*: a=rsa-sha1; c=nofws; d=gmail.com; s=gamma; h=from:to:subject:date:user-agent:mime-version:content-type :content-transfer-encoding:content-disposition:message-id; b=m6Lcbtp7B1iH0+H7gwoWfvDikyjzLwrA4mYTB630g9/f7JFZxtF20H4n1ClLqvMDvq mtZkRGOVyBnmzmnHPG2kvunsdd3dflskZ2vyN2vqorQqsigeAyayrvG9oGInBWAUA/6r I0QSfZYrG3dzmTTt3+UKzk6r6U8VSOWTrl0VQ=

I did some work on QR module to do some specialized version for small matrices. My idea is to introduce a class qr_compute, which does all the work (QR decomposition, matrixQ, matrixR, rank, solve,...) whit respect to the choosen algorithm (Householder, Gram-Schmidt or Givens). The QR class were only a front-end to this class. In this way, we don't need to specialize the QR, but only the backend qr_compute class. enum { Algorithm = rows < 3 ? 1 : ((rows < 15 && cols < 15) ? 2 : 3) }; }; - matrixQ()
- matrixR()
- rank()
- solve()
The others public members (i.e.: dimensioOfKernel(), isInjective(), isSurjective(), isInvertible()) don't need to use the compute_ member. If you mean that this work is useful, I can begin to code the three specialized versions of qr_compute, i.e.: write a simply (hard-coded) givens rotation for 2x* matrices, re-code the householder orthogonalization by move the code from QR to qr_compute and write a Gram-Schmidt orthogonalization. |

**Follow-Ups**:**Re: [eigen] Specialized QR***From:*Benoit Jacob

**Messages sorted by:**[ date | thread ]- Prev by Date:
**[eigen] git** - Next by Date:
**Re: [eigen] Specialized QR** - Previous by thread:
**Re: [eigen] git** - Next by thread:
**Re: [eigen] Specialized QR**

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