added fixtures

Author: vivekmaurya001

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

@@ -77,20 +77,26 @@ var polyvalP = require( './polyval_p.js' );
 * @param {number} hi - upper bound
 * @param {number} lo - lower bound
 * @param {integer} k - power of 2
+* @param {integer} x - power of e
 * @returns {number} function value
 */
 function expmulti( hi, lo, k ) {
 	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 );
-
-	return ldexpf( y, k );
+    if (k >= -125) {
+	    return ldexpf( y , k );
+	} else {
+		return ldexpf( y, k+100 ) * twom100 * twom100;
+	}
 }
 
 

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

@@ -68,20 +68,23 @@ var isnan = require( '@stdlib/math/base/assert/is-nan' );
 var truncf = require( '@stdlib/math/base/special/truncf' );
 var NINF = require( '@stdlib/constants/float32/ninf' );
 var PINF = require( '@stdlib/constants/float32/pinf' );
+var absf = require( '@stdlib/math/base/special/absf' );
 var expmulti = require( './expmulti.js' );
 
 
 // VARIABLES //
 
-var LN2_HI = 0.69314718; 
-var LN2_LO = 1.90821484e-10;  
-var LOG2_E = 1.44269504; 
-var OVERFLOW = 88.3762626647949; 
-var UNDERFLOW = -103.972084045410; 
+var LN2_HI = [6.9314575195e-01, -6.9314575195e-01]; 
+var LN2_LO = [1.4286067653e-06, -1.4286067653e-06]; 
+var halF = [0.5, -0.5]; 
+var INVLN2 = 1.4426950216e+00;
+var OVERFLOW = 8.8721679688e+01; 
+var UNDERFLOW = -1.0397208405e+02;
 var NEARZERO = 1.0 / (1 << 14); 
 var NEG_NEARZERO = -NEARZERO;
 
 
+
 // MAIN //
 
 /**
@@ -165,13 +168,6 @@ var NEG_NEARZERO = -NEARZERO;
 *
 * ## Notes
 *
-* -   According to an error analysis, the error is always less than \\(1\\) ulp (unit in the last place).
-*
-* -   For an IEEE double,
-*
-*     -   if \\(x > 7.09782712893383973096\mbox{e+}02\\), then \\(e^{x}\\) overflows
-*     -   if \\(x < -7.45133219101941108420\mbox{e+}02\\), then \\(e^{x}\\) underflows
-*
 * -   The hexadecimal values included in the source code are the intended ones for the used constants. Decimal values may be used, provided that the compiler will convert from decimal to binary accurately enough to produce the intended hexadecimal values.
 *
 * @param {number} x - input value
@@ -197,6 +193,9 @@ function expf( x ) {
 	var hi;
 	var lo;
 	var k;
+	var xsb;
+
+	xsb = ( x < 0.0 ) ? 1 : 0;
 
 	if ( isnan( x ) || x === PINF ) {
 		return x;
@@ -216,14 +215,11 @@ function expf( x ) {
 	) {
 		return 1.0 + x;
 	}
-	// Reduce and compute `r = hi - lo` for extra precision...
-	if ( x < 0.0 ) {
-		k = truncf( (LOG2_E*x) - 0.5 );
-	} else {
-		k = truncf( (LOG2_E*x) + 0.5 );
-	}
-	hi = x - (k*LN2_HI);
-	lo = k * LN2_LO;
+	// Argument reduction
+      
+    k = truncf(INVLN2 * x + halF[xsb]);
+    hi = x - k * LN2_HI[0];
+    lo = k * LN2_LO[0];
 
 	return expmulti( hi, lo, k );
 }

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

@@ -35,10 +35,12 @@
 * @returns {number} evaluated polynomial
 */
 function evalpoly( x ) {
+	var P1 = 1.6666625440e-1;
+    var P2 = -2.7667332906e-3;
 	if ( x === 0.0 ) {
-		return 0.16666667; 
+		return P1; 
 	}
-	return 0.16666667 + (x * (-0.0027777778 + (x * (0.00006613757 + (x * (-0.00000165339 + (x * 4.1381368e-8))))))); // eslint-disable-line max-len
+	return P1 + x*P2 // eslint-disable-line max-len
 }
 
 // EXPORTS //

lib/node_modules/@stdlib/math/base/special/expf/test/fixtures/julia/REQUIRE

@@ -0,0 +1,2 @@
+julia 1.5
+JSON 0.21

lib/node_modules/@stdlib/math/base/special/expf/test/fixtures/julia/medium_negative_float32.json

@@ -0,0 +1,80 @@
+#!/usr/bin/env julia
+#
+# @license Apache-2.0
+#
+# Copyright (c) 2018 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.
+
+import JSON
+
+"""
+    gen( domain, filepath )
+
+Generate fixture data and write to file.
+
+# Arguments
+
+* `domain`: domain
+* `filepath::AbstractString`: filepath of the output file
+
+# Examples
+
+``` julia
+julia> x = range( -100, stop = 100, length = 2001 );
+julia> gen( x, \"./data.json\" );
+```
+"""
+function gen( domain, filepath )
+	x = collect( domain );
+	y = exp.( x );
+	data = Dict([
+		("x", x),
+		("expected", y)
+	]);
+	outfile = open( filepath, "w" );
+	write( outfile, JSON.json(data) );
+	write( outfile, "\n" );
+	close( outfile );
+end
+
+# Get the filename:
+file = @__FILE__;
+
+# Extract the directory in which this file resides:
+dir = dirname( file );
+
+# Medium negative values:
+x = Float32.(range( -88.721, stop = -0.3535, length = 1000 ));
+out = joinpath( dir, "./medium_negative_float32.json" );
+gen( x, out );
+
+# Medium positive_float32 values:
+x = Float32.(range( 0.3535, stop = 88.721, length = 1000 ));
+out = joinpath( dir, "./medium_positive_float32.json" );
+gen( x, out );
+
+# Small negative_float32 values:
+x = Float32.(range( -0.3535, stop = -2.0^-14, length = 1000 ));
+out = joinpath( dir, "./small_negative_float32.json" );
+gen( x, out );
+
+# Small positive_float32 values:
+x = Float32.(range( 2.0^-14, stop = 0.3535, length = 1000 ));
+out = joinpath( dir, "./small_positive_float32.json" );
+gen( x, out );
+
+# Tiny values:
+x = Float32.(range( -2.0^-14, stop = 2.0^-14, length = 1000 ));
+out = joinpath( dir, "./tiny_float32.json" );
+gen( x, out );