eknecronzontas writes:
> The stepsize for finite-differencing in
> gsl-1.8/multiroots/fdjacobian.c is specified by the
> lines
>
> > double xj = gsl_vector_get (x, j);
> > double dx = epsrel * fabs (xj);
>
> This is, of course mathematically correct.
> Nevertheless, the behavior is less
Since I can't quite decide if this is a bug, I'll post
it here first...
The stepsize for finite-differencing in
gsl-1.8/multiroots/fdjacobian.c is specified by the
lines
> double xj = gsl_vector_get (x, j);
> double dx = epsrel * fabs (xj);
This is, of course mathematically correct.
Nevertheless