// this program is a JavaScript version of Mersenne Twister, with concealment and encapsulation in class, // an almost straight conversion from the original program, mt19937ar.c, // translated by y. okada on July 17, 2006. // Changes by Jesse Ruderman: added "var" keyword in a few spots; added export_mta etc; pasted into fuzz.js. // in this program, procedure descriptions and comments of original source code were not removed. // lines commented with //c// were originally descriptions of c procedure. and a few following lines are appropriate JavaScript descriptions. // lines commented with /* and */ are original comments. // lines commented with // are additional comments in this JavaScript version. // before using this version, create at least one instance of MersenneTwister19937 class, and initialize the each state, given below in c comments, of all the instances. /* AC-programforMT19937,withinitializationimproved2002/1/26. CodedbyTakujiNishimuraandMakotoMatsumoto.
function MersenneTwister19937()
{ /* Period parameters */ //c//#define N 624 //c//#define M 397 //c//#define MATRIX_A 0x9908b0dfUL /* constant vector a */ //c//#define UPPER_MASK 0x80000000UL /* most significant w-r bits */ //c//#define LOWER_MASK 0x7fffffffUL /* least significant r bits */ var N = 624; var M = 397; var MATRIX_A = 0x9908b0df; /* constant vector a */ var UPPER_MASK = 0x80000000; /* most significant w-r bits */ var LOWER_MASK = 0x7fffffff; /* least significant r bits */ //c//static unsigned long mt[N]; /* the array for the state vector */ //c//static int mti=N+1; /* mti==N+1 means mt[N] is not initialized */ var mt = new Array(N); /* the array for the state vector */ var mti = N+1; /* mti==N+1 means mt[N] is not initialized */
function unsigned32 (n1) // returns a 32-bits unsiged integer from an operand to which applied a bit operator.
{ return n1 < 0 ? (n1 ^ UPPER_MASK) + UPPER_MASK : n1;
}
function subtraction32 (n1, n2) // emulates lowerflow of a c 32-bits unsiged integer variable, instead of the operator -. these both arguments must be non-negative integers expressible using unsigned 32 bits.
{ return n1 < n2 ? unsigned32((0x100000000 - (n2 - n1)) & 0xffffffff) : n1 - n2;
}
function addition32 (n1, n2) // emulates overflow of a c 32-bits unsiged integer variable, instead of the operator +. these both arguments must be non-negative integers expressible using unsigned 32 bits.
{ return unsigned32((n1 + n2) & 0xffffffff)
}
function multiplication32 (n1, n2) // emulates overflow of a c 32-bits unsiged integer variable, instead of the operator *. these both arguments must be non-negative integers expressible using unsigned 32 bits.
{ var sum = 0; for (var i = 0; i < 32; ++i){ if ((n1 >>> i) & 0x1){
sum = addition32(sum, unsigned32(n2 << i));
}
} return sum;
}
/* initializes mt[N] with a seed */ //c//void init_genrand(unsigned long s) this.init_genrand = function (s)
{ //c//mt[0]= s & 0xffffffff;
mt[0]= unsigned32(s & 0xffffffff); for (mti=1; mti<N; mti++) {
mt[mti] = //c//(1812433253 * (mt[mti-1] ^ (mt[mti-1] >> 30)) + mti);
addition32(multiplication32(1812433253, unsigned32(mt[mti-1] ^ (mt[mti-1] >>> 30))), mti); /* See Knuth TAOCP Vol2. 3rd Ed. P.106 for multiplier. */ /* In the previous versions, MSBs of the seed affect */ /* only MSBs of the array mt[]. */ /* 2002/01/09 modified by Makoto Matsumoto */ //c//mt[mti] &= 0xffffffff;
mt[mti] = unsigned32(mt[mti] & 0xffffffff); /* for >32 bit machines */
}
}
/* initialize by an array with array-length */ /* init_key is the array for initializing keys */ /* key_length is its length */ /* slight change for C++, 2004/2/26 */ //c//void init_by_array(unsigned long init_key[], int key_length) this.init_by_array = function (init_key, key_length)
{ //c//int i, j, k; var i, j, k; //c//init_genrand(19650218); this.init_genrand(19650218);
i=1; j=0;
k = (N>key_length ? N : key_length); for (; k; k--) { //c//mt[i] = (mt[i] ^ ((mt[i-1] ^ (mt[i-1] >> 30)) * 1664525)) //c// + init_key[j] + j; /* non linear */
mt[i] = addition32(addition32(unsigned32(mt[i] ^ multiplication32(unsigned32(mt[i-1] ^ (mt[i-1] >>> 30)), 1664525)), init_key[j]), j);
mt[i] = //c//mt[i] &= 0xffffffff; /* for WORDSIZE > 32 machines */
unsigned32(mt[i] & 0xffffffff);
i++; j++; if (i>=N) { mt[0] = mt[N-1]; i=1; } if (j>=key_length) j=0;
} for (k=N-1; k; k--) { //c//mt[i] = (mt[i] ^ ((mt[i-1] ^ (mt[i-1] >> 30)) * 1566083941)) //c//- i; /* non linear */
mt[i] = subtraction32(unsigned32((dbg=mt[i]) ^ multiplication32(unsigned32(mt[i-1] ^ (mt[i-1] >>> 30)), 1566083941)), i); //c//mt[i] &= 0xffffffff; /* for WORDSIZE > 32 machines */
mt[i] = unsigned32(mt[i] & 0xffffffff);
i++; if (i>=N) { mt[0] = mt[N-1]; i=1; }
}
mt[0] = 0x80000000; /* MSB is 1; assuring non-zero initial array */
}
/* generates a random number on [0,0xffffffff]-interval */ //c//unsigned long genrand_int32(void) this.genrand_int32 = function ()
{ //c//unsigned long y; //c//static unsigned long mag01[2]={0x0UL, MATRIX_A}; var y; var mag01 = new Array(0x0, MATRIX_A); /* mag01[x] = x * MATRIX_A for x=0,1 */
if (mti >= N) { /* generate N words at one time */ //c//int kk; var kk;
if (mti == N+1) /* if init_genrand() has not been called, */ //c//init_genrand(5489); /* a default initial seed is used */ this.init_genrand(5489); /* a default initial seed is used */
/* generates a random number on [0,0x7fffffff]-interval */ //c//long genrand_int31(void) this.genrand_int31 = function ()
{ //c//return (genrand_int32()>>1); return (this.genrand_int32()>>>1);
}
/* generates a random number on [0,1]-real-interval */ //c//double genrand_real1(void) this.genrand_real1 = function ()
{ //c//return genrand_int32()*(1.0/4294967295.0); returnthis.genrand_int32()*(1.0/4294967295.0); /* divided by 2^32-1 */
}
/* generates a random number on [0,1)-real-interval */ //c//double genrand_real2(void) this.genrand_real2 = function ()
{ //c//return genrand_int32()*(1.0/4294967296.0); returnthis.genrand_int32()*(1.0/4294967296.0); /* divided by 2^32 */
}
/* generates a random number on (0,1)-real-interval */ //c//double genrand_real3(void) this.genrand_real3 = function ()
{ //c//return ((genrand_int32()) + 0.5)*(1.0/4294967296.0); return ((this.genrand_int32()) + 0.5)*(1.0/4294967296.0); /* divided by 2^32 */
}
/* generates a random number on [0,1) with 53-bit resolution*/ //c//double genrand_res53(void) this.genrand_res53 = function ()
{ //c//unsigned long a=genrand_int32()>>5, b=genrand_int32()>>6; var a=this.genrand_int32()>>>5, b=this.genrand_int32()>>>6; return(a*67108864.0+b)*(1.0/9007199254740992.0);
} /* These real versions are due to Isaku Wada, 2002/01/09 added */
}
function initRnd() { var fuzzMT = new MersenneTwister19937; var fuzzSeed = 53;
fuzzMT.init_genrand(fuzzSeed);
rnd = function (n) { var v = Math.floor(fuzzMT.genrand_real2() * n); return v; };
rnd.rndReal = function() { return fuzzMT.genrand_real2(); };
rnd.fuzzMT = fuzzMT;
}
Die Informationen auf dieser Webseite wurden
nach bestem Wissen sorgfältig zusammengestellt. Es wird jedoch weder Vollständigkeit, noch Richtigkeit,
noch Qualität der bereit gestellten Informationen zugesichert.
Bemerkung:
Die farbliche Syntaxdarstellung und die Messung sind noch experimentell.