diff --git a/lib/node_modules/@stdlib/stats/base/dists/cosine/pdf/test/test.factory.js b/lib/node_modules/@stdlib/stats/base/dists/cosine/pdf/test/test.factory.js
index 002d5759b617..524d8c341fac 100644
--- a/lib/node_modules/@stdlib/stats/base/dists/cosine/pdf/test/test.factory.js
+++ b/lib/node_modules/@stdlib/stats/base/dists/cosine/pdf/test/test.factory.js
@@ -22,7 +22,7 @@
 
 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' );
@@ -169,9 +169,7 @@ tape( 'if `s` equals `0`, the created function evaluates a degenerate distributi
 
 tape( 'the created function evaluates the pdf for `x` given positive `mu`', function test( t ) {
 	var expected;
-	var delta;
 	var pdf;
-	var tol;
 	var mu;
 	var s;
 	var x;
@@ -185,14 +183,8 @@ tape( 'the created function evaluates the pdf for `x` given positive `mu`', func
 	for ( i = 0; i < x.length; i++ ) {
 		pdf = factory( mu[i], s[i] );
 		y = pdf( x[i] );
-		if ( expected[i] !== null && !(expected[i] === 0.0 && y < EPS ) ) {
-			if ( y === expected[i] ) {
-				t.strictEqual( y, expected[i], 'x: '+x[i]+', mu: '+mu[i]+', s: '+s[i]+', y: '+y+', expected: '+expected[i] );
-			} else {
-				delta = abs( y - expected[ i ] );
-				tol = 3.0 * EPS * abs( expected[ i ] );
-				t.ok( delta <= tol, 'within tolerance. x: '+x[ i ]+'. mu: '+mu[i]+'. s: '+s[i]+'. y: '+y+'. E: '+expected[ i ]+'. Δ: '+delta+'. tol: '+tol+'.' );
-			}
+		if ( expected[i] !== null && !( expected[i] === 0.0 && y < EPS ) ) {
+			t.strictEqual( isAlmostSameValue( y, expected[ i ], 0 ), true, 'returns expected value' );
 		}
 	}
 	t.end();
@@ -200,9 +192,7 @@ tape( 'the created function evaluates the pdf for `x` given positive `mu`', func
 
 tape( 'the created function evaluates the pdf for `x` given negative `mu`', function test( t ) {
 	var expected;
-	var delta;
 	var pdf;
-	var tol;
 	var mu;
 	var s;
 	var x;
@@ -216,14 +206,8 @@ tape( 'the created function evaluates the pdf for `x` given negative `mu`', func
 	for ( i = 0; i < x.length; i++ ) {
 		pdf = factory( mu[i], s[i] );
 		y = pdf( x[i] );
-		if ( expected[i] !== null && !(expected[i] === 0.0 && y < EPS ) ) {
-			if ( y === expected[i] ) {
-				t.strictEqual( y, expected[i], 'x: '+x[i]+', mu: '+mu[i]+', s: '+s[i]+', y: '+y+', expected: '+expected[i] );
-			} else {
-				delta = abs( y - expected[ i ] );
-				tol = 4.0 * EPS * abs( expected[ i ] );
-				t.ok( delta <= tol, 'within tolerance. x: '+x[ i ]+'. mu: '+mu[i]+'. s: '+s[i]+'. y: '+y+'. E: '+expected[ i ]+'. Δ: '+delta+'. tol: '+tol+'.' );
-			}
+		if ( expected[i] !== null && !( expected[i] === 0.0 && y < EPS ) ) {
+			t.strictEqual( isAlmostSameValue( y, expected[ i ], 4 ), true, 'returns expected value' );
 		}
 	}
 	t.end();
@@ -231,9 +215,7 @@ tape( 'the created function evaluates the pdf for `x` given negative `mu`', func
 
 tape( 'the created function evaluates the pdf for `x` given large variance ( = large `s`)', function test( t ) {
 	var expected;
-	var delta;
 	var pdf;
-	var tol;
 	var mu;
 	var s;
 	var x;
@@ -248,13 +230,7 @@ tape( 'the created function evaluates the pdf for `x` given large variance ( = l
 		pdf = factory( mu[i], s[i] );
 		y = pdf( x[i] );
 		if ( expected[i] !== null ) {
-			if ( y === expected[i] ) {
-				t.strictEqual( y, expected[i], 'x: '+x[i]+', mu: '+mu[i]+', s: '+s[i]+', y: '+y+', expected: '+expected[i] );
-			} else {
-				delta = abs( y - expected[ i ] );
-				tol = 3.0 * EPS * abs( expected[ i ] );
-				t.ok( delta <= tol, 'within tolerance. x: '+x[ i ]+'. mu: '+mu[i]+'. s: '+s[i]+'. y: '+y+'. E: '+expected[ i ]+'. Δ: '+delta+'. tol: '+tol+'.' );
-			}
+			t.strictEqual( isAlmostSameValue( y, expected[ i ], 2 ), true, 'returns expected value' );
 		}
 	}
 	t.end();
diff --git a/lib/node_modules/@stdlib/stats/base/dists/cosine/pdf/test/test.native.js b/lib/node_modules/@stdlib/stats/base/dists/cosine/pdf/test/test.native.js
index 6cb879c58b96..d29b03682819 100644
--- a/lib/node_modules/@stdlib/stats/base/dists/cosine/pdf/test/test.native.js
+++ b/lib/node_modules/@stdlib/stats/base/dists/cosine/pdf/test/test.native.js
@@ -23,7 +23,7 @@
 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 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' );
@@ -136,8 +136,6 @@ tape( 'if provided `s` equal to `0`, the function evaluates a degenerate distrib
 
 tape( 'the function evaluates the pdf for `x` given positive `mu`', opts, function test( t ) {
 	var expected;
-	var delta;
-	var tol;
 	var mu;
 	var x;
 	var s;
@@ -150,14 +148,8 @@ tape( 'the function evaluates the pdf for `x` given positive `mu`', opts, functi
 	s = positiveMean.s;
 	for ( i = 0; i < x.length; i++ ) {
 		y = pdf( x[i], mu[i], s[i] );
-		if ( expected[i] !== null && !(expected[i] === 0.0 && y < EPS ) ) {
-			if ( y === expected[i] ) {
-				t.strictEqual( y, expected[i], 'x: '+x[i]+', mu: '+mu[i]+', s: '+s[i]+', y: '+y+', expected: '+expected[i] );
-			} else {
-				delta = abs( y - expected[ i ] );
-				tol = 3.0 * EPS * abs( expected[ i ] );
-				t.ok( delta <= tol, 'within tolerance. x: '+x[ i ]+'. mu: '+mu[i]+'. s: '+s[i]+'. y: '+y+'. E: '+expected[ i ]+'. Δ: '+delta+'. tol: '+tol+'.' );
-			}
+		if ( expected[i] !== null && !( expected[i] === 0.0 && y < EPS ) ) {
+			t.strictEqual( isAlmostSameValue( y, expected[ i ], 0 ), true, 'returns expected value' );
 		}
 	}
 	t.end();
@@ -165,8 +157,6 @@ tape( 'the function evaluates the pdf for `x` given positive `mu`', opts, functi
 
 tape( 'the function evaluates the pdf for `x` given negative `mu`', opts, function test( t ) {
 	var expected;
-	var delta;
-	var tol;
 	var mu;
 	var x;
 	var s;
@@ -179,14 +169,8 @@ tape( 'the function evaluates the pdf for `x` given negative `mu`', opts, functi
 	s = negativeMean.s;
 	for ( i = 0; i < x.length; i++ ) {
 		y = pdf( x[i], mu[i], s[i] );
-		if ( expected[i] !== null && !(expected[i] === 0.0 && y < EPS ) ) {
-			if ( y === expected[i] ) {
-				t.strictEqual( y, expected[i], 'x: '+x[i]+', mu: '+mu[i]+', s: '+s[i]+', y: '+y+', expected: '+expected[i] );
-			} else {
-				delta = abs( y - expected[ i ] );
-				tol = 4.0 * EPS * abs( expected[ i ] );
-				t.ok( delta <= tol, 'within tolerance. x: '+x[ i ]+'. mu: '+mu[i]+'. s: '+s[i]+'. y: '+y+'. E: '+expected[ i ]+'. Δ: '+delta+'. tol: '+tol+'.' );
-			}
+		if ( expected[i] !== null && !( expected[i] === 0.0 && y < EPS ) ) {
+			t.strictEqual( isAlmostSameValue( y, expected[ i ], 4 ), true, 'returns expected value' );
 		}
 	}
 	t.end();
@@ -194,8 +178,6 @@ tape( 'the function evaluates the pdf for `x` given negative `mu`', opts, functi
 
 tape( 'the function evaluates the pdf for `x` given large variance ( = large `s` )', opts, function test( t ) {
 	var expected;
-	var delta;
-	var tol;
 	var mu;
 	var x;
 	var s;
@@ -209,13 +191,7 @@ tape( 'the function evaluates the pdf for `x` given large variance ( = large `s`
 	for ( i = 0; i < x.length; i++ ) {
 		y = pdf( x[i], mu[i], s[i] );
 		if ( expected[i] !== null ) {
-			if ( y === expected[i] ) {
-				t.strictEqual( y, expected[i], 'x: '+x[i]+', mu: '+mu[i]+', s: '+s[i]+', y: '+y+', expected: '+expected[i] );
-			} else {
-				delta = abs( y - expected[ i ] );
-				tol = 3.0 * EPS * abs( expected[ i ] );
-				t.ok( delta <= tol, 'within tolerance. x: '+x[ i ]+'. mu: '+mu[i]+'. s: '+s[i]+'. y: '+y+'. E: '+expected[ i ]+'. Δ: '+delta+'. tol: '+tol+'.' );
-			}
+			t.strictEqual( isAlmostSameValue( y, expected[ i ], 2 ), true, 'returns expected value' );
 		}
 	}
 	t.end();
diff --git a/lib/node_modules/@stdlib/stats/base/dists/cosine/pdf/test/test.pdf.js b/lib/node_modules/@stdlib/stats/base/dists/cosine/pdf/test/test.pdf.js
index 9cfcb04122fa..5c47459ae45f 100644
--- a/lib/node_modules/@stdlib/stats/base/dists/cosine/pdf/test/test.pdf.js
+++ b/lib/node_modules/@stdlib/stats/base/dists/cosine/pdf/test/test.pdf.js
@@ -22,7 +22,7 @@
 
 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' );
@@ -127,8 +127,6 @@ tape( 'if provided `s` equal to `0`, the function evaluates a degenerate distrib
 
 tape( 'the function evaluates the pdf for `x` given positive `mu`', function test( t ) {
 	var expected;
-	var delta;
-	var tol;
 	var mu;
 	var x;
 	var s;
@@ -141,14 +139,8 @@ tape( 'the function evaluates the pdf for `x` given positive `mu`', function tes
 	s = positiveMean.s;
 	for ( i = 0; i < x.length; i++ ) {
 		y = pdf( x[i], mu[i], s[i] );
-		if ( expected[i] !== null && !(expected[i] === 0.0 && y < EPS ) ) {
-			if ( y === expected[i] ) {
-				t.strictEqual( y, expected[i], 'x: '+x[i]+', mu: '+mu[i]+', s: '+s[i]+', y: '+y+', expected: '+expected[i] );
-			} else {
-				delta = abs( y - expected[ i ] );
-				tol = 3.0 * EPS * abs( expected[ i ] );
-				t.ok( delta <= tol, 'within tolerance. x: '+x[ i ]+'. mu: '+mu[i]+'. s: '+s[i]+'. y: '+y+'. E: '+expected[ i ]+'. Δ: '+delta+'. tol: '+tol+'.' );
-			}
+		if ( expected[i] !== null && !( expected[i] === 0.0 && y < EPS ) ) {
+			t.strictEqual( isAlmostSameValue( y, expected[ i ], 0 ), true, 'returns expected value' );
 		}
 	}
 	t.end();
@@ -156,8 +148,6 @@ tape( 'the function evaluates the pdf for `x` given positive `mu`', function tes
 
 tape( 'the function evaluates the pdf for `x` given negative `mu`', function test( t ) {
 	var expected;
-	var delta;
-	var tol;
 	var mu;
 	var x;
 	var s;
@@ -170,14 +160,8 @@ tape( 'the function evaluates the pdf for `x` given negative `mu`', function tes
 	s = negativeMean.s;
 	for ( i = 0; i < x.length; i++ ) {
 		y = pdf( x[i], mu[i], s[i] );
-		if ( expected[i] !== null && !(expected[i] === 0.0 && y < EPS ) ) {
-			if ( y === expected[i] ) {
-				t.strictEqual( y, expected[i], 'x: '+x[i]+', mu: '+mu[i]+', s: '+s[i]+', y: '+y+', expected: '+expected[i] );
-			} else {
-				delta = abs( y - expected[ i ] );
-				tol = 4.0 * EPS * abs( expected[ i ] );
-				t.ok( delta <= tol, 'within tolerance. x: '+x[ i ]+'. mu: '+mu[i]+'. s: '+s[i]+'. y: '+y+'. E: '+expected[ i ]+'. Δ: '+delta+'. tol: '+tol+'.' );
-			}
+		if ( expected[i] !== null && !( expected[i] === 0.0 && y < EPS ) ) {
+			t.strictEqual( isAlmostSameValue( y, expected[ i ], 4 ), true, 'returns expected value' );
 		}
 	}
 	t.end();
@@ -185,8 +169,6 @@ tape( 'the function evaluates the pdf for `x` given negative `mu`', function tes
 
 tape( 'the function evaluates the pdf for `x` given large variance ( = large `s` )', function test( t ) {
 	var expected;
-	var delta;
-	var tol;
 	var mu;
 	var x;
 	var s;
@@ -200,13 +182,7 @@ tape( 'the function evaluates the pdf for `x` given large variance ( = large `s`
 	for ( i = 0; i < x.length; i++ ) {
 		y = pdf( x[i], mu[i], s[i] );
 		if ( expected[i] !== null ) {
-			if ( y === expected[i] ) {
-				t.strictEqual( y, expected[i], 'x: '+x[i]+', mu: '+mu[i]+', s: '+s[i]+', y: '+y+', expected: '+expected[i] );
-			} else {
-				delta = abs( y - expected[ i ] );
-				tol = 3.0 * EPS * abs( expected[ i ] );
-				t.ok( delta <= tol, 'within tolerance. x: '+x[ i ]+'. mu: '+mu[i]+'. s: '+s[i]+'. y: '+y+'. E: '+expected[ i ]+'. Δ: '+delta+'. tol: '+tol+'.' );
-			}
+			t.strictEqual( isAlmostSameValue( y, expected[ i ], 2 ), true, 'returns expected value' );
 		}
 	}
 	t.end();