feat: updated benchmarks with precomputation

Author: aayush0325

lib/node_modules/@stdlib/math/base/special/cinvf/benchmark/c/benchmark.c

@@ -90,20 +90,22 @@ static float rand_float( void ) {
 */
 static double benchmark( void ) {
 	double elapsed;
-	float re;
-	float im;
+	float x[ 100 ];
+	float y[ 100 ];
 	double t;
 	int i;
 
 	float complex z1;
 	float complex z2;
 
+	for ( i = 0; i < 100; i++ ) {
+		x[ i ] = ( 20.0f * rand_float() ) - 10.0f;
+		y[ i ] = ( 2048.0f * rand_float() ) - 1024.0f;
+	}
+
 	t = tic();
 	for ( i = 0; i < ITERATIONS; i++ ) {
-		re = ( 1000.0f * rand_float() ) - 500.0f;
-		im = ( 1000.0f * rand_float() ) - 500.0f;
-		z1 = re + im*I;
-
+		z1 = x[ i % 100 ] + y[ i % 100 ]*I;
 		z2 = 1.0f / z1;
 		if ( z2 != z2 ) {
 			printf( "should not return NaN\n" );

lib/node_modules/@stdlib/math/base/special/cinvf/benchmark/c/native/benchmark.c

@@ -96,18 +96,23 @@ static double benchmark( void ) {
 	float re;
 	float im;
 	double t;
-	float v;
+	float x[ 100 ];
+	float y[ 100 ];
 	int i;
 
-	stdlib_complex64_t x;
-	stdlib_complex64_t y;
+	stdlib_complex64_t a;
+	stdlib_complex64_t b;
+
+	for ( i = 0; i < 100; i++ ) {
+		x[ i ] = ( 20.0f * rand_float() ) - 10.0f;
+		y[ i ] = ( 2048.0f * rand_float() ) - 1024.0f;
+	}
 
 	t = tic();
 	for ( i = 0; i < ITERATIONS; i++ ) {
-		v = ( 1000.0f * rand_float() ) - 500.0f;
-		x = stdlib_complex64( v, v );
-		y = stdlib_base_cinvf( x );
-		stdlib_complex64_reim( y, &re, &im );
+		a = stdlib_complex64( x[ i % 100 ], y[ i % 100 ] );
+		b = stdlib_base_cinvf( a );
+		stdlib_complex64_reim( b, &re, &im );
 		if ( re != re ) {
 			printf( "unexpected result\n" );
 			break;