Skip to content

test: migrate math/base/special/betainc to ULP-based testing #11713

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Open
MANDEep22332 wants to merge 3 commits into
base: develop
Choose a base branch
from
+38 −29
Open

test: migrate math/base/special/betainc to ULP-based testing #11713

Changes from all commits
Commits
Show all changes
3 commits
Select commit Hold shift + click to select a range
bbc7228
ulp-based test for the betainc
MANDEep22332 Apr 22, 2026
bba279b
chore: trigger CI
MANDEep22332 Apr 22, 2026
879456a
chore: trigger CI
MANDEep22332 Apr 22, 2026
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
67 changes: 38 additions & 29 deletions lib/node_modules/@stdlib/math/base/special/betainc/test/test.js
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -23,7 +23,7 @@
var tape = require( 'tape' );
var isnan = require( '@stdlib/math/base/assert/is-nan' );
var abs = require( '@stdlib/math/base/special/abs' );
var EPS = require( '@stdlib/constants/float64/eps' );
var isAlmostSameValue = require( '@stdlib/assert/is-almost-same-value' );
var betainc = require( './../lib' );


Expand Down Expand Up @@ -71,83 +71,92 @@ tape( 'the function returns `NaN` negative `a` or `b`', function test( t ) {

tape( 'the function evaluates the lower regularized incomplete beta function', function test( t ) {
var expected;
var delta;
var tol;
var i;
var y;

expected = lowerRegularized;
for ( i = 0; i < x.length; i++ ) {
y = betainc( x[i], a[i], b[i] );
if ( y === expected[i] ) {
t.strictEqual( y, expected[i], 'x: '+x[i]+', y: '+y+'. a: '+a[i]+'. b: '+b[i]+', expected: '+expected[i] );
if ( isnan( expected[ i ] ) ) {
t.strictEqual( isnan( y ), true, 'returns expected value' );
} else if ( expected[ i ] === 0.0 ) {
t.strictEqual( y, 0.0, 'returns expected value' );
} else if ( abs( expected[ i ] ) < 1e-300 ) {
t.ok( abs( y ) < 1e-300, 'returns expected value' );
} else {
delta = abs( y - expected[ i ] );
tol = 50.0 * EPS * abs( expected[ i ] );
t.ok( delta <= tol, 'within tolerance. x: '+x[i]+'. a: '+a[i]+'. b: '+b[i]+'. y: '+y+'. E: '+expected[i]+'. Δ: '+delta+'. tol: '+tol );
t.strictEqual(isAlmostSameValue( y, expected[ i ], 120 ), true, 'returns expected value');
}
}
t.end();
});

tape( 'the function evaluates the lower unregularized incomplete beta function', function test( t ) {
var expected;
var delta;
var tol;
var i;
var y;

expected = lowerUnregularized;
for ( i = 0; i < x.length; i++ ) {
y = betainc( x[i], a[i], b[i], false );
if ( y === expected[i] ) {
t.strictEqual( y, expected[i], 'x: '+x[i]+', y: '+y+'. a: '+a[i]+'. b: '+b[i]+', expected: '+expected[i] );
if ( isnan( expected[ i ] ) ) {
t.strictEqual( isnan( y ), true, 'returns expected value' );
} else if ( abs( expected[ i ] ) < 1e-300 ) {
t.ok( abs( y ) < 1e-300, 'returns expected value' );
} else if (
expected[ i ] === Infinity ||
expected[ i ] === -Infinity
) {
t.strictEqual( y, expected[ i ], 'returns expected value' );
} else {
delta = abs( y - expected[ i ] );
tol = 80.0 * EPS * abs( expected[ i ] );
t.ok( delta <= tol, 'within tolerance. x: '+x[i]+'. a: '+a[i]+'. b: '+b[i]+'. y: '+y+'. E: '+expected[i]+'. Δ: '+delta+'. tol: '+tol );
t.strictEqual(isAlmostSameValue( y, expected[ i ], 128 ), true, 'returns expected value');
}
}
t.end();
});

tape( 'the function evaluates the upper regularized incomplete beta function', function test( t ) {
var expected;
var delta;
var tol;
var i;
var y;

expected = upperRegularized;
for ( i = 0; i < x.length; i++ ) {
y = betainc( x[i], a[i], b[i], true, true );
if ( y === expected[i] ) {
t.strictEqual( y, expected[i], 'x: '+x[i]+', y: '+y+'. a: '+a[i]+'. b: '+b[i]+', expected: '+expected[i] );
if ( isnan( expected[ i ] ) ) {
t.strictEqual( isnan( y ), true, 'returns expected value' );
} else if ( abs( expected[ i ] ) < 1e-300 ) {
t.ok( abs( y ) < 1e-300, 'returns expected value' );
} else if (
expected[ i ] === Infinity ||
expected[ i ] === -Infinity
) {
t.strictEqual( y, expected[ i ], 'returns expected value' );
} else {
delta = abs( y - expected[ i ] );
tol = 50.0 * EPS * abs( expected[ i ] );
t.ok( delta <= tol, 'within tolerance. x: '+x[i]+'. a: '+a[i]+'. b: '+b[i]+'. y: '+y+'. E: '+expected[i]+'. Δ: '+delta+'. tol: '+tol );
t.strictEqual(isAlmostSameValue( y, expected[ i ], 80 ), true, 'returns expected value');
}
}
t.end();
});

tape( 'the function evaluates the upper unregularized incomplete beta function', function test( t ) {
var expected;
var delta;
var tol;
var i;
var y;

expected = upperUnregularized;
for ( i = 0; i < x.length; i++ ) {
y = betainc( x[i], a[i], b[i], false, true );
if ( y === expected[i] ) {
t.strictEqual( y, expected[i], 'x: '+x[i]+', y: '+y+'. a: '+a[i]+'. b: '+b[i]+', expected: '+expected[i] );
if ( isnan( expected[ i ] ) ) {
t.strictEqual( isnan( y ), true, 'returns expected value' );
} else if ( abs( expected[ i ] ) < 1e-300 ) {
t.ok( abs( y ) < 1e-300, 'returns expected value' );
} else if (
expected[ i ] === Infinity ||
expected[ i ] === -Infinity
) {
t.strictEqual( y, expected[ i ], 'returns expected value' );
} else {
delta = abs( y - expected[ i ] );
tol = 80.0 * EPS * abs( expected[ i ] );
t.ok( delta <= tol, 'within tolerance. x: '+x[i]+'. a: '+a[i]+'. b: '+b[i]+'. y: '+y+'. E: '+expected[i]+'. Δ: '+delta+'. tol: '+tol );
t.strictEqual(isAlmostSameValue( y, expected[ i ], 120 ), true, 'returns expected value');
}
}
t.end();
Expand Down