blob: 9c6cfcda2acda78ca75f2ef3908f5d226b384f02 (
plain)
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
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
|
#include <stdlib.h>
static randbuf rand48buf;
#define A_0 0xE66D
#define A_1 0xDEEC
#define A_2 0x5
#define C 0xB
static randbuf a = { A_0, A_1, A_2 };
static unsigned short c = C;
static void calc_next(randbuf buf) {
randbuf tmp;
long t;
t = buf[0] * a[0] + c;
tmp[0] = t & 0xffff;
tmp[1] = (t >> 16) & 0xffff;
t = buf[1] * a[0] + buf[0] * a[1] + tmp[1];
tmp[1] = t & 0xffff;
tmp[2] = (t >> 16) & 0xffff;
t = buf[2] * a[0] + buf[1] * a[1] + buf[0] * a[2] + tmp[2];
tmp[2] = t & 0xffff;
buf[0] = tmp[0];
buf[1] = tmp[1];
buf[2] = tmp[2];
}
double drand48(void) {
return erand48(rand48buf);
}
long lrand48(void) {
return nrand48(rand48buf);
}
long mrand48(void) {
return jrand48(rand48buf);
}
void srand48(long seed) {
rand48buf[1] = (seed >> 16) & 0xffff;
rand48buf[2] = seed & 0xffff;
rand48buf[0] = 0x330e;
a[0] = A_0;
a[1] = A_1;
a[2] = A_2;
c = C;
}
unsigned short *seed48(randbuf buf) {
static randbuf oldx;
int i;
for (i = 0; i < 3; i++) {
oldx[i] = rand48buf[i];
rand48buf[i] = buf[i];
}
a[0] = A_0;
a[1] = A_1;
a[2] = A_2;
c = C;
return (unsigned short *)&oldx;
}
void lcong48(unsigned short param[7]) {
int i;
for (i = 0; i < 3; i++) {
rand48buf[i] = param[i];
a[i] = param[i + 3];
}
c = param[6];
}
long jrand48(randbuf buf) {
long ret;
ret = buf[2] << 16 | buf[1];
calc_next(buf);
return ret;
}
long nrand48(randbuf buf) {
return jrand48(buf) & 0x7FFFFFFFL;
}
double erand48(randbuf buf) {
double ret;
ret = ((buf[0] / 65536.0 + buf[1]) / 65536.0 + buf[2]) / 65536.0;
calc_next(buf);
return ret;
}
|