21
21
// MODULES //
22
22
23
23
var bench = require ( '@stdlib/bench' ) ;
24
- var randu = require ( '@stdlib/random/base/randu' ) ;
24
+ var uniform = require ( '@stdlib/random/base/uniform' ) ;
25
+ var Float64Array = require ( '@stdlib/array/float64' ) ;
25
26
var isnan = require ( '@stdlib/math/base/assert/is-nan' ) ;
26
27
var EPS = require ( '@stdlib/constants/float64/eps' ) ;
27
28
var pkg = require ( './../package.json' ) . name ;
@@ -32,17 +33,25 @@ var cdf = require( './../lib' );
32
33
33
34
bench ( pkg , function benchmark ( b ) {
34
35
var sigma ;
36
+ var len ;
35
37
var mu ;
36
38
var x ;
37
39
var y ;
38
40
var i ;
39
41
42
+ len = 100 ;
43
+ x = new Float64Array ( len ) ;
44
+ mu = new Float64Array ( len ) ;
45
+ sigma = new Float64Array ( len ) ;
46
+ for ( i = 0 ; i < len ; i ++ ) {
47
+ x [ i ] = uniform ( - 100.0 , 0.0 ) ;
48
+ mu [ i ] = uniform ( - 50.0 , 50.0 ) ;
49
+ sigma [ i ] = uniform ( EPS , 20.0 ) ;
50
+ }
51
+
40
52
b . tic ( ) ;
41
53
for ( i = 0 ; i < b . iterations ; i ++ ) {
42
- x = ( randu ( ) * 100.0 ) - 100 ;
43
- mu = ( randu ( ) * 100.0 ) - 50.0 ;
44
- sigma = ( randu ( ) * 20.0 ) + EPS ;
45
- y = cdf ( x , mu , sigma ) ;
54
+ y = cdf ( x [ i % len ] , mu [ i % len ] , sigma [ i % len ] ) ;
46
55
if ( isnan ( y ) ) {
47
56
b . fail ( 'should not return NaN' ) ;
48
57
}
@@ -58,6 +67,7 @@ bench( pkg, function benchmark( b ) {
58
67
bench ( pkg + ':factory' , function benchmark ( b ) {
59
68
var mycdf ;
60
69
var sigma ;
70
+ var len ;
61
71
var mu ;
62
72
var x ;
63
73
var y ;
@@ -66,11 +76,15 @@ bench( pkg+':factory', function benchmark( b ) {
66
76
mu = 10.0 ;
67
77
sigma = 4.0 ;
68
78
mycdf = cdf . factory ( mu , sigma ) ;
79
+ len = 100 ;
80
+ x = new Float64Array ( len ) ;
81
+ for ( i = 0 ; i < len ; i ++ ) {
82
+ x [ i ] = uniform ( 0.0 , 50.0 ) ;
83
+ }
69
84
70
85
b . tic ( ) ;
71
86
for ( i = 0 ; i < b . iterations ; i ++ ) {
72
- x = randu ( ) * 50.0 ;
73
- y = mycdf ( x ) ;
87
+ y = mycdf ( x [ i % len ] ) ;
74
88
if ( isnan ( y ) ) {
75
89
b . fail ( 'should not return NaN' ) ;
76
90
}
0 commit comments