Re: [eigen] LU decomposition with partial pivoting ? |

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

*To*: eigen@xxxxxxxxxxxxxxxxxxx*Subject*: Re: [eigen] LU decomposition with partial pivoting ?*From*: Benoit Jacob <jacob.benoit.1@xxxxxxxxx>*Date*: Tue, 11 May 2010 11:51:30 -0400*Dkim-signature*: v=1; a=rsa-sha256; c=relaxed/relaxed; d=gmail.com; s=gamma; h=domainkey-signature:mime-version:received:received:in-reply-to :references:date:message-id:subject:from:to:content-type; bh=ks1EBgM/TBcPGW2lRUJgAMx9yfakSk0Xp8JjwSY9l/E=; b=ogMy45fXbVdWwfeY8vg0HtFKA2hKV6X0p/WPtb9E6GPkaIv09t/k5O6KCPm4MBQHBX XxMycaQW7Qc1+cyXPQqDCvYVBctYqm9uXEG58AAkYVfd/BHwYX0Gnxf5P2GLAsf1aU0M 0qdq4Nv2xf7EQdvow52J/soRtKgp3CdmI9agE=*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; b=Icjys08N/eE60knritl3fSA4gEm/DMOBG81lFKvRFmcI5VJ29abs3rVVKismu8bh3L gVsyHHRKg5Zcjuy1s2fRTEpH4yB357GpBqWUfBx98k/ErswD/2q6V6QK24OxuTkT4cWL pTYa9XfCey8KilWYRbV6ZAvRGs9W/zM/6j5Rc=

partial pivoting LU for non-square is not a big priority for a couple of reasons. It needs to be done differently depending on whether rows>cols or cols>rows. LAPACK doesn't do it (afaicr). Then the use cases are very few. Indeed, partial pivoting LU is mostly used for solving, but exact solving on rectangular matrices is rather uncommon, i.e. with rectangular matrices on typically uses least squares solving. Benoit 2010/5/11 <vincent.lejeune@xxxxxxxxxx>: > > Hi, > > > > Eigen tells me that it does not support partial pivoting for non square > > Matrix...Is this a definitive decision from dev team, or will it change in > > a future release ? > > > > On Wed, 5 May 2010 16:52:52 +0200, Gael Guennebaud > > <gael.guennebaud@xxxxxxxxx> wrote: > >> Hi, > >> > >> in the devel branch we already have partial pivoting LU with a quite > >> efficient blocked implementation. > >> > >> I strongly encourage you to switch to the devel branch. Perhaps, the > > only > >> annoying thing with the devel branch is that the doc is not up to date > >> everywhere... > >> > >> gael. > >> > >> > >> > >> On Wed, May 5, 2010 at 4:41 PM, <vincent.lejeune@xxxxxxxxxx> wrote: > >> > >>> > >>> Hi, > >>> > >>> > >>> > >>> it seems that the LU decomposition routine provided by Eigen are doing > >>> > >>> full pivoting. > >>> > >>> However, LU decomposition in Lapack is done with only partial pivoting > >>> > >>> (see DGETRF). > >>> > >>> > >>> > >>> I was wondering if it was planned to expose a routine that do only > >>> partial > >>> > >>> pivoting (ie swap of rows only, not columns). > >>> > >>> Although full pivoting give better result in term of precisions, it is > >>> > >>> usually slower, because of cache unfriendlyness. On the other hand, > >>> partial > >>> > >>> pivoting allow to "tile" the LU algorithm, allowing more > > parallelisation > >>> > >>> for very big matrix...It would be easier to use gpu computing > >>> capabilities > >>> > >>> if the algorithm do partial pivoting only. > >>> > >>> > >>> > >>> Regards, > >>> > >>> Vincent. > >>> > >>> > >>> > > >

**References**:**[eigen] LU decomposition with partial pivoting ?***From:*vincent.lejeune

**Re: [eigen] LU decomposition with partial pivoting ?***From:*Gael Guennebaud

**Re: [eigen] LU decomposition with partial pivoting ?***From:*vincent.lejeune

**Messages sorted by:**[ date | thread ]- Prev by Date:
**Re: [eigen] Indexes: why signed instead of unsigned?** - Next by Date:
**Re: [eigen] Indexes: why signed instead of unsigned?** - Previous by thread:
**Re: [eigen] LU decomposition with partial pivoting ?** - Next by thread:
**[eigen] Array as the Scalar type**

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