diff --git a/lib/node_modules/@stdlib/math/base/special/gamma/test/test.boost.js b/lib/node_modules/@stdlib/math/base/special/gamma/test/test.boost.js
index 98919f47c99c..c450593d7ccd 100644
--- a/lib/node_modules/@stdlib/math/base/special/gamma/test/test.boost.js
+++ b/lib/node_modules/@stdlib/math/base/special/gamma/test/test.boost.js
@@ -23,10 +23,9 @@
 var tape = require( 'tape' );
 var isnan = require( '@stdlib/math/base/assert/is-nan' );
 var incrspace = require( '@stdlib/array/base/incrspace' );
-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 gamma = require( './../lib/boost/gamma.js' );
 
 
@@ -107,8 +106,6 @@ tape( 'if `x <= -170.65`, the function returns zero', function test( t ) {
 
 tape( 'the function evaluates the gamma function (positive integers)', function test( t ) {
 	var expected;
-	var delta;
-	var tol;
 	var x;
 	var v;
 	var i;
@@ -118,17 +115,13 @@ tape( 'the function evaluates the gamma function (positive integers)', function
 
 	for ( i = 0; i < x.length; i++ ) {
 		v = gamma( x[ i ] );
-		delta = abs( v - expected[ i ] );
-		tol = 10.0 * EPS * abs( expected[ i ] );
-		t.ok( delta <= tol, 'within tolerance. x: '+x[i]+'. y: '+v+'. E: '+ expected[i]+'. Δ: '+delta+'. tol: '+tol+'.' );
+		t.strictEqual( isAlmostSameValue( v, expected[ i ], 3 ), true, 'returns expected value' );
 	}
 	t.end();
 });
 
 tape( 'the function evaluates the gamma function (decimal values)', function test( t ) {
 	var expected;
-	var delta;
-	var tol;
 	var x;
 	var v;
 	var i;
@@ -138,17 +131,13 @@ tape( 'the function evaluates the gamma function (decimal values)', function tes
 
 	for ( i = 0; i < x.length; i++ ) {
 		v = gamma( x[ i ] );
-		delta = abs( v - expected[ i ] );
-		tol = 10.0 * EPS * abs( expected[ i ] );
-		t.ok( delta <= tol, 'within tolerance. x: '+x[i]+'. y: '+v+'. E: '+ expected[i]+'. Δ: '+delta+'. tol: '+tol+'.' );
+		t.strictEqual( isAlmostSameValue( v, expected[ i ], 11 ), true, 'returns expected value' );
 	}
 	t.end();
 });
 
 tape( 'the function evaluates the gamma function (values near one)', function test( t ) {
 	var expected;
-	var delta;
-	var tol;
 	var x;
 	var v;
 	var i;
@@ -158,9 +147,7 @@ tape( 'the function evaluates the gamma function (values near one)', function te
 
 	for ( i = 0; i < x.length; i++ ) {
 		v = gamma( x[ i ] );
-		delta = abs( v - expected[ i ] );
-		tol = 10.0 * EPS * abs( expected[ i ] );
-		t.ok( delta <= tol, 'within tolerance. x: '+x[i]+'. y: '+v+'. E: '+ expected[i]+'. Δ: '+delta+'. tol: '+tol+'.' );
+		t.strictEqual( isAlmostSameValue( v, expected[ i ], 13 ), true, 'returns expected value' );
 	}
 	t.end();
 });
diff --git a/lib/node_modules/@stdlib/math/base/special/gamma/test/test.js b/lib/node_modules/@stdlib/math/base/special/gamma/test/test.js
index c1f9b53ca04c..6509ae05d6f3 100644
--- a/lib/node_modules/@stdlib/math/base/special/gamma/test/test.js
+++ b/lib/node_modules/@stdlib/math/base/special/gamma/test/test.js
@@ -23,8 +23,7 @@
 var tape = require( 'tape' );
 var isnan = require( '@stdlib/math/base/assert/is-nan' );
 var incrspace = require( '@stdlib/array/base/incrspace' );
-var abs = require( '@stdlib/math/base/special/abs' );
-var max = require( '@stdlib/math/base/special/maxn' );
+var isAlmostSameValue = require( '@stdlib/assert/is-almost-same-value' );
 var PINF = require( '@stdlib/constants/float64/pinf' );
 var NINF = require( '@stdlib/constants/float64/ninf' );
 var gamma = require( './../lib' );
@@ -107,31 +106,23 @@ tape( 'if `x < -170.56749...`, the function returns zero', function test( t ) {
 });
 
 tape( 'the function evaluates the gamma function (positive integers)', function test( t ) {
-	var delta;
-	var tol;
 	var v;
 	var i;
 
 	for ( i = 0; i < data1.length; i++ ) {
 		v = gamma( data1[ i ] );
-		delta = abs( v - expected1[ i ] );
-		tol = 2.75e-12 * max( 1.0, abs( v ), abs( expected1[ i ] ) );
-		t.ok( delta <= tol, 'within tolerance. x: ' + data1[ i ] + '. Value: ' + v + '. Expected: ' + expected1[ i ] + '. Tolerance: ' + tol + '.' );
+		t.strictEqual( isAlmostSameValue( v, expected1[ i ], 1009 ), true, 'returns expected value' );
 	}
 	t.end();
 });
 
 tape( 'the function evaluates the gamma function (decimal values)', function test( t ) {
-	var delta;
-	var tol;
 	var v;
 	var i;
 
 	for ( i = 0; i < data2.length; i++ ) {
 		v = gamma( data2[ i ] );
-		delta = abs( v - expected2[ i ] );
-		tol = 2.75e-12 * max( 1.0, abs( v ), abs( expected2[ i ] ) );
-		t.ok( delta <= tol, 'within tolerance. x: ' + data2[ i ] + '. Value: ' + v + '. Expected: ' + expected2[ i ] + '. Tolerance: ' + tol + '.' );
+		t.strictEqual( isAlmostSameValue( v, expected2[ i ], 1942 ), true, 'returns expected value' );
 	}
 	t.end();
 });
@@ -145,8 +136,6 @@ tape( 'if provided a positive integer, the function returns the factorial of (n-
 
 tape( 'the function uses a small value approximation near left boundaries (abs(x)<33)', function test( t ) {
 	var expected;
-	var delta;
-	var tol;
 	var x;
 	var v;
 
@@ -155,18 +144,13 @@ tape( 'the function uses a small value approximation near left boundaries (abs(x
 
 	expected = -1.8364310433257593e-6;
 
-	delta = abs( v - expected );
-	tol = 1.0e-14;
-
-	t.ok( delta < tol, 'within tolerance. x: ' + x + '. Value: ' + v + '. Expected: ' + expected + '. Tolerance: ' + tol + '.' );
+	t.strictEqual( isAlmostSameValue( v, expected, 1 ), true, 'returns expected value' );
 
 	t.end();
 });
 
 tape( 'the function uses a small value approximation near right boundaries (abs(x)<33)', function test( t ) {
 	var expected;
-	var delta;
-	var tol;
 	var x;
 	var v;
 
@@ -175,10 +159,7 @@ tape( 'the function uses a small value approximation near right boundaries (abs(
 
 	expected = -3.856505190983844e-5;
 
-	delta = abs( v - expected );
-	tol = 1.0e-14;
-
-	t.ok( delta < tol, 'within tolerance. x: ' + x + '. Value: ' + v + '. Expected: ' + expected + '. Tolerance: ' + tol + '.' );
+	t.strictEqual( isAlmostSameValue( v, expected, 1 ), true, 'returns expected value' );
 
 	t.end();
 });
diff --git a/lib/node_modules/@stdlib/math/base/special/gamma/test/test.native.js b/lib/node_modules/@stdlib/math/base/special/gamma/test/test.native.js
index a26ba1c11772..1c5cc61fe808 100644
--- a/lib/node_modules/@stdlib/math/base/special/gamma/test/test.native.js
+++ b/lib/node_modules/@stdlib/math/base/special/gamma/test/test.native.js
@@ -24,8 +24,7 @@ var resolve = require( 'path' ).resolve;
 var tape = require( 'tape' );
 var isnan = require( '@stdlib/math/base/assert/is-nan' );
 var incrspace = require( '@stdlib/array/base/incrspace' );
-var abs = require( '@stdlib/math/base/special/abs' );
-var max = require( '@stdlib/math/base/special/maxn' );
+var isAlmostSameValue = require( '@stdlib/assert/is-almost-same-value' );
 var PINF = require( '@stdlib/constants/float64/pinf' );
 var NINF = require( '@stdlib/constants/float64/ninf' );
 var tryRequire = require( '@stdlib/utils/try-require' );
@@ -116,31 +115,23 @@ tape( 'if `x < -170.56749...`, the function returns zero', opts, function test(
 });
 
 tape( 'the function evaluates the gamma function (positive integers)', opts, function test( t ) {
-	var delta;
-	var tol;
 	var v;
 	var i;
 
 	for ( i = 0; i < data1.length; i++ ) {
 		v = gamma( data1[ i ] );
-		delta = abs( v - expected1[ i ] );
-		tol = 2.75e-12 * max( 1.0, abs( v ), abs( expected1[ i ] ) );
-		t.ok( delta <= tol, 'within tolerance. x: ' + data1[ i ] + '. Value: ' + v + '. Expected: ' + expected1[ i ] + '. Tolerance: ' + tol + '.' );
+		t.strictEqual( isAlmostSameValue( v, expected1[ i ], 1009 ), true, 'returns expected value' );
 	}
 	t.end();
 });
 
 tape( 'the function evaluates the gamma function (decimal values)', opts, function test( t ) {
-	var delta;
-	var tol;
 	var v;
 	var i;
 
 	for ( i = 0; i < data2.length; i++ ) {
 		v = gamma( data2[ i ] );
-		delta = abs( v - expected2[ i ] );
-		tol = 2.75e-12 * max( 1.0, abs( v ), abs( expected2[ i ] ) );
-		t.ok( delta <= tol, 'within tolerance. x: ' + data2[ i ] + '. Value: ' + v + '. Expected: ' + expected2[ i ] + '. Tolerance: ' + tol + '.' );
+		t.strictEqual( isAlmostSameValue( v, expected2[ i ], 1942 ), true, 'returns expected value' );
 	}
 	t.end();
 });
@@ -154,8 +145,6 @@ tape( 'if provided a positive integer, the function returns the factorial of (n-
 
 tape( 'the function uses a small value approximation near left boundaries (abs(x)<33)', opts, function test( t ) {
 	var expected;
-	var delta;
-	var tol;
 	var x;
 	var v;
 
@@ -164,18 +153,13 @@ tape( 'the function uses a small value approximation near left boundaries (abs(x
 
 	expected = -1.8364310433257593e-6;
 
-	delta = abs( v - expected );
-	tol = 1.0e-14;
-
-	t.ok( delta < tol, 'within tolerance. x: ' + x + '. Value: ' + v + '. Expected: ' + expected + '. Tolerance: ' + tol + '.' );
+	t.strictEqual( isAlmostSameValue( v, expected, 1 ), true, 'returns expected value' );
 
 	t.end();
 });
 
 tape( 'the function uses a small value approximation near right boundaries (abs(x)<33)', opts, function test( t ) {
 	var expected;
-	var delta;
-	var tol;
 	var x;
 	var v;
 
@@ -184,10 +168,7 @@ tape( 'the function uses a small value approximation near right boundaries (abs(
 
 	expected = -3.856505190983844e-5;
 
-	delta = abs( v - expected );
-	tol = 1.0e-14;
-
-	t.ok( delta < tol, 'within tolerance. x: ' + x + '. Value: ' + v + '. Expected: ' + expected + '. Tolerance: ' + tol + '.' );
+	t.strictEqual( isAlmostSameValue( v, expected, 1 ), true, 'returns expected value' );
 
 	t.end();
 });