gsl_multimin question
Rajarshi Guha
rajarshi@presidency.com
Tue Jun 3 18:33:00 GMT 2003
On Tue, 2003-06-03 at 13:13, Fabrice Rossi wrote:
> Hi.
>
> Could you show the code? Are you sure the derivatives are correct?
I dont have the actual code that was giving the problem, but I modified
the example program given in the manual to use a function similar to the
one in my original code. The same error occurs.
I've attached the C program - essentially I'm trying to minimize a sum
of squares error function. In the attached code the form is
E = sqrt( \sum x_{i} - X_{i} )
where X_{i} is my target value (observed value).
I've selected the initial argument vector randomly, but the values are
close to the target vector. I expected that the values printed in each
iteration should converge to my target vector - but they dont and after
a few iterations all I see is NaN's.
As far as I can see the derivatives are properly defined. Since you
mentioned thgat your collegue gets it to work I must be doing something
wrong somewhere but I cant see it :(
Any pointers would be appreciated.
-------------------------------------------------------------------
Rajarshi Guha <rajarshi@presidency.com> <http://jijo.cjb.net>
GPG Fingerprint: 0CCA 8EE2 2EEB 25E2 AB04 06F7 1BB9 E634 9B87 56EE
-------------------------------------------------------------------
A meeting is an event at which the minutes are kept and the hours are
lost.
-------------- next part --------------
A non-text attachment was scrubbed...
Name: bfgs.c
Type: text/x-c
Size: 2229 bytes
Desc: not available
URL: <http://sourceware.org/pipermail/gsl-discuss/attachments/20030603/007c14a3/attachment.bin>
More information about the Gsl-discuss
mailing list