Robust/src/Benchmarks/SingleTM/common/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 long genrand_int32() {
  52     long y;
  53     long[] mag01= new long[2];
  54     mag01[0] = 0x0L;
  55     mag01[1] = MATRIX_A;
  56     int mti = this.mti;
  57
  58     /* mag01[x] = x * MATRIX_A  for x=0,1 */
  59
  60     if (mti >= N) { /* generate N words at one time */
  61       int kk;
  62
  63       if (mti == N+1)   /* if init_genrand() has not been called, */
  64         init_genrand(5489L); /* a default initial seed is used */
  65
  66       for (kk=0;kk<N-M;kk++) {
  67         y = (mt[kk]&UPPER_MASK)|(mt[kk+1]&LOWER_MASK);
  68         mt[kk] = mt[kk+M] ^ (y >> 1) ^ mag01[(int)(y & 0x1L)];
  69       }
  70       for (;kk<N-1;kk++) {
  71         y = (mt[kk]&UPPER_MASK)|(mt[kk+1]&LOWER_MASK);
  72         mt[kk] = mt[kk+(M-N)] ^ (y >> 1) ^ mag01[(int)(y & 0x1L)];
  73       }
  74       y = (mt[N-1]&UPPER_MASK)|(mt[0]&LOWER_MASK);
  75       mt[N-1] = mt[M-1] ^ (y >> 1) ^ mag01[(int)(y & 0x1L)];
  76
  77       mti = 0;
  78     }
  79
  80     y = mt[mti++];
  81
  82     /* Tempering */
  83     y ^= (y >> 11);
  84     y ^= (y << 7) & 0x9d2c5680L;
  85     y ^= (y << 15) & 0xefc60000L;
  86     y ^= (y >> 18);
  87
  88     this.mti = mti;
  89
  90     return y;
  91   }
  92 }