From 6a7e5a35423f72366588aa822016b2478ab7c90c Mon Sep 17 00:00:00 2001 From: Philipp Burckhardt Date: Wed, 10 Jun 2026 02:17:33 -0500 Subject: [PATCH] test: migrate `math/base/special/rsqrtf` to ULP-based testing --- 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_pkg_readmes status: na - task: lint_markdown_docs status: na - task: lint_markdown status: na - 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: passed - 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: na - task: lint_c_benchmarks status: na - task: lint_c_tests_fixtures status: na - task: lint_shell status: na - task: lint_typescript_declarations status: passed - task: lint_typescript_tests status: na - task: lint_license_headers status: passed --- --- .../math/base/special/rsqrtf/test/test.js | 100 +++++------------- .../base/special/rsqrtf/test/test.native.js | 100 +++++------------- 2 files changed, 52 insertions(+), 148 deletions(-) diff --git a/lib/node_modules/@stdlib/math/base/special/rsqrtf/test/test.js b/lib/node_modules/@stdlib/math/base/special/rsqrtf/test/test.js index ab536c6e8d2b..c736b2fdad16 100644 --- a/lib/node_modules/@stdlib/math/base/special/rsqrtf/test/test.js +++ b/lib/node_modules/@stdlib/math/base/special/rsqrtf/test/test.js @@ -24,8 +24,8 @@ var tape = require( 'tape' ); var isnanf = require( '@stdlib/math/base/assert/is-nanf' ); var PINF = require( '@stdlib/constants/float32/pinf' ); var NINF = require( '@stdlib/constants/float32/ninf' ); -var EPS = require( '@stdlib/constants/float32/eps' ); -var abs = require( '@stdlib/math/base/special/abs' ); +var f32 = require( '@stdlib/number/float64/base/to-float32' ); +var isAlmostSameValue = require( '@stdlib/number/float32/base/assert/is-almost-same-value' ); var isPositiveZero = require( '@stdlib/math/base/assert/is-positive-zero' ); var rsqrtf = require( './../lib' ); @@ -52,184 +52,136 @@ tape( 'main export is a function', function test( t ) { tape( 'the function evaluates the reciprocal square root of `x` on the interval `[50,500]`', function test( t ) { var expected; - var delta; - var tol; var x; var y; var i; + var e; expected = veryLargePositive.expected; x = veryLargePositive.x; for ( i = 0; i < x.length; i++ ) { y = rsqrtf( x[i] ); - if ( y === expected[i] ) { - t.strictEqual( y, expected[i], 'x: '+x[i]+', y: '+y+', expected: '+expected[i] ); - } else { - delta = abs( y - expected[i] ); - tol = EPS * abs( expected[i] ); - t.strictEqual( delta <= tol, true, 'within tolerance. x: '+x[i]+'. y: '+y+'. E: '+expected[i]+'. Δ: '+delta+'. Tolerance: '+tol+'.' ); - } + e = f32( expected[ i ] ); + t.strictEqual( isAlmostSameValue( y, e, 1 ), true, 'returns expected value' ); } t.end(); }); tape( 'the function evaluates the reciprocal square root of `x` on the interval `[20,50]`', function test( t ) { var expected; - var delta; - var tol; var x; var y; var i; + var e; expected = largePositive.expected; x = largePositive.x; for ( i = 0; i < x.length; i++ ) { y = rsqrtf( x[i] ); - if ( y === expected[i] ) { - t.strictEqual( y, expected[i], 'x: '+x[i]+', y: '+y+', expected: '+expected[i] ); - } else { - delta = abs( y - expected[i] ); - tol = EPS * abs( expected[i] ); - t.strictEqual( delta <= tol, true, 'within tolerance. x: '+x[i]+'. y: '+y+'. E: '+expected[i]+'. Δ: '+delta+'. Tolerance: '+tol+'.' ); - } + e = f32( expected[ i ] ); + t.strictEqual( isAlmostSameValue( y, e, 1 ), true, 'returns expected value' ); } t.end(); }); tape( 'the function evaluates the reciprocal square root of `x` on the interval `[3,20]`', function test( t ) { var expected; - var delta; - var tol; var x; var y; var i; + var e; expected = mediumPositive.expected; x = mediumPositive.x; for ( i = 0; i < x.length; i++ ) { y = rsqrtf( x[i] ); - if ( y === expected[i] ) { - t.strictEqual( y, expected[i], 'x: '+x[i]+', y: '+y+', expected: '+expected[i] ); - } else { - delta = abs( y - expected[i] ); - tol = EPS * abs( expected[i] ); - t.strictEqual( delta <= tol, true, 'within tolerance. x: '+x[i]+'. y: '+y+'. E: '+expected[i]+'. Δ: '+delta+'. Tolerance: '+tol+'.' ); - } + e = f32( expected[ i ] ); + t.strictEqual( isAlmostSameValue( y, e, 1 ), true, 'returns expected value' ); } t.end(); }); tape( 'the function evaluates the reciprocal square root of `x` on the interval `[0.8,3]`', function test( t ) { var expected; - var delta; - var tol; var x; var y; var i; + var e; expected = smallPositive.expected; x = smallPositive.x; for ( i = 0; i < x.length; i++ ) { y = rsqrtf( x[i] ); - if ( y === expected[i] ) { - t.strictEqual( y, expected[i], 'x: '+x[i]+', y: '+y+', expected: '+expected[i] ); - } else { - delta = abs( y - expected[i] ); - tol = EPS * abs( expected[i] ); - t.strictEqual( delta <= tol, true, 'within tolerance. x: '+x[i]+'. y: '+y+'. E: '+expected[i]+'. Δ: '+delta+'. Tolerance: '+tol+'.' ); - } + e = f32( expected[ i ] ); + t.strictEqual( isAlmostSameValue( y, e, 1 ), true, 'returns expected value' ); } t.end(); }); tape( 'the function evaluates the reciprocal square root of `x` on the interval `[0.0,0.8]`', function test( t ) { var expected; - var delta; - var tol; var x; var y; var i; + var e; expected = smaller.expected; x = smaller.x; for ( i = 0; i < x.length; i++ ) { y = rsqrtf( x[i] ); - if ( y === expected[i] ) { - t.strictEqual( y, expected[i], 'x: '+x[i]+', y: '+y+', expected: '+expected[i] ); - } else { - delta = abs( y - expected[i] ); - tol = EPS * abs( expected[i] ); - t.strictEqual( delta <= tol, true, 'within tolerance. x: '+x[i]+'. y: '+y+'. E: '+expected[i]+'. Δ: '+delta+'. Tolerance: '+tol+'.' ); - } + e = f32( expected[ i ] ); + t.strictEqual( isAlmostSameValue( y, e, 1 ), true, 'returns expected value' ); } t.end(); }); tape( 'the function evaluates the reciprocal square root of `x` on the interval `[1e-300,1e-308]`', function test( t ) { var expected; - var delta; - var tol; var x; var y; var i; + var e; expected = tinyPositive.expected; x = tinyPositive.x; for ( i = 0; i < x.length; i++ ) { y = rsqrtf( x[i] ); - if ( y === expected[i] ) { - t.strictEqual( y, expected[i], 'x: '+x[i]+', y: '+y+', expected: '+expected[i] ); - } else { - delta = abs( y - expected[i] ); - tol = EPS * abs( expected[i] ); - t.strictEqual( delta <= tol, true, 'within tolerance. x: '+x[i]+'. y: '+y+'. E: '+expected[i]+'. Δ: '+delta+'. Tolerance: '+tol+'.' ); - } + e = f32( expected[ i ] ); + t.strictEqual( isAlmostSameValue( y, e, 1 ), true, 'returns expected value' ); } t.end(); }); tape( 'the function evaluates the reciprocal square root of subnormal `x`', function test( t ) { var expected; - var delta; - var tol; var x; var y; var i; + var e; expected = subnormal.expected; x = subnormal.x; for ( i = 0; i < x.length; i++ ) { y = rsqrtf( x[i] ); - if ( y === expected[i] ) { - t.strictEqual( y, expected[i], 'x: '+x[i]+', y: '+y+', expected: '+expected[i] ); - } else { - delta = abs( y - expected[i] ); - tol = EPS * abs( expected[i] ); - t.strictEqual( delta <= tol, true, 'within tolerance. x: '+x[i]+'. y: '+y+'. E: '+expected[i]+'. Δ: '+delta+'. Tolerance: '+tol+'.' ); - } + e = f32( expected[ i ] ); + t.strictEqual( isAlmostSameValue( y, e, 1 ), true, 'returns expected value' ); } t.end(); }); tape( 'the function evaluates the reciprocal square root of `x` (huge positive)', function test( t ) { var expected; - var delta; - var tol; var x; var y; var i; + var e; expected = hugePositive.expected; x = hugePositive.x; for ( i = 0; i < x.length; i++ ) { y = rsqrtf( x[i] ); - if ( y === expected[i] ) { - t.strictEqual( y, expected[i], 'x: '+x[i]+', y: '+y+', expected: '+expected[i] ); - } else { - delta = abs( y - expected[i] ); - tol = EPS * abs( expected[i] ); - t.strictEqual( delta <= tol, true, 'within tolerance. x: '+x[i]+'. y: '+y+'. E: '+expected[i]+'. Δ: '+delta+'. Tolerance: '+tol+'.' ); - } + e = f32( expected[ i ] ); + t.strictEqual( isAlmostSameValue( y, e, 1 ), true, 'returns expected value' ); } t.end(); }); diff --git a/lib/node_modules/@stdlib/math/base/special/rsqrtf/test/test.native.js b/lib/node_modules/@stdlib/math/base/special/rsqrtf/test/test.native.js index 1ac1848a3c9a..651c172e9d3b 100644 --- a/lib/node_modules/@stdlib/math/base/special/rsqrtf/test/test.native.js +++ b/lib/node_modules/@stdlib/math/base/special/rsqrtf/test/test.native.js @@ -25,8 +25,8 @@ var tape = require( 'tape' ); var isnanf = require( '@stdlib/math/base/assert/is-nanf' ); var PINF = require( '@stdlib/constants/float32/pinf' ); var NINF = require( '@stdlib/constants/float32/ninf' ); -var EPS = require( '@stdlib/constants/float32/eps' ); -var abs = require( '@stdlib/math/base/special/abs' ); +var f32 = require( '@stdlib/number/float64/base/to-float32' ); +var isAlmostSameValue = require( '@stdlib/number/float32/base/assert/is-almost-same-value' ); var isPositiveZero = require( '@stdlib/math/base/assert/is-positive-zero' ); var tryRequire = require( '@stdlib/utils/try-require' ); @@ -61,184 +61,136 @@ tape( 'main export is a function', opts, function test( t ) { tape( 'the function evaluates the reciprocal square root of `x` on the interval `[50,500]`', opts, function test( t ) { var expected; - var delta; - var tol; var x; var y; var i; + var e; expected = veryLargePositive.expected; x = veryLargePositive.x; for ( i = 0; i < x.length; i++ ) { y = rsqrtf( x[i] ); - if ( y === expected[i] ) { - t.strictEqual( y, expected[i], 'x: '+x[i]+', y: '+y+', expected: '+expected[i] ); - } else { - delta = abs( y - expected[i] ); - tol = EPS * abs( expected[i] ); - t.strictEqual( delta <= tol, true, 'within tolerance. x: '+x[i]+'. y: '+y+'. E: '+expected[i]+'. Δ: '+delta+'. Tolerance: '+tol+'.' ); - } + e = f32( expected[ i ] ); + t.strictEqual( isAlmostSameValue( y, e, 1 ), true, 'returns expected value' ); } t.end(); }); tape( 'the function evaluates the reciprocal square root of `x` on the interval `[20,50]`', opts, function test( t ) { var expected; - var delta; - var tol; var x; var y; var i; + var e; expected = largePositive.expected; x = largePositive.x; for ( i = 0; i < x.length; i++ ) { y = rsqrtf( x[i] ); - if ( y === expected[i] ) { - t.strictEqual( y, expected[i], 'x: '+x[i]+', y: '+y+', expected: '+expected[i] ); - } else { - delta = abs( y - expected[i] ); - tol = EPS * abs( expected[i] ); - t.strictEqual( delta <= tol, true, 'within tolerance. x: '+x[i]+'. y: '+y+'. E: '+expected[i]+'. Δ: '+delta+'. Tolerance: '+tol+'.' ); - } + e = f32( expected[ i ] ); + t.strictEqual( isAlmostSameValue( y, e, 1 ), true, 'returns expected value' ); } t.end(); }); tape( 'the function evaluates the reciprocal square root of `x` on the interval `[3,20]`', opts, function test( t ) { var expected; - var delta; - var tol; var x; var y; var i; + var e; expected = mediumPositive.expected; x = mediumPositive.x; for ( i = 0; i < x.length; i++ ) { y = rsqrtf( x[i] ); - if ( y === expected[i] ) { - t.strictEqual( y, expected[i], 'x: '+x[i]+', y: '+y+', expected: '+expected[i] ); - } else { - delta = abs( y - expected[i] ); - tol = EPS * abs( expected[i] ); - t.strictEqual( delta <= tol, true, 'within tolerance. x: '+x[i]+'. y: '+y+'. E: '+expected[i]+'. Δ: '+delta+'. Tolerance: '+tol+'.' ); - } + e = f32( expected[ i ] ); + t.strictEqual( isAlmostSameValue( y, e, 1 ), true, 'returns expected value' ); } t.end(); }); tape( 'the function evaluates the reciprocal square root of `x` on the interval `[0.8,3]`', opts, function test( t ) { var expected; - var delta; - var tol; var x; var y; var i; + var e; expected = smallPositive.expected; x = smallPositive.x; for ( i = 0; i < x.length; i++ ) { y = rsqrtf( x[i] ); - if ( y === expected[i] ) { - t.strictEqual( y, expected[i], 'x: '+x[i]+', y: '+y+', expected: '+expected[i] ); - } else { - delta = abs( y - expected[i] ); - tol = EPS * abs( expected[i] ); - t.strictEqual( delta <= tol, true, 'within tolerance. x: '+x[i]+'. y: '+y+'. E: '+expected[i]+'. Δ: '+delta+'. Tolerance: '+tol+'.' ); - } + e = f32( expected[ i ] ); + t.strictEqual( isAlmostSameValue( y, e, 1 ), true, 'returns expected value' ); } t.end(); }); tape( 'the function evaluates the reciprocal square root of `x` on the interval `[0.0,0.8]`', opts, function test( t ) { var expected; - var delta; - var tol; var x; var y; var i; + var e; expected = smaller.expected; x = smaller.x; for ( i = 0; i < x.length; i++ ) { y = rsqrtf( x[i] ); - if ( y === expected[i] ) { - t.strictEqual( y, expected[i], 'x: '+x[i]+', y: '+y+', expected: '+expected[i] ); - } else { - delta = abs( y - expected[i] ); - tol = EPS * abs( expected[i] ); - t.strictEqual( delta <= tol, true, 'within tolerance. x: '+x[i]+'. y: '+y+'. E: '+expected[i]+'. Δ: '+delta+'. Tolerance: '+tol+'.' ); - } + e = f32( expected[ i ] ); + t.strictEqual( isAlmostSameValue( y, e, 1 ), true, 'returns expected value' ); } t.end(); }); tape( 'the function evaluates the reciprocal square root of `x` on the interval `[1e-300,1e-308]`', opts, function test( t ) { var expected; - var delta; - var tol; var x; var y; var i; + var e; expected = tinyPositive.expected; x = tinyPositive.x; for ( i = 0; i < x.length; i++ ) { y = rsqrtf( x[i] ); - if ( y === expected[i] ) { - t.strictEqual( y, expected[i], 'x: '+x[i]+', y: '+y+', expected: '+expected[i] ); - } else { - delta = abs( y - expected[i] ); - tol = EPS * abs( expected[i] ); - t.strictEqual( delta <= tol, true, 'within tolerance. x: '+x[i]+'. y: '+y+'. E: '+expected[i]+'. Δ: '+delta+'. Tolerance: '+tol+'.' ); - } + e = f32( expected[ i ] ); + t.strictEqual( isAlmostSameValue( y, e, 1 ), true, 'returns expected value' ); } t.end(); }); tape( 'the function evaluates the reciprocal square root of subnormal `x`', opts, function test( t ) { var expected; - var delta; - var tol; var x; var y; var i; + var e; expected = subnormal.expected; x = subnormal.x; for ( i = 0; i < x.length; i++ ) { y = rsqrtf( x[i] ); - if ( y === expected[i] ) { - t.strictEqual( y, expected[i], 'x: '+x[i]+', y: '+y+', expected: '+expected[i] ); - } else { - delta = abs( y - expected[i] ); - tol = EPS * abs( expected[i] ); - t.strictEqual( delta <= tol, true, 'within tolerance. x: '+x[i]+'. y: '+y+'. E: '+expected[i]+'. Δ: '+delta+'. Tolerance: '+tol+'.' ); - } + e = f32( expected[ i ] ); + t.strictEqual( isAlmostSameValue( y, e, 1 ), true, 'returns expected value' ); } t.end(); }); tape( 'the function evaluates the reciprocal square root of `x` (huge positive)', opts, function test( t ) { var expected; - var delta; - var tol; var x; var y; var i; + var e; expected = hugePositive.expected; x = hugePositive.x; for ( i = 0; i < x.length; i++ ) { y = rsqrtf( x[i] ); - if ( y === expected[i] ) { - t.strictEqual( y, expected[i], 'x: '+x[i]+', y: '+y+', expected: '+expected[i] ); - } else { - delta = abs( y - expected[i] ); - tol = EPS * abs( expected[i] ); - t.strictEqual( delta <= tol, true, 'within tolerance. x: '+x[i]+'. y: '+y+'. E: '+expected[i]+'. Δ: '+delta+'. Tolerance: '+tol+'.' ); - } + e = f32( expected[ i ] ); + t.strictEqual( isAlmostSameValue( y, e, 1 ), true, 'returns expected value' ); } t.end(); });