diff --git a/lib/node_modules/@stdlib/stats/base/dists/normal/mgf/test/test.factory.js b/lib/node_modules/@stdlib/stats/base/dists/normal/mgf/test/test.factory.js
index 8523f3f26130..87f6b358fef2 100644
--- a/lib/node_modules/@stdlib/stats/base/dists/normal/mgf/test/test.factory.js
+++ b/lib/node_modules/@stdlib/stats/base/dists/normal/mgf/test/test.factory.js
@@ -21,11 +21,10 @@
 // MODULES //
 
 var tape = require( 'tape' );
+var isAlmostSameValue = require( '@stdlib/assert/is-almost-same-value' );
 var isnan = require( '@stdlib/math/base/assert/is-nan' );
-var abs = require( '@stdlib/math/base/special/abs' );
 var PINF = require( '@stdlib/constants/float64/pinf' );
 var NINF = require( '@stdlib/constants/float64/ninf' );
-var EPS = require( '@stdlib/constants/float64/eps' );
 var factory = require( './../lib/factory.js' );
 
 
@@ -113,10 +112,8 @@ tape( 'if provided a nonpositive `sigma`, the created function always returns `N
 
 tape( 'the created function evaluates the MGF for `x` given positive `mu`', function test( t ) {
 	var expected;
-	var delta;
 	var sigma;
 	var mgf;
-	var tol;
 	var mu;
 	var x;
 	var y;
@@ -130,13 +127,7 @@ tape( 'the created function evaluates the MGF for `x` given positive `mu`', func
 		mgf = factory( mu[i], sigma[i] );
 		y = mgf( x[i] );
 		if ( expected[i] !== null ) {
-			if ( y === expected[i] ) {
-				t.strictEqual( y, expected[i], 'x: '+x[i]+', mu: '+mu[i]+', sigma: '+sigma[i]+', y: '+y+', expected: '+expected[i] );
-			} else {
-				delta = abs( y - expected[ i ] );
-				tol = 1200.0 * EPS * abs( expected[ i ] );
-				t.ok( delta <= tol, 'within tolerance. x: '+x[ i ]+'. mu: '+mu[i]+'. sigma: '+sigma[i]+'. y: '+y+'. E: '+expected[ i ]+'. Δ: '+delta+'. tol: '+tol+'.' );
-			}
+			t.strictEqual( isAlmostSameValue( y, expected[ i ], 1867 ), true, 'returns expected value' );
 		}
 	}
 	t.end();
@@ -144,10 +135,8 @@ tape( 'the created function evaluates the MGF for `x` given positive `mu`', func
 
 tape( 'the created function evaluates the MGF for `x` given negative `mu`', function test( t ) {
 	var expected;
-	var delta;
 	var sigma;
 	var mgf;
-	var tol;
 	var mu;
 	var x;
 	var y;
@@ -161,13 +150,7 @@ tape( 'the created function evaluates the MGF for `x` given negative `mu`', func
 		mgf = factory( mu[i], sigma[i] );
 		y = mgf( x[i] );
 		if ( expected[i] !== null ) {
-			if ( y === expected[i] ) {
-				t.strictEqual( y, expected[i], 'x: '+x[i]+', mu:'+mu[i]+', sigma: '+sigma[i]+', y: '+y+', expected: '+expected[i] );
-			} else {
-				delta = abs( y - expected[ i ] );
-				tol = 1050.0 * EPS * abs( expected[ i ] );
-				t.ok( delta <= tol, 'within tolerance. x: '+x[ i ]+'. mu: '+mu[i]+'. sigma: '+sigma[i]+'. y: '+y+'. E: '+expected[ i ]+'. Δ: '+delta+'. tol: '+tol+'.' );
-			}
+			t.strictEqual( isAlmostSameValue( y, expected[ i ], 2039 ), true, 'returns expected value' );
 		}
 	}
 	t.end();
@@ -175,10 +158,8 @@ tape( 'the created function evaluates the MGF for `x` given negative `mu`', func
 
 tape( 'the created function evaluates the MGF for `x` given large variance ( = large `sigma`)', function test( t ) {
 	var expected;
-	var delta;
 	var sigma;
 	var mgf;
-	var tol;
 	var mu;
 	var x;
 	var y;
@@ -192,13 +173,7 @@ tape( 'the created function evaluates the MGF for `x` given large variance ( = l
 		mgf = factory( mu[i], sigma[i] );
 		y = mgf( x[i] );
 		if ( expected[i] !== null ) {
-			if ( y === expected[i] ) {
-				t.strictEqual( y, expected[i], 'x: '+x[i]+', mu: '+mu[i]+', sigma: '+sigma[i]+', y: '+y+', expected: '+expected[i] );
-			} else {
-				delta = abs( y - expected[ i ] );
-				tol = 700.0 * EPS * abs( expected[ i ] );
-				t.ok( delta <= tol, 'within tolerance. x: '+x[ i ]+'. mu: '+mu[i]+'. sigma: '+sigma[i]+'. y: '+y+'. E: '+expected[ i ]+'. Δ: '+delta+'. tol: '+tol+'.' );
-			}
+			t.strictEqual( isAlmostSameValue( y, expected[ i ], 969 ), true, 'returns expected value' );
 		}
 	}
 	t.end();
diff --git a/lib/node_modules/@stdlib/stats/base/dists/normal/mgf/test/test.mgf.js b/lib/node_modules/@stdlib/stats/base/dists/normal/mgf/test/test.mgf.js
index d96b3f13c450..bd50e4a1d353 100644
--- a/lib/node_modules/@stdlib/stats/base/dists/normal/mgf/test/test.mgf.js
+++ b/lib/node_modules/@stdlib/stats/base/dists/normal/mgf/test/test.mgf.js
@@ -21,11 +21,10 @@
 // MODULES //
 
 var tape = require( 'tape' );
+var isAlmostSameValue = require( '@stdlib/assert/is-almost-same-value' );
 var isnan = require( '@stdlib/math/base/assert/is-nan' );
-var abs = require( '@stdlib/math/base/special/abs' );
 var PINF = require( '@stdlib/constants/float64/pinf' );
 var NINF = require( '@stdlib/constants/float64/ninf' );
-var EPS = require( '@stdlib/constants/float64/eps' );
 var mgf = require( './../lib' );
 
 
@@ -83,9 +82,7 @@ tape( 'if provided a nonpositive `sigma`, the function returns `NaN`', function
 
 tape( 'the function evaluates the MGF for `x` given positive `mu`', function test( t ) {
 	var expected;
-	var delta;
 	var sigma;
-	var tol;
 	var mu;
 	var x;
 	var y;
@@ -98,13 +95,7 @@ tape( 'the function evaluates the MGF for `x` given positive `mu`', function tes
 	for ( i = 0; i < x.length; i++ ) {
 		y = mgf( x[i], mu[i], sigma[i] );
 		if ( expected[i] !== null ) {
-			if ( y === expected[i] ) {
-				t.strictEqual( y, expected[i], 'x: '+x[i]+', mu:'+mu[i]+', sigma: '+sigma[i]+', y: '+y+', expected: '+expected[i] );
-			} else {
-				delta = abs( y - expected[ i ] );
-				tol = 1200.0 * EPS * abs( expected[ i ] );
-				t.ok( delta <= tol, 'within tolerance. x: '+x[ i ]+'. mu: '+mu[i]+'. sigma: '+sigma[i]+'. y: '+y+'. E: '+expected[ i ]+'. Δ: '+delta+'. tol: '+tol+'.' );
-			}
+			t.strictEqual( isAlmostSameValue( y, expected[ i ], 1867 ), true, 'returns expected value' );
 		}
 	}
 	t.end();
@@ -112,9 +103,7 @@ tape( 'the function evaluates the MGF for `x` given positive `mu`', function tes
 
 tape( 'the function evaluates the MGF for `x` given negative `mu`', function test( t ) {
 	var expected;
-	var delta;
 	var sigma;
-	var tol;
 	var mu;
 	var x;
 	var y;
@@ -127,13 +116,7 @@ tape( 'the function evaluates the MGF for `x` given negative `mu`', function tes
 	for ( i = 0; i < x.length; i++ ) {
 		y = mgf( x[i], mu[i], sigma[i] );
 		if ( expected[i] !== null ) {
-			if ( y === expected[i] ) {
-				t.strictEqual( y, expected[i], 'x: '+x[i]+', mu:'+mu[i]+', sigma: '+sigma[i]+', y: '+y+', expected: '+expected[i] );
-			} else {
-				delta = abs( y - expected[ i ] );
-				tol = 1050.0 * EPS * abs( expected[ i ] );
-				t.ok( delta <= tol, 'within tolerance. x: '+x[ i ]+'. mu: '+mu[i]+'. sigma: '+sigma[i]+'. y: '+y+'. E: '+expected[ i ]+'. Δ: '+delta+'. tol: '+tol+'.' );
-			}
+			t.strictEqual( isAlmostSameValue( y, expected[ i ], 2039 ), true, 'returns expected value' );
 		}
 	}
 	t.end();
@@ -141,9 +124,7 @@ tape( 'the function evaluates the MGF for `x` given negative `mu`', function tes
 
 tape( 'the function evaluates the MGF for `x` given large variance ( = large `sigma` )', function test( t ) {
 	var expected;
-	var delta;
 	var sigma;
-	var tol;
 	var mu;
 	var x;
 	var y;
@@ -156,13 +137,7 @@ tape( 'the function evaluates the MGF for `x` given large variance ( = large `si
 	for ( i = 0; i < x.length; i++ ) {
 		y = mgf( x[i], mu[i], sigma[i] );
 		if ( expected[i] !== null ) {
-			if ( y === expected[i] ) {
-				t.strictEqual( y, expected[i], 'x: '+x[i]+', mu:'+mu[i]+', sigma: '+sigma[i]+', y: '+y+', expected: '+expected[i] );
-			} else {
-				delta = abs( y - expected[ i ] );
-				tol = 700.0 * EPS * abs( expected[ i ] );
-				t.ok( delta <= tol, 'within tolerance. x: '+x[ i ]+'. mu: '+mu[i]+'. sigma: '+sigma[i]+'. y: '+y+'. E: '+expected[ i ]+'. Δ: '+delta+'. tol: '+tol+'.' );
-			}
+			t.strictEqual( isAlmostSameValue( y, expected[ i ], 969 ), true, 'returns expected value' );
 		}
 	}
 	t.end();
diff --git a/lib/node_modules/@stdlib/stats/base/dists/normal/mgf/test/test.native.js b/lib/node_modules/@stdlib/stats/base/dists/normal/mgf/test/test.native.js
index 602b963a6f52..4eabc85e2531 100644
--- a/lib/node_modules/@stdlib/stats/base/dists/normal/mgf/test/test.native.js
+++ b/lib/node_modules/@stdlib/stats/base/dists/normal/mgf/test/test.native.js
@@ -22,12 +22,11 @@
 
 var resolve = require( 'path' ).resolve;
 var tape = require( 'tape' );
+var isAlmostSameValue = require( '@stdlib/assert/is-almost-same-value' );
 var tryRequire = require( '@stdlib/utils/try-require' );
 var isnan = require( '@stdlib/math/base/assert/is-nan' );
-var abs = require( '@stdlib/math/base/special/abs' );
 var PINF = require( '@stdlib/constants/float64/pinf' );
 var NINF = require( '@stdlib/constants/float64/ninf' );
-var EPS = require( '@stdlib/constants/float64/eps' );
 
 
 // FIXTURES //
@@ -92,9 +91,7 @@ tape( 'if provided a nonpositive `sigma`, the function returns `NaN`', opts, fun
 
 tape( 'the function evaluates the MGF for `x` given positive `mu`', opts, function test( t ) {
 	var expected;
-	var delta;
 	var sigma;
-	var tol;
 	var mu;
 	var x;
 	var y;
@@ -107,13 +104,7 @@ tape( 'the function evaluates the MGF for `x` given positive `mu`', opts, functi
 	for ( i = 0; i < x.length; i++ ) {
 		y = mgf( x[i], mu[i], sigma[i] );
 		if ( expected[i] !== null ) {
-			if ( y === expected[i] ) {
-				t.strictEqual( y, expected[i], 'x: '+x[i]+', mu:'+mu[i]+', sigma: '+sigma[i]+', y: '+y+', expected: '+expected[i] );
-			} else {
-				delta = abs( y - expected[ i ] );
-				tol = 1200.0 * EPS * abs( expected[ i ] );
-				t.ok( delta <= tol, 'within tolerance. x: '+x[ i ]+'. mu: '+mu[i]+'. sigma: '+sigma[i]+'. y: '+y+'. E: '+expected[ i ]+'. Δ: '+delta+'. tol: '+tol+'.' );
-			}
+			t.strictEqual( isAlmostSameValue( y, expected[ i ], 1867 ), true, 'returns expected value' );
 		}
 	}
 	t.end();
@@ -121,9 +112,7 @@ tape( 'the function evaluates the MGF for `x` given positive `mu`', opts, functi
 
 tape( 'the function evaluates the MGF for `x` given negative `mu`', opts, function test( t ) {
 	var expected;
-	var delta;
 	var sigma;
-	var tol;
 	var mu;
 	var x;
 	var y;
@@ -136,13 +125,7 @@ tape( 'the function evaluates the MGF for `x` given negative `mu`', opts, functi
 	for ( i = 0; i < x.length; i++ ) {
 		y = mgf( x[i], mu[i], sigma[i] );
 		if ( expected[i] !== null ) {
-			if ( y === expected[i] ) {
-				t.strictEqual( y, expected[i], 'x: '+x[i]+', mu:'+mu[i]+', sigma: '+sigma[i]+', y: '+y+', expected: '+expected[i] );
-			} else {
-				delta = abs( y - expected[ i ] );
-				tol = 1050.0 * EPS * abs( expected[ i ] );
-				t.ok( delta <= tol, 'within tolerance. x: '+x[ i ]+'. mu: '+mu[i]+'. sigma: '+sigma[i]+'. y: '+y+'. E: '+expected[ i ]+'. Δ: '+delta+'. tol: '+tol+'.' );
-			}
+			t.strictEqual( isAlmostSameValue( y, expected[ i ], 2039 ), true, 'returns expected value' );
 		}
 	}
 	t.end();
@@ -150,9 +133,7 @@ tape( 'the function evaluates the MGF for `x` given negative `mu`', opts, functi
 
 tape( 'the function evaluates the MGF for `x` given large variance ( = large `sigma` )', opts, function test( t ) {
 	var expected;
-	var delta;
 	var sigma;
-	var tol;
 	var mu;
 	var x;
 	var y;
@@ -165,13 +146,7 @@ tape( 'the function evaluates the MGF for `x` given large variance ( = large `si
 	for ( i = 0; i < x.length; i++ ) {
 		y = mgf( x[i], mu[i], sigma[i] );
 		if ( expected[i] !== null ) {
-			if ( y === expected[i] ) {
-				t.strictEqual( y, expected[i], 'x: '+x[i]+', mu:'+mu[i]+', sigma: '+sigma[i]+', y: '+y+', expected: '+expected[i] );
-			} else {
-				delta = abs( y - expected[ i ] );
-				tol = 700.0 * EPS * abs( expected[ i ] );
-				t.ok( delta <= tol, 'within tolerance. x: '+x[ i ]+'. mu: '+mu[i]+'. sigma: '+sigma[i]+'. y: '+y+'. E: '+expected[ i ]+'. Δ: '+delta+'. tol: '+tol+'.' );
-			}
+			t.strictEqual( isAlmostSameValue( y, expected[ i ], 969 ), true, 'returns expected value' );
 		}
 	}
 	t.end();