This is the mail archive of the gsl-discuss@sourceware.org mailing list for the GSL project.
Index Nav: | [Date Index] [Subject Index] [Author Index] [Thread Index] | |
---|---|---|
Message Nav: | [Date Prev] [Date Next] | [Thread Prev] [Thread Next] |
Other format: | [Raw text] |
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.
Greetings, Jochen -- Einigkeit und Recht und Freiheit http://www.Jochen-Kuepper.de Liberté, Égalité, Fraternité GnuPG key: CC1B0B4D Sex, drugs and rock-n-roll
Attachment:
PGP.sig
Description: This is a digitally signed message part
Index Nav: | [Date Index] [Subject Index] [Author Index] [Thread Index] | |
---|---|---|
Message Nav: | [Date Prev] [Date Next] | [Thread Prev] [Thread Next] |