Follow-up Comment #1, bug #42219 (project gsl):
The problem seems to be in this line of gnewton_set() in
multiroots/gnewton.c:
122 state->phi = enorm(f);
In the example program, the initial guess x is the solution vector, so the
F_DF call provides a vector f = 0.
Then in the first gnewton_iterate() call, on line
177 double theta = phi1 / phi0;
we divide by 0.
_______________________________________________________
Reply to this item at:
<http://savannah.gnu.org/bugs/?42219>
_______________________________________________
Message sent via/by Savannah
http://savannah.gnu.org/
