time-to-botec

Benchmark sampling in different programming languages
Log | Files | Refs | README

main.js (6804B)


      1 /**
      2 * @license Apache-2.0
      3 *
      4 * Copyright (c) 2018 The Stdlib Authors.
      5 *
      6 * Licensed under the Apache License, Version 2.0 (the "License");
      7 * you may not use this file except in compliance with the License.
      8 * You may obtain a copy of the License at
      9 *
     10 *    http://www.apache.org/licenses/LICENSE-2.0
     11 *
     12 * Unless required by applicable law or agreed to in writing, software
     13 * distributed under the License is distributed on an "AS IS" BASIS,
     14 * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
     15 * See the License for the specific language governing permissions and
     16 * limitations under the License.
     17 */
     18 
     19 'use strict';
     20 
     21 // MODULES //
     22 
     23 var setReadOnly = require( '@stdlib/utils/define-nonenumerable-read-only-property' );
     24 var setReadOnlyAccessor = require( '@stdlib/utils/define-nonenumerable-read-only-accessor' );
     25 var setReadWriteAccessor = require( '@stdlib/utils/define-nonenumerable-read-write-accessor' );
     26 var constantFunction = require( '@stdlib/utils/constant-function' );
     27 var noop = require( '@stdlib/utils/noop' );
     28 var copy = require( '@stdlib/utils/copy' );
     29 var isPositiveNumber = require( '@stdlib/assert/is-positive-number' ).isPrimitive;
     30 var isObject = require( '@stdlib/assert/is-plain-object' );
     31 var isNonNegativeInteger = require( '@stdlib/assert/is-nonnegative-integer' ).isPrimitive;
     32 var isPositiveInteger = require( '@stdlib/assert/is-positive-integer' ).isPrimitive;
     33 var hasOwnProp = require( '@stdlib/assert/has-own-property' );
     34 var MAX_VALUE = require( '@stdlib/constants/float64/max' );
     35 var rerlang = require( './../../../base/erlang' ).factory;
     36 var iteratorSymbol = require( '@stdlib/symbol/iterator' );
     37 
     38 
     39 // MAIN //
     40 
     41 /**
     42 * Returns an iterator for generating pseudorandom numbers drawn from an Erlang distribution.
     43 *
     44 * @param {PositiveInteger} k - shape parameter
     45 * @param {PositiveNumber} lambda  - rate parameter
     46 * @param {Options} [options] - function options
     47 * @param {PRNG} [options.prng] - pseudorandom number generator which generates uniformly distributed pseudorandom numbers
     48 * @param {PRNGSeedMT19937} [options.seed] - pseudorandom number generator seed
     49 * @param {PRNGStateMT19937} [options.state] - pseudorandom number generator state
     50 * @param {boolean} [options.copy=true] - boolean indicating whether to copy a provided pseudorandom number generator state
     51 * @param {NonNegativeInteger} [options.iter] - number of iterations
     52 * @throws {TypeError} `k` must be a positive integer
     53 * @throws {TypeError} `lambda` must be a positive number
     54 * @throws {TypeError} options argument must be an object
     55 * @throws {TypeError} must provide valid options
     56 * @throws {Error} must provide a valid state
     57 * @returns {Iterator} iterator
     58 *
     59 * @example
     60 * var iter = iterator( 2, 5.0 );
     61 *
     62 * var r = iter.next().value;
     63 * // returns <number>
     64 *
     65 * r = iter.next().value;
     66 * // returns <number>
     67 *
     68 * r = iter.next().value;
     69 * // returns <number>
     70 *
     71 * // ...
     72 */
     73 function iterator( k, lambda, options ) {
     74 	var opts;
     75 	var iter;
     76 	var rand;
     77 	var FLG;
     78 	var i;
     79 	if ( !isPositiveInteger( k ) ) {
     80 		throw new TypeError( 'invalid argument. First argument must be a positive integer. Value: `'+k+'`.' );
     81 	}
     82 	if ( !isPositiveNumber( lambda ) ) {
     83 		throw new TypeError( 'invalid argument. Second argument must be a positive number. Value: `'+lambda+'`.' );
     84 	}
     85 	if ( arguments.length > 2 ) {
     86 		if ( !isObject( options ) ) {
     87 			throw new TypeError( 'invalid argument. Options argument must be an object. Value: `'+options+'`.' );
     88 		}
     89 		opts = copy( options, 1 );
     90 		if ( hasOwnProp( opts, 'iter' ) ) {
     91 			if ( !isNonNegativeInteger( opts.iter ) ) {
     92 				throw new TypeError( 'invalid option. `iter` option must be a nonnegative integer. Option: `'+opts.iter+'`.' );
     93 			}
     94 		} else {
     95 			opts.iter = MAX_VALUE;
     96 		}
     97 		rand = rerlang( k, lambda, opts );
     98 		if ( opts.prng === void 0 && opts.copy !== false ) {
     99 			opts.state = rand.state; // cache a copy of the PRNG state
    100 		}
    101 	} else {
    102 		rand = rerlang( k, lambda );
    103 		opts = {
    104 			'iter': MAX_VALUE,
    105 			'state': rand.state // cache a copy of the PRNG state
    106 		};
    107 	}
    108 	i = 0;
    109 
    110 	// Create an iterator protocol-compliant object:
    111 	iter = {};
    112 	setReadOnly( iter, 'next', next );
    113 	setReadOnly( iter, 'return', end );
    114 
    115 	if ( opts && opts.prng ) {
    116 		setReadOnly( iter, 'seed', null );
    117 		setReadOnly( iter, 'seedLength', null );
    118 		setReadWriteAccessor( iter, 'state', constantFunction( null ), noop );
    119 		setReadOnly( iter, 'stateLength', null );
    120 		setReadOnly( iter, 'byteLength', null );
    121 	} else {
    122 		setReadOnlyAccessor( iter, 'seed', getSeed );
    123 		setReadOnlyAccessor( iter, 'seedLength', getSeedLength );
    124 		setReadWriteAccessor( iter, 'state', getState, setState );
    125 		setReadOnlyAccessor( iter, 'stateLength', getStateLength );
    126 		setReadOnlyAccessor( iter, 'byteLength', getStateSize );
    127 	}
    128 	setReadOnly( iter, 'PRNG', rand.PRNG );
    129 
    130 	// If an environment supports `Symbol.iterator`, make the iterator iterable:
    131 	if ( iteratorSymbol ) {
    132 		setReadOnly( iter, iteratorSymbol, factory );
    133 	}
    134 	return iter;
    135 
    136 	/**
    137 	* Returns an iterator protocol-compliant object containing the next iterated value.
    138 	*
    139 	* @private
    140 	* @returns {Object} iterator protocol-compliant object
    141 	*/
    142 	function next() {
    143 		i += 1;
    144 		if ( FLG || i > opts.iter ) {
    145 			return {
    146 				'done': true
    147 			};
    148 		}
    149 		return {
    150 			'value': rand(),
    151 			'done': false
    152 		};
    153 	}
    154 
    155 	/**
    156 	* Finishes an iterator.
    157 	*
    158 	* @private
    159 	* @param {*} [value] - value to return
    160 	* @returns {Object} iterator protocol-compliant object
    161 	*/
    162 	function end( value ) {
    163 		FLG = true;
    164 		if ( arguments.length ) {
    165 			return {
    166 				'value': value,
    167 				'done': true
    168 			};
    169 		}
    170 		return {
    171 			'done': true
    172 		};
    173 	}
    174 
    175 	/**
    176 	* Returns a new iterator.
    177 	*
    178 	* @private
    179 	* @returns {Iterator} iterator
    180 	*/
    181 	function factory() {
    182 		return iterator( k, lambda, opts );
    183 	}
    184 
    185 	/**
    186 	* Returns the PRNG seed.
    187 	*
    188 	* @private
    189 	* @returns {PRNGSeedMT19937} seed
    190 	*/
    191 	function getSeed() {
    192 		return rand.PRNG.seed;
    193 	}
    194 
    195 	/**
    196 	* Returns the PRNG seed length.
    197 	*
    198 	* @private
    199 	* @returns {PositiveInteger} seed length
    200 	*/
    201 	function getSeedLength() {
    202 		return rand.PRNG.seedLength;
    203 	}
    204 
    205 	/**
    206 	* Returns the PRNG state length.
    207 	*
    208 	* @private
    209 	* @returns {PositiveInteger} state length
    210 	*/
    211 	function getStateLength() {
    212 		return rand.PRNG.stateLength;
    213 	}
    214 
    215 	/**
    216 	* Returns the PRNG state size (in bytes).
    217 	*
    218 	* @private
    219 	* @returns {PositiveInteger} state size (in bytes)
    220 	*/
    221 	function getStateSize() {
    222 		return rand.PRNG.byteLength;
    223 	}
    224 
    225 	/**
    226 	* Returns the current pseudorandom number generator state.
    227 	*
    228 	* @private
    229 	* @returns {PRNGStateMT19937} current state
    230 	*/
    231 	function getState() {
    232 		return rand.PRNG.state;
    233 	}
    234 
    235 	/**
    236 	* Sets the pseudorandom number generator state.
    237 	*
    238 	* @private
    239 	* @param {PRNGStateMT19937} s - generator state
    240 	* @throws {Error} must provide a valid state
    241 	*/
    242 	function setState( s ) {
    243 		rand.PRNG.state = s;
    244 	}
    245 }
    246 
    247 
    248 // EXPORTS //
    249 
    250 module.exports = iterator;