[Patch libfortran/28452] random intrinsics for large reals
Steve Kargl
sgk@troutmask.apl.washington.edu
Tue Aug 1 02:49:00 GMT 2006
On Mon, Jul 31, 2006 at 10:00:30PM -0400, Jack Howarth wrote:
> Unfortunately the long string of zeros is still present with the
> tonights gcc trunk and the random-patch-5.diff and the change...
(snip)
> As you requested I have added...
>
> print *, x4(1:3) ; print *, x8(1:3); print*, x10(1:3)
>
> ...which I placed before the bin summaries are printed out.
> The results are...
>
0.3090169
0.309016994505426
0.3090169945054261466733706525294548
0.3524981
0.352498068323475
0.3524980683234748589735270183948808
0.8038055
0.803805520479522
0.8038055204795220681406659687936414
This looks good. Of course, doing some statistics on the
last 18 digits would be interesting.
> 195974 196081 196285 196146 195718 196224 196325 196221 196183 195869
> 195974 196079 196287 196145 195717 196226 196325 196221 196182 195870
> 195974 196079 196287 196145 195717 196226 196325 196221 196182 195870
This looks good.
> 46.47314
> 46.4762618000000
> 46.47626180000000000000000000000068
It could be the way I'm computing the chi-squared value or even it's
simply definition.
Thanks.
Thomas, I think is ready to commit.
--
Steve
More information about the Fortran
mailing list