fix: c implementation bug fix

Krishna-Sharma-g · web-flow · commit 7fc3a2ae8b8b · 2025-02-27T13:36:23.000+05:30
diff --git a/lib/node_modules/@stdlib/stats/base/dists/kumaraswamy/median/README.md b/lib/node_modules/@stdlib/stats/base/dists/kumaraswamy/median/README.md
@@ -174,8 +174,8 @@ for ( i = 0; i < 10; i++ ) {
 Returns the [median][median] of a [Kumaraswamy's double bounded][kumaraswamy-distribution] distribution with first shape parameter `a` and second shape parameter `b`.
 
 ```c
-double out = stdlib_base_dists_kumaraswamy_median( 2.0, 3.0 );
-// returns ~1.817
+double out = stdlib_base_dists_kumaraswamy_median( 1.0, 1.0 );
+// returns 0.5
 ```
 
 The function accepts the following arguments:
@@ -212,7 +212,7 @@ double stdlib_base_dists_kumaraswamy_median( const double a, const double b );
 
 static double random_uniform( const double min, const double max ) {
     double v = (double)rand() / ( (double)RAND_MAX + 1.0 );
-    return min + ( v * (max - min) );
+    return min + ( v*(max-min) );
 }
 
 int main( void ) {
@@ -225,7 +225,7 @@ int main( void ) {
         a = random_uniform( 0.1, 10.0 );
         b = random_uniform( 0.1, 10.0 );
         y = stdlib_base_dists_kumaraswamy_median( a, b );
-        printf( "a: %lf, b: %lf, Median(a, b): %lf\n", a, b, y );
+        printf( "a: %lf, b: %lf, Median(X;a,b): %lf\n", a, b, y );
     }
 
     return 0;
diff --git a/lib/node_modules/@stdlib/stats/base/dists/kumaraswamy/median/benchmark/benchmark.js b/lib/node_modules/@stdlib/stats/base/dists/kumaraswamy/median/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/kumaraswamy/median/benchmark/benchmark.native.js b/lib/node_modules/@stdlib/stats/base/dists/kumaraswamy/median/benchmark/benchmark.native.js
@@ -23,7 +23,7 @@
 var resolve = require( 'path' ).resolve;
 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 EPS = require( '@stdlib/constants/float64/eps' );
 var isnan = require( '@stdlib/math/base/assert/is-nan' );
 var tryRequire = require( '@stdlib/utils/try-require' );
@@ -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 ) + 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/kumaraswamy/median/benchmark/c/benchmark.c b/lib/node_modules/@stdlib/stats/base/dists/kumaraswamy/median/benchmark/c/benchmark.c
@@ -1,3 +1,4 @@
+
 /**
 * @license Apache-2.0
 *
@@ -93,9 +94,10 @@ static double random_uniform( const double min, const double max ) {
 * @return elapsed time in seconds
 */
 static double benchmark( void ) {
-	double elapsed;
+
 	double a[ 100 ];
 	double b[ 100 ];
+	double elapsed;
 	double y;
 	double t;
 	int i;
diff --git a/lib/node_modules/@stdlib/stats/base/dists/kumaraswamy/median/examples/c/example.c b/lib/node_modules/@stdlib/stats/base/dists/kumaraswamy/median/examples/c/example.c
@@ -36,7 +36,7 @@ int main( void ) {
 		b = random_uniform( 0.1, 10.0 );
 
 		y = stdlib_base_dists_kumaraswamy_median( a, b );
-		printf( "a: %lf, b: %lf, Median(a,b): %lf\n", a, b, y );
+		printf( "a: %lf, b: %lf, Median(X;a,b): %lf\n", a, b, y );
 	}
 
 	return 0;
diff --git a/lib/node_modules/@stdlib/stats/base/dists/kumaraswamy/median/include/stdlib/stats/base/dists/kumaraswamy/median.h b/lib/node_modules/@stdlib/stats/base/dists/kumaraswamy/median/include/stdlib/stats/base/dists/kumaraswamy/median.h
@@ -27,7 +27,7 @@ extern "C" {
 #endif
 
 /**
-* Returns the median of a Kumaraswamy distribution.
+* Returns the median of a Kumaraswamy's double bounded distribution.
 */
 double stdlib_base_dists_kumaraswamy_median( const double a, const double b );
 
diff --git a/lib/node_modules/@stdlib/stats/base/dists/kumaraswamy/median/lib/native.js b/lib/node_modules/@stdlib/stats/base/dists/kumaraswamy/median/lib/native.js
@@ -1,3 +1,4 @@
+
 /**
 * @license Apache-2.0
 *
@@ -26,23 +27,39 @@ var addon = require( './../src/addon.node' );
 // MAIN //
 
 /**
-* Returns the median of a Kumaraswamy distribution.
+* Returns the median of a Kumaraswamy's double bounded distribution.
 *
 * @private
-* @param {number} a - shape parameter
-* @param {number} b - shape parameter
-* @returns {number} median
+* @param {PositiveNumber} a - first shape parameter
+* @param {PositiveNumber} b - second shape parameter
+* @returns {PositiveNumber} median
+*
+* @example
+* var v = median( 0.5, 1.0 );
+* // returns 0.25
+*
+* @example
+* var v = median( 4.0, 12.0 );
+* // returns ~0.487
+*
+* @example
+* var v = median( 12.0, 2.0 );
+* // returns ~0.903
+*
+* @example
+* var v = median( 1.0, -0.1 );
+* // returns NaN
 *
 * @example
-* var v = median( 2.0, 3.0 );
-* // returns ~0.2062994740159002
+* var v = median( -0.1, 1.0 );
+* // returns NaN
 *
 * @example
-* var v = median( 1.0, 5.0 );
-* // returns ~0.12944943670387588
+* var v = median( 2.0, NaN );
+* // returns NaN
 *
 * @example
-* var v = median( 2.0, -0.5 );
+* var v = median( NaN, 2.0 );
 * // returns NaN
 */
 function median( a, b ) {
diff --git a/lib/node_modules/@stdlib/stats/base/dists/kumaraswamy/median/manifest.json b/lib/node_modules/@stdlib/stats/base/dists/kumaraswamy/median/manifest.json
@@ -40,8 +40,8 @@
       "dependencies": [
         "@stdlib/math/base/napi/binary",
         "@stdlib/math/base/assert/is-nan",
-        "@stdlib/math/base/special/pow",
-        "@stdlib/constants/float64/eps"
+        "@stdlib/math/base/special/pow"
+
       ]
     },
     {
@@ -74,8 +74,8 @@
       "libpath": [],
       "dependencies": [
         "@stdlib/math/base/assert/is-nan",
-        "@stdlib/math/base/special/pow",
-        "@stdlib/constants/float64/eps"
+        "@stdlib/math/base/special/pow"
+
       ]
     }
   ]
diff --git a/lib/node_modules/@stdlib/stats/base/dists/kumaraswamy/median/src/main.c b/lib/node_modules/@stdlib/stats/base/dists/kumaraswamy/median/src/main.c
@@ -21,15 +21,15 @@
 #include "stdlib/math/base/special/pow.h"
 
 /**
-* Returns the median of a Kumaraswamy distribution.
+* Returns the median of a Kumaraswamy's double bounded distribution.
 *
 * @param a       first shape parameter
 * @param b       second shape parameter
 * @return        median
 *
 * @example
-* double y = stdlib_base_dists_kumaraswamy_median( 2.0, 3.0 );
-* // returns 1.0
+* double y = stdlib_base_dists_kumaraswamy_median( 0.5, 1.0 );
+* // returns 0.25
 */
 double stdlib_base_dists_kumaraswamy_median( const double a, const double b ) {
 	if (
diff --git a/lib/node_modules/@stdlib/stats/base/dists/kumaraswamy/median/test/fixtures/julia/runner.jl b/lib/node_modules/@stdlib/stats/base/dists/kumaraswamy/median/test/fixtures/julia/runner.jl
@@ -26,8 +26,8 @@ Generate fixture data and write to file.
 
 # Arguments
 
-* `a`: shape parameter a (must be positive)
-* `b`: shape parameter b (must be positive)
+* `a`: first shape parameter
+* `b`: second shape parameter
 * `name::AbstractString`: output filename
 
 # Examples
diff --git a/lib/node_modules/@stdlib/stats/base/dists/kumaraswamy/median/test/test.js b/lib/node_modules/@stdlib/stats/base/dists/kumaraswamy/median/test/test.js
@@ -44,10 +44,10 @@ tape( 'main export is a function', function test( t ) {
 
 tape( 'if provided `NaN` for any parameter, the function returns `NaN`', function test( t ) {
 	var v = median( NaN, 0.5 );
-	t.equal( isnan( v ), true, 'returns NaN' );
+	t.equal( isnan( v ), true, 'returns expected value' );
 
 	v = median( 10.0, NaN );
-	t.equal( isnan( v ), true, 'returns NaN' );
+	t.equal( isnan( v ), true, 'returns expected value' );
 
 	t.end();
 });
@@ -56,25 +56,25 @@ tape( 'if provided a nonpositive `a`, the function returns `NaN`', function test
 	var y;
 
 	y = median( 0.0, 2.0 );
-	t.equal( isnan( y ), true, 'returns NaN' );
+	t.equal( isnan( y ), true, 'returns expected value' );
 
 	y = median( -1.0, 2.0 );
-	t.equal( isnan( y ), true, 'returns NaN' );
+	t.equal( isnan( y ), true, 'returns expected value' );
 
 	y = median( -1.0, 2.0 );
-	t.equal( isnan( y ), true, 'returns NaN' );
+	t.equal( isnan( y ), true, 'returns expected value' );
 
 	y = median( NINF, 1.0 );
-	t.equal( isnan( y ), true, 'returns NaN' );
+	t.equal( isnan( y ), true, 'returns expected value' );
 
 	y = median( NINF, PINF );
-	t.equal( isnan( y ), true, 'returns NaN' );
+	t.equal( isnan( y ), true, 'returns expected value' );
 
 	y = median( NINF, NINF );
-	t.equal( isnan( y ), true, 'returns NaN' );
+	t.equal( isnan( y ), true, 'returns expected value' );
 
 	y = median( NINF, NaN );
-	t.equal( isnan( y ), true, 'returns NaN' );
+	t.equal( isnan( y ), true, 'returns expected value' );
 
 	t.end();
 });
@@ -83,25 +83,25 @@ tape( 'if provided a nonpositive `b`, the function returns `NaN`', function test
 	var y;
 
 	y = median( 2.0, 0.0 );
-	t.equal( isnan( y ), true, 'returns NaN' );
+	t.equal( isnan( y ), true, 'returns expected value' );
 
 	y = median( 2.0, -1.0 );
-	t.equal( isnan( y ), true, 'returns NaN' );
+	t.equal( isnan( y ), true, 'returns expected value' );
 
 	y = median( 2.0, -1/0 );
-	t.equal( isnan( y ), true, 'returns NaN' );
+	t.equal( isnan( y ), true, 'returns expected value' );
 
 	y = median( 1.0, NINF );
-	t.equal( isnan( y ), true, 'returns NaN' );
+	t.equal( isnan( y ), true, 'returns expected value' );
 
 	y = median( PINF, NINF );
-	t.equal( isnan( y ), true, 'returns NaN' );
+	t.equal( isnan( y ), true, 'returns expected value' );
 
 	y = median( NINF, NINF );
-	t.equal( isnan( y ), true, 'returns NaN' );
+	t.equal( isnan( y ), true, 'returns expected value' );
 
 	y = median( NaN, NINF );
-	t.equal( isnan( y ), true, 'returns NaN' );
+	t.equal( isnan( y ), true, 'returns expected value' );
 
 	t.end();
 });
diff --git a/lib/node_modules/@stdlib/stats/base/dists/kumaraswamy/median/test/test.native.js b/lib/node_modules/@stdlib/stats/base/dists/kumaraswamy/median/test/test.native.js
@@ -53,10 +53,10 @@ tape( 'main export is a function', opts, function test( t ) {
 
 tape( 'if provided `NaN` for any parameter, the function returns `NaN`', opts, function test( t ) {
 	var v = median( NaN, 0.5 );
-	t.equal( isnan( v ), true, 'returns NaN' );
+	t.equal( isnan( v ), true, 'returns expected value' );
 
 	v = median( 10.0, NaN );
-	t.equal( isnan( v ), true, 'returns NaN' );
+	t.equal( isnan( v ), true, 'returns expected value' );
 
 	t.end();
 });
@@ -65,25 +65,25 @@ tape( 'if provided a nonpositive `a`, the function returns `NaN`', opts, functio
 	var y;
 
 	y = median( 0.0, 2.0 );
-	t.equal( isnan( y ), true, 'returns NaN' );
+	t.equal( isnan( y ), true, 'returns expected value' );
 
 	y = median( -1.0, 2.0 );
-	t.equal( isnan( y ), true, 'returns NaN' );
+	t.equal( isnan( y ), true, 'returns expected value' );
 
 	y = median( -1.0, 2.0 );
-	t.equal( isnan( y ), true, 'returns NaN' );
+	t.equal( isnan( y ), true, 'returns expected value' );
 
 	y = median( NINF, 1.0 );
-	t.equal( isnan( y ), true, 'returns NaN' );
+	t.equal( isnan( y ), true, 'returns expected value' );
 
 	y = median( NINF, PINF );
-	t.equal( isnan( y ), true, 'returns NaN' );
+	t.equal( isnan( y ), true, 'returns expected value' );
 
 	y = median( NINF, NINF );
-	t.equal( isnan( y ), true, 'returns NaN' );
+	t.equal( isnan( y ), true, 'returns expected value' );
 
 	y = median( NINF, NaN );
-	t.equal( isnan( y ), true, 'returns NaN' );
+	t.equal( isnan( y ), true, 'returns expected value' );
 
 	t.end();
 });
@@ -92,25 +92,25 @@ tape( 'if provided a nonpositive `b`, the function returns `NaN`', opts, functio
 	var y;
 
 	y = median( 2.0, 0.0 );
-	t.equal( isnan( y ), true, 'returns NaN' );
+	t.equal( isnan( y ), true, 'returns expected value' );
 
 	y = median( 2.0, -1.0 );
-	t.equal( isnan( y ), true, 'returns NaN' );
+	t.equal( isnan( y ), true, 'returns expected value' );
 
 	y = median( 2.0, -1/0 );
-	t.equal( isnan( y ), true, 'returns NaN' );
+	t.equal( isnan( y ), true, 'returns expected value' );
 
 	y = median( 1.0, NINF );
-	t.equal( isnan( y ), true, 'returns NaN' );
+	t.equal( isnan( y ), true, 'returns expected value' );
 
 	y = median( PINF, NINF );
-	t.equal( isnan( y ), true, 'returns NaN' );
+	t.equal( isnan( y ), true, 'returns expected value' );
 
 	y = median( NINF, NINF );
-	t.equal( isnan( y ), true, 'returns NaN' );
+	t.equal( isnan( y ), true, 'returns expected value' );
 
 	y = median( NaN, NINF );
-	t.equal( isnan( y ), true, 'returns NaN' );
+	t.equal( isnan( y ), true, 'returns expected value' );
 
 	t.end();
 });

Original file line number	Diff line number	Diff line change
`@@ -36,7 +36,7 @@ int main( void ) {`
`36`	`36`	`b = random_uniform( 0.1, 10.0 );`
`37`	`37`
`38`	`38`	`y = stdlib_base_dists_kumaraswamy_median( a, b );`
`39`		`- printf( "a: %lf, b: %lf, Median(a,b): %lf\n", a, b, y );`
	`39`	`+ printf( "a: %lf, b: %lf, Median(X;a,b): %lf\n", a, b, y );`
`40`	`40`	`}`
`41`	`41`
`42`	`42`	`return 0;`
Original file line number	Diff line number	Diff line change
`@@ -40,8 +40,8 @@`
`40`	`40`	`"dependencies": [`
`41`	`41`	`"@stdlib/math/base/napi/binary",`
`42`	`42`	`"@stdlib/math/base/assert/is-nan",`
`43`		`- "@stdlib/math/base/special/pow",`
`44`		`- "@stdlib/constants/float64/eps"`
	`43`	`+ "@stdlib/math/base/special/pow"`
	`44`	`+`
`45`	`45`	`]`
`46`	`46`	`},`
`47`	`47`	`{`
`@@ -74,8 +74,8 @@`
`74`	`74`	`"libpath": [],`
`75`	`75`	`"dependencies": [`
`76`	`76`	`"@stdlib/math/base/assert/is-nan",`
`77`		`- "@stdlib/math/base/special/pow",`
`78`		`- "@stdlib/constants/float64/eps"`
	`77`	`+ "@stdlib/math/base/special/pow"`
	`78`	`+`
`79`	`79`	`]`
`80`	`80`	`}`
`81`	`81`	`]`