21
21
// MODULES //
22
22
23
23
var bench = require ( '@stdlib/bench' ) ;
24
- var randu = require ( '@stdlib/random/base/randu' ) ;
24
+ var Float64Array = require ( '@stdlib/array/float64' ) ;
25
+ var uniform = require ( '@stdlib/random/base/uniform' ) ;
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,19 +33,28 @@ var cdf = require( './../lib' );
32
33
33
34
bench ( pkg , function benchmark ( b ) {
34
35
var alpha ;
36
+ var len ;
35
37
var m ;
36
38
var s ;
37
39
var x ;
38
40
var y ;
39
41
var i ;
40
42
43
+ len = 100 ;
44
+ x = new Float64Array ( len ) ;
45
+ alpha = new Float64Array ( len ) ;
46
+ s = new Float64Array ( len ) ;
47
+ m = new Float64Array ( len ) ;
48
+ for ( i = 0 ; i < len ; i ++ ) {
49
+ x [ i ] = uniform ( - 50.0 , 50.0 ) ;
50
+ alpha [ i ] = uniform ( EPS , 20.0 ) ;
51
+ m [ i ] = uniform ( - 20.0 , 40.0 ) ;
52
+ s [ i ] = uniform ( EPS , 20.0 ) ;
53
+ }
54
+
41
55
b . tic ( ) ;
42
56
for ( i = 0 ; i < b . iterations ; i ++ ) {
43
- x = ( randu ( ) * 100.0 ) - 50.0 ;
44
- alpha = ( randu ( ) * 20.0 ) + EPS ;
45
- m = ( randu ( ) * 60.0 ) - 20.0 ;
46
- s = ( randu ( ) * 20.0 ) + EPS ;
47
- y = cdf ( x , alpha , s , m ) ;
57
+ y = cdf ( x [ i % len ] , alpha [ i % len ] , s [ i % len ] , m [ i % len ] ) ;
48
58
if ( isnan ( y ) ) {
49
59
b . fail ( 'should not return NaN' ) ;
50
60
}
@@ -60,6 +70,7 @@ bench( pkg, function benchmark( b ) {
60
70
bench ( pkg + ':factory' , function benchmark ( b ) {
61
71
var mycdf ;
62
72
var alpha ;
73
+ var len ;
63
74
var m ;
64
75
var s ;
65
76
var x ;
@@ -70,11 +81,15 @@ bench( pkg+':factory', function benchmark( b ) {
70
81
s = 3.0 ;
71
82
m = 4.0 ;
72
83
mycdf = cdf . factory ( alpha , s , m ) ;
84
+ len = 100 ;
85
+ x = new Float64Array ( len ) ;
86
+ for ( i = 0 ; i < len ; i ++ ) {
87
+ x [ i ] = uniform ( 0.0 , 50.0 ) ;
88
+ }
73
89
74
90
b . tic ( ) ;
75
91
for ( i = 0 ; i < b . iterations ; i ++ ) {
76
- x = randu ( ) * 50.0 ;
77
- y = mycdf ( x ) ;
92
+ y = mycdf ( x [ i % len ] ) ;
78
93
if ( isnan ( y ) ) {
79
94
b . fail ( 'should not return NaN' ) ;
80
95
}
0 commit comments