[Help-gsl] Re: Higher order derivatives (was: GSL, C++, and function pointers)

[Jordi Gutierrez Hermoso]

On 07/03/07, eknecronzontas <address@hidden> wrote:
>      If you use the function wrapper, you can write a
> second function wrapper to call the GSL derivative
> routine again to generate the second derivative.

Is this really the way the GSL is supposed to be used? It seems
awfully inefficient, making more function calls than the typical three
that are required for an O(h^2) centred difference 2nd derivative. By
my count, using the GSL as suggested above would require six function
calls, and probably would have a different order of precision.

I'm confused,
- Jordi G. H.