Skip to content

Commit 6b76332

Browse files
chore: clean up
1 parent 2039b5f commit 6b76332

File tree

11 files changed

+92
-98
lines changed

11 files changed

+92
-98
lines changed

lib/node_modules/@stdlib/stats/base/dists/betaprime/logpdf/README.md

Lines changed: 11 additions & 33 deletions
Original file line numberDiff line numberDiff line change
@@ -22,7 +22,7 @@ limitations under the License.
2222

2323
> [Beta prime][betaprime-distribution] distribution logarithm of probability density function (PDF).
2424
25-
<section class="intro">
25+
<section class="intro"
2626

2727
The [probability density function][pdf] (PDF) for a [beta prime][betaprime-distribution] random variable is
2828

@@ -194,38 +194,11 @@ for ( i = 0; i < 10; i++ ) {
194194

195195
#### stdlib_base_dists_betaprime_logpdf( x, alpha, beta )
196196

197-
Evaluates the natural logarithm of the probability density function (logPDF) for a beta prime distribution with first shape parameter `alpha` and second shape parameter `beta` at a value `x`.
197+
Evaluates the natural logarithm of the probability density function (logPDF) for a beta prime distribution with first shape parameter `alpha` and second shape parameter `beta`.
198198

199199
```c
200200
double y = stdlib_base_dists_betaprime_logpdf( 0.5, 1.0, 1.0 );
201201
// returns ~-0.811
202-
203-
y = stdlib_base_dists_betaprime_logpdf( 0.5, 2.0, 4.0 );
204-
// returns ~-0.13
205-
206-
y = stdlib_base_dists_betaprime_logpdf( 0.2, 2.0, 2.0 );
207-
// returns ~-0.547
208-
209-
y = stdlib_base_dists_betaprime_logpdf( 0.8, 4.0, 4.0 );
210-
// returns ~-0.43
211-
212-
y = stdlib_base_dists_betaprime_logpdf( -0.5, 4.0, 2.0 );
213-
// returns -Infinity
214-
215-
y = stdlib_base_dists_betaprime_logpdf( 0.5, -1.0, 0.5 );
216-
// returns NaN
217-
218-
y = stdlib_base_dists_betaprime_logpdf( 0.5, 0.5, -1.0 );
219-
// returns NaN
220-
221-
y = stdlib_base_dists_betaprime_logpdf( NaN, 1.0, 1.0 );
222-
// returns NaN
223-
224-
y = stdlib_base_dists_betaprime_logpdf( 0.5, NaN, 1.0 );
225-
// returns NaN
226-
227-
y = stdlib_base_dists_betaprime_logpdf( 0.5, 1.0, NaN );
228-
// returns Na
229202
```
230203

231204
The function accepts the following arguments:
@@ -242,7 +215,7 @@ double stdlib_base_dists_betaprime_logpdf( const double x, double alpha, const d
242215
243216
<!-- /.usage -->
244217
245-
<!-- C API usage notes. Make sure to keep an empty line after the `section`
218+
<!-- C API usage notes. Make sure to keep an empty line after the `section`
246219
element and another before the `/section` close. -->
247220
248221
<section class="notes">
@@ -263,6 +236,11 @@ element and another before the `/section` close. -->
263236
#include <stdlib.h>
264237
#include <stdio.h>
265238
239+
static double random_uniform( const double min, const double max ) {
240+
double v = (double)rand() / ( (double)RAND_MAX + 1.0 );
241+
return min + ( v*(max-min) );
242+
}
243+
266244
int main( void ) {
267245
double alpha;
268246
double beta;
@@ -271,9 +249,9 @@ int main( void ) {
271249
int i;
272250
273251
for ( i = 0; i < 10; i++ ) {
274-
x = ( (double)rand() / (double)RAND_MAX );
275-
alpha = ( ( (double)rand() / (double)RAND_MAX )*5.0 ) + STDLIB_CONSTANT_FLOAT64_EPS;
276-
beta = ( ( (double)rand() / (double)RAND_MAX )*5.0 ) + STDLIB_CONSTANT_FLOAT64_EPS;
252+
x = random_uniform( 0.0, 1.0 );
253+
alpha = random_uniform( STDLIB_CONSTANT_FLOAT64_EPS, 5.0 );
254+
beta = random_uniform( STDLIB_CONSTANT_FLOAT64_EPS, 5.0 );
277255
y = stdlib_base_dists_betaprime_logpdf( x, alpha, beta );
278256
printf( "x: %1f, α: %1f, β: %1f, ln(f(x;α,β)): %lf\n", x, alpha, beta, y );
279257
}

lib/node_modules/@stdlib/stats/base/dists/betaprime/logpdf/benchmark/benchmark.js

Lines changed: 26 additions & 7 deletions
Original file line numberDiff line numberDiff line change
@@ -21,7 +21,8 @@
2121
// MODULES //
2222

2323
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' );
2526
var isnan = require( '@stdlib/math/base/assert/is-nan' );
2627
var EPS = require( '@stdlib/constants/float64/eps' );
2728
var pkg = require( './../package.json' ).name;
@@ -33,16 +34,28 @@ var logpdf = require( './../lib' );
3334
bench( pkg, function benchmark( b ) {
3435
var alpha;
3536
var beta;
37+
var len;
3638
var x;
3739
var y;
3840
var i;
3941

42+
len = 100;
43+
alpha = new Float64Array( len );
44+
for ( i = 0; i < len; i++ ) {
45+
alpha[ i ] = uniform( EPS, 100.0 );
46+
}
47+
beta = new Float64Array( len );
48+
for ( i = 0; i < len; i++ ) {
49+
beta[ i ] = uniform( EPS, 100.0 );
50+
}
51+
x = new Float64Array( len );
52+
for ( i = 0; i < len; i++ ) {
53+
x[ i ] = uniform( EPS, 2.0 );
54+
}
55+
4056
b.tic();
4157
for ( i = 0; i < b.iterations; i++ ) {
42-
x = ( randu()*2.0 ) + EPS;
43-
alpha = ( randu()*100.0 ) + EPS;
44-
beta = ( randu()*100.0 ) + EPS;
45-
y = logpdf( x, alpha, beta );
58+
y = logpdf( x[ i%100 ], alpha[ i%100 ], beta[ i%100 ] );
4659
if ( isnan( y ) ) {
4760
b.fail( 'should not return NaN' );
4861
}
@@ -59,6 +72,7 @@ bench( pkg+':factory', function benchmark( b ) {
5972
var mylogpdf;
6073
var alpha;
6174
var beta;
75+
var len;
6276
var x;
6377
var y;
6478
var i;
@@ -67,10 +81,15 @@ bench( pkg+':factory', function benchmark( b ) {
6781
beta = 55.54321;
6882
mylogpdf = logpdf.factory( alpha, beta );
6983

84+
len = 100;
85+
x = new Float64Array( len );
86+
for ( i = 0; i < len; i++ ) {
87+
x[ i ] = uniform( EPS, 2.0 );
88+
}
89+
7090
b.tic();
7191
for ( i = 0; i < b.iterations; i++ ) {
72-
x = ( randu()*2.0 ) + EPS;
73-
y = mylogpdf( x );
92+
y = mylogpdf( x[ i%100 ] );
7493
if ( isnan( y ) ) {
7594
b.fail( 'should not return NaN' );
7695
}

lib/node_modules/@stdlib/stats/base/dists/betaprime/logpdf/benchmark/benchmark.native.js

Lines changed: 15 additions & 30 deletions
Original file line numberDiff line numberDiff line change
@@ -22,7 +22,8 @@
2222

2323
var resolve = require( 'path' ).resolve;
2424
var bench = require( '@stdlib/bench' );
25-
var randu = require( '@stdlib/random/base/randu' );
25+
var uniform = require( '@stdlib/random/base/uniform' );
26+
var Float64Array = require( '@stdlib/array/float64' );
2627
var isnan = require( '@stdlib/math/base/assert/is-nan' );
2728
var tryRequire = require( '@stdlib/utils/try-require' );
2829
var EPS = require( '@stdlib/constants/float64/eps' );
@@ -42,44 +43,28 @@ var opts = {
4243
bench( pkg+'::native', opts, function benchmark( b ) {
4344
var alpha;
4445
var beta;
46+
var len;
4547
var x;
4648
var y;
4749
var i;
4850

49-
b.tic();
50-
for ( i = 0; i < b.iterations; i++ ) {
51-
x = ( randu()*2.0 ) + EPS;
52-
alpha = ( randu()*100.0 ) + EPS;
53-
beta = ( randu()*100.0 ) + EPS;
54-
y = logpdf( x, alpha, beta );
55-
if ( isnan( y ) ) {
56-
b.fail( 'should not return NaN' );
57-
}
51+
len = 100;
52+
alpha = new Float64Array( len );
53+
for ( i = 0; i < len; i++ ) {
54+
alpha[ i ] = uniform( EPS, 100.0 );
5855
}
59-
b.toc();
60-
if ( isnan( y ) ) {
61-
b.fail( 'should not return NaN' );
56+
beta = new Float64Array( len );
57+
for ( i = 0; i < len; i++ ) {
58+
beta[ i ] = uniform( EPS, 100.0 );
59+
}
60+
x = new Float64Array( len );
61+
for ( i = 0; i < len; i++ ) {
62+
x[ i ] = uniform( EPS, 2.0 );
6263
}
63-
b.pass( 'benchmark finished' );
64-
b.end();
65-
});
66-
67-
bench( pkg+':factory', opts, function benchmark( b ) {
68-
var mylogpdf;
69-
var alpha;
70-
var beta;
71-
var x;
72-
var y;
73-
var i;
74-
75-
alpha = 100.56789;
76-
beta = 55.54321;
77-
mylogpdf = logpdf.factory( alpha, beta );
7864

7965
b.tic();
8066
for ( i = 0; i < b.iterations; i++ ) {
81-
x = ( randu()*2.0 ) + EPS;
82-
y = mylogpdf( x );
67+
y = logpdf( x[ i%100 ], alpha[ i%100 ], beta[ i%100 ] );
8368
if ( isnan( y ) ) {
8469
b.fail( 'should not return NaN' );
8570
}

lib/node_modules/@stdlib/stats/base/dists/betaprime/logpdf/benchmark/c/Makefile

Lines changed: 1 addition & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -143,4 +143,4 @@ run: $(c_targets)
143143
clean:
144144
$(QUIET) -rm -f *.o *.out
145145

146-
.PHONY: clean
146+
.PHONY: clean

lib/node_modules/@stdlib/stats/base/dists/betaprime/logpdf/benchmark/c/benchmark.c

Lines changed: 16 additions & 9 deletions
Original file line numberDiff line numberDiff line change
@@ -83,8 +83,8 @@ static double tic( void ) {
8383
* @return random number
8484
*/
8585
static double random_uniform( const double min, const double max ) {
86-
double v = (double)rand() / ( (double)RAND_MAX + 1.0 );
87-
return min + ( v*(max-min) );
86+
double v = (double)rand() / ( (double)RAND_MAX + 1.0 );
87+
return min + ( v*(max-min) );
8888
}
8989

9090
/**
@@ -94,19 +94,26 @@ static double random_uniform( const double min, const double max ) {
9494
*/
9595
static double benchmark( void ) {
9696
double elapsed;
97-
double alpha;
98-
double beta;
99-
double x;
97+
double alpha[ 100 ];
98+
double beta[ 100 ];
99+
double x[ 100 ];
100100
double y;
101101
double t;
102102
int i;
103103

104+
for ( i = 0; i < 100; i++ ) {
105+
alpha[ i ] = random_uniform( STDLIB_CONSTANT_FLOAT64_EPS, 100.0 );
106+
}
107+
for ( i = 0; i < 100; i++ ) {
108+
beta[ i ] = random_uniform( STDLIB_CONSTANT_FLOAT64_EPS, 100.0 );
109+
}
110+
for ( i = 0; i < 100; i++ ) {
111+
x[ i ] = random_uniform( STDLIB_CONSTANT_FLOAT64_EPS, 2.0 );
112+
}
113+
104114
t = tic();
105115
for ( i = 0; i < ITERATIONS; i++ ) {
106-
x = ( ( (double)rand() / (double)RAND_MAX )*2.0 ) + STDLIB_CONSTANT_FLOAT64_EPS;
107-
alpha = ( ( (double)rand() / (double)RAND_MAX )*100.0 ) + STDLIB_CONSTANT_FLOAT64_EPS;
108-
beta = ( ( (double)rand() / (double)RAND_MAX )*100.0 ) + STDLIB_CONSTANT_FLOAT64_EPS;
109-
y = stdlib_base_dists_betaprime_logpdf( x, alpha, beta );
116+
y = stdlib_base_dists_betaprime_logpdf( x[ i%100 ], alpha[ i%100 ], beta[ i%100 ] );
110117
if ( y != y ) {
111118
printf( "should not return NaN\n" );
112119
break;

lib/node_modules/@stdlib/stats/base/dists/betaprime/logpdf/binding.gyp

Lines changed: 1 addition & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -167,4 +167,4 @@
167167
], # end actions
168168
}, # end target copy_addon
169169
], # end targets
170-
}
170+
}

lib/node_modules/@stdlib/stats/base/dists/betaprime/logpdf/examples/c/Makefile

Lines changed: 1 addition & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -143,4 +143,4 @@ run: $(c_targets)
143143
clean:
144144
$(QUIET) -rm -f *.o *.out
145145

146-
.PHONY: clean
146+
.PHONY: clean

lib/node_modules/@stdlib/stats/base/dists/betaprime/logpdf/examples/c/example.c

Lines changed: 13 additions & 8 deletions
Original file line numberDiff line numberDiff line change
@@ -21,18 +21,23 @@
2121
#include <stdlib.h>
2222
#include <stdio.h>
2323

24+
static double random_uniform( const double min, const double max ) {
25+
double v = (double)rand() / ( (double)RAND_MAX + 1.0 );
26+
return min + ( v*(max-min) );
27+
}
28+
2429
int main( void ) {
25-
double alpha;
26-
double beta;
27-
double x;
30+
double alpha;
31+
double beta;
32+
double x;
2833
double y;
2934
int i;
3035

3136
for ( i = 0; i < 10; i++ ) {
32-
x = ( (double)rand() / (double)RAND_MAX );
33-
alpha = ( ( (double)rand() / (double)RAND_MAX )*5.0 ) + STDLIB_CONSTANT_FLOAT64_EPS;
34-
beta = ( ( (double)rand() / (double)RAND_MAX )*5.0 ) + STDLIB_CONSTANT_FLOAT64_EPS;
35-
y = stdlib_base_dists_betaprime_logpdf( x, alpha, beta );
36-
printf( "x: %1f, α: %1f, β: %1f, ln(f(x;α,β)): %lf\n", x, alpha, beta, y );
37+
x = random_uniform( 0.0, 1.0 );
38+
alpha = random_uniform( STDLIB_CONSTANT_FLOAT64_EPS, 5.0 );
39+
beta = random_uniform( STDLIB_CONSTANT_FLOAT64_EPS, 5.0 );
40+
y = stdlib_base_dists_betaprime_logpdf( x, alpha, beta );
41+
printf( "x: %1f, α: %1f, β: %1f, ln(f(x;α,β)): %lf\n", x, alpha, beta, y );
3742
}
3843
}

lib/node_modules/@stdlib/stats/base/dists/betaprime/logpdf/include/stdlib/stats/base/dists/betaprime/logpdf.h

Lines changed: 1 addition & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -27,7 +27,7 @@ extern "C" {
2727
#endif
2828

2929
/**
30-
* Evaluates the probability density function (PDF) for a beta prime distribution with first shape parameter `alpha` and second shape parameter `beta` at a value `x`.
30+
* Evaluates the probability density function (PDF) for a beta prime distribution with first shape parameter `alpha` and second shape parameter `beta`.
3131
*/
3232
double stdlib_base_dists_betaprime_logpdf( const double x, const double alpha, const double beta );
3333

lib/node_modules/@stdlib/stats/base/dists/betaprime/logpdf/src/Makefile

Lines changed: 1 addition & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -67,4 +67,4 @@ clean-addon:
6767
#/
6868
clean: clean-addon
6969

70-
.PHONY: clean
70+
.PHONY: clean

0 commit comments

Comments
 (0)