This is the mail archive of the fortran@gcc.gnu.org mailing list for the GNU Fortran 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]

Re: [gfortran] PATCH -- Fix segfault in RANDOM_SEED.

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
Index Nav: [Date Index] [Subject Index] [Author Index] [Thread Index]
Message Nav: [Date Prev] [Date Next] [Thread Prev] [Thread Next]