diff --git a/lib/node_modules/@stdlib/stats/base/dists/chi/logpdf/benchmark/c/benchmark.c b/lib/node_modules/@stdlib/stats/base/dists/chi/logpdf/benchmark/c/benchmark.c
index 902a1285f1f3..28074707a968 100644
--- a/lib/node_modules/@stdlib/stats/base/dists/chi/logpdf/benchmark/c/benchmark.c
+++ b/lib/node_modules/@stdlib/stats/base/dists/chi/logpdf/benchmark/c/benchmark.c
@@ -17,6 +17,7 @@
 */
 
 #include "stdlib/stats/base/dists/chi/logpdf.h"
+#include "stdlib/constants/float64/eps.h"
 #include <sys/time.h>
 #include <math.h>
 #include <stdio.h>
@@ -101,7 +102,7 @@ static double benchmark( void ) {
 
 	for ( i = 0; i < 100; i++ ) {
 		x[ i ] = random_uniform( 0.0, 10.0 );
-		k[ i ] = random_uniform( 0.1, 10.0 );
+		k[ i ] = random_uniform( STDLIB_CONSTANT_FLOAT64_EPS, 10.0 );
 	}
 
 	t = tic();
diff --git a/lib/node_modules/@stdlib/stats/base/dists/chi/pdf/benchmark/c/benchmark.c b/lib/node_modules/@stdlib/stats/base/dists/chi/pdf/benchmark/c/benchmark.c
index 26d0f5be137c..49c7b5d86a93 100644
--- a/lib/node_modules/@stdlib/stats/base/dists/chi/pdf/benchmark/c/benchmark.c
+++ b/lib/node_modules/@stdlib/stats/base/dists/chi/pdf/benchmark/c/benchmark.c
@@ -17,6 +17,7 @@
 */
 
 #include "stdlib/stats/base/dists/chi/pdf.h"
+#include "stdlib/constants/float64/eps.h"
 #include <sys/time.h>
 #include <math.h>
 #include <stdio.h>
@@ -101,7 +102,7 @@ static double benchmark( void ) {
 
 	for ( i = 0; i < 100; i++ ) {
 		x[ i ] = random_uniform( 0.0, 10.0 );
-		k[ i ] = random_uniform( 0.1, 10.0 );
+		k[ i ] = random_uniform( STDLIB_CONSTANT_FLOAT64_EPS, 10.0 );
 	}
 
 	t = tic();
diff --git a/lib/node_modules/@stdlib/stats/base/dists/chisquare/mgf/benchmark/c/benchmark.c b/lib/node_modules/@stdlib/stats/base/dists/chisquare/mgf/benchmark/c/benchmark.c
index 80bcd9e9a958..6d41a17f255b 100644
--- a/lib/node_modules/@stdlib/stats/base/dists/chisquare/mgf/benchmark/c/benchmark.c
+++ b/lib/node_modules/@stdlib/stats/base/dists/chisquare/mgf/benchmark/c/benchmark.c
@@ -17,6 +17,7 @@
 */
 
 #include "stdlib/stats/base/dists/chisquare/mgf.h"
+#include "stdlib/constants/float64/eps.h"
 #include <math.h>
 #include <stdio.h>
 #include <stdlib.h>
@@ -101,7 +102,7 @@ static double benchmark( void ) {
 
 	for ( i = 0; i < 100; i++ ) {
 		t[ i ] = random_uniform( -0.4, 0.4 );
-		k[ i ] = random_uniform( 0.1, 10.0 );
+		k[ i ] = random_uniform( STDLIB_CONSTANT_FLOAT64_EPS, 10.0 );
 	}
 
 	start = tic();
diff --git a/lib/node_modules/@stdlib/stats/base/dists/invgamma/cdf/benchmark/c/benchmark.c b/lib/node_modules/@stdlib/stats/base/dists/invgamma/cdf/benchmark/c/benchmark.c
index 139e6f8e6b4b..4f2b691ff659 100644
--- a/lib/node_modules/@stdlib/stats/base/dists/invgamma/cdf/benchmark/c/benchmark.c
+++ b/lib/node_modules/@stdlib/stats/base/dists/invgamma/cdf/benchmark/c/benchmark.c
@@ -17,6 +17,7 @@
 */
 
 #include "stdlib/stats/base/dists/invgamma/cdf.h"
+#include "stdlib/constants/float64/eps.h"
 #include <stdlib.h>
 #include <stdio.h>
 #include <math.h>
@@ -103,7 +104,7 @@ static double benchmark( void ) {
 	for ( i = 0; i < 100; i++ ) {
 		alpha[ i ] = random_uniform( 1.0, 10.0 );
 		beta[ i ] = random_uniform( 1.0, 10.0 );
-		x[ i ] = random_uniform( 0.1, 5.0 );
+		x[ i ] = random_uniform( STDLIB_CONSTANT_FLOAT64_EPS, 5.0 );
 	}
 
 	t = tic();
diff --git a/lib/node_modules/@stdlib/stats/base/dists/levy/cdf/benchmark/c/benchmark.c b/lib/node_modules/@stdlib/stats/base/dists/levy/cdf/benchmark/c/benchmark.c
index 5de81aee2ce3..3964c1273357 100644
--- a/lib/node_modules/@stdlib/stats/base/dists/levy/cdf/benchmark/c/benchmark.c
+++ b/lib/node_modules/@stdlib/stats/base/dists/levy/cdf/benchmark/c/benchmark.c
@@ -17,6 +17,7 @@
 */
 
 #include "stdlib/stats/base/dists/levy/cdf.h"
+#include "stdlib/constants/float64/eps.h"
 #include <sys/time.h>
 #include <math.h>
 #include <stdio.h>
@@ -103,7 +104,7 @@ static double benchmark( void ) {
 	for ( i = 0; i < 100; i++ ) {
 		x[ i ] = random_uniform( -10.0, 10.0 );
 		mu[ i ] = random_uniform( -20.0, 0.0 );
-		c[ i ] = random_uniform( 0.1, 10.0 );
+		c[ i ] = random_uniform( STDLIB_CONSTANT_FLOAT64_EPS, 10.0 );
 	}
 
 	t = tic();
diff --git a/lib/node_modules/@stdlib/stats/base/dists/lognormal/pdf/benchmark/c/benchmark.c b/lib/node_modules/@stdlib/stats/base/dists/lognormal/pdf/benchmark/c/benchmark.c
index 1445a3efe20d..562d3bd9c6bb 100644
--- a/lib/node_modules/@stdlib/stats/base/dists/lognormal/pdf/benchmark/c/benchmark.c
+++ b/lib/node_modules/@stdlib/stats/base/dists/lognormal/pdf/benchmark/c/benchmark.c
@@ -102,7 +102,7 @@ static double benchmark( void ) {
 	int i;
 
 	for ( i = 0; i < 100; i++ ) {
-		x[ i ] = random_uniform( 0.1, 20.0 );
+		x[ i ] = random_uniform( STDLIB_CONSTANT_FLOAT64_EPS, 20.0 );
 		mu[ i ] = random_uniform( -5.0, 5.0 );
 		sigma[ i ] = random_uniform( STDLIB_CONSTANT_FLOAT64_EPS, 5.0 );
 	}
diff --git a/lib/node_modules/@stdlib/stats/base/dists/negative-binomial/kurtosis/benchmark/c/benchmark.c b/lib/node_modules/@stdlib/stats/base/dists/negative-binomial/kurtosis/benchmark/c/benchmark.c
index 9c3e223ec6f8..25ab7b8864f2 100644
--- a/lib/node_modules/@stdlib/stats/base/dists/negative-binomial/kurtosis/benchmark/c/benchmark.c
+++ b/lib/node_modules/@stdlib/stats/base/dists/negative-binomial/kurtosis/benchmark/c/benchmark.c
@@ -17,6 +17,7 @@
 */
 
 #include "stdlib/stats/base/dists/negative-binomial/kurtosis.h"
+#include "stdlib/constants/float64/eps.h"
 #include <stdlib.h>
 #include <stdio.h>
 #include <math.h>
@@ -101,7 +102,7 @@ static double benchmark( void ) {
 
 	for ( i = 0; i < 100; i++ ) {
 		r[ i ] = random_uniform( 1, 100 );
-		p[ i ] = random_uniform( 0.1, 1.0 );
+		p[ i ] = random_uniform( STDLIB_CONSTANT_FLOAT64_EPS, 1.0 );
 	}
 
 	t = tic();
diff --git a/lib/node_modules/@stdlib/stats/base/dists/negative-binomial/mean/benchmark/c/benchmark.c b/lib/node_modules/@stdlib/stats/base/dists/negative-binomial/mean/benchmark/c/benchmark.c
index 06ca0c818352..be4d91449c5c 100644
--- a/lib/node_modules/@stdlib/stats/base/dists/negative-binomial/mean/benchmark/c/benchmark.c
+++ b/lib/node_modules/@stdlib/stats/base/dists/negative-binomial/mean/benchmark/c/benchmark.c
@@ -17,6 +17,7 @@
 */
 
 #include "stdlib/stats/base/dists/negative-binomial/mean.h"
+#include "stdlib/constants/float64/eps.h"
 #include <stdlib.h>
 #include <stdio.h>
 #include <math.h>
@@ -101,7 +102,7 @@ static double benchmark( void ) {
 
 	for ( i = 0; i < 100; i++ ) {
 		r[ i ] = random_uniform( 1, 100 );
-		p[ i ] = random_uniform( 0.1, 1.0 );
+		p[ i ] = random_uniform( STDLIB_CONSTANT_FLOAT64_EPS, 1.0 );
 	}
 
 	t = tic();
diff --git a/lib/node_modules/@stdlib/stats/base/dists/negative-binomial/mgf/benchmark/c/benchmark.c b/lib/node_modules/@stdlib/stats/base/dists/negative-binomial/mgf/benchmark/c/benchmark.c
index 051cec45b21c..4c1e0096de37 100644
--- a/lib/node_modules/@stdlib/stats/base/dists/negative-binomial/mgf/benchmark/c/benchmark.c
+++ b/lib/node_modules/@stdlib/stats/base/dists/negative-binomial/mgf/benchmark/c/benchmark.c
@@ -17,6 +17,7 @@
 */
 
 #include "stdlib/stats/base/dists/negative-binomial/mgf.h"
+#include "stdlib/constants/float64/eps.h"
 #include <stdlib.h>
 #include <stdio.h>
 #include <math.h>
@@ -103,7 +104,7 @@ static double benchmark( void ) {
 	for ( i = 0; i < 100; i++ ) {
 		t[ i ] = random_uniform( -2.0, 0.0 );
 		r[ i ] = random_uniform( 1.0, 100.0 );
-		p[ i ] = random_uniform( 0.1, 1.0 );
+		p[ i ] = random_uniform( STDLIB_CONSTANT_FLOAT64_EPS, 1.0 );
 	}
 
 	time = tic();
diff --git a/lib/node_modules/@stdlib/stats/base/dists/negative-binomial/mode/benchmark/c/benchmark.c b/lib/node_modules/@stdlib/stats/base/dists/negative-binomial/mode/benchmark/c/benchmark.c
index aa34d63b33e4..0783f089337d 100644
--- a/lib/node_modules/@stdlib/stats/base/dists/negative-binomial/mode/benchmark/c/benchmark.c
+++ b/lib/node_modules/@stdlib/stats/base/dists/negative-binomial/mode/benchmark/c/benchmark.c
@@ -17,6 +17,7 @@
 */
 
 #include "stdlib/stats/base/dists/negative-binomial/mode.h"
+#include "stdlib/constants/float64/eps.h"
 #include <stdlib.h>
 #include <stdio.h>
 #include <math.h>
@@ -101,7 +102,7 @@ static double benchmark( void ) {
 
 	for ( i = 0; i < 100; i++ ) {
 		r[ i ] = random_uniform( 1, 100 );
-		p[ i ] = random_uniform( 0.1, 1.0 );
+		p[ i ] = random_uniform( STDLIB_CONSTANT_FLOAT64_EPS, 1.0 );
 	}
 
 	t = tic();
diff --git a/lib/node_modules/@stdlib/stats/base/dists/negative-binomial/skewness/benchmark/c/benchmark.c b/lib/node_modules/@stdlib/stats/base/dists/negative-binomial/skewness/benchmark/c/benchmark.c
index 0ec6d9856fa4..9a7deab2ac60 100644
--- a/lib/node_modules/@stdlib/stats/base/dists/negative-binomial/skewness/benchmark/c/benchmark.c
+++ b/lib/node_modules/@stdlib/stats/base/dists/negative-binomial/skewness/benchmark/c/benchmark.c
@@ -17,6 +17,7 @@
 */
 
 #include "stdlib/stats/base/dists/negative-binomial/skewness.h"
+#include "stdlib/constants/float64/eps.h"
 #include <stdlib.h>
 #include <stdio.h>
 #include <math.h>
@@ -101,7 +102,7 @@ static double benchmark( void ) {
 
 	for ( i = 0; i < 100; i++ ) {
 		r[ i ] = random_uniform( 1, 100 );
-		p[ i ] = random_uniform( 0.1, 1.0 );
+		p[ i ] = random_uniform( STDLIB_CONSTANT_FLOAT64_EPS, 1.0 );
 	}
 
 	t = tic();
diff --git a/lib/node_modules/@stdlib/stats/base/dists/negative-binomial/stdev/benchmark/c/benchmark.c b/lib/node_modules/@stdlib/stats/base/dists/negative-binomial/stdev/benchmark/c/benchmark.c
index 84cbac0ccfc5..e42ed21016a6 100644
--- a/lib/node_modules/@stdlib/stats/base/dists/negative-binomial/stdev/benchmark/c/benchmark.c
+++ b/lib/node_modules/@stdlib/stats/base/dists/negative-binomial/stdev/benchmark/c/benchmark.c
@@ -17,6 +17,7 @@
 */
 
 #include "stdlib/stats/base/dists/negative-binomial/stdev.h"
+#include "stdlib/constants/float64/eps.h"
 #include <stdlib.h>
 #include <stdio.h>
 #include <math.h>
@@ -101,7 +102,7 @@ static double benchmark( void ) {
 
 	for ( i = 0; i < 100; i++ ) {
 		r[ i ] = random_uniform( 1.0, 100.0 );
-		p[ i ] = random_uniform( 0.1, 1.0 );
+		p[ i ] = random_uniform( STDLIB_CONSTANT_FLOAT64_EPS, 1.0 );
 	}
 
 	t = tic();
diff --git a/lib/node_modules/@stdlib/stats/base/dists/negative-binomial/variance/benchmark/c/benchmark.c b/lib/node_modules/@stdlib/stats/base/dists/negative-binomial/variance/benchmark/c/benchmark.c
index fd1f203e16c6..22abc39a9c62 100644
--- a/lib/node_modules/@stdlib/stats/base/dists/negative-binomial/variance/benchmark/c/benchmark.c
+++ b/lib/node_modules/@stdlib/stats/base/dists/negative-binomial/variance/benchmark/c/benchmark.c
@@ -17,6 +17,7 @@
 */
 
 #include "stdlib/stats/base/dists/negative-binomial/variance.h"
+#include "stdlib/constants/float64/eps.h"
 #include <sys/time.h>
 #include <math.h>
 #include <stdio.h>
@@ -101,7 +102,7 @@ static double benchmark( void ) {
 
 	for ( i = 0; i < 100; i++ ) {
 		r[ i ] = random_uniform( 1.0, 10.0 );
-		p[ i ] = random_uniform( 0.1, 0.9 );
+		p[ i ] = random_uniform( STDLIB_CONSTANT_FLOAT64_EPS, 0.9 );
 	}
 
 	t = tic();
diff --git a/lib/node_modules/@stdlib/stats/base/dists/normal/cdf/benchmark/c/benchmark.c b/lib/node_modules/@stdlib/stats/base/dists/normal/cdf/benchmark/c/benchmark.c
index 9d09ab1c6de7..5d1f0b3deac8 100644
--- a/lib/node_modules/@stdlib/stats/base/dists/normal/cdf/benchmark/c/benchmark.c
+++ b/lib/node_modules/@stdlib/stats/base/dists/normal/cdf/benchmark/c/benchmark.c
@@ -17,6 +17,7 @@
 */
 
 #include "stdlib/stats/base/dists/normal/cdf.h"
+#include "stdlib/constants/float64/eps.h"
 #include <stdlib.h>
 #include <stdio.h>
 #include <math.h>
@@ -104,7 +105,7 @@ static double benchmark( void ) {
 	for ( i = 0; i < LEN; i++ ) {
 		x[ i ] = random_uniform( -100.0, 100.0 );
 		mu[ i ] = random_uniform( -50.0, 50.0 );
-		sigma[ i ] = random_uniform( 0.1, 20.0 );
+		sigma[ i ] = random_uniform( STDLIB_CONSTANT_FLOAT64_EPS, 20.0 );
 	}
 
 	t = tic();
diff --git a/lib/node_modules/@stdlib/stats/base/dists/normal/entropy/benchmark/c/benchmark.c b/lib/node_modules/@stdlib/stats/base/dists/normal/entropy/benchmark/c/benchmark.c
index fdc2b9ef291b..c969deec8643 100644
--- a/lib/node_modules/@stdlib/stats/base/dists/normal/entropy/benchmark/c/benchmark.c
+++ b/lib/node_modules/@stdlib/stats/base/dists/normal/entropy/benchmark/c/benchmark.c
@@ -17,6 +17,7 @@
 */
 
 #include "stdlib/stats/base/dists/normal/entropy.h"
+#include "stdlib/constants/float64/eps.h"
 #include <stdlib.h>
 #include <stdio.h>
 #include <math.h>
@@ -101,7 +102,7 @@ static double benchmark( void ) {
 
 	for ( i = 0; i < 100; i++ ) {
 		mu[ i ] = random_uniform( -10.0, 10.0 );
-		sigma[ i ] = random_uniform( 0.1, 10.0 );
+		sigma[ i ] = random_uniform( STDLIB_CONSTANT_FLOAT64_EPS, 10.0 );
 	}
 
 	t = tic();
diff --git a/lib/node_modules/@stdlib/stats/base/dists/normal/kurtosis/benchmark/c/benchmark.c b/lib/node_modules/@stdlib/stats/base/dists/normal/kurtosis/benchmark/c/benchmark.c
index 0d17dea1d66b..d1f33fe3019c 100644
--- a/lib/node_modules/@stdlib/stats/base/dists/normal/kurtosis/benchmark/c/benchmark.c
+++ b/lib/node_modules/@stdlib/stats/base/dists/normal/kurtosis/benchmark/c/benchmark.c
@@ -17,6 +17,7 @@
 */
 
 #include "stdlib/stats/base/dists/normal/kurtosis.h"
+#include "stdlib/constants/float64/eps.h"
 #include <stdlib.h>
 #include <stdio.h>
 #include <math.h>
@@ -101,7 +102,7 @@ static double benchmark( void ) {
 
 	for ( i = 0; i < 100; i++ ) {
 		mu[ i ] = random_uniform( -10.0, 10.0 );
-		sigma[ i ] = random_uniform( 0.1, 10.0 );
+		sigma[ i ] = random_uniform( STDLIB_CONSTANT_FLOAT64_EPS, 10.0 );
 	}
 
 	t = tic();
diff --git a/lib/node_modules/@stdlib/stats/base/dists/normal/logcdf/benchmark/c/benchmark.c b/lib/node_modules/@stdlib/stats/base/dists/normal/logcdf/benchmark/c/benchmark.c
index 6f29d8d21d51..b115de6d0fed 100644
--- a/lib/node_modules/@stdlib/stats/base/dists/normal/logcdf/benchmark/c/benchmark.c
+++ b/lib/node_modules/@stdlib/stats/base/dists/normal/logcdf/benchmark/c/benchmark.c
@@ -17,6 +17,7 @@
 */
 
 #include "stdlib/stats/base/dists/normal/logcdf.h"
+#include "stdlib/constants/float64/eps.h"
 #include <stdlib.h>
 #include <stdio.h>
 #include <math.h>
@@ -104,7 +105,7 @@ static double benchmark( void ) {
 	for ( i = 0; i < LEN; i++ ) {
 		x[ i ] = random_uniform( -100.0, 100.0 );
 		mu[ i ] = random_uniform( -50.0, 50.0 );
-		sigma[ i ] = random_uniform( 0.1, 20.0 );
+		sigma[ i ] = random_uniform( STDLIB_CONSTANT_FLOAT64_EPS, 20.0 );
 	}
 
 	t = tic();
diff --git a/lib/node_modules/@stdlib/stats/base/dists/normal/mean/benchmark/c/benchmark.c b/lib/node_modules/@stdlib/stats/base/dists/normal/mean/benchmark/c/benchmark.c
index f76d58ea3625..136ef9ca288f 100644
--- a/lib/node_modules/@stdlib/stats/base/dists/normal/mean/benchmark/c/benchmark.c
+++ b/lib/node_modules/@stdlib/stats/base/dists/normal/mean/benchmark/c/benchmark.c
@@ -17,6 +17,7 @@
 */
 
 #include "stdlib/stats/base/dists/normal/mean.h"
+#include "stdlib/constants/float64/eps.h"
 #include <stdlib.h>
 #include <stdio.h>
 #include <math.h>
@@ -101,7 +102,7 @@ static double benchmark( void ) {
 
 	for ( i = 0; i < 100; i++ ) {
 		mu[ i ] = random_uniform( -10.0, 10.0 );
-		sigma[ i ] = random_uniform( 0.1, 10.0 );
+		sigma[ i ] = random_uniform( STDLIB_CONSTANT_FLOAT64_EPS, 10.0 );
 	}
 
 	t = tic();
diff --git a/lib/node_modules/@stdlib/stats/base/dists/normal/median/benchmark/c/benchmark.c b/lib/node_modules/@stdlib/stats/base/dists/normal/median/benchmark/c/benchmark.c
index 0e6c968dc343..af4faf94ea70 100644
--- a/lib/node_modules/@stdlib/stats/base/dists/normal/median/benchmark/c/benchmark.c
+++ b/lib/node_modules/@stdlib/stats/base/dists/normal/median/benchmark/c/benchmark.c
@@ -17,6 +17,7 @@
 */
 
 #include "stdlib/stats/base/dists/normal/median.h"
+#include "stdlib/constants/float64/eps.h"
 #include <stdlib.h>
 #include <stdio.h>
 #include <math.h>
@@ -101,7 +102,7 @@ static double benchmark( void ) {
 
 	for ( i = 0; i < 100; i++ ) {
 		mu[ i ] = random_uniform( -10.0, 10.0 );
-		sigma[ i ] = random_uniform( 0.1, 10.0 );
+		sigma[ i ] = random_uniform( STDLIB_CONSTANT_FLOAT64_EPS, 10.0 );
 	}
 
 	t = tic();
diff --git a/lib/node_modules/@stdlib/stats/base/dists/normal/mode/benchmark/c/benchmark.c b/lib/node_modules/@stdlib/stats/base/dists/normal/mode/benchmark/c/benchmark.c
index 888e09a4f8eb..236ebe3c1494 100644
--- a/lib/node_modules/@stdlib/stats/base/dists/normal/mode/benchmark/c/benchmark.c
+++ b/lib/node_modules/@stdlib/stats/base/dists/normal/mode/benchmark/c/benchmark.c
@@ -17,6 +17,7 @@
 */
 
 #include "stdlib/stats/base/dists/normal/mode.h"
+#include "stdlib/constants/float64/eps.h"
 #include <stdlib.h>
 #include <stdio.h>
 #include <math.h>
@@ -101,7 +102,7 @@ static double benchmark( void ) {
 
 	for ( i = 0; i < 100; i++ ) {
 		mu[ i ] = random_uniform( -10.0, 10.0 );
-		sigma[ i ] = random_uniform( 0.1, 10.0 );
+		sigma[ i ] = random_uniform( STDLIB_CONSTANT_FLOAT64_EPS, 10.0 );
 	}
 
 	t = tic();
diff --git a/lib/node_modules/@stdlib/stats/base/dists/normal/skewness/benchmark/c/benchmark.c b/lib/node_modules/@stdlib/stats/base/dists/normal/skewness/benchmark/c/benchmark.c
index ab9fe5dbf2dc..5173b38403ad 100644
--- a/lib/node_modules/@stdlib/stats/base/dists/normal/skewness/benchmark/c/benchmark.c
+++ b/lib/node_modules/@stdlib/stats/base/dists/normal/skewness/benchmark/c/benchmark.c
@@ -17,6 +17,7 @@
 */
 
 #include "stdlib/stats/base/dists/normal/skewness.h"
+#include "stdlib/constants/float64/eps.h"
 #include <stdlib.h>
 #include <stdio.h>
 #include <math.h>
@@ -101,7 +102,7 @@ static double benchmark( void ) {
 
 	for ( i = 0; i < 100; i++ ) {
 		mu[ i ] = random_uniform( -10.0, 10.0 );
-		sigma[ i ] = random_uniform( 0.1, 10.0 );
+		sigma[ i ] = random_uniform( STDLIB_CONSTANT_FLOAT64_EPS, 10.0 );
 	}
 
 	t = tic();
diff --git a/lib/node_modules/@stdlib/stats/base/dists/normal/stdev/benchmark/c/benchmark.c b/lib/node_modules/@stdlib/stats/base/dists/normal/stdev/benchmark/c/benchmark.c
index ab739d2301c5..abedc3ae01b6 100644
--- a/lib/node_modules/@stdlib/stats/base/dists/normal/stdev/benchmark/c/benchmark.c
+++ b/lib/node_modules/@stdlib/stats/base/dists/normal/stdev/benchmark/c/benchmark.c
@@ -17,6 +17,7 @@
 */
 
 #include "stdlib/stats/base/dists/normal/stdev.h"
+#include "stdlib/constants/float64/eps.h"
 #include <stdlib.h>
 #include <stdio.h>
 #include <math.h>
@@ -101,7 +102,7 @@ static double benchmark( void ) {
 
 	for ( i = 0; i < 100; i++ ) {
 		mu[ i ] = random_uniform( -10.0, 10.0 );
-		sigma[ i ] = random_uniform( 0.1, 10.0 );
+		sigma[ i ] = random_uniform( STDLIB_CONSTANT_FLOAT64_EPS, 10.0 );
 	}
 
 	t = tic();
diff --git a/lib/node_modules/@stdlib/stats/base/dists/pareto-type1/cdf/benchmark/c/benchmark.c b/lib/node_modules/@stdlib/stats/base/dists/pareto-type1/cdf/benchmark/c/benchmark.c
index b92a9f6e0a0d..8f849e7b0c56 100644
--- a/lib/node_modules/@stdlib/stats/base/dists/pareto-type1/cdf/benchmark/c/benchmark.c
+++ b/lib/node_modules/@stdlib/stats/base/dists/pareto-type1/cdf/benchmark/c/benchmark.c
@@ -17,6 +17,7 @@
 */
 
 #include "stdlib/stats/base/dists/pareto-type1/cdf.h"
+#include "stdlib/constants/float64/eps.h"
 #include <stdlib.h>
 #include <stdio.h>
 #include <math.h>
@@ -102,8 +103,8 @@ static double benchmark( void ) {
 
 	for ( i = 0; i < 100; i++ ) {
 		x[ i ] = random_uniform( 0.0, 20.0 );
-		alpha[ i ] = random_uniform( 0.1, 100.0 );
-		beta[ i ] = random_uniform( 0.1, 100.0 );
+		alpha[ i ] = random_uniform( STDLIB_CONSTANT_FLOAT64_EPS, 100.0 );
+		beta[ i ] = random_uniform( STDLIB_CONSTANT_FLOAT64_EPS, 100.0 );
 	}
 
 	t = tic();
diff --git a/lib/node_modules/@stdlib/stats/base/dists/pareto-type1/entropy/benchmark/c/benchmark.c b/lib/node_modules/@stdlib/stats/base/dists/pareto-type1/entropy/benchmark/c/benchmark.c
index 6129693b5318..5627e9164bd8 100644
--- a/lib/node_modules/@stdlib/stats/base/dists/pareto-type1/entropy/benchmark/c/benchmark.c
+++ b/lib/node_modules/@stdlib/stats/base/dists/pareto-type1/entropy/benchmark/c/benchmark.c
@@ -17,6 +17,7 @@
 */
 
 #include "stdlib/stats/base/dists/pareto-type1/entropy.h"
+#include "stdlib/constants/float64/eps.h"
 #include <stdlib.h>
 #include <stdio.h>
 #include <math.h>
@@ -100,8 +101,8 @@ static double benchmark( void ) {
 	int i;
 
 	for ( i = 0; i < 100; i++ ) {
-		alpha[ i ] = random_uniform( 0.1, 10.0 );
-		beta[ i ] = random_uniform( 0.1, 10.0 );
+		alpha[ i ] = random_uniform( STDLIB_CONSTANT_FLOAT64_EPS, 10.0 );
+		beta[ i ] = random_uniform( STDLIB_CONSTANT_FLOAT64_EPS, 10.0 );
 	}
 
 	t = tic();
diff --git a/lib/node_modules/@stdlib/stats/base/dists/pareto-type1/logcdf/benchmark/c/benchmark.c b/lib/node_modules/@stdlib/stats/base/dists/pareto-type1/logcdf/benchmark/c/benchmark.c
index 5eababa80bd9..d0e90c7bc336 100644
--- a/lib/node_modules/@stdlib/stats/base/dists/pareto-type1/logcdf/benchmark/c/benchmark.c
+++ b/lib/node_modules/@stdlib/stats/base/dists/pareto-type1/logcdf/benchmark/c/benchmark.c
@@ -17,6 +17,7 @@
 */
 
 #include "stdlib/stats/base/dists/pareto-type1/logcdf.h"
+#include "stdlib/constants/float64/eps.h"
 #include <stdlib.h>
 #include <stdio.h>
 #include <math.h>
@@ -102,8 +103,8 @@ static double benchmark( void ) {
 
 	for ( i = 0; i < 100; i++ ) {
 		x[ i ] = random_uniform( 0.0, 8.0 );
-		alpha[ i ] = random_uniform( 0.1, 5.0 );
-		beta[ i ] = random_uniform( 0.1, 5.0 );
+		alpha[ i ] = random_uniform( STDLIB_CONSTANT_FLOAT64_EPS, 5.0 );
+		beta[ i ] = random_uniform( STDLIB_CONSTANT_FLOAT64_EPS, 5.0 );
 	}
 
 	t = tic();
diff --git a/lib/node_modules/@stdlib/stats/base/dists/pareto-type1/mode/benchmark/c/benchmark.c b/lib/node_modules/@stdlib/stats/base/dists/pareto-type1/mode/benchmark/c/benchmark.c
index 428715dddc4f..2e6ff42d1c09 100644
--- a/lib/node_modules/@stdlib/stats/base/dists/pareto-type1/mode/benchmark/c/benchmark.c
+++ b/lib/node_modules/@stdlib/stats/base/dists/pareto-type1/mode/benchmark/c/benchmark.c
@@ -17,6 +17,7 @@
 */
 
 #include "stdlib/stats/base/dists/pareto-type1/mode.h"
+#include "stdlib/constants/float64/eps.h"
 #include <stdlib.h>
 #include <stdio.h>
 #include <math.h>
@@ -100,8 +101,8 @@ static double benchmark( void ) {
 	int i;
 
 	for ( i = 0; i < 100; i++ ) {
-		alpha[ i ] = random_uniform( 0.1, 10.0 );
-		beta[ i ] = random_uniform( 0.1, 10.0 );
+		alpha[ i ] = random_uniform( STDLIB_CONSTANT_FLOAT64_EPS, 10.0 );
+		beta[ i ] = random_uniform( STDLIB_CONSTANT_FLOAT64_EPS, 10.0 );
 	}
 
 	t = tic();
diff --git a/lib/node_modules/@stdlib/stats/base/dists/pareto-type1/pdf/benchmark/c/benchmark.c b/lib/node_modules/@stdlib/stats/base/dists/pareto-type1/pdf/benchmark/c/benchmark.c
index 430b385e36d5..63330753d1c4 100644
--- a/lib/node_modules/@stdlib/stats/base/dists/pareto-type1/pdf/benchmark/c/benchmark.c
+++ b/lib/node_modules/@stdlib/stats/base/dists/pareto-type1/pdf/benchmark/c/benchmark.c
@@ -17,6 +17,7 @@
 */
 
 #include "stdlib/stats/base/dists/pareto-type1/pdf.h"
+#include "stdlib/constants/float64/eps.h"
 #include <stdlib.h>
 #include <stdio.h>
 #include <math.h>
@@ -102,7 +103,7 @@ static double benchmark( void ) {
 
 	for ( i = 0; i < 100; i++ ) {
 		x[ i ] = random_uniform( 0.0, 10.0 );
-		alpha[ i ] = random_uniform( 0.1, 10.0 );
+		alpha[ i ] = random_uniform( STDLIB_CONSTANT_FLOAT64_EPS, 10.0 );
 		beta[ i ] = random_uniform( 2.1, 10.0 );
 	}
 
diff --git a/lib/node_modules/@stdlib/stats/base/dists/pareto-type1/variance/benchmark/c/benchmark.c b/lib/node_modules/@stdlib/stats/base/dists/pareto-type1/variance/benchmark/c/benchmark.c
index e7f03f762ec1..79b6d8a63c30 100644
--- a/lib/node_modules/@stdlib/stats/base/dists/pareto-type1/variance/benchmark/c/benchmark.c
+++ b/lib/node_modules/@stdlib/stats/base/dists/pareto-type1/variance/benchmark/c/benchmark.c
@@ -17,6 +17,7 @@
 */
 
 #include "stdlib/stats/base/dists/pareto-type1/variance.h"
+#include "stdlib/constants/float64/eps.h"
 #include <stdlib.h>
 #include <stdio.h>
 #include <math.h>
@@ -100,8 +101,8 @@ static double benchmark( void ) {
 	int i;
 
 	for ( i = 0; i < 100; i++ ) {
-		alpha[ i ] = random_uniform( 0.1, 10.0 );
-		beta[ i ] = random_uniform( 0.1, 10.0 );
+		alpha[ i ] = random_uniform( STDLIB_CONSTANT_FLOAT64_EPS, 10.0 );
+		beta[ i ] = random_uniform( STDLIB_CONSTANT_FLOAT64_EPS, 10.0 );
 	}
 
 	t = tic();
diff --git a/lib/node_modules/@stdlib/stats/base/dists/planck/cdf/benchmark/c/benchmark.c b/lib/node_modules/@stdlib/stats/base/dists/planck/cdf/benchmark/c/benchmark.c
index d5dface76af1..3aff7889ad05 100644
--- a/lib/node_modules/@stdlib/stats/base/dists/planck/cdf/benchmark/c/benchmark.c
+++ b/lib/node_modules/@stdlib/stats/base/dists/planck/cdf/benchmark/c/benchmark.c
@@ -17,6 +17,7 @@
 */
 
 #include "stdlib/stats/base/dists/planck/cdf.h"
+#include "stdlib/constants/float64/eps.h"
 #include <stdlib.h>
 #include <stdio.h>
 #include <math.h>
@@ -101,7 +102,7 @@ static double benchmark( void ) {
 
 	for ( i = 0; i < 100; i++ ) {
 		x[ i ] = random_uniform( 0.0, 40.0 );
-		lambda[ i ] = random_uniform( 0.1, 10.0 );
+		lambda[ i ] = random_uniform( STDLIB_CONSTANT_FLOAT64_EPS, 10.0 );
 	}
 
 	t = tic();
diff --git a/lib/node_modules/@stdlib/stats/base/dists/planck/entropy/benchmark/c/benchmark.c b/lib/node_modules/@stdlib/stats/base/dists/planck/entropy/benchmark/c/benchmark.c
index 9b27b5a12d02..58a2504f51dd 100644
--- a/lib/node_modules/@stdlib/stats/base/dists/planck/entropy/benchmark/c/benchmark.c
+++ b/lib/node_modules/@stdlib/stats/base/dists/planck/entropy/benchmark/c/benchmark.c
@@ -17,6 +17,7 @@
 */
 
 #include "stdlib/stats/base/dists/planck/entropy.h"
+#include "stdlib/constants/float64/eps.h"
 #include <stdlib.h>
 #include <stdio.h>
 #include <math.h>
@@ -99,7 +100,7 @@ static double benchmark( void ) {
 	int i;
 
 	for ( i = 0; i < 100; i++ ) {
-		lambda[ i ] = random_uniform( 0.1, 10.0 );
+		lambda[ i ] = random_uniform( STDLIB_CONSTANT_FLOAT64_EPS, 10.0 );
 	}
 
 	t = tic();
@@ -135,5 +136,5 @@ int main( void ) {
 		printf( "ok %d benchmark finished\n", i+1 );
 	}
 	print_summary( REPEATS, REPEATS );
+	return 0;
 }
-
diff --git a/lib/node_modules/@stdlib/stats/base/dists/planck/kurtosis/benchmark/c/benchmark.c b/lib/node_modules/@stdlib/stats/base/dists/planck/kurtosis/benchmark/c/benchmark.c
index 4bd2fcc83fc4..26117e743c98 100644
--- a/lib/node_modules/@stdlib/stats/base/dists/planck/kurtosis/benchmark/c/benchmark.c
+++ b/lib/node_modules/@stdlib/stats/base/dists/planck/kurtosis/benchmark/c/benchmark.c
@@ -17,6 +17,7 @@
 */
 
 #include "stdlib/stats/base/dists/planck/kurtosis.h"
+#include "stdlib/constants/float64/eps.h"
 #include <stdlib.h>
 #include <stdio.h>
 #include <math.h>
@@ -99,7 +100,7 @@ static double benchmark( void ) {
 	int i;
 
 	for ( i = 0; i < 100; i++ ) {
-		lambda[ i ] = random_uniform( 0.1, 10.0 );
+		lambda[ i ] = random_uniform( STDLIB_CONSTANT_FLOAT64_EPS, 10.0 );
 	}
 
 	t = tic();
diff --git a/lib/node_modules/@stdlib/stats/base/dists/planck/logcdf/benchmark/c/benchmark.c b/lib/node_modules/@stdlib/stats/base/dists/planck/logcdf/benchmark/c/benchmark.c
index 36ae6bfe2bab..a3c7ccfce771 100644
--- a/lib/node_modules/@stdlib/stats/base/dists/planck/logcdf/benchmark/c/benchmark.c
+++ b/lib/node_modules/@stdlib/stats/base/dists/planck/logcdf/benchmark/c/benchmark.c
@@ -17,6 +17,7 @@
 */
 
 #include "stdlib/stats/base/dists/planck/logcdf.h"
+#include "stdlib/constants/float64/eps.h"
 #include <stdlib.h>
 #include <stdio.h>
 #include <math.h>
@@ -101,7 +102,7 @@ static double benchmark( void ) {
 
 	for ( i = 0; i < 100; i++ ) {
 		x[ i ] = random_uniform( 0.0, 40.0 );
-		lambda[ i ] = random_uniform( 0.1, 10.0 );
+		lambda[ i ] = random_uniform( STDLIB_CONSTANT_FLOAT64_EPS, 10.0 );
 	}
 
 	t = tic();
diff --git a/lib/node_modules/@stdlib/stats/base/dists/planck/logpmf/benchmark/c/benchmark.c b/lib/node_modules/@stdlib/stats/base/dists/planck/logpmf/benchmark/c/benchmark.c
index 70cd1d00dfd0..91809ed80b38 100644
--- a/lib/node_modules/@stdlib/stats/base/dists/planck/logpmf/benchmark/c/benchmark.c
+++ b/lib/node_modules/@stdlib/stats/base/dists/planck/logpmf/benchmark/c/benchmark.c
@@ -17,6 +17,7 @@
 */
 
 #include "stdlib/stats/base/dists/planck/logpmf.h"
+#include "stdlib/constants/float64/eps.h"
 #include "stdlib/math/base/special/round.h"
 #include <stdlib.h>
 #include <stdio.h>
@@ -102,7 +103,7 @@ static double benchmark( void ) {
 
 	for ( i = 0; i < 100; i++ ) {
 		x[ i ] = stdlib_base_round( random_uniform( 0.0, 40.0 ) );
-		lambda[ i ] = random_uniform( 0.1, 10.0 );
+		lambda[ i ] = random_uniform( STDLIB_CONSTANT_FLOAT64_EPS, 10.0 );
 	}
 
 	t = tic();
diff --git a/lib/node_modules/@stdlib/stats/base/dists/planck/mean/benchmark/c/benchmark.c b/lib/node_modules/@stdlib/stats/base/dists/planck/mean/benchmark/c/benchmark.c
index 66d5afed9081..b394ce474a90 100644
--- a/lib/node_modules/@stdlib/stats/base/dists/planck/mean/benchmark/c/benchmark.c
+++ b/lib/node_modules/@stdlib/stats/base/dists/planck/mean/benchmark/c/benchmark.c
@@ -17,6 +17,7 @@
 */
 
 #include "stdlib/stats/base/dists/planck/mean.h"
+#include "stdlib/constants/float64/eps.h"
 #include <stdlib.h>
 #include <stdio.h>
 #include <math.h>
@@ -99,7 +100,7 @@ static double benchmark( void ) {
 	int i;
 
 	for ( i = 0; i < 100; i++ ) {
-		lambda[ i ] = random_uniform( 0.1, 10.0 );
+		lambda[ i ] = random_uniform( STDLIB_CONSTANT_FLOAT64_EPS, 10.0 );
 	}
 
 	t = tic();
diff --git a/lib/node_modules/@stdlib/stats/base/dists/planck/median/benchmark/c/benchmark.c b/lib/node_modules/@stdlib/stats/base/dists/planck/median/benchmark/c/benchmark.c
index 61bf57af2c0d..e8008a69b549 100644
--- a/lib/node_modules/@stdlib/stats/base/dists/planck/median/benchmark/c/benchmark.c
+++ b/lib/node_modules/@stdlib/stats/base/dists/planck/median/benchmark/c/benchmark.c
@@ -17,6 +17,7 @@
 */
 
 #include "stdlib/stats/base/dists/planck/median.h"
+#include "stdlib/constants/float64/eps.h"
 #include <stdlib.h>
 #include <stdio.h>
 #include <math.h>
@@ -99,7 +100,7 @@ static double benchmark( void ) {
 	int i;
 
 	for ( i = 0; i < 100; i++ ) {
-		lambda[ i ] = random_uniform( 0.1, 11.0 );
+		lambda[ i ] = random_uniform( STDLIB_CONSTANT_FLOAT64_EPS, 11.0 );
 	}
 
 	t = tic();
diff --git a/lib/node_modules/@stdlib/stats/base/dists/planck/mgf/benchmark/c/benchmark.c b/lib/node_modules/@stdlib/stats/base/dists/planck/mgf/benchmark/c/benchmark.c
index 50c27e510d53..908997b05d25 100644
--- a/lib/node_modules/@stdlib/stats/base/dists/planck/mgf/benchmark/c/benchmark.c
+++ b/lib/node_modules/@stdlib/stats/base/dists/planck/mgf/benchmark/c/benchmark.c
@@ -17,6 +17,7 @@
 */
 
 #include "stdlib/stats/base/dists/planck/mgf.h"
+#include "stdlib/constants/float64/eps.h"
 #include <stdlib.h>
 #include <stdio.h>
 #include <math.h>
@@ -101,7 +102,7 @@ static double benchmark( void ) {
 
 	for ( i = 0; i < 100; i++ ) {
 		t[ i ] = random_uniform( 0.0, 10.0 );
-		lambda[ i ] = random_uniform( 0.1, 10.0 );
+		lambda[ i ] = random_uniform( STDLIB_CONSTANT_FLOAT64_EPS, 10.0 );
 	}
 
 	start = tic();
diff --git a/lib/node_modules/@stdlib/stats/base/dists/planck/pmf/benchmark/c/benchmark.c b/lib/node_modules/@stdlib/stats/base/dists/planck/pmf/benchmark/c/benchmark.c
index 39d5b907e7e5..fbe47da55be0 100644
--- a/lib/node_modules/@stdlib/stats/base/dists/planck/pmf/benchmark/c/benchmark.c
+++ b/lib/node_modules/@stdlib/stats/base/dists/planck/pmf/benchmark/c/benchmark.c
@@ -17,6 +17,7 @@
 */
 
 #include "stdlib/stats/base/dists/planck/pmf.h"
+#include "stdlib/constants/float64/eps.h"
 #include "stdlib/math/base/special/round.h"
 #include <stdlib.h>
 #include <stdio.h>
@@ -102,7 +103,7 @@ static double benchmark( void ) {
 
 	for ( i = 0; i < 100; i++ ) {
 		x[ i ] = stdlib_base_round( random_uniform( 0.0, 40.0 ) );
-		lambda[ i ] = random_uniform( 0.1, 10.0 );
+		lambda[ i ] = random_uniform( STDLIB_CONSTANT_FLOAT64_EPS, 10.0 );
 	}
 
 	t = tic();
diff --git a/lib/node_modules/@stdlib/stats/base/dists/planck/quantile/benchmark/c/benchmark.c b/lib/node_modules/@stdlib/stats/base/dists/planck/quantile/benchmark/c/benchmark.c
index a9463422009c..326736676116 100644
--- a/lib/node_modules/@stdlib/stats/base/dists/planck/quantile/benchmark/c/benchmark.c
+++ b/lib/node_modules/@stdlib/stats/base/dists/planck/quantile/benchmark/c/benchmark.c
@@ -17,6 +17,7 @@
 */
 
 #include "stdlib/stats/base/dists/planck/quantile.h"
+#include "stdlib/constants/float64/eps.h"
 #include <stdlib.h>
 #include <stdio.h>
 #include <math.h>
@@ -101,7 +102,7 @@ static double benchmark( void ) {
 
 	for ( i = 0; i < 100; i++ ) {
 		p[ i ] = random_uniform( 0.0, 1.0 );
-		lambda[ i ] = random_uniform( 0.1, 10.0 );
+		lambda[ i ] = random_uniform( STDLIB_CONSTANT_FLOAT64_EPS, 10.0 );
 	}
 
 	t = tic();
diff --git a/lib/node_modules/@stdlib/stats/base/dists/planck/skewness/benchmark/c/benchmark.c b/lib/node_modules/@stdlib/stats/base/dists/planck/skewness/benchmark/c/benchmark.c
index 38595c985a86..4ce9e93e90b6 100644
--- a/lib/node_modules/@stdlib/stats/base/dists/planck/skewness/benchmark/c/benchmark.c
+++ b/lib/node_modules/@stdlib/stats/base/dists/planck/skewness/benchmark/c/benchmark.c
@@ -17,6 +17,7 @@
 */
 
 #include "stdlib/stats/base/dists/planck/skewness.h"
+#include "stdlib/constants/float64/eps.h"
 #include <stdlib.h>
 #include <stdio.h>
 #include <math.h>
@@ -99,7 +100,7 @@ static double benchmark( void ) {
 	int i;
 
 	for ( i = 0; i < 100; i++ ) {
-		lambda[ i ] = random_uniform( 0.1, 10.0 );
+		lambda[ i ] = random_uniform( STDLIB_CONSTANT_FLOAT64_EPS, 10.0 );
 	}
 
 	t = tic();
diff --git a/lib/node_modules/@stdlib/stats/base/dists/planck/stdev/benchmark/c/benchmark.c b/lib/node_modules/@stdlib/stats/base/dists/planck/stdev/benchmark/c/benchmark.c
index 6a27adafedfe..1fb54e7ad652 100644
--- a/lib/node_modules/@stdlib/stats/base/dists/planck/stdev/benchmark/c/benchmark.c
+++ b/lib/node_modules/@stdlib/stats/base/dists/planck/stdev/benchmark/c/benchmark.c
@@ -17,6 +17,7 @@
 */
 
 #include "stdlib/stats/base/dists/planck/stdev.h"
+#include "stdlib/constants/float64/eps.h"
 #include <stdlib.h>
 #include <stdio.h>
 #include <math.h>
@@ -99,7 +100,7 @@ static double benchmark( void ) {
 	int i;
 
 	for ( i = 0; i < 100; i++ ) {
-		lambda[i] = random_uniform( 0.1, 10.0 );
+		lambda[i] = random_uniform( STDLIB_CONSTANT_FLOAT64_EPS, 10.0 );
 	}
 
 	t = tic();
diff --git a/lib/node_modules/@stdlib/stats/base/dists/planck/variance/benchmark/c/benchmark.c b/lib/node_modules/@stdlib/stats/base/dists/planck/variance/benchmark/c/benchmark.c
index 34678773ccc3..7bd0ff8badc2 100644
--- a/lib/node_modules/@stdlib/stats/base/dists/planck/variance/benchmark/c/benchmark.c
+++ b/lib/node_modules/@stdlib/stats/base/dists/planck/variance/benchmark/c/benchmark.c
@@ -17,6 +17,7 @@
 */
 
 #include "stdlib/stats/base/dists/planck/variance.h"
+#include "stdlib/constants/float64/eps.h"
 #include <stdlib.h>
 #include <stdio.h>
 #include <math.h>
@@ -99,7 +100,7 @@ static double benchmark( void ) {
 	int i;
 
 	for ( i = 0; i < 100; i++ ) {
-		lambda[ i ] = random_uniform( 0.1, 10.0 );
+		lambda[ i ] = random_uniform( STDLIB_CONSTANT_FLOAT64_EPS, 10.0 );
 	}
 
 	t = tic();
diff --git a/lib/node_modules/@stdlib/stats/base/dists/poisson/cdf/benchmark/c/benchmark.c b/lib/node_modules/@stdlib/stats/base/dists/poisson/cdf/benchmark/c/benchmark.c
index be0c402e3278..ca79e7fdb95c 100644
--- a/lib/node_modules/@stdlib/stats/base/dists/poisson/cdf/benchmark/c/benchmark.c
+++ b/lib/node_modules/@stdlib/stats/base/dists/poisson/cdf/benchmark/c/benchmark.c
@@ -17,6 +17,7 @@
 */
 
 #include "stdlib/stats/base/dists/poisson/cdf.h"
+#include "stdlib/constants/float64/eps.h"
 #include <stdlib.h>
 #include <stdio.h>
 #include <math.h>
@@ -100,7 +101,7 @@ static double benchmark( void ) {
 	int i;
 
 	for ( i = 0; i < 100; i++ ) {
-		lambda[ i ] = random_uniform( 0.1, 20.0 );
+		lambda[ i ] = random_uniform( STDLIB_CONSTANT_FLOAT64_EPS, 20.0 );
 		x[ i ] = random_uniform( 0.0, 10.0 );
 	}
 
diff --git a/lib/node_modules/@stdlib/stats/base/dists/poisson/mgf/benchmark/c/benchmark.c b/lib/node_modules/@stdlib/stats/base/dists/poisson/mgf/benchmark/c/benchmark.c
index fe7864f2bd7f..2edc8ebba4d4 100644
--- a/lib/node_modules/@stdlib/stats/base/dists/poisson/mgf/benchmark/c/benchmark.c
+++ b/lib/node_modules/@stdlib/stats/base/dists/poisson/mgf/benchmark/c/benchmark.c
@@ -17,6 +17,7 @@
 */
 
 #include "stdlib/stats/base/dists/poisson/mgf.h"
+#include "stdlib/constants/float64/eps.h"
 #include <sys/time.h>
 #include <math.h>
 #include <stdio.h>
@@ -101,7 +102,7 @@ static double benchmark( void ) {
 
 	for ( i = 0; i < 100; i++ ) {
 		t[ i ] = random_uniform( -2.0, 2.0 );
-		lambda[ i ] = random_uniform( 0.1, 10.0 );
+		lambda[ i ] = random_uniform( STDLIB_CONSTANT_FLOAT64_EPS, 10.0 );
 	}
 
 	start = tic();
diff --git a/lib/node_modules/@stdlib/stats/base/dists/t/pdf/benchmark/c/benchmark.c b/lib/node_modules/@stdlib/stats/base/dists/t/pdf/benchmark/c/benchmark.c
index 09578a34d005..c95b578c16da 100644
--- a/lib/node_modules/@stdlib/stats/base/dists/t/pdf/benchmark/c/benchmark.c
+++ b/lib/node_modules/@stdlib/stats/base/dists/t/pdf/benchmark/c/benchmark.c
@@ -17,6 +17,7 @@
 */
 
 #include "stdlib/stats/base/dists/t/pdf.h"
+#include "stdlib/constants/float64/eps.h"
 #include <math.h>
 #include <stdio.h>
 #include <stdlib.h>
@@ -101,7 +102,7 @@ static double benchmark( void ) {
 
 	for ( i = 0; i < 100; i++ ) {
 		x[ i ] = random_uniform( -10.0, 10.0 );
-		v[ i ] = random_uniform( 0.1, 10.0 );
+		v[ i ] = random_uniform( STDLIB_CONSTANT_FLOAT64_EPS, 10.0 );
 	}
 
 	t = tic();
diff --git a/lib/node_modules/@stdlib/stats/base/dists/weibull/entropy/benchmark/c/benchmark.c b/lib/node_modules/@stdlib/stats/base/dists/weibull/entropy/benchmark/c/benchmark.c
index 314185c3cc34..3aa2c221cbbf 100644
--- a/lib/node_modules/@stdlib/stats/base/dists/weibull/entropy/benchmark/c/benchmark.c
+++ b/lib/node_modules/@stdlib/stats/base/dists/weibull/entropy/benchmark/c/benchmark.c
@@ -17,6 +17,7 @@
 */
 
 #include "stdlib/stats/base/dists/weibull/entropy.h"
+#include "stdlib/constants/float64/eps.h"
 #include <stdlib.h>
 #include <stdio.h>
 #include <math.h>
@@ -100,8 +101,8 @@ static double benchmark( void ) {
 	int i;
 
 	for ( i = 0; i < 100; i++ ) {
-		k[ i ] = random_uniform( 0.1, 10.0 );
-		lambda[ i ] = random_uniform( 0.1, 10.0 );
+		k[ i ] = random_uniform( STDLIB_CONSTANT_FLOAT64_EPS, 10.0 );
+		lambda[ i ] = random_uniform( STDLIB_CONSTANT_FLOAT64_EPS, 10.0 );
 	}
 
 	t = tic();
diff --git a/lib/node_modules/@stdlib/stats/base/dists/weibull/quantile/benchmark/c/benchmark.c b/lib/node_modules/@stdlib/stats/base/dists/weibull/quantile/benchmark/c/benchmark.c
index 036291d5edc5..593f6264cc19 100644
--- a/lib/node_modules/@stdlib/stats/base/dists/weibull/quantile/benchmark/c/benchmark.c
+++ b/lib/node_modules/@stdlib/stats/base/dists/weibull/quantile/benchmark/c/benchmark.c
@@ -17,6 +17,7 @@
 */
 
 #include "stdlib/stats/base/dists/weibull/quantile.h"
+#include "stdlib/constants/float64/eps.h"
 #include <sys/time.h>
 #include <math.h>
 #include <stdio.h>
@@ -102,8 +103,8 @@ static double benchmark( void ) {
 
 	for ( i = 0; i < 100; i++ ) {
 		p[ i ] = random_uniform( 0.0, 1.0 );
-		k[ i ] = random_uniform( 0.1, 5.0 );
-		lambda[ i ] = random_uniform( 0.1, 5.0 );
+		k[ i ] = random_uniform( STDLIB_CONSTANT_FLOAT64_EPS, 5.0 );
+		lambda[ i ] = random_uniform( STDLIB_CONSTANT_FLOAT64_EPS, 5.0 );
 	}
 
 	t = tic();
diff --git a/lib/node_modules/@stdlib/stats/base/dists/weibull/stdev/benchmark/c/benchmark.c b/lib/node_modules/@stdlib/stats/base/dists/weibull/stdev/benchmark/c/benchmark.c
index 441db7fc3081..630205e2e501 100644
--- a/lib/node_modules/@stdlib/stats/base/dists/weibull/stdev/benchmark/c/benchmark.c
+++ b/lib/node_modules/@stdlib/stats/base/dists/weibull/stdev/benchmark/c/benchmark.c
@@ -17,6 +17,7 @@
 */
 
 #include "stdlib/stats/base/dists/weibull/stdev.h"
+#include "stdlib/constants/float64/eps.h"
 #include <stdlib.h>
 #include <stdio.h>
 #include <math.h>
@@ -100,8 +101,8 @@ static double benchmark( void ) {
 	int i;
 
 	for ( i = 0; i < 100; i++ ) {
-		k[ i ] = random_uniform( 0.1, 10.0 );
-		lambda[ i ] = random_uniform( 0.1, 10.0 );
+		k[ i ] = random_uniform( STDLIB_CONSTANT_FLOAT64_EPS, 10.0 );
+		lambda[ i ] = random_uniform( STDLIB_CONSTANT_FLOAT64_EPS, 10.0 );
 	}
 
 	t = tic();