This is the mail archive of the
gsl-discuss@sources.redhat.com
mailing list for the GSL project.
Re: GSL on Hitachi SR8000
- From: "Jochen Küpper" <jochen at unc dot edu>
- To: Reinhold Bader <Reinhold dot Bader at lrz-muenchen dot de>
- Cc: Brian Gough <bjg at network-theory dot co dot uk>, gsl-discuss at sources dot redhat dot com
- Date: 24 Aug 2002 14:07:02 -0400
- Subject: Re: GSL on Hitachi SR8000
- Organization: University of North Carolina
- References: <3D66010C.6020002@lrz.de> <15718.36054.573517.309284@debian><3D67C620.4050306@lrz.de>
- Reply-to: jochen at jochen-kuepper dot de
On Sat, 24 Aug 2002 19:45:04 +0200 Reinhold Bader wrote:
Reinhold> printf("failure of range for i=%d: %20.16f %20.16f\n",i,h->range[i],abs(h->range[i] - hh->range[i]));
Reinhold> and receive failures for h->range. Oddly enough, the absolute differences are
Reinhold> 0:
Reinhold> failure of range for i=1: 0.0025188916876574 0.0000000000000000
Reinhold> failure of range for i=2: 0.0050377833753149 0.0000000000000000
Reinhold> failure of range for i=3: 0.0075566750629723 0.0000000000000000
Reinhold> failure of range for i=4: 0.0100755667506297 0.0000000000000000
Reinhold> failure of range for i=5: 0.0125944584382872 0.0000000000000000
Reinhold> perhaps the machine differentiates 0.0 and -0.0? Generally, comparing
Reinhold> floats etc. to equality is generally not a good idea (and you do
Reinhold> read in %lg's in gsl_histogram_fscanf)...
Hmm, what's the output if you use %g instead of %f?
Greetings,
Jochen
--
University of North Carolina phone: +1-919-962-4403
Department of Chemistry phone: +1-919-962-1579
Venable Hall CB#3290 (Kenan C148) fax: +1-919-843-6041
Chapel Hill, NC 27599, USA GnuPG key: 44BCCD8E