Resolved coeffiecients in evalpoly

Author: vivekmaurya001

lib/node_modules/@stdlib/math/base/special/expf/benchmark/benchmark.js

@@ -58,11 +58,11 @@ bench( pkg+'::built-in', function benchmark( b ) {
 	b.tic();
 	for ( i = 0; i < b.iterations; i++ ) {
 		x = ( randu()*100.0 ) - 50.0;
-		y = Math.expf( x ); // eslint-disable-line stdlib/no-builtin-math
+		y = Math.exp( x ); // eslint-disable-line stdlib/no-builtin-math
 		if ( isnan( y ) ) {
 			b.fail( 'should not return NaN' );
 		}
-	}+
+	}
 	b.toc();
 	if ( isnan( y ) ) {
 		b.fail( 'should not return NaN' );

lib/node_modules/@stdlib/math/base/special/expf/lib/expmulti.js

@@ -51,22 +51,23 @@ var polyvalP = require( './polyval_p.js' );
 * @returns {number} function value
 */
 function expmulti( hi, lo, k ) {
+	var twom100;
 	var r;
 	var t;
 	var c;
 	var y;
-	var twom100;
 
 	twom100 = 7.8886090522e-31;
 	r = hi - lo;
 	t = r * r;
-    c = r - ( t*polyvalP( t ) );
+	c = r - ( t*polyvalP( t ) );
 	y = 1.0 - ( lo - ( (r*c)/(2.0-c) ) - hi );
-    if (k >= -125) {
-	    return ldexpf( y , k );
-	} else {
-		return ldexpf( y, k+100 ) * twom100 * twom100;
+
+	if ( k >= -125 ) {
+		return ldexpf( y, k );
 	}
+
+	return ldexpf( y, k + 100 ) * twom100 * twom100;
 }
 
 

lib/node_modules/@stdlib/math/base/special/expf/lib/main.js

@@ -44,17 +44,16 @@ var expmulti = require( './expmulti.js' );
 
 // VARIABLES //
 
-var LN2_HI = 6.9314575195e-01; 
-var LN2_LO = 1.4286067653e-06; 
-var halF = [0.5, -0.5]; 
+var LN2_HI = 6.9314575195e-01;
+var LN2_LO = 1.4286067653e-06;
+var halF = [ 0.5, -0.5 ];
 var INVLN2 = 1.4426950216e+00;
-var OVERFLOW = 8.8721679688e+01; 
+var OVERFLOW = 8.8721679688e+01;
 var UNDERFLOW = -1.0397208405e+02;
-var NEARZERO = 1.0 / (1 << 14); 
+var NEARZERO = 1.0 / (1 << 14);
 var NEG_NEARZERO = -NEARZERO;
 
 
-
 // MAIN //
 
 /**
@@ -145,7 +144,7 @@ var NEG_NEARZERO = -NEARZERO;
 *
 * @example
 * var v = expf( 4.0 );
-* // returns ~54.5982
+* // returns ~54.598
 *
 * @example
 * var v = expf( -9.0 );
@@ -160,10 +159,10 @@ var NEG_NEARZERO = -NEARZERO;
 * // returns NaN
 */
 function expf( x ) {
+	var xsb;
 	var hi;
 	var lo;
 	var k;
-	var xsb;
 
 	xsb = ( x < 0.0 ) ? 1 : 0;
 
@@ -186,10 +185,10 @@ function expf( x ) {
 		return 1.0 + x;
 	}
 	// Argument reduction
-      
-    k = truncf(INVLN2 * x + halF[xsb]);
-    hi = x - k * LN2_HI;
-    lo = k * LN2_LO;
+
+	k = truncf( (INVLN2 * x) + halF[ xsb ] );
+	hi = x - (k * LN2_HI);
+	lo = k * LN2_LO;
 
 	return expmulti( hi, lo, k );
 }

lib/node_modules/@stdlib/math/base/special/expf/lib/polyval_p.js

@@ -1,7 +1,7 @@
 /**
 * @license Apache-2.0
 *
-* Copyright (c) 2022 The Stdlib Authors.
+* 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.
@@ -35,14 +35,13 @@
 * @returns {number} evaluated polynomial
 */
 function evalpoly( x ) {
-	var P1 = 1.6666625440e-1;
-    var P2 = -2.7667332906e-3;
 	if ( x === 0.0 ) {
-		return P1; 
+		return 0.1666662544;
 	}
-	return P1 + x*P2 // eslint-disable-line max-len
+	return 0.1666662544 + (x * -0.0027667332906);
 }
 
+
 // EXPORTS //
 
 module.exports = evalpoly;

lib/node_modules/@stdlib/math/base/special/expf/scripts/evalpoly.js

@@ -38,12 +38,9 @@ var compileC = require( '@stdlib/math/base/tools/evalpoly-compile-c' );
 // VARIABLES //
 
 // Polynomial coefficients ordered in ascending degree...
-const P = [
-    1.6666667163e-01,  // Approximation of 1/6!
-    -2.7777778450e-03, // Approximation of 1/45!
-    6.6137559770e-05,  // Approximation of 1/5040!
-    -1.6533901999e-06, // Approximation of 1/362880!
-    4.1381369442e-08   // Approximation of 1/39916800!
+var P = [
+	1.6666625440e-1,
+	-2.7667332906e-3
 ];
 
 // Header to add to output files:

lib/node_modules/@stdlib/math/base/special/expf/src/main.c

@@ -53,7 +53,7 @@ static const float NEG_NEARZERO = -1.0 / ( 1 << 14 );
 // BEGIN: polyval_p
 
 /**
-* Evaluates a polynomial for single-precision.
+* Evaluates a polynomial.
 *
 * ## Notes
 *
@@ -65,12 +65,7 @@ static const float NEG_NEARZERO = -1.0 / ( 1 << 14 );
 * @return     evaluated polynomial
 */
 static float polyval_p( const float x ) {
-	static const float P1 = 1.6666625440e-1;
-	static const float P2 = -2.7667332906e-3;
-	if ( x == 0.0 ) {
-		return P1;
-	}
-	return P1 + x * P2;
+	return 0.1666662544f + (x * -0.0027667332906f);
 }
 
 // END: polyval_p