feat(add c implementation): add stats/base/dists/bernoulli/entropy

Author: vivekmaurya001

lib/node_modules/@stdlib/stats/base/dists/bernoulli/entropy/README.md

@@ -167,7 +167,7 @@ double out = stdlib_base_dists_bernoulli_entropy( 0.1 );
 
 The function accepts the following arguments:
 
--   **p**: `[in] double` maximum support
+-   **p**: `[in] double` success probability
 
 ```c
 double stdlib_base_dists_bernoulli_entropy( const double p );
@@ -201,8 +201,8 @@ int main( void ) {
     int i;
 
     for ( i = 0; i < 25; i++ ) {
-        p = ( (double)rand() / (double)RAND_MAX );
-        y = stdlib_base_dists_arcsine_entropy( p );
+        p = (double)rand() / ( (double)RAND_MAX + 1.0 );
+        y = stdlib_base_dists_bernoulli_entropy( p );
         printf( "x: %lf, H(X;p): %lf\n", p, y );
     }
 }

lib/node_modules/@stdlib/stats/base/dists/bernoulli/entropy/benchmark/benchmark.js

@@ -21,6 +21,7 @@
 // MODULES //
 
 var bench = require( '@stdlib/bench' );
+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;
@@ -30,14 +31,20 @@ var entropy = require( './../lib' );
 // MAIN //
 
 bench( pkg, function benchmark( b ) {
+	var len;
 	var p;
 	var y;
 	var i;
 
+	len = 100;
+	p = new Float64Array( len );
+	for ( i = 0; i < len; i++ ) {
+		p[ i ] = randu();
+	}
+
 	b.tic();
 	for ( i = 0; i < b.iterations; i++ ) {
-		p = randu();
-		y = entropy( p );
+		y = entropy( p[ i % len ] );
 		if ( isnan( y ) ) {
 			b.fail( 'should not return NaN' );
 		}

lib/node_modules/@stdlib/stats/base/dists/bernoulli/entropy/benchmark/benchmark.native.js

@@ -0,0 +1,67 @@
+/**
+* @license Apache-2.0
+*
+* Copyright (c) 2024 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.
+* You may obtain a copy of the License at
+*
+*    http://www.apache.org/licenses/LICENSE-2.0
+*
+* Unless required by applicable law or agreed to in writing, software
+* distributed under the License is distributed on an "AS IS" BASIS,
+* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
+* See the License for the specific language governing permissions and
+* limitations under the License.
+*/
+
+'use strict';
+
+// MODULES //
+
+var resolve = require( 'path' ).resolve;
+var bench = require( '@stdlib/bench' );
+var Float64Array = require( '@stdlib/array/float64' );
+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;
+
+
+// VARIABLES //
+
+var entropy = tryRequire( resolve( __dirname, './../lib/native.js' ) );
+var opts = {
+	'skip': ( entropy instanceof Error )
+};
+
+
+// MAIN //
+
+bench( pkg, opts, function benchmark( b ) {
+	var len;
+	var p;
+	var y;
+	var i;
+
+	len = 100;
+	p = new Float64Array( len );
+	for ( i = 0; i < len; i++ ) {
+		p[ i ] = randu();
+	}
+
+	b.tic();
+	for ( i = 0; i < b.iterations; i++ ) {
+		y = entropy( p[ i % len ] );
+		if ( isnan( y ) ) {
+			b.fail( 'should not return NaN' );
+		}
+	}
+	b.toc();
+	if ( isnan( y ) ) {
+		b.fail( 'should not return NaN' );
+	}
+	b.pass( 'benchmark finished' );
+	b.end();
+});

lib/node_modules/@stdlib/stats/base/dists/bernoulli/entropy/benchmark/c/benchmark.c

@@ -16,12 +16,12 @@
 * limitations under the License.
 */
 
-#include <sys/time.h>
 #include "stdlib/stats/base/dists/bernoulli/entropy.h"
 #include <math.h>
 #include <stdio.h>
 #include <stdlib.h>
 #include <time.h>
+#include <sys/time.h>
 
 #define NAME "bernoulli-entropy"
 #define ITERATIONS 1000000
@@ -91,15 +91,18 @@ static double rand_double( void ) {
 */
 static double benchmark( void ) {
 	double elapsed;
-	double p;
+	double p[ 100 ];
 	double y;
 	double t;
 	int i;
 
+	for ( i = 0; i < 100; i++ ) {
+		p[ i ] = rand_double();
+	}
+
 	t = tic();
 	for ( i = 0; i < ITERATIONS; i++ ) {
-		p = ( (double)rand() / (double)RAND_MAX );
-		y = stdlib_base_dists_bernoulli_entropy( p );
+		y = stdlib_base_dists_bernoulli_entropy( p[ i%100 ] );
 		if ( y != y ) {
 			printf( "should not return NaN\n" );
 			break;

lib/node_modules/@stdlib/stats/base/dists/bernoulli/entropy/examples/c/Makefile

@@ -143,4 +143,4 @@ run: $(c_targets)
 clean:
 	$(QUIET) -rm -f *.o *.out
 
-.PHONY: clean
+.PHONY: clean

lib/node_modules/@stdlib/stats/base/dists/bernoulli/entropy/examples/c/example.c

@@ -26,7 +26,7 @@ int main( void ) {
 	int i;
 
 	for ( i = 0; i < 25; i++ ) {
-		p = ( (double)rand() / (double)RAND_MAX );
+		p = (double)rand() / ( (double)RAND_MAX + 1.0 );
 		y = stdlib_base_dists_bernoulli_entropy( p );
 		printf( "x: %lf , H(X;p): %lf\n", p, y );
 	}

lib/node_modules/@stdlib/stats/base/dists/bernoulli/entropy/include.gypi

@@ -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/bernoulli/entropy/lib/native.js

@@ -58,4 +58,4 @@ function entropy( p ) {
 
 // EXPORTS //
 
-module.exports = entropy;
+module.exports = entropy;

lib/node_modules/@stdlib/stats/base/dists/bernoulli/entropy/manifest.json

@@ -38,14 +38,9 @@
         "libraries": [],
         "libpath": [],
         "dependencies": [
-          "@stdlib/math/base/napi/ternary",
+          "@stdlib/math/base/napi/unary",
           "@stdlib/math/base/assert/is-nan",
-          "@stdlib/math/base/special/ln",
-          "@stdlib/constants/float64/ln-pi",
-          "@stdlib/constants/float64/ln-two",
-          "@stdlib/constants/float64/ninf",
-          "@stdlib/math/base/special/asin",
-          "@stdlib/math/base/special/sqrt"
+          "@stdlib/math/base/special/ln"
         ]
       },
       {
@@ -61,12 +56,7 @@
         "libpath": [],
         "dependencies": [
           "@stdlib/math/base/assert/is-nan",
-          "@stdlib/math/base/special/ln",
-          "@stdlib/constants/float64/ln-pi",
-          "@stdlib/constants/float64/ln-two",
-          "@stdlib/constants/float64/ninf",
-          "@stdlib/math/base/special/asin",
-          "@stdlib/math/base/special/sqrt"
+          "@stdlib/math/base/special/ln"
         ]
       },
       {
@@ -82,12 +72,7 @@
         "libpath": [],
         "dependencies": [
           "@stdlib/math/base/assert/is-nan",
-          "@stdlib/math/base/special/ln",
-          "@stdlib/constants/float64/ln-pi",
-          "@stdlib/constants/float64/ln-two",
-          "@stdlib/constants/float64/ninf",
-          "@stdlib/math/base/special/asin",
-          "@stdlib/math/base/special/sqrt"
+          "@stdlib/math/base/special/ln"
         ]
       }
     ]

lib/node_modules/@stdlib/stats/base/dists/bernoulli/entropy/src/Makefile

@@ -67,4 +67,4 @@ clean-addon:
 #/
 clean: clean-addon
 
-.PHONY: clean
+.PHONY: clean