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
90
91
92
93
94
95
96
97
98
|
/* File: randsubs.c
* Author: Jean Thierry-Mieg (mieg@mrc-lmb.cam.ac.uk)
* Copyright (C) J Thierry-Mieg and R Durbin, 1992
*-------------------------------------------------------------------
* This file is part of the ACEDB genome database package, written by
* Richard Durbin (MRC LMB, UK) rd@mrc-lmb.cam.ac.uk, and
* Jean Thierry-Mieg (CRBM du CNRS, France) mieg@kaa.cnrs-mop.fr
*
* Description:
* Exported functions:
* HISTORY:
* Last edited: Dec 4 15:06 1998 (fw)
* Created: Mon Jun 15 14:44:56 1992 (mieg)
*-------------------------------------------------------------------
*/
/* $Id: randsubs.c,v 1.1 2002/11/14 20:00:06 lstein Exp $ */
#ifdef ALPHA
int random(void); /* in libc.a */
#endif /* ALPHA */
static int xrand = 18721 ;
static int yrand = 37264 ; /* original value 67571 */
static int zrand = 28737 ;
/*******************************/
double randfloat (void)
{double x ;
xrand = 171*xrand % 30269;
yrand = 172*yrand % 30307;
zrand = 170*zrand % 30323;
x = xrand/30269.0 + yrand/30307.0 + zrand/30323.0;
return (x-(int)x);
}
/*******************************/
#ifdef ALPHA
int randint (void)
{
return random() ;
}
#else
int randint (void)
{xrand = 171*xrand % 30269;
yrand = 172*yrand % 30307;
zrand = 170*zrand % 30323;
return (zrand) ;
}
#endif
/*******************************/
double randgauss (void)
{double sum ;
static double fac = 3.0/90899.0 ;
xrand = 171*xrand % 30269;
yrand = 172*yrand % 30307;
zrand = 170*zrand % 30323;
sum = xrand + yrand + zrand ;
xrand = 171*xrand % 30269;
yrand = 172*yrand % 30307;
zrand = 170*zrand % 30323;
sum += xrand + yrand + zrand ;
xrand = 171*xrand % 30269;
yrand = 172*yrand % 30307;
zrand = 170*zrand % 30323;
sum += xrand + yrand + zrand ;
xrand = 171*xrand % 30269;
yrand = 172*yrand % 30307;
zrand = 170*zrand % 30323;
sum += xrand + yrand + zrand ;
return (sum*fac - 6.0) ;
}
/*******************************/
void randsave (int *arr)
{arr[0] = xrand ;
arr[1] = yrand ;
arr[2] = zrand ;
}
/*********************************/
void randrestore (int *arr)
{xrand = arr[0] ;
yrand = arr[1] ;
zrand = arr[2] ;
}
|
