[Top][All Lists]
[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
Re: [Bug-gsl] diff.c (numerical differentiation)
From: |
Brian Gough |
Subject: |
Re: [Bug-gsl] diff.c (numerical differentiation) |
Date: |
Tue, 08 Apr 2008 21:55:31 +0100 |
User-agent: |
Wanderlust/2.14.0 (Africa) Emacs/22.1 Mule/5.0 (SAKAKI) |
At Sun, 6 Apr 2008 13:29:41 -0700,
spasmous wrote:
> h = sqrt (GSL_SQRT_DBL_EPSILON / (2.0 * a2));
> ...
> h = (x+h)-x; // ensure difference between (x+h) and x is exactly h
> *result = (GSL_FN_EVAL (f, x + h) - GSL_FN_EVAL (f, x)) / h;
Thanks for the suggestion.
--
Brian Gough