This is the mail archive of the
gsl-discuss@sourceware.org
mailing list for the GSL project.
Re: Bug report and fix for adaptive step size control in ODE solving
- From: Frank Reininghaus <frank78ac at googlemail dot com>
- To: Jochen Küpper <jochen at fhi-berlin dot mpg dot de>
- Cc: gsl-discuss at sourceware dot org
- Date: Wed, 02 Jan 2008 19:46:06 +0100
- Subject: Re: Bug report and fix for adaptive step size control in ODE solving
- References: <4772D239.9070206@googlemail.com> <0E55B479-624E-4F88-B511-CDFF93435D9E@fhi-berlin.mpg.de>
Hi Jochen,
Jochen Küpper wrote:
However, why don't you test what you really describe your comment?
This would look like the following (untested):
h0 = GSL_MAX_DBL(GSL_DBL_MIN, (t0 != t0+h0) ? h0 : GSL_DBL_EPSILON*t0);
Actually, I had been thinking of something like that first, but then I
was afraid that an optimising compiler might transform 't0 != t0+h0'
into '0 != h0', making the patch useless for the 'h0 <
GSL_DBL_EPSILON*t0' problem.
Regards,
Frank