|
From: | Brian Gladman |
Subject: | [Bug-gsl] [bug #39713] roots/secant.c "derivative value is not finite" for a good guess |
Date: | Sat, 05 Oct 2013 10:05:56 +0000 |
User-agent: | Mozilla/5.0 (Windows NT 6.2; WOW64) AppleWebKit/537.36 (KHTML, like Gecko) Chrome/30.0.1599.69 Safari/537.36 |
Follow-up Comment #15, bug #39713 (project gsl): With a small correction to the value of r (= M_PI / M_E), your modified test program (gsl-secant64.c) gives a correct result on Windows x64 using the version of secant.c now in the GSL repository: # i x_i-r x_i - x_p status 0 -7.509e-003 4.822e-002 -2 1 0.000e+000 7.509e-003 -2 2 0.000e+000 0.000e+000 0 # Converged # f(x_i) = 2.2204e-019 Have you tested this revised version of secant.c? _______________________________________________________ Reply to this item at: <http://savannah.gnu.org/bugs/?39713> _______________________________________________ Message sent via/by Savannah http://savannah.gnu.org/
[Prev in Thread] | Current Thread | [Next in Thread] |