fix: correct implementation

---
type: pre_commit_static_analysis_report
description: Results of running static analysis checks when committing changes.
report:
  - task: lint_filenames
    status: passed
  - task: lint_editorconfig
    status: passed
  - task: lint_markdown
    status: na
  - task: lint_package_json
    status: na
  - task: lint_repl_help
    status: na
  - task: lint_javascript_src
    status: passed
  - task: lint_javascript_cli
    status: na
  - task: lint_javascript_examples
    status: na
  - task: lint_javascript_tests
    status: na
  - task: lint_javascript_benchmarks
    status: na
  - task: lint_python
    status: na
  - task: lint_r
    status: na
  - task: lint_c_src
    status: na
  - task: lint_c_examples
    status: na
  - task: lint_c_benchmarks
    status: na
  - task: lint_c_tests_fixtures
    status: na
  - task: lint_shell
    status: na
  - task: lint_typescript_declarations
    status: na
  - task: lint_typescript_tests
    status: na
  - task: lint_license_headers
    status: passed
---

---
type: pre_push_report
description: Results of running various checks prior to pushing changes.
report:
  - task: run_javascript_examples
    status: na
  - task: run_c_examples
    status: passed
  - task: run_cpp_examples
    status: na
  - task: run_javascript_readme_examples
    status: passed
  - task: run_c_benchmarks
    status: passed
  - task: run_cpp_benchmarks
    status: na
  - task: run_fortran_benchmarks
    status: na
  - task: run_javascript_benchmarks
    status: passed
  - task: run_julia_benchmarks
    status: na
  - task: run_python_benchmarks
    status: na
  - task: run_r_benchmarks
    status: na
  - task: run_javascript_tests
    status: passed
---

Author: anandkaranubc

lib/node_modules/@stdlib/stats/base/dists/hypergeometric/skewness/README.md

@@ -184,12 +184,12 @@ double out = stdlib_base_dists_hypergeometric_skewness( 16, 11, 4 );
 
 The function accepts the following arguments:
 
--   **N**: `[in] double` population size.
--   **K**: `[in] double` subpopulation size.
--   **n**: `[in] double` number of draws.
+-   **N**: `[in] int32_t` population size.
+-   **K**: `[in] int32_t` subpopulation size.
+-   **n**: `[in] int32_t` number of draws.
 
 ```c
-double stdlib_base_dists_hypergeometric_skewness ( const double N, const double K, const double n );
+double stdlib_base_dists_hypergeometric_skewness ( const int32_t N, const int32_t K, const int32_t n );
 ```
 
 </section>
@@ -211,29 +211,30 @@ double stdlib_base_dists_hypergeometric_skewness ( const double N, const double
 ### Examples
 
 ```c
-#include "stdlib/math/base/special/round.h"
 #include "stdlib/stats/base/dists/hypergeometric/skewness.h"
-#include <stdio.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 ) {
-    double skew;
-    double N;
-    double K;
-    double n;
+    int32_t N;
+    int32_t K;
+    int32_t n;
+    double y;
     int i;
 
     for ( i = 0; i < 10; i++ ) {
-        N = stdlib_base_round( random_uniform( 0.0, 20.0 ) );
-        K = stdlib_base_round( random_uniform( 0.0, N ) );
-        n = stdlib_base_round( random_uniform( 0.0, K ) );
-        skew = stdlib_base_dists_hypergeometric_skewness( N, K, n );
-        printf( "N: %lf, K: %lf, n: %lf, skew(X;N,K,n): %lf\n", N, K, n, skew );
+        N = stdlib_base_ceil( random_uniform( 10.0, 100.0 ) );
+        K = stdlib_base_ceil( random_uniform( 0.0, N - 1 ) );
+        n = stdlib_base_ceil( random_uniform( 0.0, N - 1 ) );
+        y = stdlib_base_dists_hypergeometric_skewness( N, K, n );
+        printf( "N: %d, K: %d, n: %d, skew(X;N,K,n): %lf\n", N, K, n, y );
     }
 }
 ```

lib/node_modules/@stdlib/stats/base/dists/hypergeometric/skewness/benchmark/benchmark.js

@@ -21,9 +21,8 @@
 // MODULES //
 
 var bench = require( '@stdlib/bench' );
-var uniform = require( '@stdlib/random/base/uniform' );
-var randu = require( '@stdlib/random/base/randu' );
-var ceil = require( '@stdlib/math/base/special/ceil' );
+var Float64Array = require( '@stdlib/array/float64' );
+var discreteUniform = require( '@stdlib/random/base/discrete-uniform' );
 var isnan = require( '@stdlib/math/base/assert/is-nan' );
 var pkg = require( './../package.json' ).name;
 var skewness = require( './../lib' );
@@ -32,18 +31,26 @@ var skewness = require( './../lib' );
 // MAIN //
 
 bench( pkg, function benchmark( b ) {
+	var len;
 	var N;
 	var K;
 	var n;
 	var y;
 	var i;
 
+	len = 100;
+	N = new Float64Array( len );
+	K = new Float64Array( len );
+	n = new Float64Array( len );
+	for ( i = 0; i < len; i++ ) {
+		N[ i ] = discreteUniform( 10, 100 );
+		K[ i ] = discreteUniform( 1, N[ i ] - 1 );
+		n[ i ] = discreteUniform( 1, N[ i ] - 1 );
+	}
+
 	b.tic();
 	for ( i = 0; i < b.iterations; i++ ) {
-		N = ceil( uniform( 10.0, 100.0 ) );
-		K = ceil( randu()*(N-1) );
-		n = ceil( randu()*(N-1) );
-		y = skewness( N, K, n );
+		y = skewness( N[ i % len ], K[ i % len ], n[ i % len ] );
 		if ( isnan( y ) ) {
 			b.fail( 'should not return NaN' );
 		}

lib/node_modules/@stdlib/stats/base/dists/hypergeometric/skewness/benchmark/benchmark.native.js

@@ -23,11 +23,9 @@
 var resolve = require( 'path' ).resolve;
 var bench = require( '@stdlib/bench' );
 var Float64Array = require( '@stdlib/array/float64' );
-var uniform = require( '@stdlib/random/base/uniform' );
-var randu = require( '@stdlib/random/base/randu' );
+var discreteUniform = require( '@stdlib/random/base/discrete-uniform' );
 var isnan = require( '@stdlib/math/base/assert/is-nan' );
 var tryRequire = require( '@stdlib/utils/try-require' );
-var round = require( '@stdlib/math/base/special/round' );
 var pkg = require( './../package.json' ).name;
 
 
@@ -54,9 +52,9 @@ bench( pkg+'::native', opts, function benchmark( b ) {
 	K = new Float64Array( len );
 	n = new Float64Array( len );
 	for ( i = 0; i < len; i++ ) {
-		N[ i ] = round( uniform( 10.0, 100.0 ) );
-		K[ i ] = round( randu() * N[ i ] );
-		n[ i ] = round( randu() * N[ i ] );
+		N[ i ] = discreteUniform( 10, 100 );
+		K[ i ] = discreteUniform( 1, N[ i ] - 1 );
+		n[ i ] = discreteUniform( 1, N[ i ] - 1 );
 	}
 
 	b.tic();

lib/node_modules/@stdlib/stats/base/dists/hypergeometric/skewness/benchmark/c/benchmark.c

@@ -16,11 +16,11 @@
 * limitations under the License.
 */
 
-#include "stdlib/math/base/special/round.h"
 #include "stdlib/stats/base/dists/hypergeometric/skewness.h"
-#include <math.h>
-#include <stdio.h>
+#include "stdlib/math/base/special/ceil.h"
 #include <stdlib.h>
+#include <stdint.h>
+#include <stdio.h>
 #include <time.h>
 #include <sys/time.h>
 
@@ -94,29 +94,29 @@ static double random_uniform( const double min, const double max ) {
 */
 static double benchmark( void ) {
 	double elapsed;
-	double N[ 100 ];
-	double K[ 100 ];
-	double n[ 100 ];
-	double skew;
+	int32_t N[ 100 ];
+	int32_t K[ 100 ];
+	int32_t n[ 100 ];
+	double y;
 	double t;
 	int i;
 
 	for ( i = 0; i < 100; i++ ) {
-		N[ i ] = stdlib_base_round( random_uniform( 1.0, 100.0 ) );
-		K[ i ] = stdlib_base_round( random_uniform( 0.0, N[ i ] ) );
-		n[ i ] = stdlib_base_round( random_uniform( 0.0, N[ i ] ) );
+		N[ i ] = stdlib_base_ceil( random_uniform( 10.0, 100.0 ) );
+		K[ i ] = stdlib_base_ceil( random_uniform( 0.0, N[ i ] - 1 ) );
+		n[ i ] = stdlib_base_ceil( random_uniform( 0.0, N[ i ] - 1 ) );
 	}
 
 	t = tic();
 	for ( i = 0; i < ITERATIONS; i++ ) {
-		skew = stdlib_base_dists_hypergeometric_skewness( N[ i % 100 ], K[ i % 100 ], n[ i % 100 ] );
-		if ( skew != skew ) { // Check for NaN
+		y = stdlib_base_dists_hypergeometric_skewness( N[ i % 100 ], K[ i % 100 ], n[ i % 100 ] );
+		if ( y != y ) {
 			printf( "should not return NaN\n" );
 			break;
 		}
 	}
 	elapsed = tic() - t;
-	if ( skew != skew ) {
+	if ( y != y ) {
 		printf( "should not return NaN\n" );
 	}
 	return elapsed;

lib/node_modules/@stdlib/stats/base/dists/hypergeometric/skewness/examples/c/example.c

@@ -16,28 +16,29 @@
 * limitations under the License.
 */
 
-#include "stdlib/math/base/special/round.h"
 #include "stdlib/stats/base/dists/hypergeometric/skewness.h"
-#include <stdio.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 ) {
-	double skew;
-	double N;
-	double K;
-	double n;
+	int32_t N;
+	int32_t K;
+	int32_t n;
+	double y;
 	int i;
 
 	for ( i = 0; i < 10; i++ ) {
-		N = stdlib_base_round( random_uniform( 0.0, 20.0 ) );
-		K = stdlib_base_round( random_uniform( 0.0, N ) );
-		n = stdlib_base_round( random_uniform( 0.0, K ) );
-		skew = stdlib_base_dists_hypergeometric_skewness( N, K, n );
-		printf( "N: %lf, K: %lf, n: %lf, skew(X;N,K,n): %lf\n", N, K, n, skew );
+		N = stdlib_base_ceil( random_uniform( 10.0, 100.0 ) );
+		K = stdlib_base_ceil( random_uniform( 0.0, N - 1 ) );
+		n = stdlib_base_ceil( random_uniform( 0.0, N - 1 ) );
+		y = stdlib_base_dists_hypergeometric_skewness( N, K, n );
+		printf( "N: %d, K: %d, n: %d, skew(X;N,K,n): %lf\n", N, K, n, y );
 	}
 }

lib/node_modules/@stdlib/stats/base/dists/hypergeometric/skewness/include/stdlib/stats/base/dists/hypergeometric/skewness.h

@@ -19,6 +19,8 @@
 #ifndef STDLIB_STATS_BASE_DISTS_HYPERGEOMETRIC_SKEWNESS_H
 #define STDLIB_STATS_BASE_DISTS_HYPERGEOMETRIC_SKEWNESS_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,14 +29,9 @@ extern "C" {
 #endif
 
 /**
-* Returns the skewness of a hypergeometric distribution with population size `N`, subpopulation size `K`, and number of draws `n`.
-*
-* @param N    population size
-* @param K    subpopulation size
-* @param n    number of draws
-* @return     skewness of the distribution
+* Returns the skewness of a hypergeometric distribution.
 */
-double stdlib_base_dists_hypergeometric_skewness( const double N, const double K, const double n );
+double stdlib_base_dists_hypergeometric_skewness( const int32_t N, const int32_t K, const int32_t n );
 
 #ifdef __cplusplus
 }

lib/node_modules/@stdlib/stats/base/dists/hypergeometric/skewness/lib/native.js

@@ -26,48 +26,32 @@ var addon = require( './../src/addon.node' );
 // MAIN //
 
 /**
-* Returns the skewness of a hypergeometric distribution with population size `N`, subpopulation size `K`, and number of draws `n`.
+* Returns the skewness of a hypergeometric distribution.
 *
 * @private
-* @param {number} N - population size
-* @param {number} K - subpopulation size
-* @param {number} n - number of draws
+* @param {NonNegativeInteger} N - population size
+* @param {NonNegativeInteger} K - subpopulation size
+* @param {NonNegativeInteger} n - number of draws
 * @returns {number} skewness
 *
 * @example
-* var skew = skewness( 16, 11, 4 );
+* var v = skewness( 16, 11, 4 );
 * // returns ~-0.258
 *
 * @example
-* var skew = skewness( 4, 2, 2 );
+* var v = skewness( 4, 2, 2 );
 * // returns 0.0
 *
 * @example
-* var skew = skewness( 10, 5, 12 );
+* var v = skewness( 10, 5, 12 );
 * // returns NaN
 *
 * @example
-* var skew = skewness( 10.3, 10, 4 );
+* var v = skewness( 20, -2, 10 );
 * // returns NaN
 *
 * @example
-* var skew = skewness( 10, 5.5, 4 );
-* // returns NaN
-*
-* @example
-* var skew = skewness( 10, 5, 4.5 );
-* // returns NaN
-*
-* @example
-* var skew = skewness( NaN, 10, 4 );
-* // returns NaN
-*
-* @example
-* var skew = skewness( 20, NaN, 4 );
-* // returns NaN
-*
-* @example
-* var skew = skewness( 20, 10, NaN );
+* var v = skewness( -2, 4, 2 );
 * // returns NaN
 */
 function skewness( N, K, n ) {

lib/node_modules/@stdlib/stats/base/dists/hypergeometric/skewness/manifest.json

@@ -39,7 +39,6 @@
       "libpath": [],
       "dependencies": [
         "@stdlib/math/base/napi/ternary",
-        "@stdlib/math/base/assert/is-nonnegative-integer",
         "@stdlib/constants/float64/pinf",
         "@stdlib/math/base/special/sqrt"
       ]
@@ -56,10 +55,9 @@
       "libraries": [],
       "libpath": [],
       "dependencies": [
-        "@stdlib/math/base/assert/is-nonnegative-integer",
         "@stdlib/constants/float64/pinf",
         "@stdlib/math/base/special/sqrt",
-        "@stdlib/math/base/special/round"
+        "@stdlib/math/base/special/ceil"
       ]
     },
     {
@@ -74,10 +72,9 @@
       "libraries": [],
       "libpath": [],
       "dependencies": [
-        "@stdlib/math/base/assert/is-nonnegative-integer",
-        "@stdlib/math/base/special/round",
+        "@stdlib/math/base/special/sqrt",
         "@stdlib/constants/float64/pinf",
-        "@stdlib/math/base/special/sqrt"
+        "@stdlib/math/base/special/ceil"
       ]
     }
   ]

lib/node_modules/@stdlib/stats/base/dists/hypergeometric/skewness/src/addon.c

@@ -16,8 +16,8 @@
 * limitations under the License.
 */
 
-#include "stdlib/math/base/napi/ternary.h"
 #include "stdlib/stats/base/dists/hypergeometric/skewness.h"
+#include "stdlib/math/base/napi/ternary.h"
 
 // cppcheck-suppress shadowFunction
-STDLIB_MATH_BASE_NAPI_MODULE_DDD_D( stdlib_base_dists_hypergeometric_skewness )
+STDLIB_MATH_BASE_NAPI_MODULE_III_D( stdlib_base_dists_hypergeometric_skewness )

lib/node_modules/@stdlib/stats/base/dists/hypergeometric/skewness/src/main.c

@@ -16,10 +16,9 @@
 * limitations under the License.
 */
 
-#include "stdlib/constants/float64/pinf.h"
-#include "stdlib/math/base/assert/is_nonnegative_integer.h"
-#include "stdlib/math/base/special/sqrt.h"
 #include "stdlib/stats/base/dists/hypergeometric/skewness.h"
+#include "stdlib/math/base/special/sqrt.h"
+#include "stdlib/constants/float64/pinf.h"
 
 /**
 * Returns the skewness of a hypergeometric distribution.
@@ -33,13 +32,20 @@
 * double skew = stdlib_base_dists_hypergeometric_skewness( 16, 11, 4 );
 * // returns ~-0.258
 */
-double stdlib_base_dists_hypergeometric_skewness( const double N, const double K, const double n ) {
-	if ( !stdlib_base_is_nonnegative_integer( N ) || !stdlib_base_is_nonnegative_integer( K ) || !stdlib_base_is_nonnegative_integer( n ) || N == STDLIB_CONSTANT_FLOAT64_PINF || K == STDLIB_CONSTANT_FLOAT64_PINF || K > N || n > N ) {
-		return NAN;
-	}
-
-	double p = ( N - ( 2 * K ) ) * stdlib_base_sqrt( N - 1 ) * ( N - ( 2 * n ) );
-	double q = stdlib_base_sqrt( n * K * ( N - K ) * ( N - n ) ) * ( N - 2 );
+double stdlib_base_dists_hypergeometric_skewness( const int32_t N, const int32_t K, const int32_t n ) {
+	double N_d;
+	double K_d;
+	double n_d;
+	double p;
+	double q;
 
+	if (  N < 0 || K < 0 || n < 0 || N == STDLIB_CONSTANT_FLOAT64_PINF || K == STDLIB_CONSTANT_FLOAT64_PINF || K > N || n > N ) {
+		return 0.0/0.0; // NaN
+	}
+	N_d = (double)N;
+	K_d = (double)K;
+	n_d = (double)n;
+	p = ( N_d - ( 2.0 * K_d ) ) * stdlib_base_sqrt( N_d - 1.0 ) * ( N_d - ( 2.0 * n_d ) );
+	q = stdlib_base_sqrt( n_d * K_d * ( N_d - K_d ) * ( N_d - n_d ) ) * ( N_d - 2.0 );
 	return p / q;
 }