diff --git a/lib/node_modules/@stdlib/stats/base/dists/gumbel/quantile/test/test.factory.js b/lib/node_modules/@stdlib/stats/base/dists/gumbel/quantile/test/test.factory.js
index 5dee12fa987d..6c60afc97171 100644
--- a/lib/node_modules/@stdlib/stats/base/dists/gumbel/quantile/test/test.factory.js
+++ b/lib/node_modules/@stdlib/stats/base/dists/gumbel/quantile/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' );
 
 
@@ -133,9 +132,7 @@ tape( 'if provided a nonpositive `beta`, the created function always returns `Na
 tape( 'the created function evaluates the quantile function at `p` given positive `mu`', function test( t ) {
 	var quantile;
 	var expected;
-	var delta;
 	var beta;
-	var tol;
 	var mu;
 	var p;
 	var y;
@@ -149,13 +146,7 @@ tape( 'the created function evaluates the quantile function at `p` given positiv
 		quantile = factory( mu[i], beta[i] );
 		y = quantile( p[i] );
 		if ( expected[i] !== null ) {
-			if ( y === expected[i] ) {
-				t.strictEqual( y, expected[i], 'p: '+p[i]+', mu: '+mu[i]+', beta: '+beta[i]+', y: '+y+', expected: '+expected[i] );
-			} else {
-				delta = abs( y - expected[ i ] );
-				tol = 1.0 * EPS * abs( expected[ i ] );
-				t.ok( delta <= tol, 'within tolerance. p: '+p[ i ]+'. mu: '+mu[i]+'. beta: '+beta[i]+'. y: '+y+'. E: '+expected[ i ]+'. Δ: '+delta+'. tol: '+tol+'.' );
-			}
+			t.strictEqual( isAlmostSameValue( y, expected[ i ], 0 ), true, 'returns expected value' );
 		}
 	}
 	t.end();
@@ -164,9 +155,7 @@ tape( 'the created function evaluates the quantile function at `p` given positiv
 tape( 'the created function evaluates the quantile function at `p` given negative `mu`', function test( t ) {
 	var quantile;
 	var expected;
-	var delta;
 	var beta;
-	var tol;
 	var mu;
 	var p;
 	var y;
@@ -180,13 +169,7 @@ tape( 'the created function evaluates the quantile function at `p` given negativ
 		quantile = factory( mu[i], beta[i] );
 		y = quantile( p[i] );
 		if ( expected[i] !== null ) {
-			if ( y === expected[i] ) {
-				t.strictEqual( y, expected[i], 'p: '+p[i]+', mu:'+mu[i]+', beta: '+beta[i]+', y: '+y+', expected: '+expected[i] );
-			} else {
-				delta = abs( y - expected[ i ] );
-				tol = 1.0 * EPS * abs( expected[ i ] );
-				t.ok( delta <= tol, 'within tolerance. p: '+p[ i ]+'. mu: '+mu[i]+'. beta: '+beta[i]+'. y: '+y+'. E: '+expected[ i ]+'. Δ: '+delta+'. tol: '+tol+'.' );
-			}
+			t.strictEqual( isAlmostSameValue( y, expected[ i ], 0 ), true, 'returns expected value' );
 		}
 	}
 	t.end();
@@ -195,9 +178,7 @@ tape( 'the created function evaluates the quantile function at `p` given negativ
 tape( 'the created function evaluates the quantile function at `p` given large variance ( = large `beta`)', function test( t ) {
 	var quantile;
 	var expected;
-	var delta;
 	var beta;
-	var tol;
 	var mu;
 	var p;
 	var y;
@@ -211,13 +192,7 @@ tape( 'the created function evaluates the quantile function at `p` given large v
 		quantile = factory( mu[i], beta[i] );
 		y = quantile( p[i] );
 		if ( expected[i] !== null ) {
-			if ( y === expected[i] ) {
-				t.strictEqual( y, expected[i], 'p: '+p[i]+', mu: '+mu[i]+', beta: '+beta[i]+', y: '+y+', expected: '+expected[i] );
-			} else {
-				delta = abs( y - expected[ i ] );
-				tol = 1.0 * EPS * abs( expected[ i ] );
-				t.ok( delta <= tol, 'within tolerance. p: '+p[ i ]+'. mu: '+mu[i]+'. beta: '+beta[i]+'. y: '+y+'. E: '+expected[ i ]+'. Δ: '+delta+'. tol: '+tol+'.' );
-			}
+			t.strictEqual( isAlmostSameValue( y, expected[ i ], 0 ), true, 'returns expected value' );
 		}
 	}
 	t.end();
diff --git a/lib/node_modules/@stdlib/stats/base/dists/gumbel/quantile/test/test.native.js b/lib/node_modules/@stdlib/stats/base/dists/gumbel/quantile/test/test.native.js
index 3ce57a509a08..41fe23b019f0 100644
--- a/lib/node_modules/@stdlib/stats/base/dists/gumbel/quantile/test/test.native.js
+++ b/lib/node_modules/@stdlib/stats/base/dists/gumbel/quantile/test/test.native.js
@@ -23,11 +23,10 @@
 var resolve = require( 'path' ).resolve;
 var tape = require( 'tape' );
 var tryRequire = require( '@stdlib/utils/try-require' );
+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' );
 
 
 // VARIABLES //
@@ -103,9 +102,7 @@ tape( 'if provided a nonpositive `beta`, the function returns `NaN`', opts, func
 
 tape( 'the function evaluates the quantile for `x` given positive mean `mu`', opts, function test( t ) {
 	var expected;
-	var delta;
 	var beta;
-	var tol;
 	var mu;
 	var p;
 	var y;
@@ -118,13 +115,7 @@ tape( 'the function evaluates the quantile for `x` given positive mean `mu`', op
 	for ( i = 0; i < p.length; i++ ) {
 		y = quantile( p[i], mu[i], beta[i] );
 		if ( expected[i] !== null ) {
-			if ( y === expected[i] ) {
-				t.strictEqual( y, expected[i], 'p: '+p[i]+', mu:'+mu[i]+', beta: '+beta[i]+', y: '+y+', expected: '+expected[i] );
-			} else {
-				delta = abs( y - expected[ i ] );
-				tol = 50.0 * EPS * abs( expected[ i ] );
-				t.ok( delta <= tol, 'within tolerance. p: '+p[ i ]+'. mu: '+mu[i]+'. beta: '+beta[i]+'. y: '+y+'. E: '+expected[ i ]+'. Δ: '+delta+'. tol: '+tol+'.' );
-			}
+			t.strictEqual( isAlmostSameValue( y, expected[ i ], 0 ), true, 'returns expected value' );
 		}
 	}
 	t.end();
@@ -132,9 +123,7 @@ tape( 'the function evaluates the quantile for `x` given positive mean `mu`', op
 
 tape( 'the function evaluates the quantile for `x` given negative mean `mu`', opts, function test( t ) {
 	var expected;
-	var delta;
 	var beta;
-	var tol;
 	var mu;
 	var p;
 	var y;
@@ -147,13 +136,7 @@ tape( 'the function evaluates the quantile for `x` given negative mean `mu`', op
 	for ( i = 0; i < p.length; i++ ) {
 		y = quantile( p[i], mu[i], beta[i] );
 		if ( expected[i] !== null ) {
-			if ( y === expected[i] ) {
-				t.strictEqual( y, expected[i], 'p: '+p[i]+', mu:'+mu[i]+', beta: '+beta[i]+', y: '+y+', expected: '+expected[i] );
-			} else {
-				delta = abs( y - expected[ i ] );
-				tol = 50.0 * EPS * abs( expected[ i ] );
-				t.ok( delta <= tol, 'within tolerance. p: '+p[ i ]+'. mu: '+mu[i]+'. beta: '+beta[i]+'. y: '+y+'. E: '+expected[ i ]+'. Δ: '+delta+'. tol: '+tol+'.' );
-			}
+			t.strictEqual( isAlmostSameValue( y, expected[ i ], 0 ), true, 'returns expected value' );
 		}
 	}
 	t.end();
@@ -161,9 +144,7 @@ tape( 'the function evaluates the quantile for `x` given negative mean `mu`', op
 
 tape( 'the function evaluates the quantile function at `p` given large variance ( = large `beta` )', opts, function test( t ) {
 	var expected;
-	var delta;
 	var beta;
-	var tol;
 	var mu;
 	var p;
 	var y;
@@ -176,13 +157,7 @@ tape( 'the function evaluates the quantile function at `p` given large variance
 	for ( i = 0; i < p.length; i++ ) {
 		y = quantile( p[i], mu[i], beta[i] );
 		if ( expected[i] !== null ) {
-			if ( y === expected[i] ) {
-				t.strictEqual( y, expected[i], 'p: '+p[i]+', mu:'+mu[i]+', beta: '+beta[i]+', y: '+y+', expected: '+expected[i] );
-			} else {
-				delta = abs( y - expected[ i ] );
-				tol = 50.0 * EPS * abs( expected[ i ] );
-				t.ok( delta <= tol, 'within tolerance. p: '+p[ i ]+'. mu: '+mu[i]+'. beta: '+beta[i]+'. y: '+y+'. E: '+expected[ i ]+'. Δ: '+delta+'. tol: '+tol+'.' );
-			}
+			t.strictEqual( isAlmostSameValue( y, expected[ i ], 0 ), true, 'returns expected value' );
 		}
 	}
 	t.end();
diff --git a/lib/node_modules/@stdlib/stats/base/dists/gumbel/quantile/test/test.quantile.js b/lib/node_modules/@stdlib/stats/base/dists/gumbel/quantile/test/test.quantile.js
index eaa8380db4d9..efbdd0ca6448 100644
--- a/lib/node_modules/@stdlib/stats/base/dists/gumbel/quantile/test/test.quantile.js
+++ b/lib/node_modules/@stdlib/stats/base/dists/gumbel/quantile/test/test.quantile.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 quantile = require( './../lib' );
 
 
@@ -94,9 +93,7 @@ tape( 'if provided a nonpositive `beta`, the function returns `NaN`', function t
 
 tape( 'the function evaluates the quantile function at `p` given positive `mu`', function test( t ) {
 	var expected;
-	var delta;
 	var beta;
-	var tol;
 	var mu;
 	var p;
 	var y;
@@ -109,13 +106,7 @@ tape( 'the function evaluates the quantile function at `p` given positive `mu`',
 	for ( i = 0; i < p.length; i++ ) {
 		y = quantile( p[i], mu[i], beta[i] );
 		if ( expected[i] !== null ) {
-			if ( y === expected[i] ) {
-				t.strictEqual( y, expected[i], 'p: '+p[i]+', mu:'+mu[i]+', beta: '+beta[i]+', y: '+y+', expected: '+expected[i] );
-			} else {
-				delta = abs( y - expected[ i ] );
-				tol = 1.0 * EPS * abs( expected[ i ] );
-				t.ok( delta <= tol, 'within tolerance. p: '+p[ i ]+'. mu: '+mu[i]+'. beta: '+beta[i]+'. y: '+y+'. E: '+expected[ i ]+'. Δ: '+delta+'. tol: '+tol+'.' );
-			}
+			t.strictEqual( isAlmostSameValue( y, expected[ i ], 0 ), true, 'returns expected value' );
 		}
 	}
 	t.end();
@@ -123,9 +114,7 @@ tape( 'the function evaluates the quantile function at `p` given positive `mu`',
 
 tape( 'the function evaluates the quantile function at `p` given negative `mu`', function test( t ) {
 	var expected;
-	var delta;
 	var beta;
-	var tol;
 	var mu;
 	var p;
 	var y;
@@ -138,13 +127,7 @@ tape( 'the function evaluates the quantile function at `p` given negative `mu`',
 	for ( i = 0; i < p.length; i++ ) {
 		y = quantile( p[i], mu[i], beta[i] );
 		if ( expected[i] !== null ) {
-			if ( y === expected[i] ) {
-				t.strictEqual( y, expected[i], 'p: '+p[i]+', mu:'+mu[i]+', beta: '+beta[i]+', y: '+y+', expected: '+expected[i] );
-			} else {
-				delta = abs( y - expected[ i ] );
-				tol = 1.0 * EPS * abs( expected[ i ] );
-				t.ok( delta <= tol, 'within tolerance. p: '+p[ i ]+'. mu: '+mu[i]+'. beta: '+beta[i]+'. y: '+y+'. E: '+expected[ i ]+'. Δ: '+delta+'. tol: '+tol+'.' );
-			}
+			t.strictEqual( isAlmostSameValue( y, expected[ i ], 0 ), true, 'returns expected value' );
 		}
 	}
 	t.end();
@@ -152,9 +135,7 @@ tape( 'the function evaluates the quantile function at `p` given negative `mu`',
 
 tape( 'the function evaluates the quantile function at `p` given large variance ( = large `beta` )', function test( t ) {
 	var expected;
-	var delta;
 	var beta;
-	var tol;
 	var mu;
 	var p;
 	var y;
@@ -167,13 +148,7 @@ tape( 'the function evaluates the quantile function at `p` given large variance
 	for ( i = 0; i < p.length; i++ ) {
 		y = quantile( p[i], mu[i], beta[i] );
 		if ( expected[i] !== null ) {
-			if ( y === expected[i] ) {
-				t.strictEqual( y, expected[i], 'p: '+p[i]+', mu:'+mu[i]+', beta: '+beta[i]+', y: '+y+', expected: '+expected[i] );
-			} else {
-				delta = abs( y - expected[ i ] );
-				tol = 1.0 * EPS * abs( expected[ i ] );
-				t.ok( delta <= tol, 'within tolerance. p: '+p[ i ]+'. mu: '+mu[i]+'. beta: '+beta[i]+'. y: '+y+'. E: '+expected[ i ]+'. Δ: '+delta+'. tol: '+tol+'.' );
-			}
+			t.strictEqual( isAlmostSameValue( y, expected[ i ], 0 ), true, 'returns expected value' );
 		}
 	}
 	t.end();