1*1da177e4SLinus Torvalds /* 2*1da177e4SLinus Torvalds * Generic address resultion entity 3*1da177e4SLinus Torvalds * 4*1da177e4SLinus Torvalds * Authors: 5*1da177e4SLinus Torvalds * net_random Alan Cox 6*1da177e4SLinus Torvalds * net_ratelimit Andy Kleen 7*1da177e4SLinus Torvalds * 8*1da177e4SLinus Torvalds * Created by Alexey Kuznetsov <[email protected]> 9*1da177e4SLinus Torvalds * 10*1da177e4SLinus Torvalds * This program is free software; you can redistribute it and/or 11*1da177e4SLinus Torvalds * modify it under the terms of the GNU General Public License 12*1da177e4SLinus Torvalds * as published by the Free Software Foundation; either version 13*1da177e4SLinus Torvalds * 2 of the License, or (at your option) any later version. 14*1da177e4SLinus Torvalds */ 15*1da177e4SLinus Torvalds 16*1da177e4SLinus Torvalds #include <linux/module.h> 17*1da177e4SLinus Torvalds #include <linux/jiffies.h> 18*1da177e4SLinus Torvalds #include <linux/kernel.h> 19*1da177e4SLinus Torvalds #include <linux/mm.h> 20*1da177e4SLinus Torvalds #include <linux/string.h> 21*1da177e4SLinus Torvalds #include <linux/types.h> 22*1da177e4SLinus Torvalds #include <linux/random.h> 23*1da177e4SLinus Torvalds #include <linux/percpu.h> 24*1da177e4SLinus Torvalds #include <linux/init.h> 25*1da177e4SLinus Torvalds 26*1da177e4SLinus Torvalds #include <asm/system.h> 27*1da177e4SLinus Torvalds #include <asm/uaccess.h> 28*1da177e4SLinus Torvalds 29*1da177e4SLinus Torvalds 30*1da177e4SLinus Torvalds /* 31*1da177e4SLinus Torvalds This is a maximally equidistributed combined Tausworthe generator 32*1da177e4SLinus Torvalds based on code from GNU Scientific Library 1.5 (30 Jun 2004) 33*1da177e4SLinus Torvalds 34*1da177e4SLinus Torvalds x_n = (s1_n ^ s2_n ^ s3_n) 35*1da177e4SLinus Torvalds 36*1da177e4SLinus Torvalds s1_{n+1} = (((s1_n & 4294967294) <<12) ^ (((s1_n <<13) ^ s1_n) >>19)) 37*1da177e4SLinus Torvalds s2_{n+1} = (((s2_n & 4294967288) << 4) ^ (((s2_n << 2) ^ s2_n) >>25)) 38*1da177e4SLinus Torvalds s3_{n+1} = (((s3_n & 4294967280) <<17) ^ (((s3_n << 3) ^ s3_n) >>11)) 39*1da177e4SLinus Torvalds 40*1da177e4SLinus Torvalds The period of this generator is about 2^88. 41*1da177e4SLinus Torvalds 42*1da177e4SLinus Torvalds From: P. L'Ecuyer, "Maximally Equidistributed Combined Tausworthe 43*1da177e4SLinus Torvalds Generators", Mathematics of Computation, 65, 213 (1996), 203--213. 44*1da177e4SLinus Torvalds 45*1da177e4SLinus Torvalds This is available on the net from L'Ecuyer's home page, 46*1da177e4SLinus Torvalds 47*1da177e4SLinus Torvalds http://www.iro.umontreal.ca/~lecuyer/myftp/papers/tausme.ps 48*1da177e4SLinus Torvalds ftp://ftp.iro.umontreal.ca/pub/simulation/lecuyer/papers/tausme.ps 49*1da177e4SLinus Torvalds 50*1da177e4SLinus Torvalds There is an erratum in the paper "Tables of Maximally 51*1da177e4SLinus Torvalds Equidistributed Combined LFSR Generators", Mathematics of 52*1da177e4SLinus Torvalds Computation, 68, 225 (1999), 261--269: 53*1da177e4SLinus Torvalds http://www.iro.umontreal.ca/~lecuyer/myftp/papers/tausme2.ps 54*1da177e4SLinus Torvalds 55*1da177e4SLinus Torvalds ... the k_j most significant bits of z_j must be non- 56*1da177e4SLinus Torvalds zero, for each j. (Note: this restriction also applies to the 57*1da177e4SLinus Torvalds computer code given in [4], but was mistakenly not mentioned in 58*1da177e4SLinus Torvalds that paper.) 59*1da177e4SLinus Torvalds 60*1da177e4SLinus Torvalds This affects the seeding procedure by imposing the requirement 61*1da177e4SLinus Torvalds s1 > 1, s2 > 7, s3 > 15. 62*1da177e4SLinus Torvalds 63*1da177e4SLinus Torvalds */ 64*1da177e4SLinus Torvalds struct nrnd_state { 65*1da177e4SLinus Torvalds u32 s1, s2, s3; 66*1da177e4SLinus Torvalds }; 67*1da177e4SLinus Torvalds 68*1da177e4SLinus Torvalds static DEFINE_PER_CPU(struct nrnd_state, net_rand_state); 69*1da177e4SLinus Torvalds 70*1da177e4SLinus Torvalds static u32 __net_random(struct nrnd_state *state) 71*1da177e4SLinus Torvalds { 72*1da177e4SLinus Torvalds #define TAUSWORTHE(s,a,b,c,d) ((s&c)<<d) ^ (((s <<a) ^ s)>>b) 73*1da177e4SLinus Torvalds 74*1da177e4SLinus Torvalds state->s1 = TAUSWORTHE(state->s1, 13, 19, 4294967294UL, 12); 75*1da177e4SLinus Torvalds state->s2 = TAUSWORTHE(state->s2, 2, 25, 4294967288UL, 4); 76*1da177e4SLinus Torvalds state->s3 = TAUSWORTHE(state->s3, 3, 11, 4294967280UL, 17); 77*1da177e4SLinus Torvalds 78*1da177e4SLinus Torvalds return (state->s1 ^ state->s2 ^ state->s3); 79*1da177e4SLinus Torvalds } 80*1da177e4SLinus Torvalds 81*1da177e4SLinus Torvalds static void __net_srandom(struct nrnd_state *state, unsigned long s) 82*1da177e4SLinus Torvalds { 83*1da177e4SLinus Torvalds if (s == 0) 84*1da177e4SLinus Torvalds s = 1; /* default seed is 1 */ 85*1da177e4SLinus Torvalds 86*1da177e4SLinus Torvalds #define LCG(n) (69069 * n) 87*1da177e4SLinus Torvalds state->s1 = LCG(s); 88*1da177e4SLinus Torvalds state->s2 = LCG(state->s1); 89*1da177e4SLinus Torvalds state->s3 = LCG(state->s2); 90*1da177e4SLinus Torvalds 91*1da177e4SLinus Torvalds /* "warm it up" */ 92*1da177e4SLinus Torvalds __net_random(state); 93*1da177e4SLinus Torvalds __net_random(state); 94*1da177e4SLinus Torvalds __net_random(state); 95*1da177e4SLinus Torvalds __net_random(state); 96*1da177e4SLinus Torvalds __net_random(state); 97*1da177e4SLinus Torvalds __net_random(state); 98*1da177e4SLinus Torvalds } 99*1da177e4SLinus Torvalds 100*1da177e4SLinus Torvalds 101*1da177e4SLinus Torvalds unsigned long net_random(void) 102*1da177e4SLinus Torvalds { 103*1da177e4SLinus Torvalds unsigned long r; 104*1da177e4SLinus Torvalds struct nrnd_state *state = &get_cpu_var(net_rand_state); 105*1da177e4SLinus Torvalds r = __net_random(state); 106*1da177e4SLinus Torvalds put_cpu_var(state); 107*1da177e4SLinus Torvalds return r; 108*1da177e4SLinus Torvalds } 109*1da177e4SLinus Torvalds 110*1da177e4SLinus Torvalds 111*1da177e4SLinus Torvalds void net_srandom(unsigned long entropy) 112*1da177e4SLinus Torvalds { 113*1da177e4SLinus Torvalds struct nrnd_state *state = &get_cpu_var(net_rand_state); 114*1da177e4SLinus Torvalds __net_srandom(state, state->s1^entropy); 115*1da177e4SLinus Torvalds put_cpu_var(state); 116*1da177e4SLinus Torvalds } 117*1da177e4SLinus Torvalds 118*1da177e4SLinus Torvalds void __init net_random_init(void) 119*1da177e4SLinus Torvalds { 120*1da177e4SLinus Torvalds int i; 121*1da177e4SLinus Torvalds 122*1da177e4SLinus Torvalds for (i = 0; i < NR_CPUS; i++) { 123*1da177e4SLinus Torvalds struct nrnd_state *state = &per_cpu(net_rand_state,i); 124*1da177e4SLinus Torvalds __net_srandom(state, i+jiffies); 125*1da177e4SLinus Torvalds } 126*1da177e4SLinus Torvalds } 127*1da177e4SLinus Torvalds 128*1da177e4SLinus Torvalds static int net_random_reseed(void) 129*1da177e4SLinus Torvalds { 130*1da177e4SLinus Torvalds int i; 131*1da177e4SLinus Torvalds unsigned long seed[NR_CPUS]; 132*1da177e4SLinus Torvalds 133*1da177e4SLinus Torvalds get_random_bytes(seed, sizeof(seed)); 134*1da177e4SLinus Torvalds for (i = 0; i < NR_CPUS; i++) { 135*1da177e4SLinus Torvalds struct nrnd_state *state = &per_cpu(net_rand_state,i); 136*1da177e4SLinus Torvalds __net_srandom(state, seed[i]); 137*1da177e4SLinus Torvalds } 138*1da177e4SLinus Torvalds return 0; 139*1da177e4SLinus Torvalds } 140*1da177e4SLinus Torvalds late_initcall(net_random_reseed); 141*1da177e4SLinus Torvalds 142*1da177e4SLinus Torvalds int net_msg_cost = 5*HZ; 143*1da177e4SLinus Torvalds int net_msg_burst = 10; 144*1da177e4SLinus Torvalds 145*1da177e4SLinus Torvalds /* 146*1da177e4SLinus Torvalds * All net warning printk()s should be guarded by this function. 147*1da177e4SLinus Torvalds */ 148*1da177e4SLinus Torvalds int net_ratelimit(void) 149*1da177e4SLinus Torvalds { 150*1da177e4SLinus Torvalds return __printk_ratelimit(net_msg_cost, net_msg_burst); 151*1da177e4SLinus Torvalds } 152*1da177e4SLinus Torvalds 153*1da177e4SLinus Torvalds EXPORT_SYMBOL(net_random); 154*1da177e4SLinus Torvalds EXPORT_SYMBOL(net_ratelimit); 155*1da177e4SLinus Torvalds EXPORT_SYMBOL(net_srandom); 156