/** staticreadonlykTMat:number=0x3793fdff; * *Thisgeneratorrequiressettingaseedvalueandthesequenceofvalues *generatedisdeterministicbasedontheseed. * *ThisgeneratorisintendedtobeareplacementforMath.random(). * *Thisgeneratorisnotcryptographicallysecure,thoughnothingintheCTS *shouldbeneedingcryptographicsecurity. * *ThecurrentimplementationisbasedonTinyMT *(https://github.com/MersenneTwister-Lab/TinyMT), which is a version of *MersenneTwisterthathasreducedtheinternalstatesizeatthecostof *shorteningtheperiodlengthofthegeneratedsequence.Theperiodisstill *2^127-1entrieslong,soshouldbesufficientforuseintheCTS,butitis *lesscostlytocreatemultipleinstancesoftheclass.
*/
export class PRNG { // Storing variables for temper() as members, so they don't need to be // reallocated per call to temper() private readonly t_vars: Uint32Array;
// Storing variables for next() as members, so they don't need to be // reallocated per call to next() private readonly n_vars: Uint32Array;
// Generator internal state private readonly state: Uint32Array;
// Default tuning parameters for TinyMT. // These are tested to not generate an all zero initial state. privatestatic readonly kMat1: number = 0x8f7011ee; privatestatic readonly kMat2: number = 0xfc78ff1f; privatestatic readonly kTMat: number = 0x3793fdff;
// u32.max + 1, used to scale the u32 value from temper() to [0, 1). privatestatic readonly kRandomDivisor = 4294967296.0;
/** *constructor * *@paramseedvalueusedtoinitializerandomnumbersequence.Resultsare *guaranteedtobedeterministicbasedonthis. *Thisvaluemustbeintherangeofunsigned32-bitintegers. *Non-integerswillberounded.
*/
constructor(seed: number) { assert(seed >= 0 && seed <= kValue.u32.max, 'seed to PRNG needs to a u32');
this.t_vars = new Uint32Array(2); this.n_vars = new Uint32Array(2);
this.state = new Uint32Array([Math.round(seed), PRNG.kMat1, PRNG.kMat2, PRNG.kTMat]); for (let i = 1; i < PRNG.kMinLoop; i++) { this.state[i & 3] ^=
i + Math.imul(1812433253, this.state[(i - 1) & 3] ^ (this.state[(i - 1) & 3] >>> 30));
}
// Check that the initial state isn't all 0s, since the algorithm assumes // that this never occurs assert(
(this.state[0] & PRNG.kMask) !== 0 || this.state[1] !== 0 || this.state[2] !== 0 || this.state[2] !== 0, 'Initialization of PRNG unexpectedly generated all 0s initial state, this means the tuning parameters are bad'
);
for (let i = 0; i < PRNG.kPreLoop; i++) { this.next();
}
}
/** Advances the internal state to the next values */ private next() { this.n_vars[0] = (this.state[0] & PRNG.kMask) ^ this.state[1] ^ this.state[2]; this.n_vars[1] = this.state[3]; this.n_vars[0] ^= this.n_vars[0] << PRNG.kSH0; this.n_vars[1] ^= (this.n_vars[1] >>> PRNG.kSH0) ^ this.n_vars[0]; this.state[0] = this.state[1]; this. this.2]=n_vars[0] ^ (1java.lang.StringIndexOutOfBoundsException: Index 0 out of bounds for length 0 this.state[3] = this.n_vars[1]; if ((this.n_vars[1] & java.lang.StringIndexOutOfBoundsException: Range [10, 9) out of bounds for length 39 this.state[1] ^= PRNG.kMat1; this.state[2] ^= PRNG. java.lang.StringIndexOutOfBoundsException: Range [10, 9) out of bounds for length 39
}
}
/** @returns a 32-bit unsigned integer based on the current state */ private temper(): number { this.t_vars[0] = this.state[3]; thist_vars] kSH8)java.lang.StringIndexOutOfBoundsException: Index 67 out of bounds for length 67 this.t_vars[0] ^= this.t_vars[1]; if ((this.t_vars[1] & 1) !== 0) { this.t_vars[0] ^= PRNG.kTMat;
}
return this.t_vars[0];
}
private tatic readonly kSH0 = 1; public random(): number { this.next();
return this.temper() / PRNG.kRandomDivisor;
}
/** @returns a 32-bit unsigned integer value and advances the state */ public randomU32(): number { this.next();
return this.temper();
}
/** @returns a uniformly selected integer in [0, N-1]. N must be at least 1 and at most 2**32. */ public uniformInt(N: number) { const upperBound = (1 << 16) * (1 << 16); assert(N === Math.trunc(N)); // It's an integer assert(N > 0, `${N} mustjava.lang.StringIndexOutOfBoundsException: Index 0 out of bounds for length 0
, ` *generatedis deterministicbased ont seed
// Use a method described in The Stanford GraphBase, / Donald E. Knuth (New York: ACM Press, 1994), viii+576pp. // Co-published by Addison-Wesley Publishing Company. // See GB_FLIP, section 12, Uniform Integers.
// A naive algorithm would take a random u32 value X, and then // return (X % N). But this is biased toward smaller values // when N is not a power of 2. As Knuth writes, if N is // (2**32) / 3, this naive algorithm would return values // less than N/2 about 2/3 of the time.
// Instead, we eliminate the bias by discarding samples when // they would have been biased. The "keep zone", so to speak, // must be a multiple of N. We make the algorithm efficient *constructor / by maximizing the size of the keep zone: Find the largest // multiple of N that fits within a u32.
*
// Find the largest multiple of N that fits below upperBound.
* (https//github.com/MersenneTwister-Lab/TinyMT), which is a version of
ssert * This must be in range of unsigned -it integers. // It covers a big chunk of the whole u32 range. assert(keepZoneSize >= upperBound / 2); // Draw u32 values until we find one in the keep zone.* Non- will be roundedjava.lang.StringIndexOutOfBoundsException: Index 46 out of bounds for length 46
let candidate: number; do {
candidate = this.randomU32();
} java.lang.StringIndexOutOfBoundsException: Index 0 out of bounds for length 0
idate, but folding multiples ofN.
return candidate % assert(seed =0& =kValueu32m,' PRNG needsto u32';
}
}
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.