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/
