Re: GSL on Hitachi SR8000

[Reinhold Bader <Reinhold dot Bader at lrz-muenchen dot de>]

Brian Gough wrote:

> Reinhold Bader writes:
> >            printf("failure of range for i=%d: %20.16f %20.16f\n",i,h->range[i],abs(h->range[i] - hh->range[i]));
>
> abs() returns an integer value, hence the strange floating point value 
> O(1e-318).
>
> Try this instead:
>
> printf("failure i=%d: %.19e %.19e %.19e\n",i,h->range[i], hh->range[i], hh->range[i] - h->range[i])
>  
failure i=1: 2.5188916876574307000e-03 2.5188916876574298000e-03 
-8.6736173798840355000e-19
failure i=2: 5.0377833753148613000e-03 5.0377833753148596000e-03 
-1.7347234759768071000e-18
failure i=3: 7.5566750629722920000e-03 7.5566750629722903000e-03 
-1.7347234759768071000e-18
failure i=4: 1.0075566750629723000e-02 1.0075566750629719000e-02 
-3.4694469519536142000e-18
failure i=5: 1.2594458438287154000e-02 1.2594458438287151000e-02 
-3.4694469519536142000e-18
failure i=6: 1.5113350125944584000e-02 1.5113350125944581000e-02 
-3.4694469519536142000e-18
etc. (up to i=34)



--
Dr. Reinhold Bader

Leibniz-Rechenzentrum, Abt. Benutzerbetreuung	| Tel. +4989 289 28825
Barerstr. 21, 80333 Muenchen			| email Bader@lrz.de