diff --git a/lib/node_modules/@stdlib/stats/base/dists/gamma/cdf/test/test.cdf.js b/lib/node_modules/@stdlib/stats/base/dists/gamma/cdf/test/test.cdf.js
index fa3c0dd0f2f1..3694ecbaec68 100644
--- a/lib/node_modules/@stdlib/stats/base/dists/gamma/cdf/test/test.cdf.js
+++ b/lib/node_modules/@stdlib/stats/base/dists/gamma/cdf/test/test.cdf.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 isAlmostSameValue = require( '@stdlib/assert/is-almost-same-value' );
 var PINF = require( '@stdlib/constants/float64/pinf' );
 var NINF = require( '@stdlib/constants/float64/ninf' );
-var EPS = require( '@stdlib/constants/float64/eps' );
 var cdf = require( './../lib' );
 
 
@@ -144,9 +143,7 @@ tape( 'if provided an `alpha` equal to `0` and `beta` is positive, the function
 tape( 'the function evaluates the cdf for `x` given large `alpha` and `beta`', function test( t ) {
 	var expected;
 	var alpha;
-	var delta;
 	var beta;
-	var tol;
 	var x;
 	var y;
 	var i;
@@ -157,13 +154,7 @@ tape( 'the function evaluates the cdf for `x` given large `alpha` and `beta`', f
 	beta = bothLarge.beta;
 	for ( i = 0; i < x.length; i++ ) {
 		y = cdf( 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 = 140.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 ], 149 ), true, 'returns expected value' );
 	}
 	t.end();
 });
@@ -171,9 +162,7 @@ tape( 'the function evaluates the cdf for `x` given large `alpha` and `beta`', f
 tape( 'the function evaluates the cdf for `x` given large shape parameter `alpha`', function test( t ) {
 	var expected;
 	var alpha;
-	var delta;
 	var beta;
-	var tol;
 	var i;
 	var x;
 	var y;
@@ -184,13 +173,7 @@ tape( 'the function evaluates the cdf for `x` given large shape parameter `alpha
 	beta = largeShape.beta;
 	for ( i = 0; i < x.length; i++ ) {
 		y = cdf( 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 = 120.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 ], 131 ), true, 'returns expected value' );
 	}
 	t.end();
 });
@@ -198,9 +181,7 @@ tape( 'the function evaluates the cdf for `x` given large shape parameter `alpha
 tape( 'the function evaluates the cdf for `x` given large rate parameter `beta`', function test( t ) {
 	var expected;
 	var alpha;
-	var delta;
 	var beta;
-	var tol;
 	var i;
 	var x;
 	var y;
@@ -211,13 +192,7 @@ tape( 'the function evaluates the cdf for `x` given large rate parameter `beta`'
 	beta = largeRate.beta;
 	for ( i = 0; i < x.length; i++ ) {
 		y = cdf( 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 = 30.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 ], 32 ), true, 'returns expected value' );
 	}
 	t.end();
 });
diff --git a/lib/node_modules/@stdlib/stats/base/dists/gamma/cdf/test/test.factory.js b/lib/node_modules/@stdlib/stats/base/dists/gamma/cdf/test/test.factory.js
index 4ef7daaa761b..80e336de0b49 100644
--- a/lib/node_modules/@stdlib/stats/base/dists/gamma/cdf/test/test.factory.js
+++ b/lib/node_modules/@stdlib/stats/base/dists/gamma/cdf/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 isAlmostSameValue = require( '@stdlib/assert/is-almost-same-value' );
 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' );
 
 
@@ -199,10 +198,8 @@ tape( 'if `alpha` equals `0` and `beta` is positive, the created function evalua
 tape( 'the created function evaluates the cdf for `x` given large `alpha` and `beta`', function test( t ) {
 	var expected;
 	var alpha;
-	var delta;
 	var beta;
 	var cdf;
-	var tol;
 	var i;
 	var x;
 	var y;
@@ -214,13 +211,7 @@ tape( 'the created function evaluates the cdf for `x` given large `alpha` and `b
 	for ( i = 0; i < x.length; i++ ) {
 		cdf = factory( alpha[i], beta[i] );
 		y = cdf( 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 = 140.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 ], 149 ), true, 'returns expected value' );
 	}
 	t.end();
 });
@@ -228,10 +219,8 @@ tape( 'the created function evaluates the cdf for `x` given large `alpha` and `b
 tape( 'the created function evaluates the cdf for `x` given large shape parameter `alpha`', function test( t ) {
 	var expected;
 	var alpha;
-	var delta;
 	var beta;
 	var cdf;
-	var tol;
 	var i;
 	var x;
 	var y;
@@ -243,13 +232,7 @@ tape( 'the created function evaluates the cdf for `x` given large shape paramete
 	for ( i = 0; i < x.length; i++ ) {
 		cdf = factory( alpha[i], beta[i] );
 		y = cdf( 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 = 120.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 ], 131 ), true, 'returns expected value' );
 	}
 	t.end();
 });
@@ -257,10 +240,8 @@ tape( 'the created function evaluates the cdf for `x` given large shape paramete
 tape( 'the created function evaluates the cdf for `x` given large rate parameter `beta`', function test( t ) {
 	var expected;
 	var alpha;
-	var delta;
 	var beta;
 	var cdf;
-	var tol;
 	var i;
 	var x;
 	var y;
@@ -272,13 +253,7 @@ tape( 'the created function evaluates the cdf for `x` given large rate parameter
 	for ( i = 0; i < x.length; i++ ) {
 		cdf = factory( alpha[i], beta[i] );
 		y = cdf( 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 = 30.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 ], 32 ), true, 'returns expected value' );
 	}
 	t.end();
 });
diff --git a/lib/node_modules/@stdlib/stats/base/dists/gamma/cdf/test/test.native.js b/lib/node_modules/@stdlib/stats/base/dists/gamma/cdf/test/test.native.js
index b308101e7d3e..52df92df96d1 100644
--- a/lib/node_modules/@stdlib/stats/base/dists/gamma/cdf/test/test.native.js
+++ b/lib/node_modules/@stdlib/stats/base/dists/gamma/cdf/test/test.native.js
@@ -22,12 +22,11 @@
 
 var resolve = require( 'path' ).resolve;
 var tape = require( 'tape' );
-var abs = require( '@stdlib/math/base/special/abs' );
+var isAlmostSameValue = require( '@stdlib/assert/is-almost-same-value' );
 var isnan = require( '@stdlib/math/base/assert/is-nan' );
 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' );
 
 
 // FIXTURES //
@@ -121,10 +120,8 @@ tape( 'if provided a nonpositive `beta`, the function returns `NaN`', opts, func
 
 tape( 'the function evaluates the cdf for `x` given large `alpha`', opts, function test( t ) {
 	var expected;
-	var delta;
 	var alpha;
 	var beta;
-	var tol;
 	var x;
 	var y;
 	var i;
@@ -135,23 +132,15 @@ tape( 'the function evaluates the cdf for `x` given large `alpha`', opts, functi
 	beta = largeShape.beta;
 	for ( i = 0; i < x.length; i++ ) {
 		y = cdf( 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 = 120.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 ], 131 ), true, 'returns expected value' );
 	}
 	t.end();
 });
 
 tape( 'the function evaluates the cdf for `x` given large `beta`', opts, function test( t ) {
 	var expected;
-	var delta;
 	var alpha;
 	var beta;
-	var tol;
 	var x;
 	var y;
 	var i;
@@ -162,23 +151,15 @@ tape( 'the function evaluates the cdf for `x` given large `beta`', opts, functio
 	beta = largeRate.beta;
 	for ( i = 0; i < x.length; i++ ) {
 		y = cdf( 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 = 30.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 ], 32 ), true, 'returns expected value' );
 	}
 	t.end();
 });
 
 tape( 'the function evaluates the cdf 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;
@@ -189,13 +170,7 @@ tape( 'the function evaluates the cdf for `x` given large `alpha` and `beta`', o
 	beta = bothLarge.beta;
 	for ( i = 0; i < x.length; i++ ) {
 		y = cdf( 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 = 140.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 ], 149 ), true, 'returns expected value' );
 	}
 	t.end();
 });