[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
Equivalent of LSQR in matlab
From: |
Jean-Baptiste Poullet |
Subject: |
Equivalent of LSQR in matlab |
Date: |
Thu, 13 Sep 2007 09:12:39 +0200 |
Thanks for your answers! I'm also wondering whether the matlab function
LSQR (based on Conjugate Gradients on the Normal Equations) exists in
octave. I'm dealing with the least squares problem Ax=b, A being a MxN
Vandermonde matrix (complex values), x and b being vectors. Is there any
function in Octave for solving this problem, taking into account the
special structure of A (x=A\b is not sufficient given A is large (M>>N)
and ill conditioned)?
Thx.
JB
On Tue, 2007-09-11 at 22:54 +0200, David Bateman wrote:
> John W. Eaton wrote:
> > On 11-Sep-2007, David Bateman wrote:
> >
> > | David Bateman wrote:
> > | >
> > | > I've often felt that Octave should change to use the same solver as
> > | > Matlab, especially now as the sparse matrices in Octave use a QR solver
> > | > with pivoting and don't give exactly the same results as the dense
> > solvers..
> > |
> > | If you want another reason to change consider Figure 3.3 on the page
> > |
> > | http://www.netlib.org/lapack/lug/node71.html
> > |
> > | for N = 500 dgelss is about 35 times slower than dgelsy
> >
> > Switching is fine with me if the results are approximately as good. I
> > chose dgelss because it was available and I knew about it.
>
> If we want to stay with an SVD solver we should at least change to
> dgelsd that is significantly faster than dgelss and gives essentially
> the same results.
>
> D.
> _______________________________________________
> Help-octave mailing list
> address@hidden
> https://www.cae.wisc.edu/mailman/listinfo/help-octave
___________________________________________________________________________
Yahoo! Mail réinvente le mail ! Découvrez le nouveau Yahoo! Mail et son
interface révolutionnaire.
http://fr.mail.yahoo.com
- Using LAPACK in octave, Jean-Baptiste Poullet, 2007/09/11
- Re: Using LAPACK in octave, Quentin Spencer, 2007/09/11
- Re: Using LAPACK in octave, David Bateman, 2007/09/11
- Re: Using LAPACK in octave, David Bateman, 2007/09/11
- Re: Using LAPACK in octave, John W. Eaton, 2007/09/11
- Re: Using LAPACK in octave, David Bateman, 2007/09/11
- Equivalent of LSQR in matlab,
Jean-Baptiste Poullet <=
- Re: Equivalent of LSQR in matlab, Jordi Gutiérrez Hermoso, 2007/09/13
- Re: Equivalent of LSQR in matlab, John W. Eaton, 2007/09/13
- Re: Using LAPACK in octave, David Bateman, 2007/09/25
- Re: Using LAPACK in octave, John W. Eaton, 2007/09/26
- Re: Using LAPACK in octave, David Bateman, 2007/09/26
Message not available