[Top][All Lists]

[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

[Bug-gsl] [bug #42219] Division by zero in "gnewton" when "f" and "fdf"

From: Patrick Alken
Subject: [Bug-gsl] [bug #42219] Division by zero in "gnewton" when "f" and "fdf" differ
Date: Mon, 28 Apr 2014 16:10:42 +0000
User-agent: Mozilla/5.0 (X11; Linux i686 on x86_64; rv:28.0) Gecko/20100101 Firefox/28.0

Follow-up Comment #1, bug #42219 (project gsl):

The problem seems to be in this line of gnewton_set() in

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:


  Message sent via/by Savannah

reply via email to

[Prev in Thread] Current Thread [Next in Thread]