[Top][All Lists]
[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
[Bug-gsl] Multidimensional Root-Finding (gnewton)
From: |
Allan Leal |
Subject: |
[Bug-gsl] Multidimensional Root-Finding (gnewton) |
Date: |
Fri, 24 Jun 2011 14:05:13 +0100 |
User-agent: |
Mozilla/5.0 (X11; U; Linux x86_64; en-US; rv:1.9.2.17) Gecko/20110424 Thunderbird/3.1.10 |
Hello GSL,
I have detected that the solver *gnewton* does not call the *fdf*
function. It only calls *f* and *df*.
This does not happen in the *newton* solver, which only uses the *fdf*
function.
If you go to the function *gnewton_iterate* on file *gnewton.c*, you
will only see the following calls:
int status = GSL_MULTIROOT_FN_EVAL_F (fdf, state->x_trial, f);
int status = GSL_MULTIROOT_FN_EVAL_DF (fdf, x, J);
However, if you go to the function *newton_iterate* on file *newton.c*,
you will only see the following call instead:
int status = GSL_MULTIROOT_FN_EVAL_F_DF(fdf, x, f, J);
Note that in the above call, the function that defines both the function
and jacobian are called, which may be used to improve performance.
Best regards,
Allan
[Prev in Thread] |
Current Thread |
[Next in Thread] |
- [Bug-gsl] Multidimensional Root-Finding (gnewton),
Allan Leal <=