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;