[gfortran] PATCH -- Fix segfault in RANDOM_SEED.

Steve Kargl sgk@troutmask.apl.washington.edu
Tue May 25 05:03:00 GMT 2004


On Mon, May 24, 2004 at 11:58:31AM +0200, Lars Segerlund wrote:
> 
>  Actually I skipped some priming, I'd better put it back in ...
> I will also look at the put sequence if I have forgotten something.
> 
>  Basicly, the only 'really' bad sequence should be a string of zeroe's
> ... everything else should run fine, I will look at it again.
> 

I really don't understand your random_generate() function, so
I'm of limited help here.  I wrote an equipartition test for
PRNG and tested your MT.  If I seed the array with 

seed(1:624) = (/ (i,i=1,624) /) or 10 * (/ (i,i=1,624) /)

I have to burn 100000 random number before I start to get
"good" random numbers with respect to a chi square test.

I'll send my KISS implementation as soon as I receive email
from Marsaglia with the okay.

-- 
Steve



More information about the Fortran mailing list