diff --git a/lib/node_modules/@stdlib/stats/base/dists/gamma/mgf/test/test.factory.js b/lib/node_modules/@stdlib/stats/base/dists/gamma/mgf/test/test.factory.js
index b31edb8afbe1..98a153d18631 100644
--- a/lib/node_modules/@stdlib/stats/base/dists/gamma/mgf/test/test.factory.js
+++ b/lib/node_modules/@stdlib/stats/base/dists/gamma/mgf/test/test.factory.js
@@ -22,10 +22,9 @@
 
 var tape = require( 'tape' );
 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 isAlmostSameValue = require( '@stdlib/assert/is-almost-same-value' );
 var factory = require( './../lib/factory.js' );
 
 
@@ -141,11 +140,9 @@ tape( 'if provided `alpha < 0`, the created function always returns `NaN`', func
 
 tape( 'the created function evaluates the mgf for `x` given large `alpha` and `beta`', function test( t ) {
 	var expected;
-	var delta;
 	var alpha;
 	var beta;
 	var mgf;
-	var tol;
 	var x;
 	var y;
 	var i;
@@ -157,24 +154,16 @@ tape( 'the created function evaluates the mgf for `x` given large `alpha` and `b
 	for ( i = 0; i < x.length; i++ ) {
 		mgf = factory( alpha[i], beta[i] );
 		y = mgf( x[i] );
-		if ( y === expected[i] ) {
-			t.strictEqual( y, expected[i], 'x: '+x[i]+', alpha: '+alpha[i]+', beta: '+beta[i]+', y: '+y+', expected: '+expected[i] );
-		} else {
-			delta = abs( y - expected[ i ] );
-			tol = 3000.0 * EPS * abs( expected[ i ] );
-			t.ok( delta <= tol, 'within tolerance. x: '+x[ i ]+'. alpha: '+alpha[i]+'. beta: '+beta[i]+'. y: '+y+'. E: '+expected[ i ]+'. Δ: '+delta+'. tol: '+tol+'.' );
-		}
+		t.strictEqual( isAlmostSameValue( y, expected[ i ], 3259 ), true, 'returns expected value' );
 	}
 	t.end();
 });
 
 tape( 'the created function evaluates the mgf for `x` given a large shape parameter `alpha`', function test( t ) {
 	var expected;
-	var delta;
 	var alpha;
 	var beta;
 	var mgf;
-	var tol;
 	var x;
 	var y;
 	var i;
@@ -186,24 +175,16 @@ tape( 'the created function evaluates the mgf for `x` given a large shape parame
 	for ( i = 0; i < x.length; i++ ) {
 		mgf = factory( alpha[i], beta[i] );
 		y = mgf( x[i] );
-		if ( y === expected[i] ) {
-			t.strictEqual( y, expected[i], 'x: '+x[i]+', alpha: '+alpha[i]+', beta: '+beta[i]+', y: '+y+', expected: '+expected[i] );
-		} else {
-			delta = abs( y - expected[ i ] );
-			tol = 6500.0 * EPS * abs( expected[ i ] );
-			t.ok( delta <= tol, 'within tolerance. x: '+x[ i ]+'. alpha: '+alpha[i]+'. beta: '+beta[i]+'. y: '+y+'. E: '+expected[ i ]+'. Δ: '+delta+'. tol: '+tol+'.' );
-		}
+		t.strictEqual( isAlmostSameValue( y, expected[ i ], 10882 ), true, 'returns expected value' );
 	}
 	t.end();
 });
 
 tape( 'the created function evaluates the mgf for `x` given a large rate parameter `beta`', function test( t ) {
 	var expected;
-	var delta;
 	var alpha;
 	var beta;
 	var mgf;
-	var tol;
 	var x;
 	var y;
 	var i;
@@ -215,13 +196,7 @@ tape( 'the created function evaluates the mgf for `x` given a large rate paramet
 	for ( i = 0; i < x.length; i++ ) {
 		mgf = factory( alpha[i], beta[i] );
 		y = mgf( x[i] );
-		if ( y === expected[i] ) {
-			t.strictEqual( y, expected[i], 'x: '+x[i]+', alpha: '+alpha[i]+', beta: '+beta[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 ]+'. alpha: '+alpha[i]+'. beta: '+beta[i]+'. y: '+y+'. E: '+expected[ i ]+'. Δ: '+delta+'. tol: '+tol+'.' );
-		}
+		t.strictEqual( isAlmostSameValue( y, expected[ i ], 2062 ), true, 'returns expected value' );
 	}
 	t.end();
 });
diff --git a/lib/node_modules/@stdlib/stats/base/dists/gamma/mgf/test/test.mgf.js b/lib/node_modules/@stdlib/stats/base/dists/gamma/mgf/test/test.mgf.js
index 6d6d0286138a..a54bc6148540 100644
--- a/lib/node_modules/@stdlib/stats/base/dists/gamma/mgf/test/test.mgf.js
+++ b/lib/node_modules/@stdlib/stats/base/dists/gamma/mgf/test/test.mgf.js
@@ -22,10 +22,9 @@
 
 var tape = require( 'tape' );
 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 isAlmostSameValue = require( '@stdlib/assert/is-almost-same-value' );
 var mgf = require( './../lib' );
 
 
@@ -104,10 +103,8 @@ tape( 'if provided a nonpositive `beta`, the function returns `NaN`', function t
 
 tape( 'the function evaluates the MGF for `x` given large `alpha` and `beta`', function test( t ) {
 	var expected;
-	var delta;
 	var alpha;
 	var beta;
-	var tol;
 	var x;
 	var y;
 	var i;
@@ -118,23 +115,15 @@ tape( 'the function evaluates the MGF for `x` given large `alpha` and `beta`', f
 	beta = bothLarge.beta;
 	for ( i = 0; i < x.length; i++ ) {
 		y = mgf( x[i], alpha[i], beta[i] );
-		if ( y === expected[i] ) {
-			t.strictEqual( y, expected[i], 'x: '+x[i]+', alpha: '+alpha[i]+', beta: '+beta[i]+', y: '+y+', expected: '+expected[i] );
-		} else {
-			delta = abs( y - expected[ i ] );
-			tol = 3000.0 * EPS * abs( expected[ i ] );
-			t.ok( delta <= tol, 'within tolerance. x: '+x[ i ]+'. alpha: '+alpha[i]+'. beta: '+beta[i]+'. y: '+y+'. E: '+expected[ i ]+'. Δ: '+delta+'. tol: '+tol+'.' );
-		}
+		t.strictEqual( isAlmostSameValue( y, expected[ i ], 3259 ), true, 'returns expected value' );
 	}
 	t.end();
 });
 
 tape( 'the function evaluates the MGF for `x` given large shape parameter `alpha`', function test( t ) {
 	var expected;
-	var delta;
 	var alpha;
 	var beta;
-	var tol;
 	var x;
 	var y;
 	var i;
@@ -145,23 +134,15 @@ tape( 'the function evaluates the MGF for `x` given large shape parameter `alpha
 	beta = largeShape.beta;
 	for ( i = 0; i < x.length; i++ ) {
 		y = mgf( x[i], alpha[i], beta[i] );
-		if ( y === expected[i] ) {
-			t.strictEqual( y, expected[i], 'x: '+x[i]+', alpha: '+alpha[i]+', beta: '+beta[i]+', y: '+y+', expected: '+expected[i] );
-		} else {
-			delta = abs( y - expected[ i ] );
-			tol = 6500.0 * EPS * abs( expected[ i ] );
-			t.ok( delta <= tol, 'within tolerance. x: '+x[ i ]+'. alpha: '+alpha[i]+'. beta: '+beta[i]+'. y: '+y+'. E: '+expected[ i ]+'. Δ: '+delta+'. tol: '+tol+'.' );
-		}
+		t.strictEqual( isAlmostSameValue( y, expected[ i ], 10882 ), true, 'returns expected value' );
 	}
 	t.end();
 });
 
 tape( 'the function evaluates the MGF for `x` given large rate parameter `beta`', function test( t ) {
 	var expected;
-	var delta;
 	var alpha;
 	var beta;
-	var tol;
 	var x;
 	var y;
 	var i;
@@ -172,13 +153,7 @@ tape( 'the function evaluates the MGF for `x` given large rate parameter `beta`'
 	beta = largeRate.beta;
 	for ( i = 0; i < x.length; i++ ) {
 		y = mgf( x[i], alpha[i], beta[i] );
-		if ( y === expected[i] ) {
-			t.strictEqual( y, expected[i], 'x: '+x[i]+', alpha: '+alpha[i]+', beta: '+beta[i]+', y: '+y+', expected: '+expected[i] );
-		} else {
-			delta = abs( y - expected[ i ] );
-			tol = 2200.0 * EPS * abs( expected[ i ] );
-			t.ok( delta <= tol, 'within tolerance. x: '+x[ i ]+'. alpha: '+alpha[i]+'. beta: '+beta[i]+'. y: '+y+'. E: '+expected[ i ]+'. Δ: '+delta+'. tol: '+tol+'.' );
-		}
+		t.strictEqual( isAlmostSameValue( y, expected[ i ], 2062 ), true, 'returns expected value' );
 	}
 	t.end();
 });
diff --git a/lib/node_modules/@stdlib/stats/base/dists/gamma/mgf/test/test.native.js b/lib/node_modules/@stdlib/stats/base/dists/gamma/mgf/test/test.native.js
index 04624a430e24..de716e59eacc 100644
--- a/lib/node_modules/@stdlib/stats/base/dists/gamma/mgf/test/test.native.js
+++ b/lib/node_modules/@stdlib/stats/base/dists/gamma/mgf/test/test.native.js
@@ -23,11 +23,10 @@
 var resolve = require( 'path' ).resolve;
 var tape = require( 'tape' );
 var isnan = require( '@stdlib/math/base/assert/is-nan' );
-var abs = require( '@stdlib/math/base/special/abs' );
 var tryRequire = require( '@stdlib/utils/try-require' );
 var PINF = require( '@stdlib/constants/float64/pinf' );
 var NINF = require( '@stdlib/constants/float64/ninf' );
-var EPS = require( '@stdlib/constants/float64/eps' );
+var isAlmostSameValue = require( '@stdlib/assert/is-almost-same-value' );
 
 
 // FIXTURES //
@@ -113,10 +112,8 @@ tape( 'if provided a nonpositive `beta`, the function returns `NaN`', opts, func
 
 tape( 'the function evaluates the MGF for `x` given large `alpha` and `beta`', opts, function test( t ) {
 	var expected;
-	var delta;
 	var alpha;
 	var beta;
-	var tol;
 	var x;
 	var y;
 	var i;
@@ -127,23 +124,15 @@ tape( 'the function evaluates the MGF for `x` given large `alpha` and `beta`', o
 	beta = bothLarge.beta;
 	for ( i = 0; i < x.length; i++ ) {
 		y = mgf( x[i], alpha[i], beta[i] );
-		if ( y === expected[i] ) {
-			t.strictEqual( y, expected[i], 'x: '+x[i]+', alpha: '+alpha[i]+', beta: '+beta[i]+', y: '+y+', expected: '+expected[i] );
-		} else {
-			delta = abs( y - expected[ i ] );
-			tol = 3000.0 * EPS * abs( expected[ i ] );
-			t.ok( delta <= tol, 'within tolerance. x: '+x[ i ]+'. alpha: '+alpha[i]+'. beta: '+beta[i]+'. y: '+y+'. E: '+expected[ i ]+'. Δ: '+delta+'. tol: '+tol+'.' );
-		}
+		t.strictEqual( isAlmostSameValue( y, expected[ i ], 3259 ), true, 'returns expected value' );
 	}
 	t.end();
 });
 
 tape( 'the function evaluates the MGF for `x` given large shape parameter `alpha`', opts, function test( t ) {
 	var expected;
-	var delta;
 	var alpha;
 	var beta;
-	var tol;
 	var x;
 	var y;
 	var i;
@@ -154,23 +143,15 @@ tape( 'the function evaluates the MGF for `x` given large shape parameter `alpha
 	beta = largeShape.beta;
 	for ( i = 0; i < x.length; i++ ) {
 		y = mgf( x[i], alpha[i], beta[i] );
-		if ( y === expected[i] ) {
-			t.strictEqual( y, expected[i], 'x: '+x[i]+', alpha: '+alpha[i]+', beta: '+beta[i]+', y: '+y+', expected: '+expected[i] );
-		} else {
-			delta = abs( y - expected[ i ] );
-			tol = 6500.0 * EPS * abs( expected[ i ] );
-			t.ok( delta <= tol, 'within tolerance. x: '+x[ i ]+'. alpha: '+alpha[i]+'. beta: '+beta[i]+'. y: '+y+'. E: '+expected[ i ]+'. Δ: '+delta+'. tol: '+tol+'.' );
-		}
+		t.strictEqual( isAlmostSameValue( y, expected[ i ], 10882 ), true, 'returns expected value' );
 	}
 	t.end();
 });
 
 tape( 'the function evaluates the MGF for `x` given large rate parameter `beta`', opts, function test( t ) {
 	var expected;
-	var delta;
 	var alpha;
 	var beta;
-	var tol;
 	var x;
 	var y;
 	var i;
@@ -181,13 +162,7 @@ tape( 'the function evaluates the MGF for `x` given large rate parameter `beta`'
 	beta = largeRate.beta;
 	for ( i = 0; i < x.length; i++ ) {
 		y = mgf( x[i], alpha[i], beta[i] );
-		if ( y === expected[i] ) {
-			t.strictEqual( y, expected[i], 'x: '+x[i]+', alpha: '+alpha[i]+', beta: '+beta[i]+', y: '+y+', expected: '+expected[i] );
-		} else {
-			delta = abs( y - expected[ i ] );
-			tol = 2200.0 * EPS * abs( expected[ i ] );
-			t.ok( delta <= tol, 'within tolerance. x: '+x[ i ]+'. alpha: '+alpha[i]+'. beta: '+beta[i]+'. y: '+y+'. E: '+expected[ i ]+'. Δ: '+delta+'. tol: '+tol+'.' );
-		}
+		t.strictEqual( isAlmostSameValue( y, expected[ i ], 2062 ), true, 'returns expected value' );
 	}
 	t.end();
 });