feat: added new implementation and other minor updates

Author: 0PrashantYadav0

lib/node_modules/@stdlib/stats/base/dists/binomial/median/README.md

@@ -169,20 +169,20 @@ for ( i = 0; i < 10; i++ ) {
 
 #### stdlib_base_dists_binomial_median( n, p )
 
-Returns the [median][median] of a [Binomial][binomial-distribution] distribution with number of trials \( n \) and success probability \( p \).
+Returns the [median][median] of a [binomial][binomial-distribution] distribution with number of trials `n` and success probability `p`.
 
 ```c
-double out = stdlib_base_dists_binomial_median( 10, 0.5 );
-// returns 5
+double out = stdlib_base_dists_binomial_median( 100, 0.1 );
+// returns 10
 ```
 
 The function accepts the following arguments:
 
--   **n**: `[in] int` number of trials.
+-   **n**: `[in] int32_t` number of trials.
 -   **p**: `[in] double` success probability.
 
 ```c
-double stdlib_base_dists_binomial_median( const int n, const double p );
+double stdlib_base_dists_binomial_median( const int32_t n, const double p );
 ```
 
 </section>
@@ -205,21 +205,30 @@ double stdlib_base_dists_binomial_median( const int n, const double p );
 
 ```c
 #include "stdlib/stats/base/dists/binomial/median.h"
+#include "stdlib/math/base/special/ceil.h"
 #include <stdlib.h>
+#include <stdint.h>
 #include <stdio.h>
 
+static double random_uniform( const double min, const double max ) {
+    double v = (double)rand() / ( (double)RAND_MAX + 1.0 );
+    return min + ( v * (max - min) );
+}
+
 int main( void ) {
-    int n;
+    int32_t n;
     double p;
     double y;
     int i;
 
     for ( i = 0; i < 25; i++ ) {
-        n = rand() % 20 + 1; // random number of trials between 1 and 20
-        p = (double)rand() / ( (double)RAND_MAX + 1.0 );
+        n = stdlib_base_ceil( random_uniform( 0.0, 100.0 ) );
+        p = random_uniform( 0.0, 1.0 );
         y = stdlib_base_dists_binomial_median( n, p );
         printf( "n: %d, p: %lf, Median(X;n,p): %lf\n", n, p, y );
     }
+
+    return 0;
 }
 ```
 

lib/node_modules/@stdlib/stats/base/dists/binomial/median/benchmark/benchmark.js

@@ -22,8 +22,9 @@
 
 var bench = require( '@stdlib/bench' );
 var ceil = require( '@stdlib/math/base/special/ceil' );
-var randu = require( '@stdlib/random/base/randu' );
+var Int32Array = require( '@stdlib/array/int32' );
 var Float64Array = require( '@stdlib/array/float64' );
+var randu = require( '@stdlib/random/base/randu' );
 var isnan = require( '@stdlib/math/base/assert/is-nan' );
 var pkg = require( './../package.json' ).name;
 var median = require( './../lib' );
@@ -39,17 +40,16 @@ bench( pkg, function benchmark( b ) {
 	var i;
 
 	len = 100;
-	n = new Float64Array(len);
-	p = new Float64Array(len);
-
-	for (i = 0; i < len; i++) {
-		n[i] = ceil( randu()*100.0 );
-		p[i] = randu();
+	n = new Int32Array( len );
+	p = new Float64Array( len );
+	for ( i = 0; i < len; i++ ) {
+		n[ i ] = ceil( randu() * 100.0 );
+		p[ i ] = randu();
 	}
 
 	b.tic();
 	for ( i = 0; i < b.iterations; i++ ) {
-		y = median(n[i % len], p[i % len]);
+		y = median( n[ i%len ], p[ i%len ] );
 		if ( isnan( y ) ) {
 			b.fail( 'should not return NaN' );
 		}

lib/node_modules/@stdlib/stats/base/dists/binomial/median/benchmark/benchmark.native.js

@@ -1,7 +1,7 @@
 /**
 * @license Apache-2.0
 *
-* Copyright (c) 2024 The Stdlib Authors.
+* Copyright (c) 2025 The Stdlib Authors.
 *
 * Licensed under the Apache License, Version 2.0 (the "License");
 * you may not use this file except in compliance with the License.
@@ -22,11 +22,12 @@
 
 var resolve = require( 'path' ).resolve;
 var bench = require( '@stdlib/bench' );
+var Int32Array = require( '@stdlib/array/int32' );
 var Float64Array = require( '@stdlib/array/float64' );
-var randint = require( '@stdlib/random/base/discrete-uniform' );
+var tryRequire = require( '@stdlib/utils/try-require' );
+var ceil = require( '@stdlib/math/base/special/ceil' );
 var randu = require( '@stdlib/random/base/randu' );
 var isnan = require( '@stdlib/math/base/assert/is-nan' );
-var tryRequire = require( '@stdlib/utils/try-require' );
 var pkg = require( './../package.json' ).name;
 
 
@@ -48,10 +49,10 @@ bench( pkg+'::native', opts, function benchmark( b ) {
 	var i;
 
 	len = 100;
-	n = new Float64Array( len );
+	n = new Int32Array( len );
 	p = new Float64Array( len );
 	for ( i = 0; i < len; i++ ) {
-		n[ i ] = randint( 1, 100 );
+		n[ i ] = ceil( randu() * 100.0 );
 		p[ i ] = randu();
 	}
 

lib/node_modules/@stdlib/stats/base/dists/binomial/median/benchmark/c/Makefile

@@ -1,7 +1,7 @@
 #/
 # @license Apache-2.0
 #
-# Copyright (c) 2024 The Stdlib Authors.
+# Copyright (c) 2025 The Stdlib Authors.
 #
 # Licensed under the Apache License, Version 2.0 (the "License");
 # you may not use this file except in compliance with the License.

lib/node_modules/@stdlib/stats/base/dists/binomial/median/benchmark/c/benchmark.c

@@ -1,7 +1,7 @@
 /**
 * @license Apache-2.0
 *
-* Copyright (c) 2024 The Stdlib Authors.
+* Copyright (c) 2025 The Stdlib Authors.
 *
 * Licensed under the Apache License, Version 2.0 (the "License");
 * you may not use this file except in compliance with the License.
@@ -17,7 +17,9 @@
 */
 
 #include "stdlib/stats/base/dists/binomial/median.h"
+#include "stdlib/math/base/special/ceil.h"
 #include <stdlib.h>
+#include <stdint.h>
 #include <stdio.h>
 #include <math.h>
 #include <time.h>
@@ -37,8 +39,8 @@ static void print_version( void ) {
 /**
 * Prints the TAP summary.
 *
-* @param total		total number of tests
-* @param passing	total number of passing tests
+* @param total    total number of tests
+* @param passing  total number of passing tests
 */
 static void print_summary( int total, int passing ) {
 	printf( "#\n" );
@@ -52,7 +54,7 @@ static void print_summary( int total, int passing ) {
 /**
 * Prints benchmarks results.
 *
-* @param elapsed   elapsed time in seconds
+* @param elapsed    elapsed time in seconds
 */
 static void print_results( double elapsed ) {
 	double rate = (double)ITERATIONS / elapsed;
@@ -75,22 +77,15 @@ static double tic( void ) {
 }
 
 /**
-* Generates a random number on the interval [0,1).
+* Generates a random number on the interval [min,max).
 *
-* @return random number
+* @param min    minimum value (inclusive)
+* @param max    maximum value (exclusive)
+* @return       random number
 */
-static double rand_double( void ) {
-	int r = rand();
-	return (double)r / ( (double)RAND_MAX + 1.0 );
-}
-
-/**
-* Generates a random integer in the range [1, 100].
-*
-* @return random integer
-*/
-static int rand_int( void ) {
-	return (rand() % 100) + 1;
+static double random_uniform( const double min, const double max ) {
+	double v = (double)rand() / ( (double)RAND_MAX + 1.0 );
+	return min + ( v*(max-min) );
 }
 
 /**
@@ -100,20 +95,20 @@ static int rand_int( void ) {
 */
 static double benchmark( void ) {
 	double elapsed;
+	int32_t n[ 100 ];
 	double p[ 100 ];
-	int n[ 100 ];
 	double y;
 	double t;
 	int i;
 
 	for ( i = 0; i < 100; i++ ) {
-		p[ i ] = rand_double();
-		n[ i ] = rand_int();
+		n[ i ] = stdlib_base_ceil( random_uniform( 0.0, 100.0 ) );
+		p[ i ] = random_uniform( 0.0, 1.0 );
 	}
 
 	t = tic();
 	for ( i = 0; i < ITERATIONS; i++ ) {
-		y = stdlib_base_dists_binomial_median( n[ i%100 ], p[ i%100 ] );
+		y = stdlib_base_dists_binomial_median( n[ i % 100 ], p[ i % 100 ] );
 		if ( y != y ) {
 			printf( "should not return NaN\n" );
 			break;

lib/node_modules/@stdlib/stats/base/dists/binomial/median/binding.gyp

@@ -1,6 +1,6 @@
 # @license Apache-2.0
 #
-# Copyright (c) 2024 The Stdlib Authors.
+# Copyright (c) 2025 The Stdlib Authors.
 #
 # Licensed under the Apache License, Version 2.0 (the "License");
 # you may not use this file except in compliance with the License.
@@ -167,4 +167,4 @@
       ], # end actions
     }, # end target copy_addon
   ], # end targets
-}
+}

lib/node_modules/@stdlib/stats/base/dists/binomial/median/examples/c/Makefile

@@ -1,7 +1,7 @@
 #/
 # @license Apache-2.0
 #
-# Copyright (c) 2024 The Stdlib Authors.
+# Copyright (c) 2025 The Stdlib Authors.
 #
 # Licensed under the Apache License, Version 2.0 (the "License");
 # you may not use this file except in compliance with the License.

lib/node_modules/@stdlib/stats/base/dists/binomial/median/examples/c/example.c

@@ -1,7 +1,7 @@
 /**
 * @license Apache-2.0
 *
-* Copyright (c) 2024 The Stdlib Authors.
+* Copyright (c) 2025 The Stdlib Authors.
 *
 * Licensed under the Apache License, Version 2.0 (the "License");
 * you may not use this file except in compliance with the License.
@@ -17,19 +17,28 @@
 */
 
 #include "stdlib/stats/base/dists/binomial/median.h"
+#include "stdlib/math/base/special/ceil.h"
 #include <stdlib.h>
+#include <stdint.h>
 #include <stdio.h>
 
+static double random_uniform( const double min, const double max ) {
+	double v = (double)rand() / ( (double)RAND_MAX + 1.0 );
+	return min + ( v * (max - min) );
+}
+
 int main( void ) {
+	int32_t n;
 	double p;
-	double n;
 	double y;
 	int i;
 
 	for ( i = 0; i < 25; i++ ) {
-		p = (double)rand() / ( (double)RAND_MAX + 1.0 );
-		n = ( rand() % 100 ) + 1;
+		n = stdlib_base_ceil( random_uniform( 0.0, 100.0 ) );
+		p = random_uniform( 0.0, 1.0 );
 		y = stdlib_base_dists_binomial_median( n, p );
-		printf( "n: %lf , p: %lf , H(X;n,p): %lf\n", n, p, y );
+		printf( "n: %d, p: %lf, Median(X;n,p): %lf\n", n, p, y );
 	}
+
+	return 0;
 }

lib/node_modules/@stdlib/stats/base/dists/binomial/median/include.gypi

@@ -1,6 +1,6 @@
 # @license Apache-2.0
 #
-# Copyright (c) 2024 The Stdlib Authors.
+# Copyright (c) 2025 The Stdlib Authors.
 #
 # Licensed under the Apache License, Version 2.0 (the "License");
 # you may not use this file except in compliance with the License.
@@ -50,4 +50,4 @@
       '<!@(node -e "var arr = require(\'@stdlib/utils/library-manifest\')(\'./manifest.json\',{},{\'basedir\':process.cwd(),\'paths\':\'posix\'}).libpath; for ( var i = 0; i < arr.length; i++ ) { console.log( arr[ i ] ); }")',
     ],
   }, # end variables
-}
+}

lib/node_modules/@stdlib/stats/base/dists/binomial/median/include/stdlib/stats/base/dists/binomial/median.h

@@ -1,7 +1,7 @@
 /**
 * @license Apache-2.0
 *
-* Copyright (c) 2024 The Stdlib Authors.
+* Copyright (c) 2025 The Stdlib Authors.
 *
 * Licensed under the Apache License, Version 2.0 (the "License");
 * you may not use this file except in compliance with the License.
@@ -19,6 +19,8 @@
 #ifndef STDLIB_STATS_BASE_DISTS_BINOMIAL_MEDIAN_H
 #define STDLIB_STATS_BASE_DISTS_BINOMIAL_MEDIAN_H
 
+#include <stdint.h>
+
 /*
 * If C++, prevent name mangling so that the compiler emits a binary file having undecorated names, thus mirroring the behavior of a C compiler.
 */
@@ -27,9 +29,9 @@ extern "C" {
 #endif
 
 /**
-* Returns the median of a binomial distribution with `n` trials and success probability `p`.
+* Returns the median of a binomial distribution with number of trials `n` and success probability `p`.
 */
-double stdlib_base_dists_binomial_median( const double n, const double p );
+double stdlib_base_dists_binomial_median( const int32_t n, const double p );
 
 #ifdef __cplusplus
 }