Re: [eigen] qr().solve() implementation

[ Thread Index | Date Index | More Archives ]


thanks for the patch, however there were some shortcomings because you
assumed the input matrix was full rank.
here is the commit log:

[10:46] <CIA-60> ggael * r917639 eigen2/trunk/kdesupport/eigen2/ (4
files in 3 dirs):
[10:46] <CIA-60> QR: add isInjective(), isSurjective(),
[10:46] <CIA-60> mark isFullRank() deprecated,
[10:46] <CIA-60> add solve() (mix of Keir's patch and LU::solve())
[10:46] <CIA-60> => there is big problem with complex which are not working


On Tue, Jan 27, 2009 at 7:33 PM, Keir Mierle <mierle@xxxxxxxxx> wrote:
> Here is a (incomplete) implementation of qr().solve(b, &x).
> There is something wrong with the generated docs; it seems to take the
> results from LU_solve.cpp rather than QR_solve.cpp but I don't have
> time to figure this out now.
> It is probably wasteful in that it explicitly computes matrixQ() and
> matrixR(). However, it does work.
> Keir

Mail converted by MHonArc 2.6.19+