[go: up one dir, main page]
Menu ▾ ▴

[r377]: / fish_shark / ran0.c  Maximize  Restore  History

Download this file

28 lines (25 with data), 419 Bytes 

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
#define IA 16807

#define IM 2147483647

#define AM (1.0/IM)

#define IQ 127773

#define IR 2836

#define MASK 123459876



float ran0(long *idum)

{

	long k;

	float ans;



	*idum ^= MASK;

	k=(*idum)/IQ;

	*idum=IA*(*idum-k*IQ)-IR*k;

	if (*idum < 0) *idum += IM;

	ans=AM*(*idum);

	*idum ^= MASK;

	return ans;

}

#undef IA

#undef IM

#undef AM

#undef IQ

#undef IR

#undef MASK

/* (C) Copr. 1986-92 Numerical Recipes Software )&%. */