random.c revision 1.1 1 1.1 cgd /*
2 1.1 cgd * Copyright (c) 1988 The Regents of the University of California.
3 1.1 cgd * All rights reserved.
4 1.1 cgd *
5 1.1 cgd * This code is derived from software contributed to Berkeley by
6 1.1 cgd * Timothy C. Stoehr.
7 1.1 cgd *
8 1.1 cgd * Redistribution and use in source and binary forms, with or without
9 1.1 cgd * modification, are permitted provided that the following conditions
10 1.1 cgd * are met:
11 1.1 cgd * 1. Redistributions of source code must retain the above copyright
12 1.1 cgd * notice, this list of conditions and the following disclaimer.
13 1.1 cgd * 2. Redistributions in binary form must reproduce the above copyright
14 1.1 cgd * notice, this list of conditions and the following disclaimer in the
15 1.1 cgd * documentation and/or other materials provided with the distribution.
16 1.1 cgd * 3. All advertising materials mentioning features or use of this software
17 1.1 cgd * must display the following acknowledgement:
18 1.1 cgd * This product includes software developed by the University of
19 1.1 cgd * California, Berkeley and its contributors.
20 1.1 cgd * 4. Neither the name of the University nor the names of its contributors
21 1.1 cgd * may be used to endorse or promote products derived from this software
22 1.1 cgd * without specific prior written permission.
23 1.1 cgd *
24 1.1 cgd * THIS SOFTWARE IS PROVIDED BY THE REGENTS AND CONTRIBUTORS ``AS IS'' AND
25 1.1 cgd * ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
26 1.1 cgd * IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
27 1.1 cgd * ARE DISCLAIMED. IN NO EVENT SHALL THE REGENTS OR CONTRIBUTORS BE LIABLE
28 1.1 cgd * FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL
29 1.1 cgd * DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS
30 1.1 cgd * OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION)
31 1.1 cgd * HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT
32 1.1 cgd * LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY
33 1.1 cgd * OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF
34 1.1 cgd * SUCH DAMAGE.
35 1.1 cgd */
36 1.1 cgd
37 1.1 cgd #ifndef lint
38 1.1 cgd static char sccsid[] = "@(#)random.c 5.3 (Berkeley) 6/1/90";
39 1.1 cgd #endif /* not lint */
40 1.1 cgd
41 1.1 cgd /*
42 1.1 cgd * random.c
43 1.1 cgd *
44 1.1 cgd * This source herein may be modified and/or distributed by anybody who
45 1.1 cgd * so desires, with the following restrictions:
46 1.1 cgd * 1.) No portion of this notice shall be removed.
47 1.1 cgd * 2.) Credit shall not be taken for the creation of this source.
48 1.1 cgd * 3.) This code is not to be traded, sold, or used for personal
49 1.1 cgd * gain or profit.
50 1.1 cgd *
51 1.1 cgd */
52 1.1 cgd
53 1.1 cgd static long rntb[32] = {
54 1.1 cgd 3, 0x9a319039, 0x32d9c024, 0x9b663182, 0x5da1f342,
55 1.1 cgd 0xde3b81e0, 0xdf0a6fb5, 0xf103bc02, 0x48f340fb, 0x7449e56b,
56 1.1 cgd 0xbeb1dbb0, 0xab5c5918, 0x946554fd, 0x8c2e680f, 0xeb3d799f,
57 1.1 cgd 0xb11ee0b7, 0x2d436b86, 0xda672e2a, 0x1588ca88, 0xe369735d,
58 1.1 cgd 0x904f35f7, 0xd7158fd6, 0x6fa6f051, 0x616e6b96, 0xac94efdc,
59 1.1 cgd 0x36413f93, 0xc622c298, 0xf5a42ab8, 0x8a88d77b, 0xf5ad9d0e,
60 1.1 cgd 0x8999220b, 0x27fb47b9
61 1.1 cgd };
62 1.1 cgd
63 1.1 cgd static long *fptr = &rntb[4];
64 1.1 cgd static long *rptr = &rntb[1];
65 1.1 cgd static long *state = &rntb[1];
66 1.1 cgd static int rand_type = 3;
67 1.1 cgd static int rand_deg = 31;
68 1.1 cgd static int rand_sep = 3;
69 1.1 cgd static long *end_ptr = &rntb[32];
70 1.1 cgd
71 1.1 cgd srrandom(x)
72 1.1 cgd int x;
73 1.1 cgd {
74 1.1 cgd register int i;
75 1.1 cgd long rrandom();
76 1.1 cgd
77 1.1 cgd state[0] = (long) x;
78 1.1 cgd if (rand_type != 0) {
79 1.1 cgd for (i = 1; i < rand_deg; i++) {
80 1.1 cgd state[i] = 1103515245 * state[i - 1] + 12345;
81 1.1 cgd }
82 1.1 cgd fptr = &state[rand_sep];
83 1.1 cgd rptr = &state[0];
84 1.1 cgd for (i = 0; i < 10 * rand_deg; i++) {
85 1.1 cgd (void) rrandom();
86 1.1 cgd }
87 1.1 cgd }
88 1.1 cgd }
89 1.1 cgd
90 1.1 cgd long
91 1.1 cgd rrandom()
92 1.1 cgd {
93 1.1 cgd long i;
94 1.1 cgd
95 1.1 cgd if (rand_type == 0) {
96 1.1 cgd i = state[0] = (state[0]*1103515245 + 12345) & 0x7fffffff;
97 1.1 cgd } else {
98 1.1 cgd *fptr += *rptr;
99 1.1 cgd i = (*fptr >> 1) & 0x7fffffff;
100 1.1 cgd if (++fptr >= end_ptr) {
101 1.1 cgd fptr = state;
102 1.1 cgd ++rptr;
103 1.1 cgd } else {
104 1.1 cgd if (++rptr >= end_ptr) {
105 1.1 cgd rptr = state;
106 1.1 cgd }
107 1.1 cgd }
108 1.1 cgd }
109 1.1 cgd return(i);
110 1.1 cgd }
111 1.1 cgd
112 1.1 cgd get_rand(x, y)
113 1.1 cgd register int x, y;
114 1.1 cgd {
115 1.1 cgd register int r, t;
116 1.1 cgd long lr;
117 1.1 cgd
118 1.1 cgd if (x > y) {
119 1.1 cgd t = y;
120 1.1 cgd y = x;
121 1.1 cgd x = t;
122 1.1 cgd }
123 1.1 cgd lr = rrandom();
124 1.1 cgd lr &= (long) 0x00003fff;
125 1.1 cgd r = (int) lr;
126 1.1 cgd r = (r % ((y - x) + 1)) + x;
127 1.1 cgd return(r);
128 1.1 cgd }
129 1.1 cgd
130 1.1 cgd rand_percent(percentage)
131 1.1 cgd register int percentage;
132 1.1 cgd {
133 1.1 cgd return(get_rand(1, 100) <= percentage);
134 1.1 cgd }
135 1.1 cgd
136 1.1 cgd coin_toss()
137 1.1 cgd {
138 1.1 cgd
139 1.1 cgd return(((rrandom() & 01) ? 1 : 0));
140 1.1 cgd }
141