diff --git a/lib/node_modules/@stdlib/stats/base/dists/invgamma/quantile/test/test.factory.js b/lib/node_modules/@stdlib/stats/base/dists/invgamma/quantile/test/test.factory.js
index c4527ac6a72f..75e831a5e477 100644
--- a/lib/node_modules/@stdlib/stats/base/dists/invgamma/quantile/test/test.factory.js
+++ b/lib/node_modules/@stdlib/stats/base/dists/invgamma/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' );
 
 
@@ -156,10 +155,8 @@ tape( 'if provided a negative `alpha`, the created function always returns `NaN`
 tape( 'the created function evaluates the quantile for `p` given large `alpha` and `beta`', function test( t ) {
 	var expected;
 	var quantile;
-	var delta;
 	var alpha;
 	var beta;
-	var tol;
 	var i;
 	var p;
 	var y;
@@ -171,13 +168,7 @@ tape( 'the created function evaluates the quantile for `p` given large `alpha` a
 	for ( i = 0; i < p.length; i++ ) {
 		quantile = factory( alpha[ i ], beta[ i ] );
 		y = quantile( p[ i ] );
-		if ( y === expected[ i ] ) {
-			t.strictEqual( y, expected[ i ], 'p: '+p[ i ]+', alpha:'+alpha[ i ]+', beta: '+beta[ i ]+', y: '+y+', expected: '+expected[ i ] );
-		} else {
-			delta = abs( y - expected[ i ] );
-			tol = 200.0 * EPS * abs( expected[ i ] );
-			t.ok( delta <= tol, 'within tolerance. p: '+p[ i ]+'. alpha: '+alpha[ i ]+'. beta: '+beta[ i ]+'. y: '+y+'. E: '+expected[ i ]+'. Δ: '+delta+'. tol: '+tol+'.' );
-		}
+		t.strictEqual( isAlmostSameValue( y, expected[ i ], 23 ), true, 'returns expected value' );
 	}
 	t.end();
 });
@@ -185,10 +176,8 @@ tape( 'the created function evaluates the quantile for `p` given large `alpha` a
 tape( 'the created function evaluates the quantile for `p` given large shape parameter `alpha`', function test( t ) {
 	var expected;
 	var quantile;
-	var delta;
 	var alpha;
 	var beta;
-	var tol;
 	var i;
 	var p;
 	var y;
@@ -200,13 +189,7 @@ tape( 'the created function evaluates the quantile for `p` given large shape par
 	for ( i = 0; i < p.length; i++ ) {
 		quantile = factory( alpha[ i ], beta[ i ] );
 		y = quantile( p[ i ] );
-		if ( y === expected[ i ] ) {
-			t.strictEqual( y, expected[ i ], 'p: '+p[ i ]+', alpha:'+alpha[ i ]+', beta: '+beta[ i ]+', y: '+y+', expected: '+expected[ i ] );
-		} else {
-			delta = abs( y - expected[ i ] );
-			tol = 100.0 * EPS * abs( expected[ i ] );
-			t.ok( delta <= tol, 'within tolerance. p: '+p[ i ]+'. alpha: '+alpha[ i ]+'. beta: '+beta[ i ]+'. y: '+y+'. E: '+expected[ i ]+'. Δ: '+delta+'. tol: '+tol+'.' );
-		}
+		t.strictEqual( isAlmostSameValue( y, expected[ i ], 23 ), true, 'returns expected value' );
 	}
 	t.end();
 });
@@ -214,10 +197,8 @@ tape( 'the created function evaluates the quantile for `p` given large shape par
 tape( 'the created function evaluates the quantile for `p` given large scale parameter `beta`', function test( t ) {
 	var expected;
 	var quantile;
-	var delta;
 	var alpha;
 	var beta;
-	var tol;
 	var i;
 	var p;
 	var y;
@@ -229,13 +210,7 @@ tape( 'the created function evaluates the quantile for `p` given large scale par
 	for ( i = 0; i < p.length; i++ ) {
 		quantile = factory( alpha[ i ], beta[ i ] );
 		y = quantile( p[ i ] );
-		if ( y === expected[ i ] ) {
-			t.strictEqual( y, expected[ i ], 'p: '+p[ i ]+', alpha:'+alpha[ 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 ]+'. alpha: '+alpha[ i ]+'. beta: '+beta[ i ]+'. y: '+y+'. E: '+expected[ i ]+'. Δ: '+delta+'. tol: '+tol+'.' );
-		}
+		t.strictEqual( isAlmostSameValue( y, expected[ i ], 81 ), true, 'returns expected value' );
 	}
 	t.end();
 });
diff --git a/lib/node_modules/@stdlib/stats/base/dists/invgamma/quantile/test/test.native.js b/lib/node_modules/@stdlib/stats/base/dists/invgamma/quantile/test/test.native.js
index da74540e3eb5..1f1d408677a5 100644
--- a/lib/node_modules/@stdlib/stats/base/dists/invgamma/quantile/test/test.native.js
+++ b/lib/node_modules/@stdlib/stats/base/dists/invgamma/quantile/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 //
@@ -127,10 +126,8 @@ tape( 'if provided a nonpositive `beta`, the function returns `NaN`', opts, func
 
 tape( 'the function evaluates the quantile for `x` given large parameters `alpha` and `beta`', opts, function test( t ) {
 	var expected;
-	var delta;
 	var alpha;
 	var beta;
-	var tol;
 	var p;
 	var y;
 	var i;
@@ -141,23 +138,15 @@ tape( 'the function evaluates the quantile for `x` given large parameters `alpha
 	p = bothLarge.p;
 	for ( i = 0; i < p.length; i++ ) {
 		y = quantile( p[ i ], alpha[ i ], beta[ i ] );
-		if ( y === expected[ i ] ) {
-			t.strictEqual( y, expected[ i ], 'p: '+p[ i ]+', alpha: '+alpha[ i ]+', beta: '+beta[ i ]+', y: '+y+', expected: '+expected[ i ] );
-		} else {
-			delta = abs( y - expected[ i ] );
-			tol = 200.0 * EPS * abs( expected[ i ] );
-			t.ok( delta <= tol, 'within tolerance. p: '+p[ i ]+'. alpha: '+alpha[ i ]+'. beta: '+beta[ i ]+'. y: '+y+'. E: '+expected[ i ]+'. Δ: '+delta+'. tol: '+tol+'.' );
-		}
+		t.strictEqual( isAlmostSameValue( y, expected[ i ], 23 ), true, 'returns expected value' );
 	}
 	t.end();
 });
 
 tape( 'the function evaluates the quantile for `x` given large shape parameter `alpha`', opts, function test( t ) {
 	var expected;
-	var delta;
 	var alpha;
 	var beta;
-	var tol;
 	var p;
 	var y;
 	var i;
@@ -168,23 +157,15 @@ tape( 'the function evaluates the quantile for `x` given large shape parameter `
 	p = largeShape.p;
 	for ( i = 0; i < p.length; i++ ) {
 		y = quantile( p[ i ], alpha[ i ], beta[ i ] );
-		if ( y === expected[ i ] ) {
-			t.strictEqual( y, expected[ i ], 'p: '+p[ i ]+', alpha: '+alpha[ i ]+', beta: '+beta[ i ]+', y: '+y+', expected: '+expected[ i ] );
-		} else {
-			delta = abs( y - expected[ i ] );
-			tol = 100.0 * EPS * abs( expected[ i ] );
-			t.ok( delta <= tol, 'within tolerance. p: '+p[ i ]+'. alpha: '+alpha[ i ]+'. beta: '+beta[ i ]+'. y: '+y+'. E: '+expected[ i ]+'. Δ: '+delta+'. tol: '+tol+'.' );
-		}
+		t.strictEqual( isAlmostSameValue( y, expected[ i ], 23 ), true, 'returns expected value' );
 	}
 	t.end();
 });
 
 tape( 'the function evaluates the quantile for `x` given large scale parameter `beta`', opts, function test( t ) {
 	var expected;
-	var delta;
 	var alpha;
 	var beta;
-	var tol;
 	var p;
 	var y;
 	var i;
@@ -195,13 +176,7 @@ tape( 'the function evaluates the quantile for `x` given large scale parameter `
 	p = largeRate.p;
 	for ( i = 0; i < p.length; i++ ) {
 		y = quantile( p[ i ], alpha[ i ], beta[ i ] );
-		if ( y === expected[ i ] ) {
-			t.strictEqual( y, expected[ i ], 'p: '+p[ i ]+', alpha: '+alpha[ 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 ]+'. alpha: '+alpha[ i ]+'. beta: '+beta[ i ]+'. y: '+y+'. E: '+expected[ i ]+'. Δ: '+delta+'. tol: '+tol+'.' );
-		}
+		t.strictEqual( isAlmostSameValue( y, expected[ i ], 81 ), true, 'returns expected value' );
 	}
 	t.end();
 });
diff --git a/lib/node_modules/@stdlib/stats/base/dists/invgamma/quantile/test/test.quantile.js b/lib/node_modules/@stdlib/stats/base/dists/invgamma/quantile/test/test.quantile.js
index a8240882f3fb..77134dcefa19 100644
--- a/lib/node_modules/@stdlib/stats/base/dists/invgamma/quantile/test/test.quantile.js
+++ b/lib/node_modules/@stdlib/stats/base/dists/invgamma/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' );
 
 
@@ -112,10 +111,8 @@ tape( 'if provided a negative `beta`, the function returns `NaN`', function test
 
 tape( 'the function evaluates the quantile for `x` given large parameters `alpha` and `beta`', function test( t ) {
 	var expected;
-	var delta;
 	var alpha;
 	var beta;
-	var tol;
 	var i;
 	var p;
 	var y;
@@ -126,23 +123,15 @@ tape( 'the function evaluates the quantile for `x` given large parameters `alpha
 	beta = bothLarge.beta;
 	for ( i = 0; i < p.length; i++ ) {
 		y = quantile( p[ i ], alpha[ i ], beta[ i ] );
-		if ( y === expected[ i ] ) {
-			t.strictEqual( y, expected[ i ], 'p: '+p[ i ]+', alpha:'+alpha[ i ]+', beta: '+beta[ i ]+', y: '+y+', expected: '+expected[ i ] );
-		} else {
-			delta = abs( y - expected[ i ] );
-			tol = 200.0 * EPS * abs( expected[ i ] );
-			t.ok( delta <= tol, 'within tolerance. p: '+p[ i ]+'. alpha: '+alpha[ i ]+'. beta: '+beta[ i ]+'. y: '+y+'. E: '+expected[ i ]+'. Δ: '+delta+'. tol: '+tol+'.' );
-		}
+		t.strictEqual( isAlmostSameValue( y, expected[ i ], 23 ), true, 'returns expected value' );
 	}
 	t.end();
 });
 
 tape( 'the function evaluates the quantile for `x` given large shape parameter `alpha`', function test( t ) {
 	var expected;
-	var delta;
 	var alpha;
 	var beta;
-	var tol;
 	var i;
 	var p;
 	var y;
@@ -153,23 +142,15 @@ tape( 'the function evaluates the quantile for `x` given large shape parameter `
 	beta = largeShape.beta;
 	for ( i = 0; i < p.length; i++ ) {
 		y = quantile( p[ i ], alpha[ i ], beta[ i ] );
-		if ( y === expected[ i ] ) {
-			t.strictEqual( y, expected[ i ], 'p: '+p[ i ]+', alpha:'+alpha[ i ]+', beta: '+beta[ i ]+', y: '+y+', expected: '+expected[ i ] );
-		} else {
-			delta = abs( y - expected[ i ] );
-			tol = 100.0 * EPS * abs( expected[ i ] );
-			t.ok( delta <= tol, 'within tolerance. p: '+p[ i ]+'. alpha: '+alpha[ i ]+'. beta: '+beta[ i ]+'. y: '+y+'. E: '+expected[ i ]+'. Δ: '+delta+'. tol: '+tol+'.' );
-		}
+		t.strictEqual( isAlmostSameValue( y, expected[ i ], 23 ), true, 'returns expected value' );
 	}
 	t.end();
 });
 
 tape( 'the function evaluates the quantile for `x` given large scale parameter `beta`', function test( t ) {
 	var expected;
-	var delta;
 	var alpha;
 	var beta;
-	var tol;
 	var i;
 	var p;
 	var y;
@@ -180,13 +161,7 @@ tape( 'the function evaluates the quantile for `x` given large scale parameter `
 	beta = largeRate.beta;
 	for ( i = 0; i < p.length; i++ ) {
 		y = quantile( p[ i ], alpha[ i ], beta[ i ] );
-		if ( y === expected[ i ] ) {
-			t.strictEqual( y, expected[ i ], 'p: '+p[ i ]+', alpha:'+alpha[ 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 ]+'. alpha: '+alpha[ i ]+'. beta: '+beta[ i ]+'. y: '+y+'. E: '+expected[ i ]+'. Δ: '+delta+'. tol: '+tol+'.' );
-		}
+		t.strictEqual( isAlmostSameValue( y, expected[ i ], 81 ), true, 'returns expected value' );
 	}
 	t.end();
 });