Robust/src/Benchmarks/SingleTM/Bayes/Random.java

   1 public class Random {
   2   long[] mt;
   3   int mti;
   4   long RANDOM_DEFAULT_SEED;
   5   /* period parameter */
   6   int N;
   7   int M;
   8   long MATRIX_A;
   9   long UPPER_MASK;
  10   long LOWER_MASK;
  11
  12   public Random() {
  13     RANDOM_DEFAULT_SEED = 0L;
  14     N = 624;
  15     M = 397;
  16     mt = new long[N];
  17     mti = N;
  18     MATRIX_A = 0x9908b0dfL;   /* constant vector a */
  19     UPPER_MASK = 0x80000000L; /* most significant w-r bits */
  20     LOWER_MASK = 0x7fffffffL; /* least significant r bits */
  21   }
  22
  23   public void random_alloc() {
  24     init_genrand(this.RANDOM_DEFAULT_SEED);
  25   }
  26
  27   /* initializes mt[N] with a seed */
  28   public void init_genrand(long s) {
  29     int mti;
  30     mt[0]= s & 0xFFFFFFFFL;
  31     for (mti=1; mti<N; mti++) {
  32      mt[mti] = (1812433253L * (mt[mti-1] ^ (mt[mti-1] >> 30)) + mti);
  33       /* See Knuth TAOCP Vol2. 3rd Ed. P.106 for multiplier. */
  34       /* In the previous versions, MSBs of the seed affect   */
  35       /* only MSBs of the array mt[].                        */
  36       /* 2002/01/09 modified by Makoto Matsumoto             */
  37       mt[mti] &= 0xFFFFFFFFL;
  38       /* for >32 bit machines */
  39     }
  40     this.mti=mti;
  41   }
  42
  43   public void random_seed(long seed) {
  44     init_genrand(seed);
  45   }
  46
  47   public long random_generate() {
  48     return genrand_int32();
  49   }
  50
  51   //public static long genrand_int32(long[] mt, long mtiPtr) {
  52   public long genrand_int32() {
  53     long y;
  54     long[] mag01= new long[2];
  55     mag01[0] = 0x0L;
  56     mag01[1] = MATRIX_A;
  57     int mti = this.mti;
  58
  59     /* mag01[x] = x * MATRIX_A  for x=0,1 */
  60
  61     if (mti >= N) { /* generate N words at one time */
  62       int kk;
  63
  64       if (mti == N+1)   /* if init_genrand() has not been called, */
  65         init_genrand(5489L); /* a default initial seed is used */
  66
  67       for (kk=0;kk<N-M;kk++) {
  68         y = (mt[kk]&UPPER_MASK)|(mt[kk+1]&LOWER_MASK);
  69         mt[kk] = mt[kk+M] ^ (y >> 1) ^ mag01[(int)(y & 0x1L)];
  70       }
  71       for (;kk<N-1;kk++) {
  72         y = (mt[kk]&UPPER_MASK)|(mt[kk+1]&LOWER_MASK);
  73         mt[kk] = mt[kk+(M-N)] ^ (y >> 1) ^ mag01[(int)(y & 0x1L)];
  74       }
  75       y = (mt[N-1]&UPPER_MASK)|(mt[0]&LOWER_MASK);
  76       mt[N-1] = mt[M-1] ^ (y >> 1) ^ mag01[(int)(y & 0x1L)];
  77
  78       mti = 0;
  79     }
  80
  81     y = mt[mti++];
  82
  83     /* Tempering */
  84     y ^= (y >> 11);
  85     y ^= (y << 7) & 0x9d2c5680L;
  86     y ^= (y << 15) & 0xefc60000L;
  87     y ^= (y >> 18);
  88
  89     this.mti = mti;
  90
  91     return y;
  92   }
  93 }