bench: refactor C benchmark and update main.c formating

Author: Neerajpathak07

lib/node_modules/@stdlib/math/base/special/kernel-betainc/README.md

@@ -167,9 +167,9 @@ Evaluates the incomplete beta function and its first derivative.
 
 ```c
 double out;
-double deriv;
+double derivative;
 
-stdlib_base_kernel_betainc( 0.5, 2.0, 5.0, true, false, &out, &deriv );
+stdlib_base_kernel_betainc( 0.2, 1.0, 2.0, true, false, &out, &derivative );
 ```
 
 The function accepts the following arguments:
@@ -212,14 +212,13 @@ void stdlib_base_kernel_betainc( double x, double a, double b, const bool regula
 #include <stdbool.h>
 
 int main( void ) {
-
-    double out;
+    struct BasePRNGObject *obj = stdlib_base_random_randu_allocate( 0 );
     double deriv;
+    double out;
+    int32_t i;
     double x;
     double a;
     double b;
-    int32_t i;
-    struct BasePRNGObject *obj = stdlib_base_random_randu_allocate( 0 );
 
     for ( i = 0; i < 100; i++ ) {
         x = stdlib_base_random_randu( obj );

lib/node_modules/@stdlib/math/base/special/kernel-betainc/benchmark/c/native/benchmark.c

@@ -17,6 +17,7 @@
 */
 
 #include "stdlib/math/base/special/kernel_betainc.h"
+#include "stdlib/constants/float64/eps.h"
 #include <stdlib.h>
 #include <stdio.h>
 #include <stdbool.h>
@@ -77,13 +78,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 );
+static double random_uniform( const double min, const double max ) {
+	double v = (double)rand() / ( (double)RAND_MAX + 1.0 );
+	return min + ( v*(max-min) );
 }
 
 /**
@@ -93,19 +96,23 @@ static double rand_double( void ) {
 */
 static double benchmark( void ) {
 	double x[ 100 ];
+	double a[ 100 ];
+	double b[ 100 ];
 	double elapsed;
-	double out;
 	double deriv;
+	double out;
 	double t;
 	int i;
 
 	for ( i = 0; i < 100; i++ ) {
-		x[ i ] = rand_double();
+		x[ i ] = random_uniform( 0.0, 1.0 );
+		a[ i ] = random_uniform( STDLIB_FLOAT64_EPS, 100.0 );
+		b[ i ] = random_uniform( STDLIB_FLOAT64_EPS, 100.0 );
 	}
 
 	t = tic();
 	for ( i = 0; i < ITERATIONS; i++ ) {
-		stdlib_base_kernel_betainc( x[ i % 100 ], 2.0, 5.0, true, false, &out, &deriv );
+		stdlib_base_kernel_betainc( x[ i % 100 ], a[ i % 100 ], b[ i % 100 ], true, false, &out, &deriv );
 		if ( out != out || deriv != deriv ) {
 			printf( "unexpected results\n" );
 			break;

lib/node_modules/@stdlib/math/base/special/kernel-betainc/examples/c/example.c

@@ -23,14 +23,14 @@
 #include <stdbool.h>
 
 int main( void ) {
-
-	double out;
+	struct BasePRNGObject *obj = stdlib_base_random_randu_allocate( 0 );
 	double deriv;
+	double out;
 	int32_t i;
 	double x;
 	double a;
 	double b;
-	struct BasePRNGObject *obj = stdlib_base_random_randu_allocate( 0 );
+
 	for ( i = 0; i < 100; i++ ) {
 		x = stdlib_base_random_randu( obj );
 		a = stdlib_base_random_randu( obj ) * 10.0;

lib/node_modules/@stdlib/math/base/special/kernel-betainc/src/main.c

@@ -54,7 +54,6 @@
 #include <stdint.h>
 #include <stddef.h>
 
-// VARIABLES //
 static const double ONE_OVER_PI = 1.0 / STDLIB_CONSTANT_FLOAT64_PI;
 
 /**
@@ -67,6 +66,7 @@ static const double ONE_OVER_PI = 1.0 / STDLIB_CONSTANT_FLOAT64_PI;
 static double fullIGammaPrefix( const double a, const double z ) {
 	double prefix;
 	double alz;
+
 	alz = a * stdlib_base_ln( z );
 	if ( z >= 1.0 ) {
 		if ( ( alz < STDLIB_CONSTANT_FLOAT64_MAX_LN ) && ( -z > STDLIB_CONSTANT_FLOAT64_MIN_LN ) ) {
@@ -101,6 +101,7 @@ static double regularizedGammaPrefix( const double a, const double z ) {
 	double amz;
 	double sq;
 	double d;
+
 	agh = a + STDLIB_CONSTANT_FLOAT64_GAMMA_LANCZOS_G - 0.5;
 	d = ( ( z - a ) - STDLIB_CONSTANT_FLOAT64_GAMMA_LANCZOS_G + 0.5 ) / agh;
 	if ( a < 1.0 ) {
@@ -158,6 +159,7 @@ static void ibetaFraction2t( const double a, const double b, const double x, dou
 	double denom;
 	double aN;
 	double bN;
+
 	aN = ( a + ( *m ) - 1 ) * ( a + b + ( *m ) - 1 ) * ( *m ) * ( b - ( *m ) ) * x * x;
 	denom = a + ( 2.0 * ( *m ) ) - 1.0;
 	aN /= denom * denom;
@@ -177,6 +179,7 @@ static double ibetaFraction2CF( const double a, const double b, const double x,
 	double D;
 	double f;
 	double m;
+
 	m = 0.0;
 	ibetaFraction2t( a, b, x, &m, y, outTerms );
 	f = outTerms[ 1 ];
@@ -262,6 +265,7 @@ static double ibetaPowerTerms( const double a, const double b, const double x, c
 	double b2;
 	double c;
 	double l;
+
 	if ( !normalized ) {
 		// Can we do better here?
 		return stdlib_base_pow( x, a ) * stdlib_base_pow( y, b );
@@ -386,6 +390,7 @@ static double ibetaPowerTerms( const double a, const double b, const double x, c
 static double ibetaFraction2( const double a, const double b, const double x, const double y, const bool normalized, double *out ) {
 	double result;
 	double fract;
+
 	result = ibetaPowerTerms( a, b, x, y, normalized );
 	if ( out ) {
 		*out = result;
@@ -410,8 +415,8 @@ static double ibetaFraction2( const double a, const double b, const double x, co
 * @return     			function value
 */
 static double betaSmallBLargeASeries( const double a, const double b, const double x, const double y, const double s0, const double mult, const bool normalized ) {
-	double prefix;
 	double p[ 30 ];
+	double prefix;
 	double tmp1;
 	double tnp1;
 	double sum;
@@ -422,13 +427,14 @@ static double betaSmallBLargeASeries( const double a, const double b, const doub
 	double mbn;
 	double lx;
 	double t4;
-	double h;
-	double j;
 	int32_t m;
 	int32_t n;
+	double h;
+	double j;
 	double r;
 	double t;
 	double u;
+
 	// Some values we'll need later, these are Eq 9.1:
 	bm1 = b - 1.0;
 	t = a + ( bm1 / 2.0 );
@@ -511,6 +517,7 @@ static double ibetaAStep( const double a, const double b, const double x, const
 	double term;
 	double sum;
 	double i;
+
 	prefix = ibetaPowerTerms( a, b, x, y, normalized );
 	if ( out ) {
 		*out = prefix;
@@ -551,6 +558,7 @@ static double ibetaAStep( const double a, const double b, const double x, const
 static double risingFactorialRatio( const double a, const double b, const double k ) {
 	double result;
 	double i;
+
 	if ( k == 0 ) {
 		return 1.0;
 	}
@@ -572,6 +580,7 @@ static double risingFactorialRatio( const double a, const double b, const double
 */
 static double ibetaSeriesT( double *a, const double b, const double x, double *poch, int32_t *n, double *result ) {
 	double r;
+
 	r = ( *result ) / ( *a );
 	*a += 1.0;
 	*result *= ( *poch ) * x / (double)( *n );
@@ -595,6 +604,7 @@ static double binomialCCDF( const double n, const double k, const double x, cons
 	double start;
 	double term;
 	double i;
+
 	result = stdlib_base_pow( x, n );
 	if ( result > STDLIB_CONSTANT_FLOAT64_SMALLEST_NORMAL ) {
 		term = result;
@@ -646,19 +656,20 @@ static double binomialCCDF( const double n, const double k, const double x, cons
 * @return     			function value
 */
 static double ibetaSeries( const double a, const double b, const double x, const double s0, const bool normalized, double *out, const double y ) {
+	int32_t counter;
 	double result;
+	double poch;
 	double agh;
 	double bgh;
 	double cgh;
+	double sum;
 	double ac;
-	double poch;
 	int32_t n;
-	int32_t counter;
 	double l1;
 	double l2;
 	double c;
 	double r;
-	double sum;
+
 	if ( normalized ) {
 		c = a + b;
 		// Incomplete beta power term, combined with the Lanczos approximation:
@@ -721,6 +732,7 @@ static double ibetaSeries( const double a, const double b, const double x, const
 * @example
 * double out;
 * double derivative;
+*
 * stdlib_base_kernel_betainc( 0.2, 1.0, 2.0, true, false, &out, &derivative );
 */
 void stdlib_base_kernel_betainc( double x, double a, double b, const bool regularized, bool upper, double *out, double *derivative ) {
@@ -734,6 +746,7 @@ void stdlib_base_kernel_betainc( double x, double a, double b, const bool regula
 	double n;
 	double p;
 	double y;
+
 	y = 1.0 - x;
 	// Derivative not set...
 	*derivative = -1;