lustre/lvfs/prng.c

   1 /* -*- mode: c; c-basic-offset: 8; indent-tabs-mode: nil; -*-
   2  * vim:expandtab:shiftwidth=8:tabstop=8:
   3  *
   4  *   This file is part of the Lustre file system, http://www.lustre.org
   5  *   Lustre is a trademark of Cluster File Systems, Inc.
   6  *
   7  * concatenation of following two 16-bit multiply with carry generators
   8  * x(n)=a*x(n-1)+carry mod 2^16 and y(n)=b*y(n-1)+carry mod 2^16,
   9  * number and carry packed within the same 32 bit integer.
  10  * algorithm recommended by Marsaglia
  11  ******************************************************************/
  12 #ifndef EXPORT_SYMTAB
  13 # define EXPORT_SYMTAB
  14 #endif
  15
  16 #ifndef __KERNEL__
  17 #include <liblustre.h>
  18 #define get_random_bytes(val, size)     (*val) = 0
  19 #endif
  20 #include <obd_class.h>
  21 #include <linux/random.h>
  22
  23 /*
  24 From: George Marsaglia <geo@stat.fsu.edu>
  25 Newsgroups: sci.math
  26 Subject: Re: A RANDOM NUMBER GENERATOR FOR C
  27 Date: Tue, 30 Sep 1997 05:29:35 -0700
  28
  29  * You may replace the two constants 36969 and 18000 by any
  30  * pair of distinct constants from this list:
  31  * 18000 18030 18273 18513 18879 19074 19098 19164 19215 19584
  32  * 19599 19950 20088 20508 20544 20664 20814 20970 21153 21243
  33  * 21423 21723 21954 22125 22188 22293 22860 22938 22965 22974
  34  * 23109 23124 23163 23208 23508 23520 23553 23658 23865 24114
  35  * 24219 24660 24699 24864 24948 25023 25308 25443 26004 26088
  36  * 26154 26550 26679 26838 27183 27258 27753 27795 27810 27834
  37  * 27960 28320 28380 28689 28710 28794 28854 28959 28980 29013
  38  * 29379 29889 30135 30345 30459 30714 30903 30963 31059 31083
  39  * (or any other 16-bit constants k for which both k*2^16-1
  40  * and k*2^15-1 are prime) */
  41
  42 #define RANDOM_CONST_A 18030
  43 #define RANDOM_CONST_B 29013
  44
  45 static unsigned int seed_x = 521288629;
  46 static unsigned int seed_y = 362436069;
  47 unsigned int ll_rand(void)
  48 {
  49
  50         seed_x = RANDOM_CONST_A * (seed_x & 65535) + (seed_x >> 16);
  51         seed_y = RANDOM_CONST_B * (seed_y & 65535) + (seed_y >> 16);
  52
  53         return ((seed_x << 16) + (seed_y & 65535));
  54 }
  55 EXPORT_SYMBOL(ll_rand);
  56
  57 /* Note that if the input seeds are not completely random, then there is
  58  * a preferred location for the entropy in the two seeds, in order to avoid
  59  * the initial values from the PRNG to be the same each time.
  60  *
  61  * seed1 (seed_x) should have the most entropy in the low bits of the word
  62  * seed2 (seed_y) should have the most entropy in the high bits of the word */
  63 void ll_srand(unsigned int seed1, unsigned int seed2)
  64 {
  65         if (seed1)
  66                 seed_x = seed1; /* use default seeds if parameter is 0 */
  67         if (seed2)
  68                 seed_y = seed2;
  69 }
  70 EXPORT_SYMBOL(ll_srand);
  71
  72 void ll_get_random_bytes(void *buf, int size)
  73 {
  74         int *p = buf;
  75         int rem, tmp;
  76
  77         LASSERT(size >= 0);
  78
  79         rem = min((int)((unsigned long)buf & (sizeof(int) - 1)), size);
  80         if (rem) {
  81                 get_random_bytes(&tmp, sizeof(tmp));
  82                 tmp ^= ll_rand();
  83                 memcpy(buf, &tmp, rem);
  84                 p = buf + rem;
  85                 size -= rem;
  86         }
  87
  88         while (size >= sizeof(int)) {
  89                 get_random_bytes(&tmp, sizeof(tmp));
  90                 *p = ll_rand() ^ tmp;
  91                 size -= sizeof(int);
  92                 p++;
  93         }
  94         buf = p;
  95         if (size) {
  96                 get_random_bytes(&tmp, sizeof(tmp));
  97                 tmp ^= ll_rand();
  98                 memcpy(buf, &tmp, size);
  99         }
 100 }
 101 EXPORT_SYMBOL(ll_get_random_bytes);
 102
 103 void ll_generate_random_uuid(class_uuid_t uuid_out)
 104 {
 105         ll_get_random_bytes(uuid_out, sizeof(class_uuid_t));
 106 }
 107 EXPORT_SYMBOL(ll_generate_random_uuid);
 108