3 * $Id: fibrand.c,v 1.1 2000/05/21 11:28:30 mdw Exp $
7 * (c) 1999 Straylight/Edgeware
8 * (c) 2000 Mark Wooding
11 /*----- Licensing notice --------------------------------------------------*
13 * Copyright (c) 2000 Mark Wooding
14 * All rights reserved.
16 * Redistribution and use in source and binary forms, with or without
17 * modification, are permitted provided that the following conditions are
20 * 1. Redistributions of source code must retain the above copyright
21 * notice, this list of conditions and the following disclaimer.
23 * 2, Redistributions in binary form must reproduce the above copyright
24 * notice, this list of conditions and the following disclaimer in the
25 * documentation and/or other materials provided with the distribution.
27 * 3. The name of the authors may not be used to endorse or promote
28 * products derived from this software without specific prior written
31 * THIS SOFTWARE IS PROVIDED ``AS IS'' AND ANY EXPRESS OR IMPLIED
32 * WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED WARRANTIES OF
33 * MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE ARE DISCLAIMED. IN
34 * NO EVENT SHALL THE REGENTS OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT,
35 * INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES
36 * (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR
37 * SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION)
38 * HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT,
39 * STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN
40 * ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
41 * POSSIBILITY OF SUCH DAMAGE.
43 * Instead of accepting the above terms, you may redistribute and/or modify
44 * this software under the terms of either the GNU General Public License,
45 * or the GNU Library General Public License, published by the Free
46 * Software Foundation; either version 2 of the License, or (at your
47 * option) any later version.
50 /*----- Revision history --------------------------------------------------*
53 * Revision 1.1 2000/05/21 11:28:30 mdw
56 * --- Past lives (Catacomb) --- *
58 * Revision 1.1 1999/12/10 23:15:27 mdw
59 * Noncryptographic random number generator.
63 /*----- Header files ------------------------------------------------------*/
73 /*----- Main code ---------------------------------------------------------*/
75 /* --- @fibrand_step@ --- *
77 * Arguments: @fibrand *f@ = pointer to Fibonacci generator context
79 * Returns: Next output from generator.
81 * Use: Steps the generator. Returns
82 * %$x_{i - 24} + x_{i - 55} \bmod 2^{32}$%.
85 uint32
fibrand_step(fibrand
*f
)
88 unsigned j
= i
+ (FIB_SZ
- FIB_TAP
);
92 x
= f
->x
[i
] = U32(f
->x
[i
] + f
->x
[j
]);
100 /* --- @fibrand_lcseed@ --- *
102 * Arguments: @fibrand *f@ = pointer to Fibonacci generator context
103 * @uint32 seed@ = seed value
107 * Use: Initializes a Fibonacci generator using outputs from the
108 * @lcrand@ generator seeded from @seed@. This is faster than
109 * using a generic @lcrand@-based generator and @fibrand_rseed@
110 * because it uses raw outputs rather than uniformly distributed
114 void fibrand_lcseed(fibrand
*f
, uint32 seed
)
119 for (i
= 0; i
< FIB_SZ
; i
++)
120 p
|= f
->x
[i
] = seed
= lcrand(seed
);
122 i
= lcrand_range(&seed
, FIB_SZ
);
128 /* --- @fibrand_range@ --- *
130 * Arguments: @fibrand *f@ = pointer to Fibonacci generator context
133 * Returns: A uniformly distributed pseudorandom integer in the interval
137 uint32
fibrand_range(fibrand
*f
, uint32 m
)
139 uint32 r
= 0xffffffff - (0xffffffff % m
);
142 /* --- Now generate numbers until a good one comes along --- */
144 do x
= fibrand_step(f
); while (x
>= r
);
145 return (x
/ (r
/ m
));
148 /*----- That's all, folks -------------------------------------------------*/