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