diff --git a/lib/node_modules/@stdlib/stats/base/dists/gamma/cdf/benchmark/benchmark.js b/lib/node_modules/@stdlib/stats/base/dists/gamma/cdf/benchmark/benchmark.js index 507d2a6c7b33..98476110a77b 100644 --- a/lib/node_modules/@stdlib/stats/base/dists/gamma/cdf/benchmark/benchmark.js +++ b/lib/node_modules/@stdlib/stats/base/dists/gamma/cdf/benchmark/benchmark.js @@ -21,7 +21,8 @@ // MODULES // var bench = require( '@stdlib/bench' ); -var randu = require( '@stdlib/random/base/randu' ); +var Float64Array = require( '@stdlib/array/float64' ); +var uniform = require( '@stdlib/random/base/uniform' ); var isnan = require( '@stdlib/math/base/assert/is-nan' ); var EPS = require( '@stdlib/constants/float64/eps' ); var pkg = require( './../package.json' ).name; @@ -33,16 +34,24 @@ var cdf = require( './../lib' ); bench( pkg, function benchmark( b ) { var alpha; var beta; + var len; var x; var y; var i; + len = 100; + x = new Float64Array( len ); + alpha = new Float64Array( len ); + beta = new Float64Array( len ); + for ( i = 0; i < len; i++ ) { + x[ i ] = uniform( EPS, 100.0 ); + alpha[ i ] = uniform( EPS, 100.0 ); + beta[ i ] = uniform( EPS, 100.0 ); + } + b.tic(); for ( i = 0; i < b.iterations; i++ ) { - x = ( randu()*100.0 ) + EPS; - alpha = ( randu()*100.0 ) + EPS; - beta = ( randu()*100.0 ) + EPS; - y = cdf( x, alpha, beta ); + y = cdf( x[ i % len ], alpha[ i % len ], beta[ i % len ] ); if ( isnan( y ) ) { b.fail( 'should not return NaN' ); } @@ -59,6 +68,7 @@ bench( pkg+':factory', function benchmark( b ) { var mycdf; var alpha; var beta; + var len; var x; var y; var i; @@ -66,11 +76,15 @@ bench( pkg+':factory', function benchmark( b ) { alpha = 20.0; beta = 15.0; mycdf = cdf.factory( alpha, beta ); + len = 100; + x = new Float64Array( len ); + for ( i = 0; i < len; i++ ) { + x[ i ] = uniform( EPS, 50.0 ); + } b.tic(); for ( i = 0; i < b.iterations; i++ ) { - x = ( randu()*50.0 ) + EPS; - y = mycdf( x ); + y = mycdf( x[ i % len ] ); if ( isnan( y ) ) { b.fail( 'should not return NaN' ); } diff --git a/lib/node_modules/@stdlib/stats/base/dists/gamma/ctor/benchmark/benchmark.js b/lib/node_modules/@stdlib/stats/base/dists/gamma/ctor/benchmark/benchmark.js index 85ffd08542a2..3ee50ab3ecff 100644 --- a/lib/node_modules/@stdlib/stats/base/dists/gamma/ctor/benchmark/benchmark.js +++ b/lib/node_modules/@stdlib/stats/base/dists/gamma/ctor/benchmark/benchmark.js @@ -21,7 +21,8 @@ // MODULES // var bench = require( '@stdlib/bench' ); -var randu = require( '@stdlib/random/base/randu' ); +var Float64Array = require( '@stdlib/array/float64' ); +var uniform = require( '@stdlib/random/base/uniform' ); var isnan = require( '@stdlib/math/base/assert/is-nan' ); var EPS = require( '@stdlib/constants/float64/eps' ); var pkg = require( './../package.json' ).name; @@ -34,13 +35,20 @@ bench( pkg+'::instantiation', function benchmark( b ) { var alpha; var beta; var dist; + var len; var i; + len = 100; + alpha = new Float64Array( len ); + beta = new Float64Array( len ); + for ( i = 0; i < len; i++ ) { + alpha[ i ] = uniform( EPS, 10.0 ); + beta[ i ] = uniform( EPS, 10.0 ); + } + b.tic(); for ( i = 0; i < b.iterations; i++ ) { - alpha = ( randu() * 10.0 ) + EPS; - beta = ( randu() * 10.0 ) + EPS; - dist = new Gamma( alpha, beta ); + dist = new Gamma( alpha[ i % len ], beta[ i % len ] ); if ( !( dist instanceof Gamma ) ) { b.fail( 'should return a distribution instance' ); } @@ -83,18 +91,23 @@ bench( pkg+'::set:alpha', function benchmark( b ) { var alpha; var beta; var dist; + var len; var y; var i; alpha = 100.56789; beta = 55.54321; dist = new Gamma( alpha, beta ); + len = 100; + y = new Float64Array( len ); + for ( i = 0; i < len; i++ ) { + y[ i ] = uniform( EPS, 100.0 ); + } b.tic(); for ( i = 0; i < b.iterations; i++ ) { - y = ( 100.0*randu() ) + EPS; - dist.alpha = y; - if ( dist.alpha !== y ) { + dist.alpha = y[ i % len ]; + if ( dist.alpha !== y[ i % len ] ) { b.fail( 'should return set value' ); } } @@ -136,18 +149,23 @@ bench( pkg+'::set:beta', function benchmark( b ) { var alpha; var beta; var dist; + var len; var y; var i; alpha = 100.56789; beta = 55.54321; dist = new Gamma( alpha, beta ); + len = 100; + y = new Float64Array( len ); + for ( i = 0; i < len; i++ ) { + y[ i ] = uniform( EPS, 100.0 ); + } b.tic(); for ( i = 0; i < b.iterations; i++ ) { - y = ( 100.0*randu() ) + EPS; - dist.beta = y; - if ( dist.beta !== y ) { + dist.beta = y[ i % len ]; + if ( dist.beta !== y[ i % len ] ) { b.fail( 'should return set value' ); } } @@ -163,16 +181,23 @@ bench( pkg+':entropy', function benchmark( b ) { var alpha; var beta; var dist; + var len; + var x; var y; var i; alpha = 100.56789; beta = 55.54321; dist = new Gamma( alpha, beta ); + len = 100; + x = new Float64Array( len ); + for ( i = 0; i < len; i++ ) { + x[ i ] = uniform( EPS, 100.0 ); + } b.tic(); for ( i = 0; i < b.iterations; i++ ) { - dist.alpha = ( 100.0*randu() ) + EPS; + dist.alpha = x[ i % len ]; y = dist.entropy; if ( isnan( y ) ) { b.fail( 'should not return NaN' ); @@ -190,16 +215,23 @@ bench( pkg+':kurtosis', function benchmark( b ) { var alpha; var beta; var dist; + var len; + var x; var y; var i; alpha = 100.56789; beta = 55.54321; dist = new Gamma( alpha, beta ); + len = 100; + x = new Float64Array( len ); + for ( i = 0; i < len; i++ ) { + x[ i ] = uniform( EPS, 100.0 ); + } b.tic(); for ( i = 0; i < b.iterations; i++ ) { - dist.alpha = ( 100.0*randu() ) + EPS; + dist.alpha = x[ i % len ]; y = dist.kurtosis; if ( isnan( y ) ) { b.fail( 'should not return NaN' ); @@ -217,16 +249,23 @@ bench( pkg+':mean', function benchmark( b ) { var alpha; var beta; var dist; + var len; + var x; var y; var i; alpha = 100.56789; beta = 55.54321; dist = new Gamma( alpha, beta ); + len = 100; + x = new Float64Array( len ); + for ( i = 0; i < len; i++ ) { + x[ i ] = uniform( EPS, 100.0 ); + } b.tic(); for ( i = 0; i < b.iterations; i++ ) { - dist.alpha = ( 100.0*randu() ) + EPS; + dist.alpha = x[ i % len ]; y = dist.mean; if ( isnan( y ) ) { b.fail( 'should not return NaN' ); @@ -244,16 +283,23 @@ bench( pkg+':median', function benchmark( b ) { var alpha; var beta; var dist; + var len; + var x; var y; var i; alpha = 100.56789; beta = 55.54321; dist = new Gamma( alpha, beta ); + len = 100; + x = new Float64Array( len ); + for ( i = 0; i < len; i++ ) { + x[ i ] = uniform( EPS, 100.0 ); + } b.tic(); for ( i = 0; i < b.iterations; i++ ) { - dist.alpha = ( 100.0*randu() ) + EPS; + dist.alpha = x[ i % len ]; y = dist.median; if ( isnan( y ) ) { b.fail( 'should not return NaN' ); @@ -271,16 +317,23 @@ bench( pkg+':mode', function benchmark( b ) { var alpha; var beta; var dist; + var len; + var x; var y; var i; alpha = 100.56789; beta = 55.54321; dist = new Gamma( alpha, beta ); + len = 100; + x = new Float64Array( len ); + for ( i = 0; i < len; i++ ) { + x[ i ] = uniform( 1.0 + EPS, 100.0 ); + } b.tic(); for ( i = 0; i < b.iterations; i++ ) { - dist.alpha = ( 100.0*randu() ) + 1.0 + EPS; + dist.alpha = x[ i % len ]; y = dist.mode; if ( isnan( y ) ) { b.fail( 'should not return NaN' ); @@ -298,16 +351,23 @@ bench( pkg+':skewness', function benchmark( b ) { var alpha; var beta; var dist; + var len; + var x; var y; var i; alpha = 100.56789; beta = 55.54321; dist = new Gamma( alpha, beta ); + len = 100; + x = new Float64Array( len ); + for ( i = 0; i < len; i++ ) { + x[ i ] = uniform( EPS, 100.0 ); + } b.tic(); for ( i = 0; i < b.iterations; i++ ) { - dist.alpha = ( 100.0*randu() ) + EPS; + dist.alpha = x[ i % len ]; y = dist.skewness; if ( isnan( y ) ) { b.fail( 'should not return NaN' ); @@ -325,16 +385,23 @@ bench( pkg+':stdev', function benchmark( b ) { var alpha; var beta; var dist; + var len; + var x; var y; var i; alpha = 100.56789; beta = 55.54321; dist = new Gamma( alpha, beta ); + len = 100; + x = new Float64Array( len ); + for ( i = 0; i < len; i++ ) { + x[ i ] = uniform( EPS, 100.0 ); + } b.tic(); for ( i = 0; i < b.iterations; i++ ) { - dist.alpha = ( 100.0*randu() ) + EPS; + dist.alpha = x[ i % len ]; y = dist.stdev; if ( isnan( y ) ) { b.fail( 'should not return NaN' ); @@ -352,16 +419,23 @@ bench( pkg+':variance', function benchmark( b ) { var alpha; var beta; var dist; + var len; + var x; var y; var i; alpha = 100.56789; beta = 55.54321; dist = new Gamma( alpha, beta ); + len = 100; + x = new Float64Array( len ); + for ( i = 0; i < len; i++ ) { + x[ i ] = uniform( EPS, 100.0 ); + } b.tic(); for ( i = 0; i < b.iterations; i++ ) { - dist.alpha = ( 100.0*randu() ) + EPS; + dist.alpha = x[ i % len ]; y = dist.variance; if ( isnan( y ) ) { b.fail( 'should not return NaN' ); @@ -379,6 +453,7 @@ bench( pkg+':cdf', function benchmark( b ) { var alpha; var beta; var dist; + var len; var x; var y; var i; @@ -386,11 +461,15 @@ bench( pkg+':cdf', function benchmark( b ) { alpha = 100.56789; beta = 55.54321; dist = new Gamma( alpha, beta ); + len = 100; + x = new Float64Array( len ); + for ( i = 0; i < len; i++ ) { + x[ i ] = uniform( 0.0, 1.0 ); + } b.tic(); for ( i = 0; i < b.iterations; i++ ) { - x = randu(); - y = dist.cdf( x ); + y = dist.cdf( x[ i % len ] ); if ( isnan( y ) ) { b.fail( 'should not return NaN' ); } @@ -407,6 +486,7 @@ bench( pkg+':logcdf', function benchmark( b ) { var alpha; var beta; var dist; + var len; var x; var y; var i; @@ -414,11 +494,15 @@ bench( pkg+':logcdf', function benchmark( b ) { alpha = 100.56789; beta = 55.54321; dist = new Gamma( alpha, beta ); + len = 100; + x = new Float64Array( len ); + for ( i = 0; i < len; i++ ) { + x[ i ] = uniform( 0.0, 1.0 ); + } b.tic(); for ( i = 0; i < b.iterations; i++ ) { - x = randu(); - y = dist.logcdf( x ); + y = dist.logcdf( x[ i % len ] ); if ( isnan( y ) ) { b.fail( 'should not return NaN' ); } @@ -435,6 +519,7 @@ bench( pkg+':logpdf', function benchmark( b ) { var alpha; var beta; var dist; + var len; var x; var y; var i; @@ -442,11 +527,15 @@ bench( pkg+':logpdf', function benchmark( b ) { alpha = 100.56789; beta = 55.54321; dist = new Gamma( alpha, beta ); + len = 100; + x = new Float64Array( len ); + for ( i = 0; i < len; i++ ) { + x[ i ] = uniform( 0.0, 1.0 ); + } b.tic(); for ( i = 0; i < b.iterations; i++ ) { - x = randu(); - y = dist.logpdf( x ); + y = dist.logpdf( x[ i % len ] ); if ( isnan( y ) ) { b.fail( 'should not return NaN' ); } @@ -463,6 +552,7 @@ bench( pkg+':mgf', function benchmark( b ) { var alpha; var beta; var dist; + var len; var x; var y; var i; @@ -470,11 +560,15 @@ bench( pkg+':mgf', function benchmark( b ) { alpha = 100.56789; beta = 55.54321; dist = new Gamma( alpha, beta ); + len = 100; + x = new Float64Array( len ); + for ( i = 0; i < len; i++ ) { + x[ i ] = uniform( 0.0, 1.0 ); + } b.tic(); for ( i = 0; i < b.iterations; i++ ) { - x = randu(); - y = dist.mgf( x ); + y = dist.mgf( x[ i % len ] ); if ( isnan( y ) ) { b.fail( 'should not return NaN' ); } @@ -491,6 +585,7 @@ bench( pkg+':pdf', function benchmark( b ) { var alpha; var beta; var dist; + var len; var x; var y; var i; @@ -498,11 +593,15 @@ bench( pkg+':pdf', function benchmark( b ) { alpha = 100.56789; beta = 55.54321; dist = new Gamma( alpha, beta ); + len = 100; + x = new Float64Array( len ); + for ( i = 0; i < len; i++ ) { + x[ i ] = uniform( 0.0, 1.0 ); + } b.tic(); for ( i = 0; i < b.iterations; i++ ) { - x = randu(); - y = dist.pdf( x ); + y = dist.pdf( x[ i % len ] ); if ( isnan( y ) ) { b.fail( 'should not return NaN' ); } @@ -519,6 +618,7 @@ bench( pkg+':quantile', function benchmark( b ) { var alpha; var beta; var dist; + var len; var x; var y; var i; @@ -526,11 +626,15 @@ bench( pkg+':quantile', function benchmark( b ) { alpha = 100.56789; beta = 55.54321; dist = new Gamma( alpha, beta ); + len = 100; + x = new Float64Array( len ); + for ( i = 0; i < len; i++ ) { + x[ i ] = uniform( 0.0, 1.0 ); + } b.tic(); for ( i = 0; i < b.iterations; i++ ) { - x = randu(); - y = dist.quantile( x ); + y = dist.quantile( x[ i % len ] ); if ( isnan( y ) ) { b.fail( 'should not return NaN' ); } diff --git a/lib/node_modules/@stdlib/stats/base/dists/gamma/entropy/benchmark/benchmark.js b/lib/node_modules/@stdlib/stats/base/dists/gamma/entropy/benchmark/benchmark.js index a53a74eb6bab..e6546d02fb66 100644 --- a/lib/node_modules/@stdlib/stats/base/dists/gamma/entropy/benchmark/benchmark.js +++ b/lib/node_modules/@stdlib/stats/base/dists/gamma/entropy/benchmark/benchmark.js @@ -21,7 +21,8 @@ // MODULES // var bench = require( '@stdlib/bench' ); -var randu = require( '@stdlib/random/base/randu' ); +var uniform = require( '@stdlib/random/base/uniform' ); +var Float64Array = require( '@stdlib/array/float64' ); var isnan = require( '@stdlib/math/base/assert/is-nan' ); var EPS = require( '@stdlib/constants/float64/eps' ); var pkg = require( './../package.json' ).name; @@ -33,14 +34,21 @@ var entropy = require( './../lib' ); bench( pkg, function benchmark( b ) { var alpha; var beta; + var len; var y; var i; + len = 100; + alpha = new Float64Array( len ); + beta = new Float64Array( len ); + for ( i = 0; i < len; i++ ) { + alpha[ i ] = uniform( EPS, 10.0 ); + beta[ i ] = uniform( EPS, 10.0 ); + } + b.tic(); for ( i = 0; i < b.iterations; i++ ) { - alpha = ( randu()*10.0 ) + EPS; - beta = ( randu()*10.0 ) + EPS; - y = entropy( alpha, beta ); + y = entropy( alpha[ i % len ], beta[ i % len ] ); if ( isnan( y ) ) { b.fail( 'should not return NaN' ); } diff --git a/lib/node_modules/@stdlib/stats/base/dists/gamma/kurtosis/benchmark/benchmark.js b/lib/node_modules/@stdlib/stats/base/dists/gamma/kurtosis/benchmark/benchmark.js index 467e4afd5c70..d074ef57515b 100644 --- a/lib/node_modules/@stdlib/stats/base/dists/gamma/kurtosis/benchmark/benchmark.js +++ b/lib/node_modules/@stdlib/stats/base/dists/gamma/kurtosis/benchmark/benchmark.js @@ -22,7 +22,7 @@ var bench = require( '@stdlib/bench' ); var Float64Array = require( '@stdlib/array/float64' ); -var randu = require( '@stdlib/random/base/randu' ); +var uniform = require( '@stdlib/random/base/uniform' ); var isnan = require( '@stdlib/math/base/assert/is-nan' ); var EPS = require( '@stdlib/constants/float64/eps' ); var pkg = require( './../package.json' ).name; @@ -42,8 +42,8 @@ bench( pkg, function benchmark( b ) { alpha = new Float64Array( len ); beta = new Float64Array( len ); for ( i = 0; i < len; i++ ) { - alpha[ i ] = ( randu() * 20.0 ) + EPS; - beta[ i ] = ( randu() * 20.0 ) + EPS; + alpha[ i ] = uniform( EPS, 20.0 ); + beta[ i ] = uniform( EPS, 20.0 ); } b.tic(); diff --git a/lib/node_modules/@stdlib/stats/base/dists/gamma/kurtosis/benchmark/benchmark.native.js b/lib/node_modules/@stdlib/stats/base/dists/gamma/kurtosis/benchmark/benchmark.native.js index 654ad7f26cab..88160d21e836 100644 --- a/lib/node_modules/@stdlib/stats/base/dists/gamma/kurtosis/benchmark/benchmark.native.js +++ b/lib/node_modules/@stdlib/stats/base/dists/gamma/kurtosis/benchmark/benchmark.native.js @@ -24,7 +24,7 @@ var resolve = require( 'path' ).resolve; var bench = require( '@stdlib/bench' ); var Float64Array = require( '@stdlib/array/float64' ); var tryRequire = require( '@stdlib/utils/try-require' ); -var randu = require( '@stdlib/random/base/randu' ); +var uniform = require( '@stdlib/random/base/uniform' ); var isnan = require( '@stdlib/math/base/assert/is-nan' ); var EPS = require( '@stdlib/constants/float64/eps' ); var pkg = require( './../package.json' ).name; @@ -51,8 +51,8 @@ bench( pkg+'::native', opts, function benchmark( b ) { alpha = new Float64Array( len ); beta = new Float64Array( len ); for ( i = 0; i < len; i++ ) { - alpha[ i ] = ( randu() * 20.0 ) + EPS; - beta[ i ] = ( randu() * 20.0 ) + EPS; + alpha[ i ] = uniform( EPS, 20.0 ); + beta[ i ] = uniform( EPS, 20.0 ); } b.tic(); diff --git a/lib/node_modules/@stdlib/stats/base/dists/gamma/logcdf/benchmark/benchmark.js b/lib/node_modules/@stdlib/stats/base/dists/gamma/logcdf/benchmark/benchmark.js index 0021cc6767e4..4ef48333accb 100644 --- a/lib/node_modules/@stdlib/stats/base/dists/gamma/logcdf/benchmark/benchmark.js +++ b/lib/node_modules/@stdlib/stats/base/dists/gamma/logcdf/benchmark/benchmark.js @@ -21,7 +21,8 @@ // MODULES // var bench = require( '@stdlib/bench' ); -var randu = require( '@stdlib/random/base/randu' ); +var Float64Array = require( '@stdlib/array/float64' ); +var uniform = require( '@stdlib/random/base/uniform' ); var isnan = require( '@stdlib/math/base/assert/is-nan' ); var EPS = require( '@stdlib/constants/float64/eps' ); var pkg = require( './../package.json' ).name; @@ -33,16 +34,24 @@ var logcdf = require( './../lib' ); bench( pkg, function benchmark( b ) { var alpha; var beta; + var len; var x; var y; var i; + len = 100; + x = new Float64Array( len ); + alpha = new Float64Array( len ); + beta = new Float64Array( len ); + for ( i = 0; i < len; i++ ) { + x[ i ] = uniform( EPS, 100.0 ); + alpha[ i ] = uniform( EPS, 100.0 ); + beta[ i ] = uniform( EPS, 100.0 ); + } + b.tic(); for ( i = 0; i < b.iterations; i++ ) { - x = ( randu()*100.0 ) + EPS; - alpha = ( randu()*100.0 ) + EPS; - beta = ( randu()*100.0 ) + EPS; - y = logcdf( x, alpha, beta ); + y = logcdf( x[ i % len ], alpha[ i % len ], beta[ i % len ] ); if ( isnan( y ) ) { b.fail( 'should not return NaN' ); } @@ -59,6 +68,7 @@ bench( pkg+':factory', function benchmark( b ) { var mylogcdf; var alpha; var beta; + var len; var x; var y; var i; @@ -66,11 +76,15 @@ bench( pkg+':factory', function benchmark( b ) { alpha = 20.0; beta = 15.0; mylogcdf = logcdf.factory( alpha, beta ); + len = 100; + x = new Float64Array( len ); + for ( i = 0; i < len; i++ ) { + x[ i ] = uniform( EPS, 50.0 ); + } b.tic(); for ( i = 0; i < b.iterations; i++ ) { - x = ( randu()*50.0 ) + EPS; - y = mylogcdf( x ); + y = mylogcdf( x[ i % len ] ); if ( isnan( y ) ) { b.fail( 'should not return NaN' ); } diff --git a/lib/node_modules/@stdlib/stats/base/dists/gamma/logpdf/benchmark/benchmark.js b/lib/node_modules/@stdlib/stats/base/dists/gamma/logpdf/benchmark/benchmark.js index fd84a6dec988..22e5c3027a86 100644 --- a/lib/node_modules/@stdlib/stats/base/dists/gamma/logpdf/benchmark/benchmark.js +++ b/lib/node_modules/@stdlib/stats/base/dists/gamma/logpdf/benchmark/benchmark.js @@ -21,7 +21,8 @@ // MODULES // var bench = require( '@stdlib/bench' ); -var randu = require( '@stdlib/random/base/randu' ); +var Float64Array = require( '@stdlib/array/float64' ); +var uniform = require( '@stdlib/random/base/uniform' ); var isnan = require( '@stdlib/math/base/assert/is-nan' ); var EPS = require( '@stdlib/constants/float64/eps' ); var pkg = require( './../package.json' ).name; @@ -33,16 +34,24 @@ var logpdf = require( './../lib' ); bench( pkg, function benchmark( b ) { var alpha; var beta; + var len; var x; var y; var i; + len = 100; + x = new Float64Array( len ); + alpha = new Float64Array( len ); + beta = new Float64Array( len ); + for ( i = 0; i < len; i++ ) { + x[ i ] = uniform( EPS, 100.0 ); + alpha[ i ] = uniform( EPS, 100.0 ); + beta[ i ] = uniform( EPS, 100.0 ); + } + b.tic(); for ( i = 0; i < b.iterations; i++ ) { - x = ( randu()*100.0 ) + EPS; - alpha = ( randu()*100.0 ) + EPS; - beta = ( randu()*100.0 ) + EPS; - y = logpdf( x, alpha, beta ); + y = logpdf( x[ i % len ], alpha[ i % len ], beta[ i % len ] ); if ( isnan( y ) ) { b.fail( 'should not return NaN' ); } @@ -59,6 +68,7 @@ bench( pkg+':factory', function benchmark( b ) { var mylogpdf; var alpha; var beta; + var len; var x; var y; var i; @@ -66,11 +76,15 @@ bench( pkg+':factory', function benchmark( b ) { alpha = 20.0; beta = 15.0; mylogpdf = logpdf.factory( alpha, beta ); + len = 100; + x = new Float64Array( len ); + for ( i = 0; i < len; i++ ) { + x[ i ] = uniform( EPS, 50.0 ); + } b.tic(); for ( i = 0; i < b.iterations; i++ ) { - x = ( randu()*50.0 ) + EPS; - y = mylogpdf( x ); + y = mylogpdf( x[ i % len ] ); if ( isnan( y ) ) { b.fail( 'should not return NaN' ); } diff --git a/lib/node_modules/@stdlib/stats/base/dists/gamma/mean/benchmark/benchmark.js b/lib/node_modules/@stdlib/stats/base/dists/gamma/mean/benchmark/benchmark.js index e7ca2d7aa5a2..8bd0d699b1ca 100644 --- a/lib/node_modules/@stdlib/stats/base/dists/gamma/mean/benchmark/benchmark.js +++ b/lib/node_modules/@stdlib/stats/base/dists/gamma/mean/benchmark/benchmark.js @@ -22,7 +22,7 @@ var bench = require( '@stdlib/bench' ); var Float64Array = require( '@stdlib/array/float64' ); -var randu = require( '@stdlib/random/base/randu' ); +var uniform = require( '@stdlib/random/base/uniform' ); var isnan = require( '@stdlib/math/base/assert/is-nan' ); var EPS = require( '@stdlib/constants/float64/eps' ); var pkg = require( './../package.json' ).name; @@ -42,8 +42,8 @@ bench( pkg, function benchmark( b ) { alpha = new Float64Array( len ); beta = new Float64Array( len ); for ( i = 0; i < len; i++ ) { - alpha[ i ] = ( randu() * 20.0 ) + EPS; - beta[ i ] = ( randu() * 20.0 ) + EPS; + alpha[ i ] = uniform( EPS, 20.0 ); + beta[ i ] = uniform( EPS, 20.0 ); } b.tic(); diff --git a/lib/node_modules/@stdlib/stats/base/dists/gamma/mean/benchmark/benchmark.native.js b/lib/node_modules/@stdlib/stats/base/dists/gamma/mean/benchmark/benchmark.native.js index ff7a0e749af7..336ebd70f3df 100644 --- a/lib/node_modules/@stdlib/stats/base/dists/gamma/mean/benchmark/benchmark.native.js +++ b/lib/node_modules/@stdlib/stats/base/dists/gamma/mean/benchmark/benchmark.native.js @@ -24,7 +24,7 @@ var resolve = require( 'path' ).resolve; var bench = require( '@stdlib/bench' ); var Float64Array = require( '@stdlib/array/float64' ); var tryRequire = require( '@stdlib/utils/try-require' ); -var randu = require( '@stdlib/random/base/randu' ); +var uniform = require( '@stdlib/random/base/uniform' ); var isnan = require( '@stdlib/math/base/assert/is-nan' ); var EPS = require( '@stdlib/constants/float64/eps' ); var pkg = require( './../package.json' ).name; @@ -51,8 +51,8 @@ bench( pkg+'::native', opts, function benchmark( b ) { alpha = new Float64Array( len ); beta = new Float64Array( len ); for ( i = 0; i < len; i++ ) { - alpha[ i ] = ( randu() * 20.0 ) + EPS; - beta[ i ] = ( randu() * 20.0 ) + EPS; + alpha[ i ] = uniform( EPS, 20.0 ); + beta[ i ] = uniform( EPS, 20.0 ); } b.tic(); diff --git a/lib/node_modules/@stdlib/stats/base/dists/gamma/mgf/benchmark/benchmark.js b/lib/node_modules/@stdlib/stats/base/dists/gamma/mgf/benchmark/benchmark.js index 8f3d54042cda..4c9acf0818af 100644 --- a/lib/node_modules/@stdlib/stats/base/dists/gamma/mgf/benchmark/benchmark.js +++ b/lib/node_modules/@stdlib/stats/base/dists/gamma/mgf/benchmark/benchmark.js @@ -21,7 +21,8 @@ // MODULES // var bench = require( '@stdlib/bench' ); -var randu = require( '@stdlib/random/base/randu' ); +var uniform = require( '@stdlib/random/base/uniform' ); +var Float64Array = require( '@stdlib/array/float64' ); var isnan = require( '@stdlib/math/base/assert/is-nan' ); var EPS = require( '@stdlib/constants/float64/eps' ); var pkg = require( './../package.json' ).name; @@ -33,16 +34,24 @@ var mgf = require( './../lib' ); bench( pkg, function benchmark( b ) { var alpha; var beta; + var len; var t; var y; var i; + len = 100; + alpha = new Float64Array( len ); + beta = new Float64Array( len ); + t = new Float64Array( len ); + for ( i = 0; i < len; i++ ) { + alpha[ i ] = uniform( EPS, 100.0 ); + beta[ i ] = uniform( EPS, 100.0 ); + t[ i ] = uniform( 0.0, beta[ i ] ); + } + b.tic(); for ( i = 0; i < b.iterations; i++ ) { - alpha = ( randu()*100.0 ) + EPS; - beta = ( randu()*100.0 ) + EPS; - t = randu() * beta; - y = mgf( t, alpha, beta ); + y = mgf( t[ i % len ], alpha[ i % len ], beta[ i % len ] ); if ( isnan( y ) ) { b.fail( 'should not return NaN' ); } @@ -59,6 +68,7 @@ bench( pkg+':factory', function benchmark( b ) { var mymgf; var alpha; var beta; + var len; var t; var y; var i; @@ -66,11 +76,15 @@ bench( pkg+':factory', function benchmark( b ) { alpha = 20.0; beta = 15.0; mymgf = mgf.factory( alpha, beta ); + len = 100; + t = new Float64Array( len ); + for ( i = 0; i < len; i++ ) { + t[ i ] = uniform( 0.0, beta ); + } b.tic(); for ( i = 0; i < b.iterations; i++ ) { - t = randu() * beta; - y = mymgf( t ); + y = mymgf( t[ i % len ] ); if ( isnan( y ) ) { b.fail( 'should not return NaN' ); } diff --git a/lib/node_modules/@stdlib/stats/base/dists/gamma/mode/benchmark/benchmark.js b/lib/node_modules/@stdlib/stats/base/dists/gamma/mode/benchmark/benchmark.js index f0ddab5bc9fa..42e5c8bd2082 100644 --- a/lib/node_modules/@stdlib/stats/base/dists/gamma/mode/benchmark/benchmark.js +++ b/lib/node_modules/@stdlib/stats/base/dists/gamma/mode/benchmark/benchmark.js @@ -22,7 +22,7 @@ var bench = require( '@stdlib/bench' ); var Float64Array = require( '@stdlib/array/float64' ); -var randu = require( '@stdlib/random/base/randu' ); +var uniform = require( '@stdlib/random/base/uniform' ); var isnan = require( '@stdlib/math/base/assert/is-nan' ); var EPS = require( '@stdlib/constants/float64/eps' ); var pkg = require( './../package.json' ).name; @@ -42,8 +42,8 @@ bench( pkg, function benchmark( b ) { alpha = new Float64Array( len ); beta = new Float64Array( len ); for ( i = 0; i < len; i++ ) { - alpha[ i ] = ( randu() * 10.0 ) + 1.0 + EPS; - beta[ i ] = ( randu() * 10.0 ) + EPS; + alpha[ i ] = uniform( 1.0 + EPS, 10.0 ); + beta[ i ] = uniform( EPS, 10.0 ); } b.tic(); diff --git a/lib/node_modules/@stdlib/stats/base/dists/gamma/mode/benchmark/benchmark.native.js b/lib/node_modules/@stdlib/stats/base/dists/gamma/mode/benchmark/benchmark.native.js index 6d6727e2e48c..5e821e95d87d 100644 --- a/lib/node_modules/@stdlib/stats/base/dists/gamma/mode/benchmark/benchmark.native.js +++ b/lib/node_modules/@stdlib/stats/base/dists/gamma/mode/benchmark/benchmark.native.js @@ -24,7 +24,7 @@ var resolve = require( 'path' ).resolve; var bench = require( '@stdlib/bench' ); var Float64Array = require( '@stdlib/array/float64' ); var tryRequire = require( '@stdlib/utils/try-require' ); -var randu = require( '@stdlib/random/base/randu' ); +var uniform = require( '@stdlib/random/base/uniform' ); var isnan = require( '@stdlib/math/base/assert/is-nan' ); var EPS = require( '@stdlib/constants/float64/eps' ); var pkg = require( './../package.json' ).name; @@ -51,8 +51,8 @@ bench( pkg+'::native', opts, function benchmark( b ) { alpha = new Float64Array( len ); beta = new Float64Array( len ); for ( i = 0; i < len; i++ ) { - alpha[ i ] = ( randu() * 10.0 ) + 1.0 + EPS; - beta[ i ] = ( randu() * 10.0 ) + EPS; + alpha[ i ] = uniform( 1.0 + EPS, 10.0 ); + beta[ i ] = uniform( EPS, 10.0 ); } b.tic(); diff --git a/lib/node_modules/@stdlib/stats/base/dists/gamma/pdf/benchmark/benchmark.js b/lib/node_modules/@stdlib/stats/base/dists/gamma/pdf/benchmark/benchmark.js index dbe69a73f5f1..3365a4416f16 100644 --- a/lib/node_modules/@stdlib/stats/base/dists/gamma/pdf/benchmark/benchmark.js +++ b/lib/node_modules/@stdlib/stats/base/dists/gamma/pdf/benchmark/benchmark.js @@ -21,7 +21,8 @@ // MODULES // var bench = require( '@stdlib/bench' ); -var randu = require( '@stdlib/random/base/randu' ); +var Float64Array = require( '@stdlib/array/float64' ); +var uniform = require( '@stdlib/random/base/uniform' ); var isnan = require( '@stdlib/math/base/assert/is-nan' ); var EPS = require( '@stdlib/constants/float64/eps' ); var pkg = require( './../package.json' ).name; @@ -33,16 +34,24 @@ var pdf = require( './../lib' ); bench( pkg, function benchmark( b ) { var alpha; var beta; + var len; var x; var y; var i; + len = 100; + x = new Float64Array( len ); + alpha = new Float64Array( len ); + beta = new Float64Array( len ); + for ( i = 0; i < len; i++ ) { + x[ i ] = uniform( EPS, 100.0 ); + alpha[ i ] = uniform( EPS, 100.0 ); + beta[ i ] = uniform( EPS, 100.0 ); + } + b.tic(); for ( i = 0; i < b.iterations; i++ ) { - x = ( randu()*100.0 ) + EPS; - alpha = ( randu()*100.0 ) + EPS; - beta = ( randu()*100.0 ) + EPS; - y = pdf( x, alpha, beta ); + y = pdf( x[ i % len ], alpha[ i % len ], beta[ i % len ] ); if ( isnan( y ) ) { b.fail( 'should not return NaN' ); } @@ -59,6 +68,7 @@ bench( pkg+':factory', function benchmark( b ) { var mypdf; var alpha; var beta; + var len; var x; var y; var i; @@ -66,11 +76,15 @@ bench( pkg+':factory', function benchmark( b ) { alpha = 20.0; beta = 15.0; mypdf = pdf.factory( alpha, beta ); + len = 100; + x = new Float64Array( len ); + for ( i = 0; i < len; i++ ) { + x[ i ] = uniform( EPS, 50.0 ); + } b.tic(); for ( i = 0; i < b.iterations; i++ ) { - x = ( randu()*50.0 ) + EPS; - y = mypdf( x ); + y = mypdf( x[ i % len ] ); if ( isnan( y ) ) { b.fail( 'should not return NaN' ); } diff --git a/lib/node_modules/@stdlib/stats/base/dists/gamma/quantile/benchmark/benchmark.js b/lib/node_modules/@stdlib/stats/base/dists/gamma/quantile/benchmark/benchmark.js index 14ccc93ec6b0..77305a94d894 100644 --- a/lib/node_modules/@stdlib/stats/base/dists/gamma/quantile/benchmark/benchmark.js +++ b/lib/node_modules/@stdlib/stats/base/dists/gamma/quantile/benchmark/benchmark.js @@ -21,7 +21,8 @@ // MODULES // var bench = require( '@stdlib/bench' ); -var randu = require( '@stdlib/random/base/randu' ); +var Float64Array = require( '@stdlib/array/float64' ); +var uniform = require( '@stdlib/random/base/uniform' ); var isnan = require( '@stdlib/math/base/assert/is-nan' ); var EPS = require( '@stdlib/constants/float64/eps' ); var pkg = require( './../package.json' ).name; @@ -33,16 +34,24 @@ var quantile = require( './../lib' ); bench( pkg, function benchmark( b ) { var alpha; var beta; + var len; var p; var y; var i; + len = 100; + p = new Float64Array( len ); + alpha = new Float64Array( len ); + beta = new Float64Array( len ); + for ( i = 0; i < len; i++ ) { + p[ i ] = uniform( 0.0, 1.0 ); + alpha[ i ] = uniform( EPS, 100.0 ); + beta[ i ] = uniform( EPS, 100.0 ); + } + b.tic(); for ( i = 0; i < b.iterations; i++ ) { - p = randu(); - alpha = ( randu()*100.0 ) + EPS; - beta = ( randu()*100.0 ) + EPS; - y = quantile( p, alpha, beta ); + y = quantile( p[ i % len ], alpha[ i % len ], beta[ i % len ] ); if ( isnan( y ) ) { b.fail( 'should not return NaN' ); } @@ -59,6 +68,7 @@ bench( pkg+':factory', function benchmark( b ) { var myquantile; var alpha; var beta; + var len; var p; var y; var i; @@ -66,11 +76,15 @@ bench( pkg+':factory', function benchmark( b ) { alpha = 20.0; beta = 15.0; myquantile = quantile.factory( alpha, beta ); + len = 100; + p = new Float64Array( len ); + for ( i = 0; i < len; i++ ) { + p[ i ] = uniform( 0.0, 1.0 ); + } b.tic(); for ( i = 0; i < b.iterations; i++ ) { - p = randu(); - y = myquantile( p ); + y = myquantile( p[ i % len ] ); if ( isnan( y ) ) { b.fail( 'should not return NaN' ); } diff --git a/lib/node_modules/@stdlib/stats/base/dists/gamma/skewness/benchmark/benchmark.js b/lib/node_modules/@stdlib/stats/base/dists/gamma/skewness/benchmark/benchmark.js index 174f6e8fa6db..64d525c2ed13 100644 --- a/lib/node_modules/@stdlib/stats/base/dists/gamma/skewness/benchmark/benchmark.js +++ b/lib/node_modules/@stdlib/stats/base/dists/gamma/skewness/benchmark/benchmark.js @@ -22,7 +22,7 @@ var bench = require( '@stdlib/bench' ); var Float64Array = require( '@stdlib/array/float64' ); -var randu = require( '@stdlib/random/base/randu' ); +var uniform = require( '@stdlib/random/base/uniform' ); var isnan = require( '@stdlib/math/base/assert/is-nan' ); var EPS = require( '@stdlib/constants/float64/eps' ); var pkg = require( './../package.json' ).name; @@ -42,8 +42,8 @@ bench( pkg, function benchmark( b ) { alpha = new Float64Array( len ); beta = new Float64Array( len ); for ( i = 0; i < len; i++ ) { - alpha[ i ] = ( randu() * 10.0 ) + EPS; - beta[ i ] = ( randu() * 10.0 ) + EPS; + alpha[ i ] = uniform( EPS, 10.0 ); + beta[ i ] = uniform( EPS, 10.0 ); } b.tic(); diff --git a/lib/node_modules/@stdlib/stats/base/dists/gamma/skewness/benchmark/benchmark.native.js b/lib/node_modules/@stdlib/stats/base/dists/gamma/skewness/benchmark/benchmark.native.js index 8142d99e30dc..2a7a09604bc5 100644 --- a/lib/node_modules/@stdlib/stats/base/dists/gamma/skewness/benchmark/benchmark.native.js +++ b/lib/node_modules/@stdlib/stats/base/dists/gamma/skewness/benchmark/benchmark.native.js @@ -24,7 +24,7 @@ var resolve = require( 'path' ).resolve; var bench = require( '@stdlib/bench' ); var Float64Array = require( '@stdlib/array/float64' ); var tryRequire = require( '@stdlib/utils/try-require' ); -var randu = require( '@stdlib/random/base/randu' ); +var uniform = require( '@stdlib/random/base/uniform' ); var isnan = require( '@stdlib/math/base/assert/is-nan' ); var EPS = require( '@stdlib/constants/float64/eps' ); var pkg = require( './../package.json' ).name; @@ -51,8 +51,8 @@ bench( pkg+'::native', opts, function benchmark( b ) { alpha = new Float64Array( len ); beta = new Float64Array( len ); for ( i = 0; i < len; i++ ) { - alpha[ i ] = ( randu() * 20.0 ) + EPS; - beta[ i ] = ( randu() * 20.0 ) + EPS; + alpha[ i ] = uniform( EPS, 10.0 ); + beta[ i ] = uniform( EPS, 10.0 ); } b.tic(); diff --git a/lib/node_modules/@stdlib/stats/base/dists/gamma/stdev/benchmark/benchmark.js b/lib/node_modules/@stdlib/stats/base/dists/gamma/stdev/benchmark/benchmark.js index 3ca20815d248..2d9142f8845c 100644 --- a/lib/node_modules/@stdlib/stats/base/dists/gamma/stdev/benchmark/benchmark.js +++ b/lib/node_modules/@stdlib/stats/base/dists/gamma/stdev/benchmark/benchmark.js @@ -22,7 +22,7 @@ var bench = require( '@stdlib/bench' ); var Float64Array = require( '@stdlib/array/float64' ); -var randu = require( '@stdlib/random/base/randu' ); +var uniform = require( '@stdlib/random/base/uniform' ); var isnan = require( '@stdlib/math/base/assert/is-nan' ); var EPS = require( '@stdlib/constants/float64/eps' ); var pkg = require( './../package.json' ).name; @@ -42,8 +42,8 @@ bench( pkg, function benchmark( b ) { alpha = new Float64Array( len ); beta = new Float64Array( len ); for ( i = 0; i < len; i++ ) { - alpha[ i ] = ( randu() * 20.0 ) + EPS; - beta[ i ] = ( randu() * 20.0 ) + EPS; + alpha[ i ] = uniform( EPS, 20.0 ); + beta[ i ] = uniform( EPS, 20.0 ); } b.tic(); diff --git a/lib/node_modules/@stdlib/stats/base/dists/gamma/stdev/benchmark/benchmark.native.js b/lib/node_modules/@stdlib/stats/base/dists/gamma/stdev/benchmark/benchmark.native.js index 1f50463f91a0..0f63b3b12d71 100644 --- a/lib/node_modules/@stdlib/stats/base/dists/gamma/stdev/benchmark/benchmark.native.js +++ b/lib/node_modules/@stdlib/stats/base/dists/gamma/stdev/benchmark/benchmark.native.js @@ -24,7 +24,7 @@ var resolve = require( 'path' ).resolve; var bench = require( '@stdlib/bench' ); var Float64Array = require( '@stdlib/array/float64' ); var tryRequire = require( '@stdlib/utils/try-require' ); -var randu = require( '@stdlib/random/base/randu' ); +var uniform = require( '@stdlib/random/base/uniform' ); var isnan = require( '@stdlib/math/base/assert/is-nan' ); var EPS = require( '@stdlib/constants/float64/eps' ); var pkg = require( './../package.json' ).name; @@ -51,8 +51,8 @@ bench( pkg+'::native', opts, function benchmark( b ) { alpha = new Float64Array( len ); beta = new Float64Array( len ); for ( i = 0; i < len; i++ ) { - alpha[ i ] = ( randu() * 20.0 ) + EPS; - beta[ i ] = ( randu() * 20.0 ) + EPS; + alpha[ i ] = uniform( EPS, 20.0 ); + beta[ i ] = uniform( EPS, 20.0 ); } b.tic(); diff --git a/lib/node_modules/@stdlib/stats/base/dists/gamma/variance/benchmark/benchmark.js b/lib/node_modules/@stdlib/stats/base/dists/gamma/variance/benchmark/benchmark.js index fa0fb9917365..35413c3b5d8b 100644 --- a/lib/node_modules/@stdlib/stats/base/dists/gamma/variance/benchmark/benchmark.js +++ b/lib/node_modules/@stdlib/stats/base/dists/gamma/variance/benchmark/benchmark.js @@ -21,7 +21,8 @@ // MODULES // var bench = require( '@stdlib/bench' ); -var randu = require( '@stdlib/random/base/randu' ); +var Float64Array = require( '@stdlib/array/float64' ); +var uniform = require( '@stdlib/random/base/uniform' ); var isnan = require( '@stdlib/math/base/assert/is-nan' ); var EPS = require( '@stdlib/constants/float64/eps' ); var pkg = require( './../package.json' ).name; @@ -33,14 +34,21 @@ var variance = require( './../lib' ); bench( pkg, function benchmark( b ) { var alpha; var beta; + var len; var y; var i; + len = 100; + alpha = new Float64Array( len ); + beta = new Float64Array( len ); + for ( i = 0; i < len; i++ ) { + alpha[ i ] = uniform( EPS, 10.0 ); + beta[ i ] = uniform( EPS, 10.0 ); + } + b.tic(); for ( i = 0; i < b.iterations; i++ ) { - alpha = ( randu()*10.0 ) + EPS; - beta = ( randu()*10.0 ) + EPS; - y = variance( alpha, beta ); + y = variance( alpha[ i % len ], beta[ i % len ] ); if ( isnan( y ) ) { b.fail( 'should not return NaN' ); }