chore: directly draw from the desired distribution instead of adding constants

Planeshifter · Planeshifter · commit fc0ff171dab5 · 2025-01-12T17:22:40.000-05:00
---
type: pre_commit_static_analysis_report
description: Results of running static analysis checks when committing changes.
report:
  - task: lint_filenames
    status: passed
  - task: lint_editorconfig
    status: passed
  - task: lint_markdown
    status: passed
  - task: lint_package_json
    status: na
  - task: lint_repl_help
    status: na
  - task: lint_javascript_src
    status: na
  - task: lint_javascript_cli
    status: na
  - task: lint_javascript_examples
    status: na
  - task: lint_javascript_tests
    status: na
  - task: lint_javascript_benchmarks
    status: na
  - task: lint_python
    status: na
  - task: lint_r
    status: na
  - task: lint_c_src
    status: na
  - task: lint_c_examples
    status: passed
  - task: lint_c_benchmarks
    status: passed
  - task: lint_c_tests_fixtures
    status: na
  - task: lint_shell
    status: na
  - task: lint_typescript_declarations
    status: na
  - task: lint_typescript_tests
    status: na
  - task: lint_license_headers
    status: passed
---
diff --git a/lib/node_modules/@stdlib/stats/base/dists/beta/skewness/README.md b/lib/node_modules/@stdlib/stats/base/dists/beta/skewness/README.md
@@ -224,8 +224,8 @@ int main( void ) {
     int i;
 
     for ( i = 0; i < 25; i++ ) {
-        alpha = random_uniform( 0.0, 10.0 ) + STDLIB_CONSTANT_FLOAT64_EPS;
-        beta = random_uniform( 0.0, 10.0 ) + STDLIB_CONSTANT_FLOAT64_EPS;
+        alpha = random_uniform( STDLIB_CONSTANT_FLOAT64_EPS, 10.0 );
+        beta = random_uniform( STDLIB_CONSTANT_FLOAT64_EPS, 10.0 );
         y = stdlib_base_dists_beta_skewness( alpha, beta );
         printf( "α: %lf, β: %lf, skew(X;α,β): %lf\n", alpha, beta, y );
     }
diff --git a/lib/node_modules/@stdlib/stats/base/dists/beta/skewness/benchmark/c/benchmark.c b/lib/node_modules/@stdlib/stats/base/dists/beta/skewness/benchmark/c/benchmark.c
@@ -101,8 +101,8 @@ static double benchmark( void ) {
 	int i;
 
 	for ( i = 0; i < 100; i++ ) {
-		alpha[ i ] = random_uniform( 0.0, 10.0 ) + STDLIB_CONSTANT_FLOAT64_EPS;
-		beta[ i ] = random_uniform( 0.0, 10.0 ) + STDLIB_CONSTANT_FLOAT64_EPS;
+		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/beta/skewness/examples/c/example.c b/lib/node_modules/@stdlib/stats/base/dists/beta/skewness/examples/c/example.c
@@ -33,8 +33,8 @@ int main( void ) {
 	int i;
 
 	for ( i = 0; i < 25; i++ ) {
-		alpha = random_uniform( 0.0, 10.0 ) + STDLIB_CONSTANT_FLOAT64_EPS;
-		beta = random_uniform( 0.0, 10.0 ) + STDLIB_CONSTANT_FLOAT64_EPS;
+		alpha = random_uniform( STDLIB_CONSTANT_FLOAT64_EPS, 10.0 );
+		beta = random_uniform( STDLIB_CONSTANT_FLOAT64_EPS, 10.0 );
 		y = stdlib_base_dists_beta_skewness( alpha, beta );
 		printf( "α: %lf, β: %lf, skew(X;α,β): %lf\n", alpha, beta, y );
 	}
diff --git a/lib/node_modules/@stdlib/stats/base/dists/beta/stdev/README.md b/lib/node_modules/@stdlib/stats/base/dists/beta/stdev/README.md
@@ -224,8 +224,8 @@ int main( void ) {
     int i;
 
     for ( i = 0; i < 25; i++ ) {
-        alpha = random_uniform( 0.0, 10.0 ) + STDLIB_CONSTANT_FLOAT64_EPS;
-        beta = random_uniform( 0.0, 10.0 ) + STDLIB_CONSTANT_FLOAT64_EPS;
+        alpha = random_uniform( STDLIB_CONSTANT_FLOAT64_EPS, 10.0 );
+        beta = random_uniform( STDLIB_CONSTANT_FLOAT64_EPS, 10.0 );
         y = stdlib_base_dists_beta_stdev( alpha, beta );
         printf( "alpha: %lf, beta: %lf, F(X;a,b): %lf\n", alpha, beta, y );
     }
diff --git a/lib/node_modules/@stdlib/stats/base/dists/beta/stdev/benchmark/c/benchmark.c b/lib/node_modules/@stdlib/stats/base/dists/beta/stdev/benchmark/c/benchmark.c
@@ -101,8 +101,8 @@ static double benchmark( void ) {
 	int i;
 
 	for ( i = 0; i < 100; i++ ) {
-		alpha[ i ] = random_uniform( 0.0, 10.0 ) + STDLIB_CONSTANT_FLOAT64_EPS;
-		beta[ i ] = random_uniform( 0.0, 10.0 ) + STDLIB_CONSTANT_FLOAT64_EPS;
+		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/beta/stdev/examples/c/example.c b/lib/node_modules/@stdlib/stats/base/dists/beta/stdev/examples/c/example.c
@@ -33,8 +33,8 @@ int main( void ) {
 	int i;
 
 	for ( i = 0; i < 25; i++ ) {
-		beta = random_uniform( 0.0, 10.0 ) + STDLIB_CONSTANT_FLOAT64_EPS;
-		alpha = random_uniform( 0.0, 10.0 ) + STDLIB_CONSTANT_FLOAT64_EPS;
+		beta = random_uniform( STDLIB_CONSTANT_FLOAT64_EPS, 10.0 );
+		alpha = random_uniform( STDLIB_CONSTANT_FLOAT64_EPS, 10.0 );
 		y = stdlib_base_dists_beta_stdev( alpha, beta );
 		printf( "alpha: %lf, beta: %lf, F(X;alpha,beta): %lf\n", alpha, beta, y );
 	}
diff --git a/lib/node_modules/@stdlib/stats/base/dists/beta/variance/benchmark/c/benchmark.c b/lib/node_modules/@stdlib/stats/base/dists/beta/variance/benchmark/c/benchmark.c
@@ -101,8 +101,8 @@ static double benchmark( void ) {
 	int i;
 
 	for ( i = 0; i < 100; i++ ) {
-		alpha[ i ] = random_uniform( 0.0, 10.0 ) + STDLIB_CONSTANT_FLOAT64_EPS;
-		beta[ i ] = random_uniform( 0.0, 10.0 ) + STDLIB_CONSTANT_FLOAT64_EPS;
+		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/beta/variance/examples/c/example.c b/lib/node_modules/@stdlib/stats/base/dists/beta/variance/examples/c/example.c
@@ -33,8 +33,8 @@ int main( void ) {
 	int i;
 
 	for ( i = 0; i < 25; i++ ) {
-		alpha = random_uniform( 0.0, 10.0 ) + STDLIB_CONSTANT_FLOAT64_EPS;
-		beta = random_uniform( 0.0, 10.0 ) + STDLIB_CONSTANT_FLOAT64_EPS;
+		alpha = random_uniform( STDLIB_CONSTANT_FLOAT64_EPS, 10.0 );
+		beta = random_uniform( STDLIB_CONSTANT_FLOAT64_EPS, 10.0 );
 		y = stdlib_base_dists_beta_variance( alpha, beta );
 		printf( "α: %lf, β: %lf, Var(X;α,β): %lf\n", alpha, beta, y );
 	}
diff --git a/lib/node_modules/@stdlib/stats/base/dists/betaprime/kurtosis/benchmark/c/benchmark.c b/lib/node_modules/@stdlib/stats/base/dists/betaprime/kurtosis/benchmark/c/benchmark.c
@@ -101,7 +101,7 @@ static double benchmark( void ) {
 
 	for ( i = 0; i < 100; i++ ) {
 		alpha[ i ] = random_uniform( 0.0, 20.0 );
-		beta[ i ] = random_uniform( 0.0, 20.0 ) + 4.0;
+		beta[ i ] = random_uniform( 4.0, 24.0 );
 	}
 
 	t = tic();
diff --git a/lib/node_modules/@stdlib/stats/base/dists/cauchy/entropy/examples/c/example.c b/lib/node_modules/@stdlib/stats/base/dists/cauchy/entropy/examples/c/example.c
@@ -34,7 +34,7 @@ int main( void ) {
 
 	for ( i = 0; i < 25; i++ ) {
 		x0 = random_uniform( 0.0, 100.0 );
-		gamma = random_uniform( 0.0, 10.0 ) + STDLIB_CONSTANT_FLOAT64_EPS;
+		gamma = random_uniform( STDLIB_CONSTANT_FLOAT64_EPS, 10.0 );
 		y = stdlib_base_dists_cauchy_entropy( x0, gamma );
 		printf( "x0: %lf, gamma: %lf, F(x0;gamma): %lf\n", x0, gamma, y );
 	}
diff --git a/lib/node_modules/@stdlib/stats/base/dists/cauchy/mode/examples/c/example.c b/lib/node_modules/@stdlib/stats/base/dists/cauchy/mode/examples/c/example.c
@@ -34,7 +34,7 @@ int main( void ) {
 
 	for ( i = 0; i < 25; i++ ) {
 		x0 = random_uniform( 0.0, 100.0 );
-		gamma = random_uniform( 0.0, 10.0 ) + STDLIB_CONSTANT_FLOAT64_EPS;
+		gamma = random_uniform( STDLIB_CONSTANT_FLOAT64_EPS, 10.0 );
 		y = stdlib_base_dists_cauchy_mode( x0, gamma );
 		printf( "x0: %lf, gamma: %lf, mode(X;x0,gamma): %lf\n", x0, gamma, y );
 	}
diff --git a/lib/node_modules/@stdlib/stats/base/dists/cosine/mode/examples/c/example.c b/lib/node_modules/@stdlib/stats/base/dists/cosine/mode/examples/c/example.c
@@ -34,7 +34,7 @@ int main( void ) {
 
 	for ( i = 0; i < 10; i++ ) {
 		mu = random_uniform( -50.0, 50.0 );
-		s = random_uniform( 0.0, 20.0 ) + STDLIB_CONSTANT_FLOAT64_EPS;
+		s = random_uniform( STDLIB_CONSTANT_FLOAT64_EPS, 20.0 );
 		y = stdlib_base_dists_cosine_mode( mu, s );
 		printf( "µ: %lf, s: %lf, mode(X;µ,s): %lf\n", mu, s , y );
 	}
diff --git a/lib/node_modules/@stdlib/stats/base/dists/exponential/stdev/benchmark/c/benchmark.c b/lib/node_modules/@stdlib/stats/base/dists/exponential/stdev/benchmark/c/benchmark.c
@@ -100,7 +100,7 @@ static double benchmark( void ) {
 	int i;
 
 	for ( i = 0; i < 100; i++ ) {
-		lambda[ i ] = random_uniform( 0.0, 20.0 ) + STDLIB_CONSTANT_FLOAT64_EPS;
+		lambda[ i ] = random_uniform( STDLIB_CONSTANT_FLOAT64_EPS, 20.0 );
 	}
 
 	t = tic();
diff --git a/lib/node_modules/@stdlib/stats/base/dists/f/entropy/README.md b/lib/node_modules/@stdlib/stats/base/dists/f/entropy/README.md
@@ -224,8 +224,8 @@ int main( void ) {
     int i;
 
     for ( i = 0; i < 25; i++ ) {
-        d1 = random_uniform( 0.0, 10.0 ) + STDLIB_CONSTANT_FLOAT64_EPS;
-        d2 = random_uniform( 0.0, 10.0 ) + STDLIB_CONSTANT_FLOAT64_EPS;
+        d1 = random_uniform( STDLIB_CONSTANT_FLOAT64_EPS, 10.0 );
+        d2 = random_uniform( STDLIB_CONSTANT_FLOAT64_EPS, 10.0 );
         y = stdlib_base_dists_f_entropy( d1, d2 );
         printf( "d1: %lf, d2: %lf, h(X;d1,d2): %lf\n", d1, d2, y );
     }
diff --git a/lib/node_modules/@stdlib/stats/base/dists/f/entropy/examples/c/example.c b/lib/node_modules/@stdlib/stats/base/dists/f/entropy/examples/c/example.c
@@ -33,8 +33,8 @@ int main( void ) {
 	int i;
 
 	for ( i = 0; i < 25; i++ ) {
-		d1 = random_uniform( 0.0, 10.0 ) + STDLIB_CONSTANT_FLOAT64_EPS;
-		d2 = random_uniform( 0.0, 10.0 ) + STDLIB_CONSTANT_FLOAT64_EPS;
+		d1 = random_uniform( STDLIB_CONSTANT_FLOAT64_EPS, 10.0 );
+		d2 = random_uniform( STDLIB_CONSTANT_FLOAT64_EPS, 10.0 );
 		y = stdlib_base_dists_f_entropy( d1, d2 );
 		printf( "d1: %lf, d2: %lf, h(X;d1,d2): %lf\n", d1, d2, y );
 	}
diff --git a/lib/node_modules/@stdlib/stats/base/dists/f/skewness/README.md b/lib/node_modules/@stdlib/stats/base/dists/f/skewness/README.md
@@ -227,8 +227,8 @@ int main( void ) {
     int i;
 
     for ( i = 0; i < 25; i++ ) {
-        d1 = random_uniform( 0.0, 10.0 ) + STDLIB_CONSTANT_FLOAT64_EPS;
-        d2 = random_uniform( 0.0, 20.0 ) + STDLIB_CONSTANT_FLOAT64_EPS;
+        d1 = random_uniform( STDLIB_CONSTANT_FLOAT64_EPS, 10.0 );
+        d2 = random_uniform( STDLIB_CONSTANT_FLOAT64_EPS, 20.0 );
         y = stdlib_base_dists_f_skewness( d1, d2 );
         printf( "d1: %lf, d2: %lf, skew(X;d1,d2): %lf\n", d1, d2, y );
     }
diff --git a/lib/node_modules/@stdlib/stats/base/dists/f/skewness/benchmark/c/benchmark.c b/lib/node_modules/@stdlib/stats/base/dists/f/skewness/benchmark/c/benchmark.c
@@ -101,8 +101,8 @@ static double benchmark( void ) {
 	int i;
 
 	for ( i = 0; i < 100; i++ ) {
-		d1[ i ] = random_uniform( 0.0, 10.0 ) + STDLIB_CONSTANT_FLOAT64_EPS;
-		d2[ i ] = random_uniform( 0.0, 10.0 ) + 6.0 + STDLIB_CONSTANT_FLOAT64_EPS;
+		d1[ i ] = random_uniform( STDLIB_CONSTANT_FLOAT64_EPS, 10.0 );
+		d2[ i ] = random_uniform( 6.0, 16.0 );
 	}
 
 	t = tic();
diff --git a/lib/node_modules/@stdlib/stats/base/dists/f/skewness/examples/c/example.c b/lib/node_modules/@stdlib/stats/base/dists/f/skewness/examples/c/example.c
@@ -33,8 +33,8 @@ int main( void ) {
 	int i;
 
 	for ( i = 0; i < 25; i++ ) {
-		d1 = random_uniform( 0.0, 10.0 ) + STDLIB_CONSTANT_FLOAT64_EPS;
-		d2 = random_uniform( 0.0, 20.0 ) + STDLIB_CONSTANT_FLOAT64_EPS;
+		d1 = random_uniform( STDLIB_CONSTANT_FLOAT64_EPS, 10.0 );
+		d2 = random_uniform( STDLIB_CONSTANT_FLOAT64_EPS, 20.0 );
 		y = stdlib_base_dists_f_skewness( d1, d2 );
 		printf( "d1: %lf, d2: %lf, skew(X;d1,d2): %lf\n", d1, d2, y );
 	}
diff --git a/lib/node_modules/@stdlib/stats/base/dists/f/stdev/README.md b/lib/node_modules/@stdlib/stats/base/dists/f/stdev/README.md
@@ -224,8 +224,8 @@ int main( void ) {
     int i;
 
     for ( i = 0; i < 25; i++ ) {
-        d1 = random_uniform( 0.0, 10.0 ) + STDLIB_CONSTANT_FLOAT64_EPS;
-        d2 = random_uniform( 0.0, 20.0 ) + STDLIB_CONSTANT_FLOAT64_EPS;
+        d1 = random_uniform( STDLIB_CONSTANT_FLOAT64_EPS, 10.0 );
+        d2 = random_uniform( STDLIB_CONSTANT_FLOAT64_EPS, 20.0 );
         y = stdlib_base_dists_f_stdev( d1, d2 );
         printf( "d1: %lf, d2: %lf, SD(X;d1,d2): %lf\n", d1, d2, y );
     }
diff --git a/lib/node_modules/@stdlib/stats/base/dists/f/stdev/benchmark/c/benchmark.c b/lib/node_modules/@stdlib/stats/base/dists/f/stdev/benchmark/c/benchmark.c
@@ -101,8 +101,8 @@ static double benchmark( void ) {
 	int i;
 
 	for ( i = 0; i < 100; i++ ) {
-		d1[ i ] = random_uniform( 0.0, 10.0 ) + STDLIB_CONSTANT_FLOAT64_EPS;
-		d2[ i ] = random_uniform( 0.0, 10.0 ) + 4.0 + STDLIB_CONSTANT_FLOAT64_EPS;
+		d1[ i ] = random_uniform( STDLIB_CONSTANT_FLOAT64_EPS, 10.0 );
+		d2[ i ] = random_uniform( 4.0, 14.0 );
 	}
 
 	t = tic();
diff --git a/lib/node_modules/@stdlib/stats/base/dists/f/variance/README.md b/lib/node_modules/@stdlib/stats/base/dists/f/variance/README.md
@@ -224,8 +224,8 @@ int main( void ) {
     int i;
 
     for ( i = 0; i < 25; i++ ) {
-        d1 = random_uniform( 0.0, 10.0 ) + STDLIB_CONSTANT_FLOAT64_EPS;
-        d2 = random_uniform( 0.0, 20.0 ) + STDLIB_CONSTANT_FLOAT64_EPS;
+        d1 = random_uniform( STDLIB_CONSTANT_FLOAT64_EPS, 10.0 );
+        d2 = random_uniform( STDLIB_CONSTANT_FLOAT64_EPS, 20.0 );
         y = stdlib_base_dists_f_variance( d1, d2 );
         printf( "d1: %lf, d2: %lf, Var(X;d1,d2): %lf\n", d1, d2, y );
     }
diff --git a/lib/node_modules/@stdlib/stats/base/dists/f/variance/benchmark/c/benchmark.c b/lib/node_modules/@stdlib/stats/base/dists/f/variance/benchmark/c/benchmark.c
@@ -101,8 +101,8 @@ static double benchmark( void ) {
 	int i;
 
 	for ( i = 0; i < 100; i++ ) {
-		d1[ i ] = random_uniform( 0.0, 10.0 ) + STDLIB_CONSTANT_FLOAT64_EPS;
-		d2[ i ] = random_uniform( 0.0, 10.0 ) + 4.0 + STDLIB_CONSTANT_FLOAT64_EPS;
+		d1[ i ] = random_uniform( STDLIB_CONSTANT_FLOAT64_EPS, 10.0 );
+		d2[ i ] = random_uniform( 4.0, 14.0 );
 	}
 
 	t = tic();
diff --git a/lib/node_modules/@stdlib/stats/base/dists/gamma/mean/benchmark/c/benchmark.c b/lib/node_modules/@stdlib/stats/base/dists/gamma/mean/benchmark/c/benchmark.c
@@ -101,8 +101,8 @@ static double benchmark( void ) {
 	int i;
 
 	for ( i = 0; i < 100; i++ ) {
-		alpha[ i ] = random_uniform( 0.0, 20.0 ) + STDLIB_CONSTANT_FLOAT64_EPS;
-		beta[ i ] = random_uniform( 0.0, 20.0 ) + STDLIB_CONSTANT_FLOAT64_EPS;
+		alpha[ i ] = random_uniform( STDLIB_CONSTANT_FLOAT64_EPS, 20.0 );
+		beta[ i ] = random_uniform( STDLIB_CONSTANT_FLOAT64_EPS, 20.0 );
 	}
 
 	t = tic();
diff --git a/lib/node_modules/@stdlib/stats/base/dists/gamma/mode/benchmark/c/benchmark.c b/lib/node_modules/@stdlib/stats/base/dists/gamma/mode/benchmark/c/benchmark.c
@@ -101,8 +101,8 @@ static double benchmark( void ) {
 	int i;
 
 	for ( i = 0; i < 100; i++ ) {
-		alpha[ i ] = random_uniform( 0.0, 20.0 ) + STDLIB_CONSTANT_FLOAT64_EPS;
-		beta[ i ] = random_uniform( 0.0, 20.0 ) + STDLIB_CONSTANT_FLOAT64_EPS;
+		alpha[ i ] = random_uniform( STDLIB_CONSTANT_FLOAT64_EPS, 20.0 );
+		beta[ i ] = random_uniform( STDLIB_CONSTANT_FLOAT64_EPS, 20.0 );
 	}
 
 	t = tic();
diff --git a/lib/node_modules/@stdlib/stats/base/dists/gamma/skewness/benchmark/c/benchmark.c b/lib/node_modules/@stdlib/stats/base/dists/gamma/skewness/benchmark/c/benchmark.c
@@ -101,8 +101,8 @@ static double benchmark( void ) {
 	int i;
 
 	for ( i = 0; i < 100; i++ ) {
-		alpha[ i ] = random_uniform( 0.0, 20.0 ) + STDLIB_CONSTANT_FLOAT64_EPS;
-		beta[ i ] = random_uniform( 0.0, 20.0 ) + STDLIB_CONSTANT_FLOAT64_EPS;
+		alpha[ i ] = random_uniform( STDLIB_CONSTANT_FLOAT64_EPS, 20.0 );
+		beta[ i ] = random_uniform( STDLIB_CONSTANT_FLOAT64_EPS, 20.0 );
 	}
 
 	t = tic();
diff --git a/lib/node_modules/@stdlib/stats/base/dists/geometric/mean/benchmark/c/benchmark.c b/lib/node_modules/@stdlib/stats/base/dists/geometric/mean/benchmark/c/benchmark.c
@@ -100,7 +100,7 @@ static double benchmark( void ) {
 	int i;
 
 	for ( i = 0; i < 100; i++ ) {
-		p[ i ] = random_uniform( 0.0, 1.0 ) + STDLIB_CONSTANT_FLOAT64_EPS;
+		p[ i ] = random_uniform( STDLIB_CONSTANT_FLOAT64_EPS, 1.0 );
 	}
 
 	t = tic();
diff --git a/lib/node_modules/@stdlib/stats/base/dists/gumbel/variance/examples/c/example.c b/lib/node_modules/@stdlib/stats/base/dists/gumbel/variance/examples/c/example.c
@@ -32,7 +32,7 @@ int main( void ) {
 	int i;
 
 	for ( i = 0; i < 25; i++ ) {
-		mu = random_uniform( 0.0, 10.0 ) - 5.0;
+		mu = random_uniform( -5.0, 5.0 );
 		beta = random_uniform( 0.0, 20.0 );
 		y = stdlib_base_dists_gumbel_variance( mu, beta );
 		printf( "µ: %lf, β: %lf, Var(X;µ,β): %lf\n", mu, beta, y );
diff --git a/lib/node_modules/@stdlib/stats/base/dists/invgamma/mode/README.md b/lib/node_modules/@stdlib/stats/base/dists/invgamma/mode/README.md
@@ -224,8 +224,8 @@ int main( void ) {
     int i;
 
     for ( i = 0; i < 25; i++ ) {
-        a = random_uniform( 0.0, 10.0 ) + STDLIB_CONSTANT_FLOAT64_EPS;
-        b = random_uniform( 0.0, 10.0 ) + STDLIB_CONSTANT_FLOAT64_EPS;
+        a = random_uniform( STDLIB_CONSTANT_FLOAT64_EPS, 10.0 );
+        b = random_uniform( STDLIB_CONSTANT_FLOAT64_EPS, 10.0 );
         y = stdlib_base_dists_invgamma_mode( alpha, beta );
         printf( "α: %lf, β: %lf, mode(X;α,β): %lf\n", alpha, beta, y );
     }
diff --git a/lib/node_modules/@stdlib/stats/base/dists/invgamma/mode/benchmark/c/benchmark.c b/lib/node_modules/@stdlib/stats/base/dists/invgamma/mode/benchmark/c/benchmark.c
@@ -101,8 +101,8 @@ static double benchmark( void ) {
 	int i;
 
 	for ( i = 0; i < 100; i++ ) {
-		alpha[ i ] = random_uniform( 0.0, 10.0 ) + STDLIB_CONSTANT_FLOAT64_EPS;
-		beta[ i ] = random_uniform( 0.0, 10.0 ) + STDLIB_CONSTANT_FLOAT64_EPS;
+		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/invgamma/mode/examples/c/example.c b/lib/node_modules/@stdlib/stats/base/dists/invgamma/mode/examples/c/example.c
@@ -33,8 +33,8 @@ int main( void ) {
 	int i;
 
 	for ( i = 0; i < 25; i++ ) {
-		alpha = random_uniform( 0.0, 10.0 ) + STDLIB_CONSTANT_FLOAT64_EPS;
-		beta = random_uniform( 0.0, 10.0 ) + STDLIB_CONSTANT_FLOAT64_EPS;
+		alpha = random_uniform( STDLIB_CONSTANT_FLOAT64_EPS, 10.0 );
+		beta = random_uniform( STDLIB_CONSTANT_FLOAT64_EPS, 10.0 );
 		y = stdlib_base_dists_invgamma_mode( alpha, beta );
 		printf( "α: %lf, β: %lf, mode(X;α,β): %lf\n", alpha, beta, y );
 	}
diff --git a/lib/node_modules/@stdlib/stats/base/dists/invgamma/skewness/README.md b/lib/node_modules/@stdlib/stats/base/dists/invgamma/skewness/README.md
@@ -221,8 +221,8 @@ int main( void ) {
     int i;
 
     for ( i = 0; i < 25; i++ ) {
-        alpha = random_uniform( 0.0, 10.0 ) + STDLIB_CONSTANT_FLOAT64_EPS;
-        beta = random_uniform( 0.0, 10.0 ) + STDLIB_CONSTANT_FLOAT64_EPS;
+        alpha = random_uniform( STDLIB_CONSTANT_FLOAT64_EPS, 10.0 );
+        beta = random_uniform( STDLIB_CONSTANT_FLOAT64_EPS, 10.0 );
         y = stdlib_base_dists_invgamma_skewness( alpha, beta );
         printf( "α: %lf, β: %lf, skew(X;α,β): %lf\n", alpha, beta, y );
     }
diff --git a/lib/node_modules/@stdlib/stats/base/dists/invgamma/skewness/benchmark/c/benchmark.c b/lib/node_modules/@stdlib/stats/base/dists/invgamma/skewness/benchmark/c/benchmark.c
@@ -101,8 +101,8 @@ static double benchmark( void ) {
 	int i;
 
 	for ( i = 0; i < 100; i++ ) {
-		alpha[ i ] = random_uniform( 0.0, 10.0 ) + 3.0 + STDLIB_CONSTANT_FLOAT64_EPS;
-		beta[ i ] = random_uniform( 0.0, 10.0 ) + STDLIB_CONSTANT_FLOAT64_EPS;
+		alpha[ i ] = random_uniform( 3.0, 13.0 );
+		beta[ i ] = random_uniform( STDLIB_CONSTANT_FLOAT64_EPS, 10.0 );
 	}
 
 	t = tic();
diff --git a/lib/node_modules/@stdlib/stats/base/dists/invgamma/skewness/examples/c/example.c b/lib/node_modules/@stdlib/stats/base/dists/invgamma/skewness/examples/c/example.c
@@ -33,8 +33,8 @@ int main( void ) {
 	int i;
 
 	for ( i = 0; i < 25; i++ ) {
-		alpha = random_uniform( 0.0, 10.0 ) + STDLIB_CONSTANT_FLOAT64_EPS;
-		beta = random_uniform( 0.0, 10.0 ) + STDLIB_CONSTANT_FLOAT64_EPS;
+		alpha = random_uniform( STDLIB_CONSTANT_FLOAT64_EPS, 10.0 );
+		beta = random_uniform( STDLIB_CONSTANT_FLOAT64_EPS, 10.0 );
 		y = stdlib_base_dists_invgamma_skewness( alpha, beta );
 		printf( "α: %lf, β: %lf, skew(X;α,β): %lf\n", alpha, beta, y );
 	}
diff --git a/lib/node_modules/@stdlib/stats/base/dists/logistic/median/benchmark/c/benchmark.c b/lib/node_modules/@stdlib/stats/base/dists/logistic/median/benchmark/c/benchmark.c
@@ -101,8 +101,8 @@ static double benchmark( void ) {
 	int i;
 
 	for ( i = 0; i < 100; i++ ) {
-		mu[ i ] = random_uniform( 0.0, 100.0 ) - 50.0;
-		s[ i ] = random_uniform( 0.0, 20.0 ) + STDLIB_CONSTANT_FLOAT64_EPS;
+		mu[ i ] = random_uniform( -50.0, 50.0 );
+		s[ i ] = random_uniform( STDLIB_CONSTANT_FLOAT64_EPS, 20.0 );
 	}
 
 	t = tic();
diff --git a/lib/node_modules/@stdlib/stats/base/dists/logistic/mode/examples/c/example.c b/lib/node_modules/@stdlib/stats/base/dists/logistic/mode/examples/c/example.c
@@ -32,7 +32,7 @@ int main( void ) {
 	int i;
 
 	for ( i = 0; i < 25; i++ ) {
-		mu = random_uniform( 0.0, 10.0 ) - 5.0;
+		mu = random_uniform( -5.0, 5.0 );
 		s = random_uniform( 0.0, 20.0 );
 		y = stdlib_base_dists_logistic_mode( mu, s );
 		printf( "µ: %lf, s: %lf, mode(X;µ,s): %lf\n", mu, s, y );
diff --git a/lib/node_modules/@stdlib/stats/base/dists/rayleigh/mean/benchmark/c/benchmark.c b/lib/node_modules/@stdlib/stats/base/dists/rayleigh/mean/benchmark/c/benchmark.c
@@ -100,7 +100,7 @@ static double benchmark( void ) {
 	int i;
 
 	for ( i = 0; i < 100; i++ ) {
-		sigma[ i ] = random_uniform( 0.0, 20.0 ) + STDLIB_CONSTANT_FLOAT64_EPS;
+		sigma[ i ] = random_uniform( STDLIB_CONSTANT_FLOAT64_EPS, 20.0 );
 	}
 
 	t = tic();
diff --git a/lib/node_modules/@stdlib/stats/base/dists/rayleigh/median/benchmark/c/benchmark.c b/lib/node_modules/@stdlib/stats/base/dists/rayleigh/median/benchmark/c/benchmark.c
@@ -100,7 +100,7 @@ static double benchmark( void ) {
 	int i;
 
 	for ( i = 0; i < 100; i++ ) {
-		sigma[ i ] = random_uniform( 0.0, 20.0 ) + STDLIB_CONSTANT_FLOAT64_EPS;
+		sigma[ i ] = random_uniform( STDLIB_CONSTANT_FLOAT64_EPS, 20.0 );
 	}
 
 	t = tic();
diff --git a/lib/node_modules/@stdlib/stats/base/dists/rayleigh/mode/benchmark/c/benchmark.c b/lib/node_modules/@stdlib/stats/base/dists/rayleigh/mode/benchmark/c/benchmark.c
@@ -100,7 +100,7 @@ static double benchmark( void ) {
 	int i;
 
 	for ( i = 0; i < 100; i++ ) {
-		sigma[ i ] = random_uniform( 0.0, 20.0 ) + STDLIB_CONSTANT_FLOAT64_EPS;
+		sigma[ i ] = random_uniform( STDLIB_CONSTANT_FLOAT64_EPS, 20.0 );
 	}
 
 	t = tic();
diff --git a/lib/node_modules/@stdlib/stats/base/dists/rayleigh/variance/benchmark/c/benchmark.c b/lib/node_modules/@stdlib/stats/base/dists/rayleigh/variance/benchmark/c/benchmark.c

Original file line number	Diff line number	Diff line change
`@@ -101,8 +101,8 @@ static double benchmark( void ) {`
`101`	`101`	`int i;`
`102`	`102`
`103`	`103`	`for ( i = 0; i < 100; i++ ) {`
`104`		`- alpha[ i ] = random_uniform( 0.0, 10.0 ) + STDLIB_CONSTANT_FLOAT64_EPS;`
`105`		`- beta[ i ] = random_uniform( 0.0, 10.0 ) + STDLIB_CONSTANT_FLOAT64_EPS;`
	`104`	`+ alpha[ i ] = random_uniform( STDLIB_CONSTANT_FLOAT64_EPS, 10.0 );`
	`105`	`+ beta[ i ] = random_uniform( STDLIB_CONSTANT_FLOAT64_EPS, 10.0 );`
`106`	`106`	`}`
`107`	`107`
`108`	`108`	`t = tic();`
Original file line number	Diff line number	Diff line change
`@@ -101,7 +101,7 @@ static double benchmark( void ) {`
`101`	`101`
`102`	`102`	`for ( i = 0; i < 100; i++ ) {`
`103`	`103`	`alpha[ i ] = random_uniform( 0.0, 20.0 );`
`104`		`- beta[ i ] = random_uniform( 0.0, 20.0 ) + 4.0;`
	`104`	`+ beta[ i ] = random_uniform( 4.0, 24.0 );`
`105`	`105`	`}`
`106`	`106`
`107`	`107`	`t = tic();`
Original file line number	Diff line number	Diff line change
`@@ -34,7 +34,7 @@ int main( void ) {`
`34`	`34`
`35`	`35`	`for ( i = 0; i < 25; i++ ) {`
`36`	`36`	`x0 = random_uniform( 0.0, 100.0 );`
`37`		`- gamma = random_uniform( 0.0, 10.0 ) + STDLIB_CONSTANT_FLOAT64_EPS;`
	`37`	`+ gamma = random_uniform( STDLIB_CONSTANT_FLOAT64_EPS, 10.0 );`
`38`	`38`	`y = stdlib_base_dists_cauchy_entropy( x0, gamma );`
`39`	`39`	`printf( "x0: %lf, gamma: %lf, F(x0;gamma): %lf\n", x0, gamma, y );`
`40`	`40`	`}`