gx
chenyc
2025-02-12 ea42ff3ebee1eeb3fb29423aa848a249441db81c
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
// A simple smoke test and benchmark for the generators.
 
var assert = require('assert');
var xor128 = require('../lib/xor128');
var xorwow = require('../lib/xorwow');
var xs7 = require('../lib/xorshift7');
var xor4096 = require('../lib/xor4096');
var tychei = require('../lib/tychei');
var alea = require('../lib/alea');
var sr = require('../seedrandom');
 
describe("XOR-Shift generator test", function() {
 
var benchmarks = { native: { rand: Math.random, times: [] } };
 
function test(label, alg, double1, float3, int4, hc, qc, ec, e2c) {
  var fn = alg(1);
  var fn2 = alg('hello.', { state: true });
  benchmarks[label] = {rand: fn.quick, times: []};
  it("should use " + label + " correctly", function() {
    if (double1 != null) assert.equal(fn.double(), double1);
    if (float3 != null) assert.equal(fn.quick(), float3);
    if (int4 != null) assert.equal(fn.int32(), int4);
    assert(fn() > 0);
    assert(fn() < 1);
    assert(fn2() > 0);
    // Internal state is visible only if requested.
    assert(!('state' in fn));
    assert('state' in fn2);
    var ss = fn2.state();
    var rs = fn2();
    assert(rs < 1);
    var j, h = 0, q = 0, e = 0, r, p, e2 = 0;
    for (j = 0; j < 1024; ++j) {
      r = fn();
      if (r < 0.5) h += 1;
      if (r < 0.25) q += 1;
      if (r < 0.125) e += 1;
      r2 = fn2();
      if (r2 < 0.125) e2 += 1;
    }
    if (hc != null) {
      assert.equal(h, hc);
      assert.equal(q, qc);
      assert.equal(e, ec);
      assert.equal(e2, e2c);
      h = q = e = p = 0;
      for (j = 0; j < 1024; ++j) {
        r = fn.double();
        if (r < 0.5) h += 1;
        if (r < 0.25) q += 1;
        if (r < 0.125) e += 1;
        if (fn.int32() >= 0) p += 1;
      }
      // Sanity-check double() and int32.
      assert(h >= 480 && h <= 543, h);
      assert(q >= 226 && q <= 286, q);
      assert(e >= 100 && e <= 156, e);
      assert(e2 >= 100 && e2 <= 156, e2);
      assert(p >= 482 && p <= 543, p);
    }
    var fn3 = alg(0, { state: ss });
    assert.equal(fn3(), rs);
  });
}
 
test("xor128", xor128,
    0.7963797148975774, 0.22171171731315553, 317177041, 498, 236, 110, 115);
test("xorwow", xorwow,
    0.8178000247146859, 0.8407576507888734, 533150816, 519, 228, 121, 123);
test("xorshift7", xs7,
    0.21241471533241418, 0.9957620368804783, -1678071207, 510, 261, 143, 124);
test("tychei", tychei,
    0.42331440041340196, 0.9365617581643164, -884984569, 521, 242, 116, 126);
test("seedrandom", sr,
    0.1776348083296759, 0.2160690303426236, 1397712774, 526, 282, 131, 137);
test("xor4096", xor4096,
    0.1520436450538547, 0.4206166828516871, 1312695376, 496, 241, 113, 142);
test("alea", alea,
    0.5260470956849501, 0.47771977609954774, -1625913352, 494, 246, 125, 122);
 
it("runs benchmarks", function() {
  var n = 4;
  var trials = 10;
  var top = 4;
  this.timeout(200 * n * trials);
  this.slow(30 * n * trials);
  var fn, k, start, end, j, t;
  for (k in benchmarks) {
    fn = benchmarks[k].rand;
    // warmup.
    for (j = 0; j < 1e5; ++j) fn();
  }
  for (t = 0; t < trials; ++t) {
    for (k in benchmarks) {
      fn = benchmarks[k].rand;
      start = +new Date;
      // benchmark.
      for (j = 0; j < n * 1e5; ++j) fn();
      end = +new Date;
      benchmarks[k].times.push(end - start);
    }
  }
  for (k in benchmarks) {
    benchmarks[k].times.sort();
  }
  function fastest(array) {
    var sum = 0;
    for (var j = 0; j < top; ++j) {
       sum += array[j];
    }
    return sum / top;
  }
  var nativetime = fastest(benchmarks.native.times);
  for (k in benchmarks) {
    var time = fastest(benchmarks[k].times);
    console.log(k+ ': ' + time / n + ' nanoseconds per call, ' +
       (time / nativetime).toFixed(1) + 'x native random.');
  }
});
 
});