Re: Unusual behaviour in gsl_diff_central?

[Brian Gough <bjg at network-theory dot co dot uk>]

eknecronzontas writes:
 > Hello!
 > 
 >      I notice that gsl_diff_central seems to fail
 > for large values of the function argument (see code
 > below). (I am using gsl 1.3 on Redhat 6.2.) I would 
 > expect it to return a non-zero value in this case,
 > but it returns zero, in spite of trying to evaluate
 > the function at x=NaN.
 >      Am I missing something here?

Thanks for the bug report.  Quickly looking at the source it
seems that the algorithm doesn't handle the case of large
arguments.  We'll either have to note it in the documentation or
change the algorithm to be more adaptive.

-- 
Brian