From 4786cbcd24c585a043f009c607538570614bbed0 Mon Sep 17 00:00:00 2001
From: Claude <noreply@anthropic.com>
Date: Thu, 27 Aug 2026 14:31:57 +0000
Subject: [PATCH] test: migrate `stats/base/dists/beta/quantile` to ULP-based
 assertions

Ref: https://github.com/stdlib-js/stdlib/issues/11352

Co-Authored-By: Claude Opus 5 <noreply@anthropic.com>
Claude-Session: https://claude.ai/code/session_014Cv2JMqyGWr6vf5Yd2aGyH
---
 .../dists/beta/quantile/test/test.factory.js  | 33 +++----------------
 .../dists/beta/quantile/test/test.quantile.js | 33 +++----------------
 2 files changed, 8 insertions(+), 58 deletions(-)

diff --git a/lib/node_modules/@stdlib/stats/base/dists/beta/quantile/test/test.factory.js b/lib/node_modules/@stdlib/stats/base/dists/beta/quantile/test/test.factory.js
index 67692c858654..370d18f63716 100644
--- a/lib/node_modules/@stdlib/stats/base/dists/beta/quantile/test/test.factory.js
+++ b/lib/node_modules/@stdlib/stats/base/dists/beta/quantile/test/test.factory.js
@@ -21,11 +21,10 @@
 // MODULES //
 
 var tape = require( 'tape' );
+var isAlmostSameValue = require( '@stdlib/assert/is-almost-same-value' );
 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 factory = require( './../lib/factory.js' );
 
 
@@ -161,10 +160,8 @@ tape( 'if provided a nonpositive `alpha`, the created function always returns `N
 tape( 'the created function evaluates the quantile for `p` given large `alpha` and `beta`', function test( t ) {
 	var expected;
 	var quantile;
-	var delta;
 	var alpha;
 	var beta;
-	var tol;
 	var p;
 	var y;
 	var i;
@@ -176,13 +173,7 @@ tape( 'the created function evaluates the quantile for `p` given large `alpha` a
 	for ( i = 0; i < p.length; i++ ) {
 		quantile = factory( alpha[i], beta[i] );
 		y = quantile( p[i] );
-		if ( y === expected[i] ) {
-			t.strictEqual( y, expected[i], 'p: '+p[i]+', alpha: '+alpha[i]+', beta: '+beta[i]+', y: '+y+', expected: '+expected[i] );
-		} else {
-			delta = abs( y - expected[ i ] );
-			tol = 20.0 * EPS * abs( expected[ i ] );
-			t.ok( delta <= tol, 'within tolerance. p: '+p[ i ]+'. alpha: '+alpha[i]+'. beta: '+beta[i]+'. y: '+y+'. E: '+expected[ i ]+'. Δ: '+delta+'. tol: '+tol+'.' );
-		}
+		t.strictEqual( isAlmostSameValue( y, expected[ i ], 18 ), true, 'returns expected value' );
 	}
 	t.end();
 });
@@ -190,10 +181,8 @@ tape( 'the created function evaluates the quantile for `p` given large `alpha` a
 tape( 'the created function evaluates the quantile for `p` given large `alpha`', function test( t ) {
 	var expected;
 	var quantile;
-	var delta;
 	var alpha;
 	var beta;
-	var tol;
 	var p;
 	var y;
 	var i;
@@ -205,13 +194,7 @@ tape( 'the created function evaluates the quantile for `p` given large `alpha`',
 	for ( i = 0; i < p.length; i++ ) {
 		quantile = factory( alpha[i], beta[i] );
 		y = quantile( p[i] );
-		if ( y === expected[i] ) {
-			t.strictEqual( y, expected[i], 'p: '+p[i]+', alpha: '+alpha[i]+', beta: '+beta[i]+', y: '+y+', expected: '+expected[i] );
-		} else {
-			delta = abs( y - expected[ i ] );
-			tol = 20.0 * EPS * abs( expected[ i ] );
-			t.ok( delta <= tol, 'within tolerance. p: '+p[ i ]+'. alpha: '+alpha[i]+'. beta: '+beta[i]+'. y: '+y+'. E: '+expected[ i ]+'. Δ: '+delta+'. tol: '+tol+'.' );
-		}
+		t.strictEqual( isAlmostSameValue( y, expected[ i ], 9 ), true, 'returns expected value' );
 	}
 	t.end();
 });
@@ -219,10 +202,8 @@ tape( 'the created function evaluates the quantile for `p` given large `alpha`',
 tape( 'the created function evaluates the quantile for `p` given large `beta`', function test( t ) {
 	var expected;
 	var quantile;
-	var delta;
 	var alpha;
 	var beta;
-	var tol;
 	var p;
 	var y;
 	var i;
@@ -234,13 +215,7 @@ tape( 'the created function evaluates the quantile for `p` given large `beta`',
 	for ( i = 0; i < p.length; i++ ) {
 		quantile = factory( alpha[i], beta[i] );
 		y = quantile( p[i] );
-		if ( y === expected[i] ) {
-			t.strictEqual( y, expected[i], 'p: '+p[i]+', alpha: '+alpha[i]+', beta: '+beta[i]+', y: '+y+', expected: '+expected[i] );
-		} else {
-			delta = abs( y - expected[ i ] );
-			tol = 20.0 * EPS * abs( expected[ i ] );
-			t.ok( delta <= tol, 'within tolerance. p: '+p[ i ]+'. alpha: '+alpha[i]+'. beta: '+beta[i]+'. y: '+y+'. E: '+expected[ i ]+'. Δ: '+delta+'. tol: '+tol+'.' );
-		}
+		t.strictEqual( isAlmostSameValue( y, expected[ i ], 29 ), true, 'returns expected value' );
 	}
 	t.end();
 });
diff --git a/lib/node_modules/@stdlib/stats/base/dists/beta/quantile/test/test.quantile.js b/lib/node_modules/@stdlib/stats/base/dists/beta/quantile/test/test.quantile.js
index 584c71305003..6586162a72b6 100644
--- a/lib/node_modules/@stdlib/stats/base/dists/beta/quantile/test/test.quantile.js
+++ b/lib/node_modules/@stdlib/stats/base/dists/beta/quantile/test/test.quantile.js
@@ -21,11 +21,10 @@
 // MODULES //
 
 var tape = require( 'tape' );
+var isAlmostSameValue = require( '@stdlib/assert/is-almost-same-value' );
 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 quantile = require( './../lib' );
 
 
@@ -118,10 +117,8 @@ tape( 'if provided a nonpositive `beta`, the function returns `NaN`', function t
 
 tape( 'the function evaluates the quantile for `x` given large `alpha` and `beta`', function test( t ) {
 	var expected;
-	var delta;
 	var alpha;
 	var beta;
-	var tol;
 	var p;
 	var y;
 	var i;
@@ -132,13 +129,7 @@ tape( 'the function evaluates the quantile for `x` given large `alpha` and `beta
 	beta = bothLarge.beta;
 	for ( i = 0; i < p.length; i++ ) {
 		y = quantile( p[i], alpha[i], beta[i] );
-		if ( y === expected[i] ) {
-			t.strictEqual( y, expected[i], 'p: '+p[i]+', alpha: '+alpha[i]+', beta: '+beta[i]+', y: '+y+', expected: '+expected[i] );
-		} else {
-			delta = abs( y - expected[ i ] );
-			tol = 20.0 * EPS * abs( expected[ i ] );
-			t.ok( delta <= tol, 'within tolerance. p: '+p[ i ]+'. alpha: '+alpha[i]+'. beta: '+beta[i]+'. y: '+y+'. E: '+expected[ i ]+'. Δ: '+delta+'. tol: '+tol+'.' );
-		}
+		t.strictEqual( isAlmostSameValue( y, expected[ i ], 18 ), true, 'returns expected value' );
 	}
 	t.end();
 });
@@ -146,9 +137,7 @@ tape( 'the function evaluates the quantile for `x` given large `alpha` and `beta
 tape( 'the function evaluates the quantile for `x` given large `alpha`', function test( t ) {
 	var expected;
 	var alpha;
-	var delta;
 	var beta;
-	var tol;
 	var p;
 	var y;
 	var i;
@@ -159,13 +148,7 @@ tape( 'the function evaluates the quantile for `x` given large `alpha`', functio
 	beta = largeAlpha.beta;
 	for ( i = 0; i < p.length; i++ ) {
 		y = quantile( p[i], alpha[i], beta[i] );
-		if ( y === expected[i] ) {
-			t.strictEqual( y, expected[i], 'p: '+p[i]+', alpha: '+alpha[i]+', beta: '+beta[i]+', y: '+y+', expected: '+expected[i] );
-		} else {
-			delta = abs( y - expected[ i ] );
-			tol = 20.0 * EPS * abs( expected[ i ] );
-			t.ok( delta <= tol, 'within tolerance. p: '+p[ i ]+'. alpha: '+alpha[i]+'. beta: '+beta[i]+'. y: '+y+'. E: '+expected[ i ]+'. Δ: '+delta+'. tol: '+tol+'.' );
-		}
+		t.strictEqual( isAlmostSameValue( y, expected[ i ], 9 ), true, 'returns expected value' );
 	}
 	t.end();
 });
@@ -173,9 +156,7 @@ tape( 'the function evaluates the quantile for `x` given large `alpha`', functio
 tape( 'the function evaluates the quantile for `x` given large `beta`', function test( t ) {
 	var expected;
 	var alpha;
-	var delta;
 	var beta;
-	var tol;
 	var p;
 	var y;
 	var i;
@@ -186,13 +167,7 @@ tape( 'the function evaluates the quantile for `x` given large `beta`', function
 	beta = largeBeta.beta;
 	for ( i = 0; i < p.length; i++ ) {
 		y = quantile( p[i], alpha[i], beta[i] );
-		if ( y === expected[i] ) {
-			t.strictEqual( y, expected[i], 'p: '+p[i]+', alpha: '+alpha[i]+', beta: '+beta[i]+', y: '+y+', expected: '+expected[i] );
-		} else {
-			delta = abs( y - expected[ i ] );
-			tol = 20.0 * EPS * abs( expected[ i ] );
-			t.ok( delta <= tol, 'within tolerance. p: '+p[ i ]+'. alpha: '+alpha[i]+'. beta: '+beta[i]+'. y: '+y+'. E: '+expected[ i ]+'. Δ: '+delta+'. tol: '+tol+'.' );
-		}
+		t.strictEqual( isAlmostSameValue( y, expected[ i ], 29 ), true, 'returns expected value' );
 	}
 	t.end();
 });