Re: [eigen] Non linear equations solver |

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

*To*: eigen@xxxxxxxxxxxxxxxxxxx*Subject*: Re: [eigen] Non linear equations solver*From*: Koldo Ramirez <koldo.ramirez@xxxxxxxxx>*Date*: Tue, 11 Aug 2009 08:56:53 +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; bh=1ITCiUPGZPPjpA11iJ2kSYFDqK3UcHhrN2lxviKTssY=; b=gfGXVIzoC9gTvswMYOGL6KM2IuCTwJgnWMfX/G5HcamBM3TtuA1j7peWEUpCmc+jpO yY88oKVjb4dJyM8m2FH32gz9nuhwECCv/Fsn3TkeXCpwB2S/C4gZevJFf0Au7eiMbXem ZgfOKkZfpvlnMt6+JjvWWLnFzmw1SXhehWfco=*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=TNAZghGnwbJxWGFk977/7nnD7l4xJA13oTj2YgRIGtszQbCtLPi5c2lBxIwUSFMKV3 V6yxLf5D/LVJy8+c4VN9hsiIddbmju7odVPrQVLgLdyFm54+ITv1v3sJQacoy1r06Hjo 3bNQ1scXSNxiLo9Dl4StO+qMSY+SXarh2l4D4=

Hello Gael

If I were you I would be a good mathematician and programmer :-)

You and other people in Eigen's group are at least as good as Jorge Moré et al that developed Minpack in the Argonne National Laboratory years ago, a sample of robust and reliable software (http://en.wikipedia.org/wiki/MINPACK)

Something like this means many months of excelent design and programing effort from real experts.

Thomas Capricelli confirm me (and I think he is much better programmer than me) that this is too much difficult for us.

We prefer to offer Eigen's community in weeks a robust well tested library, better than to have in months or years a "best in class" library.

This way in only few weeks we can have something really useful. And we are sure we will do it. We think that re design Minpack from scratch is out of our capability.

Thak you

Koldo

2009/8/10 Gael Guennebaud <gael.guennebaud@xxxxxxxxx>

my 2 cents:

if I were you, I'd start from Keir implementation, extend it to make

the underlying linear solver a template parameter so that one can use

LU, QR, SVD, LLT, etc. or any sparse solver when using sparse

matrices. Actually, we already discussed these things in another

threads.

Then you could have a look at cminpack implementation to see the

differences, and in particular to see whether there are some good

tricks to make the algorithm better (better convergence, more

accurate, etc.)

An in parallel to that, we could discuss the API design, in

particular, how to control the algorithm, how to specify the Jacobian

(user defined, AD, finite differences), etc.

cheers,

gael.

--

On Mon, Aug 10, 2009 at 9:07 AM, Koldo Ramirez<koldo.ramirez@xxxxxxxxx> wrote:

> Hello Thomas

>

> I suppose Minpack 2 is based in version 1 as they present it as an

> improvement focusing in large scale problems and using Blas and Lapack..

>

> About the Octave reference is just to show that an important project like

> Octave has considered Minpack.

>

> Best regards

> Koldo

>

> 2009/8/8 Thomas Capricelli <orzel@xxxxxxxxxxxxxxx>

>>

>>

>> In data lunedì 27 luglio 2009 10:13:58, Koldo Ramirez ha scritto:

>> > It exists Minpack version 2. Between other advantages it uses BLAS and

>> > LAPACK. This is the presentation (

>> > http://ftp.mcs.anl.gov/pub/MINPACK-2/README or

>> > http://www-new.mcs.anl.gov/otc/minpack/summary.html)

>>

>> It doesn't seem to be related to minpack actually. only the name is the

>> same,

>> i think it's not based on the minpack code. is it ?

>>

>>

>> > ..., the old minpack code is gone. ...

>> > Essentially the same method as in minpack is used - double

>> > dogleg

>> > trust region steps with Broyden updating of the jacobian. For

>> > large

>> > enough (no. of variables > 10) square and overdetermined(!)

>> > systems, the QR factorization of the jacobian is updated rather than the

>> > jacobian itself, avoiding costly re-factorizations. Broyden updating may

>> > be

>> > switched off by the user, yielding plain levenberg-marquardt method.

>>

>> interesting to not, those are things we could do later on. I still prefer

>> considering minpack currently has it has a lot broader user test range.

>>

>>

>> Thomas

>> --

>> Thomas Capricelli <orzel@xxxxxxxxxxxxxxx>

>> http://www.freehackers.org/thomas

>>

>>

>

>

Gaël Guennebaud

Iparla - INRIA Bordeaux

(+33)5 40 00 37 95

**References**:**Re: [eigen] Non linear equations solver***From:*Thomas Capricelli

**Re: [eigen] Non linear equations solver***From:*Koldo Ramirez

**Re: [eigen] Non linear equations solver***From:*Gael Guennebaud

**Messages sorted by:**[ date | thread ]- Prev by Date:
**Re: [eigen] Non linear equations solver** - Next by Date:
**[eigen] a couple of RFCs** - Previous by thread:
**Re: [eigen] Non linear equations solver** - Next by thread:
**[eigen] a couple of RFCs**

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