Home | History | Annotate | Line # | Download | only in rogue
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