diff --git a/lib/node_modules/@stdlib/blas/base/zdrot/test/test.ndarray.js b/lib/node_modules/@stdlib/blas/base/zdrot/test/test.ndarray.js
index c940adafd196..88410c0fff54 100644
--- a/lib/node_modules/@stdlib/blas/base/zdrot/test/test.ndarray.js
+++ b/lib/node_modules/@stdlib/blas/base/zdrot/test/test.ndarray.js
@@ -23,8 +23,7 @@
 var tape = require( 'tape' );
 var Float64Array = require( '@stdlib/array/float64' );
 var Complex128Array = require( '@stdlib/array/complex128' );
-var EPS = require( '@stdlib/constants/float64/eps' );
-var abs = require( '@stdlib/math/base/special/abs' );
+var isAlmostSameValue = require( '@stdlib/assert/is-almost-same-value' );
 var zdrot = require( './../lib/ndarray.js' );
 
 
@@ -37,22 +36,14 @@ var zdrot = require( './../lib/ndarray.js' );
 * @param {Object} t - test object
 * @param {Collection} actual - actual values
 * @param {Collection} expected - expected values
-* @param {number} rtol - relative tolerance
+* @param {NonNegativeInteger} ulp - maximum allowed ULP difference
 */
-function isApprox( t, actual, expected, rtol ) {
-	var delta;
-	var tol;
+function isApprox( t, actual, expected, ulp ) {
 	var i;
 
 	t.strictEqual( actual.length, expected.length, 'returns expected value' );
 	for ( i = 0; i < expected.length; i++ ) {
-		if ( actual[ i ] === expected[ i ] ) {
-			t.strictEqual( actual[ i ], expected[ i ], 'returns expected value' );
-		} else {
-			delta = abs( actual[ i ] - expected[ i ] );
-			tol = rtol * EPS * abs( expected[ i ] );
-			t.ok( delta <= tol, 'within tolerance. actual: '+actual[ i ]+'. expected: '+expected[ i ]+'. delta: '+delta+'. tol: '+tol+'.' );
-		}
+		t.strictEqual( isAlmostSameValue( actual[ i ], expected[ i ], ulp ), true, 'returns expected value' );
 	}
 }
 
@@ -125,8 +116,8 @@ tape( 'the function applies a plane rotation', function test( t ) {
 	]);
 
 	out = zdrot( x.length, x, 1, 0, y, 1, 0, 0.8, 0.6 );
-	isApprox( t, viewX, xe, 2.0 );
-	isApprox( t, viewY, ye, 2.0 );
+	isApprox( t, viewX, xe, 1 );
+	isApprox( t, viewY, ye, 1 );
 	t.strictEqual( out, y, 'returns expected value' );
 	t.end();
 });
@@ -186,8 +177,8 @@ tape( 'the function supports an `x` stride', function test( t ) {
 	]);
 
 	out = zdrot( 2, x, 2, 0, y, 1, 0, 0.8, 0.6 );
-	isApprox( t, viewX, xe, 2.0 );
-	isApprox( t, viewY, ye, 2.0 );
+	isApprox( t, viewX, xe, 1 );
+	isApprox( t, viewY, ye, 1 );
 	t.strictEqual( out, y, 'returns expected value' );
 	t.end();
 });
@@ -247,8 +238,8 @@ tape( 'the function supports an `x` offset', function test( t ) {
 	]);
 
 	out = zdrot( 3, x, 1, 1, y, 1, 0, 0.8, 0.6 );
-	isApprox( t, viewX, xe, 2.0 );
-	isApprox( t, viewY, ye, 2.0 );
+	isApprox( t, viewX, xe, 1 );
+	isApprox( t, viewY, ye, 1 );
 	t.strictEqual( out, y, 'returns expected value' );
 	t.end();
 });
@@ -308,8 +299,8 @@ tape( 'the function supports a `y` stride', function test( t ) {
 	]);
 
 	out = zdrot( 2, x, 1, 0, y, 2, 0, 0.8, 0.6 );
-	isApprox( t, viewX, xe, 2.0 );
-	isApprox( t, viewY, ye, 2.0 );
+	isApprox( t, viewX, xe, 1 );
+	isApprox( t, viewY, ye, 1 );
 	t.strictEqual( out, y, 'returns expected value' );
 	t.end();
 });
@@ -369,8 +360,8 @@ tape( 'the function supports a `y` offset', function test( t ) {
 	]);
 
 	out = zdrot( 2, x, 2, 0, y, -1, 1, 0.8, 0.6 );
-	isApprox( t, viewX, xe, 2.0 );
-	isApprox( t, viewY, ye, 2.0 );
+	isApprox( t, viewX, xe, 1 );
+	isApprox( t, viewY, ye, 1 );
 	t.strictEqual( out, y, 'returns expected value' );
 	t.end();
 });
@@ -472,8 +463,8 @@ tape( 'the function supports negative strides', function test( t ) {
 	]);
 
 	out = zdrot( 2, x, -1, 1, y, -2, 2, 0.8, 0.6 );
-	isApprox( t, viewX, xe, 2.0 );
-	isApprox( t, viewY, ye, 2.0 );
+	isApprox( t, viewX, xe, 1 );
+	isApprox( t, viewY, ye, 1 );
 	t.strictEqual( out, y, 'returns expected value' );
 	t.end();
 });
@@ -533,8 +524,8 @@ tape( 'the function supports complex access patterns', function test( t ) {
 	]);
 
 	out = zdrot( 2, x, 2, 1, y, 1, 2, 0.8, 0.6 );
-	isApprox( t, viewX, xe, 2.0 );
-	isApprox( t, viewY, ye, 2.0 );
+	isApprox( t, viewX, xe, 1 );
+	isApprox( t, viewY, ye, 1 );
 	t.strictEqual( out, y, 'returns expected value' );
 	t.end();
 });
diff --git a/lib/node_modules/@stdlib/blas/base/zdrot/test/test.ndarray.native.js b/lib/node_modules/@stdlib/blas/base/zdrot/test/test.ndarray.native.js
index 30c5fe0f7b97..af2a09d69eb5 100644
--- a/lib/node_modules/@stdlib/blas/base/zdrot/test/test.ndarray.native.js
+++ b/lib/node_modules/@stdlib/blas/base/zdrot/test/test.ndarray.native.js
@@ -24,8 +24,7 @@ var resolve = require( 'path' ).resolve;
 var tape = require( 'tape' );
 var Float64Array = require( '@stdlib/array/float64' );
 var Complex128Array = require( '@stdlib/array/complex128' );
-var EPS = require( '@stdlib/constants/float64/eps' );
-var abs = require( '@stdlib/math/base/special/abs' );
+var isAlmostSameValue = require( '@stdlib/assert/is-almost-same-value' );
 var tryRequire = require( '@stdlib/utils/try-require' );
 
 
@@ -46,22 +45,14 @@ var opts = {
 * @param {Object} t - test object
 * @param {Collection} actual - actual values
 * @param {Collection} expected - expected values
-* @param {number} rtol - relative tolerance
+* @param {NonNegativeInteger} ulp - maximum allowed ULP difference
 */
-function isApprox( t, actual, expected, rtol ) {
-	var delta;
-	var tol;
+function isApprox( t, actual, expected, ulp ) {
 	var i;
 
 	t.strictEqual( actual.length, expected.length, 'returns expected value' );
 	for ( i = 0; i < expected.length; i++ ) {
-		if ( actual[ i ] === expected[ i ] ) {
-			t.strictEqual( actual[ i ], expected[ i ], 'returns expected value' );
-		} else {
-			delta = abs( actual[ i ] - expected[ i ] );
-			tol = rtol * EPS * abs( expected[ i ] );
-			t.ok( delta <= tol, 'within tolerance. actual: '+actual[ i ]+'. expected: '+expected[ i ]+'. delta: '+delta+'. tol: '+tol+'.' );
-		}
+		t.strictEqual( isAlmostSameValue( actual[ i ], expected[ i ], ulp ), true, 'returns expected value' );
 	}
 }
 
@@ -134,8 +125,8 @@ tape( 'the function applies a plane rotation', opts, function test( t ) {
 	]);
 
 	out = zdrot( x.length, x, 1, 0, y, 1, 0, 0.8, 0.6 );
-	isApprox( t, viewX, xe, 2.0 );
-	isApprox( t, viewY, ye, 2.0 );
+	isApprox( t, viewX, xe, 1 );
+	isApprox( t, viewY, ye, 1 );
 	t.strictEqual( out, y, 'returns expected value' );
 	t.end();
 });
@@ -195,8 +186,8 @@ tape( 'the function supports an `x` stride', opts, function test( t ) {
 	]);
 
 	out = zdrot( 2, x, 2, 0, y, 1, 0, 0.8, 0.6 );
-	isApprox( t, viewX, xe, 2.0 );
-	isApprox( t, viewY, ye, 2.0 );
+	isApprox( t, viewX, xe, 1 );
+	isApprox( t, viewY, ye, 1 );
 	t.strictEqual( out, y, 'returns expected value' );
 	t.end();
 });
@@ -256,8 +247,8 @@ tape( 'the function supports an `x` offset', opts, function test( t ) {
 	]);
 
 	out = zdrot( 3, x, 1, 1, y, 1, 0, 0.8, 0.6 );
-	isApprox( t, viewX, xe, 2.0 );
-	isApprox( t, viewY, ye, 2.0 );
+	isApprox( t, viewX, xe, 1 );
+	isApprox( t, viewY, ye, 1 );
 	t.strictEqual( out, y, 'returns expected value' );
 	t.end();
 });
@@ -317,8 +308,8 @@ tape( 'the function supports a `y` stride', opts, function test( t ) {
 	]);
 
 	out = zdrot( 2, x, 1, 0, y, 2, 0, 0.8, 0.6 );
-	isApprox( t, viewX, xe, 2.0 );
-	isApprox( t, viewY, ye, 2.0 );
+	isApprox( t, viewX, xe, 1 );
+	isApprox( t, viewY, ye, 1 );
 	t.strictEqual( out, y, 'returns expected value' );
 	t.end();
 });
@@ -378,8 +369,8 @@ tape( 'the function supports a `y` offset', opts, function test( t ) {
 	]);
 
 	out = zdrot( 2, x, 2, 0, y, -1, 1, 0.8, 0.6 );
-	isApprox( t, viewX, xe, 2.0 );
-	isApprox( t, viewY, ye, 2.0 );
+	isApprox( t, viewX, xe, 1 );
+	isApprox( t, viewY, ye, 1 );
 	t.strictEqual( out, y, 'returns expected value' );
 	t.end();
 });
@@ -481,8 +472,8 @@ tape( 'the function supports negative strides', opts, function test( t ) {
 	]);
 
 	out = zdrot( 2, x, -1, 1, y, -2, 2, 0.8, 0.6 );
-	isApprox( t, viewX, xe, 2.0 );
-	isApprox( t, viewY, ye, 2.0 );
+	isApprox( t, viewX, xe, 1 );
+	isApprox( t, viewY, ye, 1 );
 	t.strictEqual( out, y, 'returns expected value' );
 	t.end();
 });
@@ -542,8 +533,8 @@ tape( 'the function supports complex access patterns', opts, function test( t )
 	]);
 
 	out = zdrot( 2, x, 2, 1, y, 1, 2, 0.8, 0.6 );
-	isApprox( t, viewX, xe, 2.0 );
-	isApprox( t, viewY, ye, 2.0 );
+	isApprox( t, viewX, xe, 1 );
+	isApprox( t, viewY, ye, 1 );
 	t.strictEqual( out, y, 'returns expected value' );
 	t.end();
 });
diff --git a/lib/node_modules/@stdlib/blas/base/zdrot/test/test.zdrot.js b/lib/node_modules/@stdlib/blas/base/zdrot/test/test.zdrot.js
index 08dd719cc6b0..fdbd57b99b9f 100644
--- a/lib/node_modules/@stdlib/blas/base/zdrot/test/test.zdrot.js
+++ b/lib/node_modules/@stdlib/blas/base/zdrot/test/test.zdrot.js
@@ -23,8 +23,7 @@
 var tape = require( 'tape' );
 var Float64Array = require( '@stdlib/array/float64' );
 var Complex128Array = require( '@stdlib/array/complex128' );
-var EPS = require( '@stdlib/constants/float64/eps' );
-var abs = require( '@stdlib/math/base/special/abs' );
+var isAlmostSameValue = require( '@stdlib/assert/is-almost-same-value' );
 var zdrot = require( './../lib/zdrot.js' );
 
 
@@ -37,22 +36,14 @@ var zdrot = require( './../lib/zdrot.js' );
 * @param {Object} t - test object
 * @param {Collection} actual - actual values
 * @param {Collection} expected - expected values
-* @param {number} rtol - relative tolerance
+* @param {NonNegativeInteger} ulp - maximum allowed ULP difference
 */
-function isApprox( t, actual, expected, rtol ) {
-	var delta;
-	var tol;
+function isApprox( t, actual, expected, ulp ) {
 	var i;
 
 	t.strictEqual( actual.length, expected.length, 'returns expected value' );
 	for ( i = 0; i < expected.length; i++ ) {
-		if ( actual[ i ] === expected[ i ] ) {
-			t.strictEqual( actual[ i ], expected[ i ], 'returns expected value' );
-		} else {
-			delta = abs( actual[ i ] - expected[ i ] );
-			tol = rtol * EPS * abs( expected[ i ] );
-			t.ok( delta <= tol, 'within tolerance. actual: '+actual[ i ]+'. expected: '+expected[ i ]+'. delta: '+delta+'. tol: '+tol+'.' );
-		}
+		t.strictEqual( isAlmostSameValue( actual[ i ], expected[ i ], ulp ), true, 'returns expected value' );
 	}
 }
 
@@ -125,8 +116,8 @@ tape( 'the function applies a plane rotation', function test( t ) {
 	]);
 
 	out = zdrot( x.length, x, 1, y, 1, 0.8, 0.6 );
-	isApprox( t, viewX, xe, 2.0 );
-	isApprox( t, viewY, ye, 2.0 );
+	isApprox( t, viewX, xe, 1 );
+	isApprox( t, viewY, ye, 1 );
 	t.strictEqual( out, y, 'returns expected value' );
 	t.end();
 });
@@ -186,8 +177,8 @@ tape( 'the function supports an `x` stride', function test( t ) {
 	]);
 
 	out = zdrot( 2, x, 2, y, 1, 0.8, 0.6 );
-	isApprox( t, viewX, xe, 2.0 );
-	isApprox( t, viewY, ye, 2.0 );
+	isApprox( t, viewX, xe, 1 );
+	isApprox( t, viewY, ye, 1 );
 	t.strictEqual( out, y, 'returns expected value' );
 	t.end();
 });
@@ -247,8 +238,8 @@ tape( 'the function supports a `y` stride', function test( t ) {
 	]);
 
 	out = zdrot( 2, x, 1, y, 2, 0.8, 0.6 );
-	isApprox( t, viewX, xe, 2.0 );
-	isApprox( t, viewY, ye, 2.0 );
+	isApprox( t, viewX, xe, 1 );
+	isApprox( t, viewY, ye, 1 );
 	t.strictEqual( out, y, 'returns expected value' );
 	t.end();
 });
@@ -350,8 +341,8 @@ tape( 'the function supports negative strides', function test( t ) {
 	]);
 
 	out = zdrot( 2, x, -1, y, -2, 0.8, 0.6 );
-	isApprox( t, viewX, xe, 2.0 );
-	isApprox( t, viewY, ye, 2.0 );
+	isApprox( t, viewX, xe, 1 );
+	isApprox( t, viewY, ye, 1 );
 	t.strictEqual( out, y, 'returns expected value' );
 	t.end();
 });
@@ -411,8 +402,8 @@ tape( 'the function supports complex access patterns', function test( t ) {
 	]);
 
 	out = zdrot( 2, x, 2, y, -1, 0.8, 0.6 );
-	isApprox( t, viewX, xe, 2.0 );
-	isApprox( t, viewY, ye, 2.0 );
+	isApprox( t, viewX, xe, 1 );
+	isApprox( t, viewY, ye, 1 );
 	t.strictEqual( out, y, 'returns expected value' );
 	t.end();
 });
@@ -479,8 +470,8 @@ tape( 'the function supports view offsets', function test( t ) {
 	]);
 
 	out = zdrot( 2, x1, -2, y1, 1, 0.8, 0.6 );
-	isApprox( t, viewX, xe, 2.0 );
-	isApprox( t, viewY, ye, 2.0 );
+	isApprox( t, viewX, xe, 1 );
+	isApprox( t, viewY, ye, 1 );
 	t.strictEqual( out, y1, 'returns expected value' );
 	t.end();
 });
diff --git a/lib/node_modules/@stdlib/blas/base/zdrot/test/test.zdrot.native.js b/lib/node_modules/@stdlib/blas/base/zdrot/test/test.zdrot.native.js
index 4acdcfcde091..f6bbd253a12e 100644
--- a/lib/node_modules/@stdlib/blas/base/zdrot/test/test.zdrot.native.js
+++ b/lib/node_modules/@stdlib/blas/base/zdrot/test/test.zdrot.native.js
@@ -24,8 +24,7 @@ var resolve = require( 'path' ).resolve;
 var tape = require( 'tape' );
 var Float64Array = require( '@stdlib/array/float64' );
 var Complex128Array = require( '@stdlib/array/complex128' );
-var EPS = require( '@stdlib/constants/float64/eps' );
-var abs = require( '@stdlib/math/base/special/abs' );
+var isAlmostSameValue = require( '@stdlib/assert/is-almost-same-value' );
 var tryRequire = require( '@stdlib/utils/try-require' );
 
 
@@ -46,22 +45,14 @@ var opts = {
 * @param {Object} t - test object
 * @param {Collection} actual - actual values
 * @param {Collection} expected - expected values
-* @param {number} rtol - relative tolerance
+* @param {NonNegativeInteger} ulp - maximum allowed ULP difference
 */
-function isApprox( t, actual, expected, rtol ) {
-	var delta;
-	var tol;
+function isApprox( t, actual, expected, ulp ) {
 	var i;
 
 	t.strictEqual( actual.length, expected.length, 'returns expected value' );
 	for ( i = 0; i < expected.length; i++ ) {
-		if ( actual[ i ] === expected[ i ] ) {
-			t.strictEqual( actual[ i ], expected[ i ], 'returns expected value' );
-		} else {
-			delta = abs( actual[ i ] - expected[ i ] );
-			tol = rtol * EPS * abs( expected[ i ] );
-			t.ok( delta <= tol, 'within tolerance. actual: '+actual[ i ]+'. expected: '+expected[ i ]+'. delta: '+delta+'. tol: '+tol+'.' );
-		}
+		t.strictEqual( isAlmostSameValue( actual[ i ], expected[ i ], ulp ), true, 'returns expected value' );
 	}
 }
 
@@ -134,8 +125,8 @@ tape( 'the function applies a plane rotation', opts, function test( t ) {
 	]);
 
 	out = zdrot( x.length, x, 1, y, 1, 0.8, 0.6 );
-	isApprox( t, viewX, xe, 2.0 );
-	isApprox( t, viewY, ye, 2.0 );
+	isApprox( t, viewX, xe, 1 );
+	isApprox( t, viewY, ye, 1 );
 	t.strictEqual( out, y, 'returns expected value' );
 	t.end();
 });
@@ -195,8 +186,8 @@ tape( 'the function supports an `x` stride', opts, function test( t ) {
 	]);
 
 	out = zdrot( 2, x, 2, y, 1, 0.8, 0.6 );
-	isApprox( t, viewX, xe, 2.0 );
-	isApprox( t, viewY, ye, 2.0 );
+	isApprox( t, viewX, xe, 1 );
+	isApprox( t, viewY, ye, 1 );
 	t.strictEqual( out, y, 'returns expected value' );
 	t.end();
 });
@@ -256,8 +247,8 @@ tape( 'the function supports a `y` stride', opts, function test( t ) {
 	]);
 
 	out = zdrot( 2, x, 1, y, 2, 0.8, 0.6 );
-	isApprox( t, viewX, xe, 2.0 );
-	isApprox( t, viewY, ye, 2.0 );
+	isApprox( t, viewX, xe, 1 );
+	isApprox( t, viewY, ye, 1 );
 	t.strictEqual( out, y, 'returns expected value' );
 	t.end();
 });
@@ -359,8 +350,8 @@ tape( 'the function supports negative strides', opts, function test( t ) {
 	]);
 
 	out = zdrot( 2, x, -1, y, -2, 0.8, 0.6 );
-	isApprox( t, viewX, xe, 2.0 );
-	isApprox( t, viewY, ye, 2.0 );
+	isApprox( t, viewX, xe, 1 );
+	isApprox( t, viewY, ye, 1 );
 	t.strictEqual( out, y, 'returns expected value' );
 	t.end();
 });
@@ -420,8 +411,8 @@ tape( 'the function supports complex access patterns', opts, function test( t )
 	]);
 
 	out = zdrot( 2, x, 2, y, -1, 0.8, 0.6 );
-	isApprox( t, viewX, xe, 2.0 );
-	isApprox( t, viewY, ye, 2.0 );
+	isApprox( t, viewX, xe, 1 );
+	isApprox( t, viewY, ye, 1 );
 	t.strictEqual( out, y, 'returns expected value' );
 	t.end();
 });
@@ -488,8 +479,8 @@ tape( 'the function supports view offsets', opts, function test( t ) {
 	]);
 
 	out = zdrot( 2, x1, -2, y1, 1, 0.8, 0.6 );
-	isApprox( t, viewX, xe, 2.0 );
-	isApprox( t, viewY, ye, 2.0 );
+	isApprox( t, viewX, xe, 1 );
+	isApprox( t, viewY, ye, 1 );
 	t.strictEqual( out, y1, 'returns expected value' );
 	t.end();
 });