🎧 van der Corput low-discrepancy sequences
Go to file
2016-09-24 14:54:13 -04:00
bin version bump 0.2.0: cleanup 2016-01-16 13:47:42 -05:00
misc version bump 0.2.1: cleanup 2016-09-24 14:54:13 -04:00
test_files version bump 0.2.0: cleanup 2016-01-16 13:47:42 -05:00
.flowconfig version bump 0.2.0: cleanup 2016-01-16 13:47:42 -05:00
.gitignore version bump 0.2.1: cleanup 2016-09-24 14:54:13 -04:00
.jscs.json version bump 0.2.0: cleanup 2016-01-16 13:47:42 -05:00
.jshintrc version bump 0.2.0: cleanup 2016-01-16 13:47:42 -05:00
.npmignore version bump 0.2.0: cleanup 2016-01-16 13:47:42 -05:00
.travis.yml version bump 0.2.0: cleanup 2016-01-16 13:47:42 -05:00
index.html version bump 0.2.1: cleanup 2016-09-24 14:54:13 -04:00
LICENSE version bump 0.2.0: cleanup 2016-01-16 13:47:42 -05:00
Makefile version bump 0.2.1: cleanup 2016-09-24 14:54:13 -04:00
package.json version bump 0.2.1: cleanup 2016-09-24 14:54:13 -04:00
README.md version bump 0.2.1: cleanup 2016-09-24 14:54:13 -04:00
test.js version bump 0.2.0: cleanup 2016-01-16 13:47:42 -05:00
vdc.flow.js version bump 0.2.0: cleanup 2016-01-16 13:47:42 -05:00
vdc.js version bump 0.2.0: cleanup 2016-01-16 13:47:42 -05:00

VDC

Pure JS implementation of van der Corput low-discrepancy sequences.

Installation

With npm:

$ npm install vdc

In the browser:

<script src="vdc.js"></script>

The browser exposes a variable VDC

Usage

The exported function VDC accepts a opts object with the following fields:

  • b (default 2) base for the sequence.
  • n (default 0) starting index for the sequence.

Calling without arguments will default to the aforementioned values.

The object returned by VDC exposes a next() method to get the next element.

The field last holds the most recently generated value (accessing the field does not trigger a recalculation)

For example:

//var VDC = require('vdc') // uncomment this line if in node

var opts = {'n':0, 'b':2};
var generator = VDC(opts);
for(var i = 0; i != 10; ++i) console.log(generator.next());

The expected output is

0       (0/1)
0.5     (1/2)
0.25    (1/4)
0.75    (3/4)
0.125   (1/8)
0.625   (5/8)
0.375   (3/8)
0.875   (7/8)
0.0625  (1/16)
0.5625  (9/16)

Testing

make test will run the nodejs-based test.

make baseline will generate the test baselines using Mathematica by explicitly extracting and reversing the digits. The implementation is based off a tutorial http://reference.wolfram.com/mathematica/CUDALink/tutorial/Applications.html

VanDerCorput[base_][len_] := Table[
  With[{digits = Reverse@IntegerDigits[n, base]},
   Sum[2^(-ii)*digits[[ii]], {ii, Length[digits]}]
   ], {n, len}]

Notes

0 is the first value. Some sources (notably Wikipedia) start the sequence at 1/2, but others (Glasserman "Monte Carlo Methods in Financial Engineering") claim that the original sequence definition started at 0.

License

Please consult the attached LICENSE file for details. All rights not explicitly granted by the Apache 2.0 license are reserved by the Original Author.

Badges

Build Status

Coverage Status

Analytics