diff --git a/lib/node_modules/@stdlib/stats/base/dists/poisson/pmf/test/test.factory.js b/lib/node_modules/@stdlib/stats/base/dists/poisson/pmf/test/test.factory.js
index d078eb1ec5e4..c2d51827a9ad 100644
--- a/lib/node_modules/@stdlib/stats/base/dists/poisson/pmf/test/test.factory.js
+++ b/lib/node_modules/@stdlib/stats/base/dists/poisson/pmf/test/test.factory.js
@@ -21,10 +21,9 @@
 // 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 EPS = require( '@stdlib/constants/float64/eps' );
 var factory = require( './../lib/factory.js' );
 
 
@@ -147,9 +146,7 @@ tape( 'if provided `lambda` equals `0`, the created function evaluates a degener
 tape( 'the created function evaluates the pmf for `x` given small `lambda`', function test( t ) {
 	var expected;
 	var lambda;
-	var delta;
 	var pmf;
-	var tol;
 	var i;
 	var x;
 	var y;
@@ -160,13 +157,7 @@ tape( 'the created function evaluates the pmf for `x` given small `lambda`', fun
 	for ( i = 0; i < x.length; i++ ) {
 		pmf = factory( lambda[i] );
 		y = pmf( x[i] );
-		if ( y === expected[i] ) {
-			t.strictEqual( y, expected[i], 'x: '+x[i]+'. lambda: '+lambda[i]+', y: '+y+', expected: '+expected[i] );
-		} else {
-			delta = abs( y - expected[ i ] );
-			tol = 280.0 * EPS * abs( expected[ i ] );
-			t.ok( delta <= tol, 'within tolerance. x: '+x[ i ]+'. lambda: '+lambda[i]+'. y: '+y+'. E: '+expected[ i ]+'. Δ: '+delta+'. tol: '+tol+'.' );
-		}
+		t.strictEqual( isAlmostSameValue( y, expected[ i ], 448 ), true, 'returns expected value' );
 	}
 	t.end();
 });
@@ -174,9 +165,7 @@ tape( 'the created function evaluates the pmf for `x` given small `lambda`', fun
 tape( 'the function evaluates the pmf for `x` given medium `lambda`', function test( t ) {
 	var expected;
 	var lambda;
-	var delta;
 	var pmf;
-	var tol;
 	var i;
 	var x;
 	var y;
@@ -187,13 +176,7 @@ tape( 'the function evaluates the pmf for `x` given medium `lambda`', function t
 	for ( i = 0; i < x.length; i++ ) {
 		pmf = factory( lambda[i] );
 		y = pmf( x[i] );
-		if ( y === expected[i] ) {
-			t.strictEqual( y, expected[i], 'x: '+x[i]+'. lambda: '+lambda[i]+', y: '+y+', expected: '+expected[i] );
-		} else {
-			delta = abs( y - expected[ i ] );
-			tol = 160.0 * EPS * abs( expected[ i ] );
-			t.ok( delta <= tol, 'within tolerance. x: '+x[ i ]+'. lambda: '+lambda[i]+'. y: '+y+'. E: '+expected[ i ]+'. Δ: '+delta+'. tol: '+tol+'.' );
-		}
+		t.strictEqual( isAlmostSameValue( y, expected[ i ], 242 ), true, 'returns expected value' );
 	}
 	t.end();
 });
@@ -201,9 +184,7 @@ tape( 'the function evaluates the pmf for `x` given medium `lambda`', function t
 tape( 'the function evaluates the pmf for `x` given large `lambda`', function test( t ) {
 	var expected;
 	var lambda;
-	var delta;
 	var pmf;
-	var tol;
 	var i;
 	var x;
 	var y;
@@ -214,13 +195,7 @@ tape( 'the function evaluates the pmf for `x` given large `lambda`', function te
 	for ( i = 0; i < x.length; i++ ) {
 		pmf = factory( lambda[i] );
 		y = pmf( x[i] );
-		if ( y === expected[i] ) {
-			t.strictEqual( y, expected[i], 'x: '+x[i]+'. lambda: '+lambda[i]+', y: '+y+', expected: '+expected[i] );
-		} else {
-			delta = abs( y - expected[ i ] );
-			tol = 160.0 * EPS * abs( expected[ i ] );
-			t.ok( delta <= tol, 'within tolerance. x: '+x[ i ]+'. lambda: '+lambda[i]+'. y: '+y+'. E: '+expected[ i ]+'. Δ: '+delta+'. tol: '+tol+'.' );
-		}
+		t.strictEqual( isAlmostSameValue( y, expected[ i ], 205 ), true, 'returns expected value' );
 	}
 	t.end();
 });
diff --git a/lib/node_modules/@stdlib/stats/base/dists/poisson/pmf/test/test.native.js b/lib/node_modules/@stdlib/stats/base/dists/poisson/pmf/test/test.native.js
index 64d855be9f7b..a59945a5891a 100644
--- a/lib/node_modules/@stdlib/stats/base/dists/poisson/pmf/test/test.native.js
+++ b/lib/node_modules/@stdlib/stats/base/dists/poisson/pmf/test/test.native.js
@@ -22,11 +22,10 @@
 
 var resolve = require( 'path' ).resolve;
 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 tryRequire = require( '@stdlib/utils/try-require' );
 
 
@@ -126,8 +125,6 @@ tape( 'if provided `lambda` equals `0`, the function evaluates a degenerate dist
 tape( 'the function evaluates the pmf for `x` given a small parameter `lambda`', opts, function test( t ) {
 	var expected;
 	var lambda;
-	var delta;
-	var tol;
 	var i;
 	var x;
 	var y;
@@ -137,13 +134,7 @@ tape( 'the function evaluates the pmf for `x` given a small parameter `lambda`',
 	lambda = smallMean.lambda;
 	for ( i = 0; i < x.length; i++ ) {
 		y = pmf( x[ i ], lambda[ i ] );
-		if ( y === expected[ i ] ) {
-			t.strictEqual( y, expected[ i ], 'x: '+x[ i ]+'. lambda: '+lambda[ i ]+', y: '+y+', expected: '+expected[ i ] );
-		} else {
-			delta = abs( y - expected[ i ] );
-			tol = 280.0 * EPS * abs( expected[ i ] );
-			t.ok( delta <= tol, 'within tolerance. x: '+x[ i ]+'. lambda: '+lambda[ i ]+'. y: '+y+'. E: '+expected[ i ]+'. Δ: '+delta+'. tol: '+tol+'.' );
-		}
+		t.strictEqual( isAlmostSameValue( y, expected[ i ], 448 ), true, 'returns expected value' );
 	}
 	t.end();
 });
@@ -151,8 +142,6 @@ tape( 'the function evaluates the pmf for `x` given a small parameter `lambda`',
 tape( 'the function evaluates the pmf for `x` given a medium parameter `lambda`', opts, function test( t ) {
 	var expected;
 	var lambda;
-	var delta;
-	var tol;
 	var i;
 	var x;
 	var y;
@@ -162,13 +151,7 @@ tape( 'the function evaluates the pmf for `x` given a medium parameter `lambda`'
 	lambda = mediumMean.lambda;
 	for ( i = 0; i < x.length; i++ ) {
 		y = pmf( x[ i ], lambda[ i ] );
-		if ( y === expected[ i ] ) {
-			t.strictEqual( y, expected[ i ], 'x: '+x[ i ]+'. lambda: '+lambda[ i ]+', y: '+y+', expected: '+expected[ i ] );
-		} else {
-			delta = abs( y - expected[ i ] );
-			tol = 160.0 * EPS * abs( expected[ i ] );
-			t.ok( delta <= tol, 'within tolerance. x: '+x[ i ]+'. lambda: '+lambda[ i ]+'. y: '+y+'. E: '+expected[ i ]+'. Δ: '+delta+'. tol: '+tol+'.' );
-		}
+		t.strictEqual( isAlmostSameValue( y, expected[ i ], 242 ), true, 'returns expected value' );
 	}
 	t.end();
 });
@@ -176,8 +159,6 @@ tape( 'the function evaluates the pmf for `x` given a medium parameter `lambda`'
 tape( 'the function evaluates the pmf for `x` given a large parameter `lambda`', opts, function test( t ) {
 	var expected;
 	var lambda;
-	var delta;
-	var tol;
 	var i;
 	var x;
 	var y;
@@ -187,13 +168,7 @@ tape( 'the function evaluates the pmf for `x` given a large parameter `lambda`',
 	lambda = largeMean.lambda;
 	for ( i = 0; i < x.length; i++ ) {
 		y = pmf( x[ i ], lambda[ i ] );
-		if ( y === expected[ i ] ) {
-			t.strictEqual( y, expected[ i ], 'x: '+x[ i ]+'. lambda: '+lambda[ i ]+', y: '+y+', expected: '+expected[ i ] );
-		} else {
-			delta = abs( y - expected[ i ] );
-			tol = 160.0 * EPS * abs( expected[ i ] );
-			t.ok( delta <= tol, 'within tolerance. x: '+x[ i ]+'. lambda: '+lambda[ i ]+'. y: '+y+'. E: '+expected[ i ]+'. Δ: '+delta+'. tol: '+tol+'.' );
-		}
+		t.strictEqual( isAlmostSameValue( y, expected[ i ], 205 ), true, 'returns expected value' );
 	}
 	t.end();
 });
diff --git a/lib/node_modules/@stdlib/stats/base/dists/poisson/pmf/test/test.pmf.js b/lib/node_modules/@stdlib/stats/base/dists/poisson/pmf/test/test.pmf.js
index 57970ff1ccc5..346aff174e4e 100644
--- a/lib/node_modules/@stdlib/stats/base/dists/poisson/pmf/test/test.pmf.js
+++ b/lib/node_modules/@stdlib/stats/base/dists/poisson/pmf/test/test.pmf.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 pmf = require( './../lib' );
 
 
@@ -117,8 +116,6 @@ tape( 'if provided `lambda` equals `0`, the function evaluates a degenerate dist
 tape( 'the function evaluates the pmf for `x` given a small parameter `lambda`', function test( t ) {
 	var expected;
 	var lambda;
-	var delta;
-	var tol;
 	var i;
 	var x;
 	var y;
@@ -128,13 +125,7 @@ tape( 'the function evaluates the pmf for `x` given a small parameter `lambda`',
 	lambda = smallMean.lambda;
 	for ( i = 0; i < x.length; i++ ) {
 		y = pmf( x[i], lambda[i] );
-		if ( y === expected[i] ) {
-			t.strictEqual( y, expected[i], 'x: '+x[i]+'. lambda: '+lambda[i]+', y: '+y+', expected: '+expected[i] );
-		} else {
-			delta = abs( y - expected[ i ] );
-			tol = 280.0 * EPS * abs( expected[ i ] );
-			t.ok( delta <= tol, 'within tolerance. x: '+x[ i ]+'. lambda: '+lambda[i]+'. y: '+y+'. E: '+expected[ i ]+'. Δ: '+delta+'. tol: '+tol+'.' );
-		}
+		t.strictEqual( isAlmostSameValue( y, expected[ i ], 448 ), true, 'returns expected value' );
 	}
 	t.end();
 });
@@ -142,8 +133,6 @@ tape( 'the function evaluates the pmf for `x` given a small parameter `lambda`',
 tape( 'the function evaluates the pmf for `x` given a medium parameter `lambda`', function test( t ) {
 	var expected;
 	var lambda;
-	var delta;
-	var tol;
 	var i;
 	var x;
 	var y;
@@ -153,13 +142,7 @@ tape( 'the function evaluates the pmf for `x` given a medium parameter `lambda`'
 	lambda = mediumMean.lambda;
 	for ( i = 0; i < x.length; i++ ) {
 		y = pmf( x[i], lambda[i] );
-		if ( y === expected[i] ) {
-			t.strictEqual( y, expected[i], 'x: '+x[i]+'. lambda: '+lambda[i]+', y: '+y+', expected: '+expected[i] );
-		} else {
-			delta = abs( y - expected[ i ] );
-			tol = 160.0 * EPS * abs( expected[ i ] );
-			t.ok( delta <= tol, 'within tolerance. x: '+x[ i ]+'. lambda: '+lambda[i]+'. y: '+y+'. E: '+expected[ i ]+'. Δ: '+delta+'. tol: '+tol+'.' );
-		}
+		t.strictEqual( isAlmostSameValue( y, expected[ i ], 242 ), true, 'returns expected value' );
 	}
 	t.end();
 });
@@ -167,8 +150,6 @@ tape( 'the function evaluates the pmf for `x` given a medium parameter `lambda`'
 tape( 'the function evaluates the pmf for `x` given a large parameter `lambda`', function test( t ) {
 	var expected;
 	var lambda;
-	var delta;
-	var tol;
 	var i;
 	var x;
 	var y;
@@ -178,13 +159,7 @@ tape( 'the function evaluates the pmf for `x` given a large parameter `lambda`',
 	lambda = largeMean.lambda;
 	for ( i = 0; i < x.length; i++ ) {
 		y = pmf( x[i], lambda[i] );
-		if ( y === expected[i] ) {
-			t.strictEqual( y, expected[i], 'x: '+x[i]+'. lambda: '+lambda[i]+', y: '+y+', expected: '+expected[i] );
-		} else {
-			delta = abs( y - expected[ i ] );
-			tol = 160.0 * EPS * abs( expected[ i ] );
-			t.ok( delta <= tol, 'within tolerance. x: '+x[ i ]+'. lambda: '+lambda[i]+'. y: '+y+'. E: '+expected[ i ]+'. Δ: '+delta+'. tol: '+tol+'.' );
-		}
+		t.strictEqual( isAlmostSameValue( y, expected[ i ], 205 ), true, 'returns expected value' );
 	}
 	t.end();
 });