[Bug-gsl] [bug #39713] roots/secant.c "derivative value is not finite" for a good guess

[Brian Gladman]

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/