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 than ideal
Jochen Küpper writes:
I checked the source code for gsl_matrix_free and it will fail if passed
a null pointer.
Shouldn't the GSL free-routines be changed like the following patch?
My thinking on that was that for most people calling the free()
functions on a null pointer is usually an
Hello,
Is it possible to set the maximum step size in the ode (ordinary
differential equation) solvers? It's pretty important that the solver
doesn't start running away... Also, I will have large and sudden changes
in some of the parameters and I think I will need to reset the step size
to