diff --git a/lib/node_modules/@stdlib/stats/base/dists/arcsine/mode/src/addon.c b/lib/node_modules/@stdlib/stats/base/dists/arcsine/mode/src/addon.c
index e3068087b26e..bda748ea2fd2 100644
--- a/lib/node_modules/@stdlib/stats/base/dists/arcsine/mode/src/addon.c
+++ b/lib/node_modules/@stdlib/stats/base/dists/arcsine/mode/src/addon.c
@@ -16,8 +16,8 @@
 * limitations under the License.
 */
 
-#include "stdlib/math/base/napi/binary.h"
 #include "stdlib/stats/base/dists/arcsine/mode.h"
+#include "stdlib/math/base/napi/binary.h"
 
 // cppcheck-suppress shadowFunction
 STDLIB_MATH_BASE_NAPI_MODULE_DD_D( stdlib_base_dists_arcsine_mode )
diff --git a/lib/node_modules/@stdlib/stats/base/dists/arcsine/mode/src/main.c b/lib/node_modules/@stdlib/stats/base/dists/arcsine/mode/src/main.c
index ae2fa8465f13..e681597c3673 100644
--- a/lib/node_modules/@stdlib/stats/base/dists/arcsine/mode/src/main.c
+++ b/lib/node_modules/@stdlib/stats/base/dists/arcsine/mode/src/main.c
@@ -16,8 +16,8 @@
 * limitations under the License.
 */
 
-#include "stdlib/math/base/assert/is_nan.h"
 #include "stdlib/stats/base/dists/arcsine/mode.h"
+#include "stdlib/math/base/assert/is_nan.h"
 
 /**
 * Returns the mode of an arcsine distribution.
diff --git a/lib/node_modules/@stdlib/stats/base/dists/arcsine/mode/test/test.js b/lib/node_modules/@stdlib/stats/base/dists/arcsine/mode/test/test.js
index 5dc055e15ac7..83491462682e 100644
--- a/lib/node_modules/@stdlib/stats/base/dists/arcsine/mode/test/test.js
+++ b/lib/node_modules/@stdlib/stats/base/dists/arcsine/mode/test/test.js
@@ -21,11 +21,9 @@
 // MODULES //
 
 var tape = require( 'tape' );
-var abs = require( '@stdlib/math/base/special/abs' );
 var isnan = require( '@stdlib/math/base/assert/is-nan' );
 var PINF = require( '@stdlib/constants/float64/pinf' );
 var NINF = require( '@stdlib/constants/float64/ninf' );
-var EPS = require( '@stdlib/constants/float64/eps' );
 var mode = require( './../lib' );
 
 
@@ -72,8 +70,6 @@ tape( 'if provided `a >= b`, the function returns `NaN`', function test( t ) {
 
 tape( 'the function returns the mode of an arcsine distribution', function test( t ) {
 	var expected;
-	var delta;
-	var tol;
 	var a;
 	var b;
 	var i;
@@ -84,13 +80,7 @@ tape( 'the function returns the mode of an arcsine distribution', function test(
 	b = data.b;
 	for ( i = 0; i < expected.length; i++ ) {
 		y = mode( a[i], b[i] );
-		if ( y === expected[i] ) {
-			t.equal( y, expected[i], 'a: '+a[i]+', b: '+b[i]+', y: '+y+', expected: '+expected[i] );
-		} else {
-			delta = abs( y - expected[ i ] );
-			tol = 1.0 * EPS * abs( expected[ i ] );
-			t.ok( delta <= tol, 'within tolerance. a: '+a[i]+'. b: '+b[i]+'. y: '+y+'. E: '+expected[ i ]+'. Δ: '+delta+'. tol: '+tol+'.' );
-		}
+		t.equal( y, expected[i], 'a: '+a[i]+', b: '+b[i]+', y: '+y+', expected: '+expected[i] );
 	}
 	t.end();
 });
diff --git a/lib/node_modules/@stdlib/stats/base/dists/arcsine/mode/test/test.native.js b/lib/node_modules/@stdlib/stats/base/dists/arcsine/mode/test/test.native.js
index e08d8d7506c2..4cf6ad8b1565 100644
--- a/lib/node_modules/@stdlib/stats/base/dists/arcsine/mode/test/test.native.js
+++ b/lib/node_modules/@stdlib/stats/base/dists/arcsine/mode/test/test.native.js
@@ -23,11 +23,9 @@
 var resolve = require( 'path' ).resolve;
 var tape = require( 'tape' );
 var tryRequire = require( '@stdlib/utils/try-require' );
-var abs = require( '@stdlib/math/base/special/abs' );
 var isnan = require( '@stdlib/math/base/assert/is-nan' );
 var PINF = require( '@stdlib/constants/float64/pinf' );
 var NINF = require( '@stdlib/constants/float64/ninf' );
-var EPS = require( '@stdlib/constants/float64/eps' );
 
 
 // VARIABLES //
@@ -81,8 +79,6 @@ tape( 'if provided `a >= b`, the function returns `NaN`', opts, function test( t
 
 tape( 'the function returns the mode of an arcsine distribution', opts, function test( t ) {
 	var expected;
-	var delta;
-	var tol;
 	var a;
 	var b;
 	var i;
@@ -93,13 +89,7 @@ tape( 'the function returns the mode of an arcsine distribution', opts, function
 	b = data.b;
 	for ( i = 0; i < expected.length; i++ ) {
 		y = mode( a[i], b[i] );
-		if ( y === expected[i] ) {
-			t.equal( y, expected[i], 'a: '+a[i]+', b: '+b[i]+', y: '+y+', expected: '+expected[i] );
-		} else {
-			delta = abs( y - expected[ i ] );
-			tol = 1.0 * EPS * abs( expected[ i ] );
-			t.ok( delta <= tol, 'within tolerance. a: '+a[i]+'. b: '+b[i]+'. y: '+y+'. E: '+expected[ i ]+'. Δ: '+delta+'. tol: '+tol+'.' );
-		}
+		t.equal( y, expected[i], 'a: '+a[i]+', b: '+b[i]+', y: '+y+', expected: '+expected[i] );
 	}
 	t.end();
 });
diff --git a/lib/node_modules/@stdlib/stats/base/dists/arcsine/skewness/benchmark/c/benchmark.c b/lib/node_modules/@stdlib/stats/base/dists/arcsine/skewness/benchmark/c/benchmark.c
index 2b6eef8c0cf7..98837da2fd1a 100644
--- a/lib/node_modules/@stdlib/stats/base/dists/arcsine/skewness/benchmark/c/benchmark.c
+++ b/lib/node_modules/@stdlib/stats/base/dists/arcsine/skewness/benchmark/c/benchmark.c
@@ -19,6 +19,7 @@
 #include "stdlib/stats/base/dists/arcsine/skewness.h"
 #include <stdio.h>
 #include <stdlib.h>
+#include <math.h>
 #include <time.h>
 #include <sys/time.h>
 
diff --git a/lib/node_modules/@stdlib/stats/base/dists/arcsine/skewness/src/addon.c b/lib/node_modules/@stdlib/stats/base/dists/arcsine/skewness/src/addon.c
index 3d68bbb0876c..ca2920f656ec 100644
--- a/lib/node_modules/@stdlib/stats/base/dists/arcsine/skewness/src/addon.c
+++ b/lib/node_modules/@stdlib/stats/base/dists/arcsine/skewness/src/addon.c
@@ -16,8 +16,8 @@
 * limitations under the License.
 */
 
-#include "stdlib/math/base/napi/binary.h"
 #include "stdlib/stats/base/dists/arcsine/skewness.h"
+#include "stdlib/math/base/napi/binary.h"
 
 // cppcheck-suppress shadowFunction
 STDLIB_MATH_BASE_NAPI_MODULE_DD_D( stdlib_base_dists_arcsine_skewness )
diff --git a/lib/node_modules/@stdlib/stats/base/dists/arcsine/skewness/src/main.c b/lib/node_modules/@stdlib/stats/base/dists/arcsine/skewness/src/main.c
index 944b3a8973e8..522b83d33321 100644
--- a/lib/node_modules/@stdlib/stats/base/dists/arcsine/skewness/src/main.c
+++ b/lib/node_modules/@stdlib/stats/base/dists/arcsine/skewness/src/main.c
@@ -16,8 +16,8 @@
 * limitations under the License.
 */
 
-#include "stdlib/math/base/assert/is_nan.h"
 #include "stdlib/stats/base/dists/arcsine/skewness.h"
+#include "stdlib/math/base/assert/is_nan.h"
 
 /**
 * Returns the skewness of an arcsine distribution.
diff --git a/lib/node_modules/@stdlib/stats/base/dists/arcsine/stdev/src/main.c b/lib/node_modules/@stdlib/stats/base/dists/arcsine/stdev/src/main.c
index f06cb2fd387a..844087f53b4f 100644
--- a/lib/node_modules/@stdlib/stats/base/dists/arcsine/stdev/src/main.c
+++ b/lib/node_modules/@stdlib/stats/base/dists/arcsine/stdev/src/main.c
@@ -19,7 +19,7 @@
 #include "stdlib/stats/base/dists/arcsine/stdev.h"
 #include "stdlib/math/base/assert/is_nan.h"
 
-#define SQRT1OVER8 0.35355339059327373 // sqrt(1.0 / 8.0)
+static const double SQRT1OVER8 = 0.35355339059327373 // sqrt(1.0 / 8.0)
 
 /**
 * Returns the standard deviation of an arcsine distribution.
diff --git a/lib/node_modules/@stdlib/stats/base/dists/cosine/mode/manifest.json b/lib/node_modules/@stdlib/stats/base/dists/cosine/mode/manifest.json
index 93ccded754bf..86be782a5d98 100644
--- a/lib/node_modules/@stdlib/stats/base/dists/cosine/mode/manifest.json
+++ b/lib/node_modules/@stdlib/stats/base/dists/cosine/mode/manifest.json
@@ -39,8 +39,7 @@
       "libpath": [],
       "dependencies": [
         "@stdlib/math/base/napi/binary",
-        "@stdlib/math/base/assert/is-nan",
-        "@stdlib/constants/float64/eps"
+        "@stdlib/math/base/assert/is-nan"
       ]
     },
     {
diff --git a/lib/node_modules/@stdlib/stats/base/dists/cosine/mode/test/test.js b/lib/node_modules/@stdlib/stats/base/dists/cosine/mode/test/test.js
index 6ea6651625f2..64f4f5616359 100644
--- a/lib/node_modules/@stdlib/stats/base/dists/cosine/mode/test/test.js
+++ b/lib/node_modules/@stdlib/stats/base/dists/cosine/mode/test/test.js
@@ -22,10 +22,8 @@
 
 var tape = require( 'tape' );
 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 mode = require( './../lib' );
 
 
@@ -79,8 +77,6 @@ tape( 'if provided a nonpositive `s`, the function returns `NaN`', function test
 
 tape( 'the function returns the mode of a raised cosine distribution', function test( t ) {
 	var expected;
-	var delta;
-	var tol;
 	var mu;
 	var s;
 	var y;
@@ -92,13 +88,7 @@ tape( 'the function returns the mode of a raised cosine distribution', function
 	for ( i = 0; i < mu.length; i++ ) {
 		y = mode( mu[i], s[i] );
 		if ( expected[i] !== null ) {
-			if ( y === expected[i] ) {
-				t.equal( y, expected[i], 'mu:'+mu[i]+', s: '+s[i]+', y: '+y+', expected: '+expected[i] );
-			} else {
-				delta = abs( y - expected[ i ] );
-				tol = 1.0 * EPS * abs( expected[ i ] );
-				t.ok( delta <= tol, 'within tolerance. mu: '+mu[i]+'. s: '+s[i]+'. y: '+y+'. E: '+expected[ i ]+'. Δ: '+delta+'. tol: '+tol+'.' );
-			}
+			t.equal( y, expected[i], 'mu:'+mu[i]+', s: '+s[i]+', y: '+y+', expected: '+expected[i] );
 		}
 	}
 	t.end();
diff --git a/lib/node_modules/@stdlib/stats/base/dists/cosine/mode/test/test.native.js b/lib/node_modules/@stdlib/stats/base/dists/cosine/mode/test/test.native.js
index 2bf0223d7b7c..85d41162ee02 100644
--- a/lib/node_modules/@stdlib/stats/base/dists/cosine/mode/test/test.native.js
+++ b/lib/node_modules/@stdlib/stats/base/dists/cosine/mode/test/test.native.js
@@ -23,10 +23,8 @@
 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 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' );
 
 
@@ -88,8 +86,6 @@ tape( 'if provided a nonpositive `s`, the function returns `NaN`', opts, functio
 
 tape( 'the function returns the mode of a raised cosine distribution', opts, function test( t ) {
 	var expected;
-	var delta;
-	var tol;
 	var mu;
 	var s;
 	var y;
@@ -101,13 +97,7 @@ tape( 'the function returns the mode of a raised cosine distribution', opts, fun
 	for ( i = 0; i < mu.length; i++ ) {
 		y = mode( mu[i], s[i] );
 		if ( expected[i] !== null ) {
-			if ( y === expected[i] ) {
-				t.equal( y, expected[i], 'mu:'+mu[i]+', s: '+s[i]+', y: '+y+', expected: '+expected[i] );
-			} else {
-				delta = abs( y - expected[ i ] );
-				tol = 1.0 * EPS * abs( expected[ i ] );
-				t.ok( delta <= tol, 'within tolerance. mu: '+mu[i]+'. s: '+s[i]+'. y: '+y+'. E: '+expected[ i ]+'. Δ: '+delta+'. tol: '+tol+'.' );
-			}
+			t.equal( y, expected[i], 'mu:'+mu[i]+', s: '+s[i]+', y: '+y+', expected: '+expected[i] );
 		}
 	}
 	t.end();