main.js (6457B)
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 isPositive = 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 hasOwnProp = require( '@stdlib/assert/has-own-property' ); 33 var MAX_VALUE = require( '@stdlib/constants/float64/max' ); 34 var rrayleigh = require( './../../../base/rayleigh' ).factory; 35 var iteratorSymbol = require( '@stdlib/symbol/iterator' ); 36 37 38 // MAIN // 39 40 /** 41 * Returns an iterator for generating pseudorandom numbers drawn from a Rayleigh distribution. 42 * 43 * @param {PositiveNumber} sigma - scale parameter 44 * @param {Options} [options] - function options 45 * @param {PRNG} [options.prng] - pseudorandom number generator which generates uniformly distributed pseudorandom numbers 46 * @param {PRNGSeedMT19937} [options.seed] - pseudorandom number generator seed 47 * @param {PRNGStateMT19937} [options.state] - pseudorandom number generator state 48 * @param {boolean} [options.copy=true] - boolean indicating whether to copy a provided pseudorandom number generator state 49 * @param {NonNegativeInteger} [options.iter] - number of iterations 50 * @throws {TypeError} first argument must be a positive number 51 * @throws {TypeError} options argument must be an object 52 * @throws {TypeError} must provide valid options 53 * @throws {Error} must provide a valid state 54 * @returns {Iterator} iterator 55 * 56 * @example 57 * var iter = iterator( 2.0 ); 58 * 59 * var r = iter.next().value; 60 * // returns <number> 61 * 62 * r = iter.next().value; 63 * // returns <number> 64 * 65 * r = iter.next().value; 66 * // returns <number> 67 * 68 * // ... 69 */ 70 function iterator( sigma, options ) { 71 var opts; 72 var iter; 73 var rand; 74 var FLG; 75 var i; 76 if ( !isPositive( sigma ) ) { 77 throw new TypeError( 'invalid argument. First argument must be a positive number. Value: `'+sigma+'`.' ); 78 } 79 if ( arguments.length > 1 ) { 80 if ( !isObject( options ) ) { 81 throw new TypeError( 'invalid argument. Options argument must be an object. Value: `'+options+'`.' ); 82 } 83 opts = copy( options, 1 ); 84 if ( hasOwnProp( opts, 'iter' ) ) { 85 if ( !isNonNegativeInteger( opts.iter ) ) { 86 throw new TypeError( 'invalid option. `iter` option must be a nonnegative integer. Option: `'+opts.iter+'`.' ); 87 } 88 } else { 89 opts.iter = MAX_VALUE; 90 } 91 rand = rrayleigh( sigma, opts ); 92 if ( opts.prng === void 0 && opts.copy !== false ) { 93 opts.state = rand.state; // cache a copy of the PRNG state 94 } 95 } else { 96 rand = rrayleigh( sigma ); 97 opts = { 98 'iter': MAX_VALUE, 99 'state': rand.state // cache a copy of the PRNG state 100 }; 101 } 102 i = 0; 103 104 // Create an iterator protocol-compliant object: 105 iter = {}; 106 setReadOnly( iter, 'next', next ); 107 setReadOnly( iter, 'return', end ); 108 109 if ( opts && opts.prng ) { 110 setReadOnly( iter, 'seed', null ); 111 setReadOnly( iter, 'seedLength', null ); 112 setReadWriteAccessor( iter, 'state', constantFunction( null ), noop ); 113 setReadOnly( iter, 'stateLength', null ); 114 setReadOnly( iter, 'byteLength', null ); 115 } else { 116 setReadOnlyAccessor( iter, 'seed', getSeed ); 117 setReadOnlyAccessor( iter, 'seedLength', getSeedLength ); 118 setReadWriteAccessor( iter, 'state', getState, setState ); 119 setReadOnlyAccessor( iter, 'stateLength', getStateLength ); 120 setReadOnlyAccessor( iter, 'byteLength', getStateSize ); 121 } 122 setReadOnly( iter, 'PRNG', rand.PRNG ); 123 124 // If an environment supports `Symbol.iterator`, make the iterator iterable: 125 if ( iteratorSymbol ) { 126 setReadOnly( iter, iteratorSymbol, factory ); 127 } 128 return iter; 129 130 /** 131 * Returns an iterator protocol-compliant object containing the next iterated value. 132 * 133 * @private 134 * @returns {Object} iterator protocol-compliant object 135 */ 136 function next() { 137 i += 1; 138 if ( FLG || i > opts.iter ) { 139 return { 140 'done': true 141 }; 142 } 143 return { 144 'value': rand(), 145 'done': false 146 }; 147 } 148 149 /** 150 * Finishes an iterator. 151 * 152 * @private 153 * @param {*} [value] - value to return 154 * @returns {Object} iterator protocol-compliant object 155 */ 156 function end( value ) { 157 FLG = true; 158 if ( arguments.length ) { 159 return { 160 'value': value, 161 'done': true 162 }; 163 } 164 return { 165 'done': true 166 }; 167 } 168 169 /** 170 * Returns a new iterator. 171 * 172 * @private 173 * @returns {Iterator} iterator 174 */ 175 function factory() { 176 return iterator( sigma, opts ); 177 } 178 179 /** 180 * Returns the PRNG seed. 181 * 182 * @private 183 * @returns {PRNGSeedMT19937} seed 184 */ 185 function getSeed() { 186 return rand.PRNG.seed; 187 } 188 189 /** 190 * Returns the PRNG seed length. 191 * 192 * @private 193 * @returns {PositiveInteger} seed length 194 */ 195 function getSeedLength() { 196 return rand.PRNG.seedLength; 197 } 198 199 /** 200 * Returns the PRNG state length. 201 * 202 * @private 203 * @returns {PositiveInteger} state length 204 */ 205 function getStateLength() { 206 return rand.PRNG.stateLength; 207 } 208 209 /** 210 * Returns the PRNG state size (in bytes). 211 * 212 * @private 213 * @returns {PositiveInteger} state size (in bytes) 214 */ 215 function getStateSize() { 216 return rand.PRNG.byteLength; 217 } 218 219 /** 220 * Returns the current pseudorandom number generator state. 221 * 222 * @private 223 * @returns {PRNGStateMT19937} current state 224 */ 225 function getState() { 226 return rand.PRNG.state; 227 } 228 229 /** 230 * Sets the pseudorandom number generator state. 231 * 232 * @private 233 * @param {PRNGStateMT19937} s - generator state 234 * @throws {Error} must provide a valid state 235 */ 236 function setState( s ) { 237 rand.PRNG.state = s; 238 } 239 } 240 241 242 // EXPORTS // 243 244 module.exports = iterator;