diff options
Diffstat (limited to 'js/src/tests/non262/Math/pow-approx.js')
| -rw-r--r-- | js/src/tests/non262/Math/pow-approx.js | 568 |
1 files changed, 568 insertions, 0 deletions
diff --git a/js/src/tests/non262/Math/pow-approx.js b/js/src/tests/non262/Math/pow-approx.js new file mode 100644 index 0000000000..11d009158c --- /dev/null +++ b/js/src/tests/non262/Math/pow-approx.js @@ -0,0 +1,568 @@ +if (typeof fdlibm === "undefined") { + var fdlibm = SpecialPowers.Cu.getJSTestingFunctions().fdlibm; +} + +if (typeof getBuildConfiguration === "undefined") { + var getBuildConfiguration = SpecialPowers.Cu.getJSTestingFunctions().getBuildConfiguration; +} + +const f64 = new Float64Array(1); +const ui64 = new BigUint64Array(f64.buffer); + +function toBits(n) { + f64[0] = n; + return ui64[0]; +} + +function errorInULP(actual, expected) { + // Handle NaN and +0/-0. + if (Object.is(actual, expected)) { + return 0; + } + + let x = toBits(actual); + let y = toBits(expected); + return x <= y ? Number(y - x) : Number(x - y); +} + +// Test methodology: +// +// Generate test cases for inputs where the original js::powi implementation +// returns a different result than std::pow. If such inputs where found, compare +// them against fdlibm::pow to find inputs where the error is larger than 1 ULP. +// +// Compile with: +// -std=c++17 -O3 -msse -msse2 -mfpmath=sse -fno-math-errno -fno-exceptions -fno-rtti -march=native +// +// static bool test(double x, int32_t y) { +// if (std::isnan(x)) { +// return true; +// } +// +// double t = js::powi(x, y); +// double u = std::pow(x, static_cast<double>(y)); +// if (t == u) { +// return true; +// } +// +// uint64_t a; +// std::memcpy(&a, &t, sizeof(double)); +// +// uint64_t b; +// std::memcpy(&b, &u, sizeof(double)); +// +// double v = fdlibm::pow(x, y); +// +// uint64_t c; +// std::memcpy(&c, &v, sizeof(double)); +// +// double w = musl::pow(x, y); +// +// uint64_t d; +// std::memcpy(&d, &w, sizeof(double)); +// +// // Expect at most 1 ULP difference between std::pow and fdlibm::pow. +// if ((b < c && c - b > 1) || (b > c && b - c > 1)) { +// printf("!!! [fdlibm] %.53f ** %d: 0x%" PRIx64 " != 0x%" PRIx64 "\n", x, y, b, c); +// exit(1); +// } +// +// // Expect at most 1 ULP difference between std::pow and musl::pow. +// if ((b < d && d - b > 1) || (b > d && b - d > 1)) { +// printf("!!! [musl] %.53f ** %d: 0x%" PRIx64 " != 0x%" PRIx64 "\n", x, y, b, d); +// exit(1); +// } +// +// // Accept 1 ULP difference between js::powi and fdlibm::pow. +// if ((a <= c && c - a <= 1) || (a >= c && a - c <= 1)) { +// return true; +// } +// +// // Output if a larger error was found. +// printf("%.53f ** %d: 0x%" PRIx64 " != 0x%" PRIx64 " (0x%" PRIx64 ") (0x%" PRIx64 ")\n", x, y, a, b, c, d); +// return false; +// } +// +// int main() { +// // Use mt19937 for reproducible results. +// std::mt19937_64 gen64; +// std::mt19937 gen32; +// +// for (uint64_t i = 0; i < 100'000'000'000; ++i) { +// uint64_t x = gen64(); +// int32_t y = gen32(); +// +// double f; +// std::memcpy(&f, &x, sizeof(double)); +// +// test(f, y); +// } +// } + +// Raw output: +// +// 0.99998738156596089776684266325901262462139129638671875 ** 38583256: 0x140854811fb319e7 != 0x140854811fe4d778 (0x140854811fe4d778) (0x140854811fe4d778) +// -0.99843469603485224261874009243911132216453552246093750 ** 326215: 0x91dad4716de6fc4b != 0x91dad4716de5e587 (0x91dad4716de5e588) (0x91dad4716de5e587) +// 0.00003722856305626354357250426541092735988058848306537 ** -33: 0x5e47357c3582e49e != 0x5e47357c3582e4a3 (0x5e47357c3582e4a3) (0x5e47357c3582e4a3) +// -0.99996909838479330900895547529216855764389038085937500 ** 17078527: 0x9058409e5ea3b80a != 0x9058409e5eb11ef4 (0x9058409e5eb11ef4) (0x9058409e5eb11ef4) +// 0.99992690642006631929206150743993930518627166748046875 ** -6725291: 0x6c42a167a8b7c0b2 != 0x6c42a167a8b81d0e (0x6c42a167a8b81d0e) (0x6c42a167a8b81d0e) +// -0.99879181217764612110698863034485839307308197021484375 ** 485128: 0xb0d9c6f2f710d24 != 0xb0d9c6f2f71708d (0xb0d9c6f2f71708d) (0xb0d9c6f2f71708d) +// -1.00560838484317760510577954846667125821113586425781250 ** 92252: 0x6e744b727536056b != 0x6e744b72753599a4 (0x6e744b72753599a4) (0x6e744b72753599a4) +// 0.99999532655875444930870798998512327671051025390625000 ** 93511912: 0x1886c29a53ed9332 != 0x1886c29a53cba724 (0x1886c29a53cba724) (0x1886c29a53cba724) +// -0.99989751779212987514711130643263459205627441406250000 ** -2864087: 0xda664b586d48712f != 0xda664b586d437e8c (0xda664b586d437e8c) (0xda664b586d437e8c) +// -239.35307289280868303649185691028833389282226562500000000 ** -90: 0x137a8b43006c4438 != 0x137a8b43006c443e (0x137a8b43006c443e) (0x137a8b43006c443e) +// 0.96128212369452570307259975379565730690956115722656250 ** -9670: 0x625d7eb275191f6f != 0x625d7eb2751920bf (0x625d7eb2751920c0) (0x625d7eb2751920bf) +// 0.99996078564218904283222855156054720282554626464843750 ** 10583765: 0x1a829de67930f619 != 0x1a829de67951cc2d (0x1a829de67951cc2d) (0x1a829de67951cc2d) +// -953.14032530394126752071315422654151916503906250000000000 ** 22: 0x4d8a6d863703112c != 0x4d8a6d863703112e (0x4d8a6d863703112e) (0x4d8a6d863703112e) +// 0.99857985216514444370972114484175108373165130615234375 ** 335918: 0x14e345eb84f09d46 != 0x14e345eb84f036f4 (0x14e345eb84f036f4) (0x14e345eb84f036f4) +// -1.20521595553711002857255607523256912827491760253906250 ** -2760: 0x117b0064dd165101 != 0x117b0064dd16511a (0x117b0064dd16511a) (0x117b0064dd16511a) +// -1.19074911947068473594413262617308646440505981445312500 ** 3884: 0x7d132c80ed6973f6 != 0x7d132c80ed697072 (0x7d132c80ed697072) (0x7d132c80ed697072) +// -0.99999908129426284819629699995857663452625274658203125 ** -172780371: 0xce400f20e4a13b1a != 0xce400f20e3e56454 (0xce400f20e3e56454) (0xce400f20e3e56454) +// -0.00000000000000000000000000007930552628950037082519209 ** 8: 0x1142888ad3062fc1 != 0x1142888ad3062fbe (0x1142888ad3062fbe) (0x1142888ad3062fbe) +// -0.99998583604065760521706351937609724700450897216796875 ** -5861784: 0x476b83d92617a928 != 0x476b83d9261b0d4e (0x476b83d9261b0d4e) (0x476b83d9261b0d4e) +// 0.99989915564587761309667257592082023620605468750000000 ** 5468367: 0xe34d25f36eef64b != 0xe34d25f36f555ca (0xe34d25f36f555ca) (0xe34d25f36f555ca) +// 0.99977805581863743444870351595454849302768707275390625 ** -130493: 0x428ba17ba9286df6 != 0x428ba17ba9282f94 (0x428ba17ba9282f94) (0x428ba17ba9282f94) +// 29.19821057723854806909002945758402347564697265625000000 ** -20: 0x39d8ffec76e30251 != 0x39d8ffec76e3024a (0x39d8ffec76e3024a) (0x39d8ffec76e3024a) +// 0.99985373283040668290766461723251268267631530761718750 ** 2345687: 0x20ff8c2fd8e5b4e0 != 0x20ff8c2fd8e00564 (0x20ff8c2fd8e00564) (0x20ff8c2fd8e00564) +// -0.88383265987178571965188211834174580872058868408203125 ** -841: 0xc94c6878de27b17c != 0xc94c6878de27b20d (0xc94c6878de27b20d) (0xc94c6878de27b20d) +// 0.99999589815682188298495702838408760726451873779296875 ** 72449292: 0x25233af2e809c6a6 != 0x25233af2e87ddc61 (0x25233af2e87ddc61) (0x25233af2e87ddc61) +// 345736476.13618659973144531250000000000000000000000000000000000 ** -16: 0x2391db755176ac1b != 0x2391db755176ac19 (0x2391db755176ac19) (0x2391db755176ac19) +// -0.99999307321818442506611290809814818203449249267578125 ** -55045397: 0xe250f3d69f25ec86 != 0xe250f3d69f03e875 (0xe250f3d69f03e875) (0xe250f3d69f03e875) +// 1419676.56599932140670716762542724609375000000000000000000000 ** 25: 0x5fde72aa74287c2d != 0x5fde72aa74287c30 (0x5fde72aa74287c30) (0x5fde72aa74287c30) +// 0.95797249286536323431562323094112798571586608886718750 ** -11483: 0x6c63b79e88c07b6f != 0x6c63b79e88c07a3f (0x6c63b79e88c07a3f) (0x6c63b79e88c07a3f) +// 0.99998135132609855535434917328529991209506988525390625 ** 5682278: 0x3661650feb28b969 != 0x3661650feb22b7ed (0x3661650feb22b7ed) (0x3661650feb22b7ed) +// -1.02020595459010832151136582979233935475349426269531250 ** -1668: 0x3ced0e90ddfec9a3 != 0x3ced0e90ddfecabc (0x3ced0e90ddfecabc) (0x3ced0e90ddfecabc) +// 0.97281701550260646360612781791132874786853790283203125 ** 13717: 0x1dd88a88f24fc0d5 != 0x1dd88a88f24fb801 (0x1dd88a88f24fb801) (0x1dd88a88f24fb801) +// -0.88724290003841266294415390802896581590175628662109375 ** -3437: 0xe502ab8ea591420d != 0xe502ab8ea5914139 (0xe502ab8ea5914139) (0xe502ab8ea5914139) +// -0.99998630320599690701754980182158760726451873779296875 ** -11251995: 0xcdd44ff462cfbf32 != 0xcdd44ff462dbdcb2 (0xcdd44ff462dbdcb2) (0xcdd44ff462dbdcb2) +// -0.99995743703658013235013868325040675699710845947265625 ** 13995099: 0x8a38604324e009d5 != 0x8a38604324c2ec7d (0x8a38604324c2ec7d) (0x8a38604324c2ec7d) +// 0.99991090354494038816568490801728330552577972412109375 ** 7116340: 0x6c2ca56237c8161 != 0x6c2ca562366c00b (0x6c2ca562366c00b) (0x6c2ca562366c00b) +// 0.00000022955540324908999561342678487341206761129797087 ** 27: 0x1ab703277bbb112d != 0x1ab703277bbb1131 (0x1ab703277bbb1130) (0x1ab703277bbb1131) +// -1.00000041289256280663266807096078991889953613281250000 ** -365287834: 0x3255339a24caec8a != 0x3255339a26f00dc8 (0x3255339a26f00dc8) (0x3255339a26f00dc8) +// -1.38949508997780957209045027411775663495063781738281250 ** 1996: 0x7b22ad71344bbd0b != 0x7b22ad71344bbddf (0x7b22ad71344bbddf) (0x7b22ad71344bbddf) +// 0.99999867528282249118376512342365458607673645019531250 ** 164253172: 0x2c50f93fbc72a2b4 != 0x2c50f93fbb2d2fd4 (0x2c50f93fbb2d2fd4) (0x2c50f93fbb2d2fd4) +// 1.00356688770562074708436739456374198198318481445312500 ** -141698: 0x12717fb35c5fd169 != 0x12717fb35c5ff8c8 (0x12717fb35c5ff8c8) (0x12717fb35c5ff8c8) +// 368710687472107.18750000000000000000000000000000000000000000000000000 ** -20: 0x37282f0ae9be13c != 0x37282f0ae9be138 (0x37282f0ae9be138) (0x37282f0ae9be138) +// 0.99246668780181890312519499275367707014083862304687500 ** -44617: 0x5e5ad2c000333e50 != 0x5e5ad2c0003351f5 (0x5e5ad2c0003351f5) (0x5e5ad2c0003351f5) +// 1.13820783188362395499382273555966094136238098144531250 ** 1411: 0x506701df16f3a891 != 0x506701df16f3a70d (0x506701df16f3a70d) (0x506701df16f3a70d) +// -0.99671841783028414241130121808964759111404418945312500 ** 97041: 0xa32c44e6e77f8d3b != 0xa32c44e6e77f6a7a (0xa32c44e6e77f6a7a) (0xa32c44e6e77f6a7a) +// -0.57021831816264889614132016504299826920032501220703125 ** -802: 0x688ef2cc36fa60b3 != 0x688ef2cc36fa6064 (0x688ef2cc36fa6064) (0x688ef2cc36fa6064) +// -0.97423450510790443601649712945800274610519409179687500 ** 23570: 0x874c760d601ec94 != 0x874c760d601e66f (0x874c760d601e66f) (0x874c760d601e66f) +// -0.98067196425761504752216524138930253684520721435546875 ** -19882: 0x62ec606ceb9af0ae != 0x62ec606ceb9ae89c (0x62ec606ceb9ae89c) (0x62ec606ceb9ae89c) +// 0.99683039770073134100414335989626124501228332519531250 ** -29823: 0x487816b919332b03 != 0x487816b919333fe4 (0x487816b919333fe4) (0x487816b919333fe4) +// 0.99882797644578258378089685720624402165412902832031250 ** -540990: 0x792372efd5ca5ad2 != 0x792372efd5c92857 (0x792372efd5c92857) (0x792372efd5c92857) + +const testCases = [ + [0.99998738156596089776684266325901262462139129638671875 , 38583256], + [-0.99843469603485224261874009243911132216453552246093750 , 326215], + [0.00003722856305626354357250426541092735988058848306537 , -33], + [-0.99996909838479330900895547529216855764389038085937500 , 17078527], + [0.99992690642006631929206150743993930518627166748046875 , -6725291], + [-0.99879181217764612110698863034485839307308197021484375 , 485128], + [-1.00560838484317760510577954846667125821113586425781250 , 92252], + [0.99999532655875444930870798998512327671051025390625000 , 93511912], + [-0.99989751779212987514711130643263459205627441406250000 , -2864087], + [-239.35307289280868303649185691028833389282226562500000000 , -90], + [0.96128212369452570307259975379565730690956115722656250 , -9670], + [0.99996078564218904283222855156054720282554626464843750 , 10583765], + [-953.14032530394126752071315422654151916503906250000000000 , 22], + [0.99857985216514444370972114484175108373165130615234375 , 335918], + [-1.20521595553711002857255607523256912827491760253906250 , -2760], + [-1.19074911947068473594413262617308646440505981445312500 , 3884], + [-0.99999908129426284819629699995857663452625274658203125 , -172780371], + [-0.00000000000000000000000000007930552628950037082519209 , 8], + [-0.99998583604065760521706351937609724700450897216796875 , -5861784], + [0.99989915564587761309667257592082023620605468750000000 , 5468367], + [0.99977805581863743444870351595454849302768707275390625 , -130493], + [29.19821057723854806909002945758402347564697265625000000 , -20], + [0.99985373283040668290766461723251268267631530761718750 , 2345687], + [-0.88383265987178571965188211834174580872058868408203125 , -841], + [0.99999589815682188298495702838408760726451873779296875 , 72449292], + [345736476.13618659973144531250000000000000000000000000000000000 , -16], + [-0.99999307321818442506611290809814818203449249267578125 , -55045397], + [1419676.56599932140670716762542724609375000000000000000000000 , 25], + [0.95797249286536323431562323094112798571586608886718750 , -11483], + [0.99998135132609855535434917328529991209506988525390625 , 5682278], + [-1.02020595459010832151136582979233935475349426269531250 , -1668], + [0.97281701550260646360612781791132874786853790283203125 , 13717], + [-0.88724290003841266294415390802896581590175628662109375 , -3437], + [-0.99998630320599690701754980182158760726451873779296875 , -11251995], + [-0.99995743703658013235013868325040675699710845947265625 , 13995099], + [0.99991090354494038816568490801728330552577972412109375 , 7116340], + [0.00000022955540324908999561342678487341206761129797087 , 27], + [-1.00000041289256280663266807096078991889953613281250000 , -365287834], + [-1.38949508997780957209045027411775663495063781738281250 , 1996], + [0.99999867528282249118376512342365458607673645019531250 , 164253172], + [1.00356688770562074708436739456374198198318481445312500 , -141698], + [368710687472107.18750000000000000000000000000000000000000000000000000 , -20], + [0.99246668780181890312519499275367707014083862304687500 , -44617], + [1.13820783188362395499382273555966094136238098144531250 , 1411], + [-0.99671841783028414241130121808964759111404418945312500 , 97041], + [-0.57021831816264889614132016504299826920032501220703125 , -802], + [-0.97423450510790443601649712945800274610519409179687500 , 23570], + [-0.98067196425761504752216524138930253684520721435546875 , -19882], + [0.99683039770073134100414335989626124501228332519531250 , -29823], + [0.99882797644578258378089685720624402165412902832031250 , -540990], +]; + +// Test program modified to avoid bases with |abs(x) < 1| and large exponents. +// +// ```cpp +// // Skip over likely denormals. +// if (-1 < f && f < 0) { +// f -= 1; +// } else if (0 < f && f < 1) { +// f += 1; +// } +// +// // Keep the power small. +// y &= 63; +// ``` +// +// 7.86990183266223297664510027971118688583374023437500000 ** 54: 0x49fa67548289784a != 0x49fa675482897851 (0x49fa675482897850) (0x49fa675482897851) +// -1.00000018751738117828153917798772454261779785156250000 ** 25: 0xbff00004ea6921f6 != 0xbff00004ea6921fc (0xbff00004ea6921fc) (0xbff00004ea6921fc) +// 1.19908234423429393977755808009533211588859558105468750 ** 58: 0x40e246fe7b30c6ec != 0x40e246fe7b30c6e6 (0x40e246fe7b30c6e6) (0x40e246fe7b30c6e6) +// 1.00000649317438283780745678086532279849052429199218750 ** 42: 0x3ff0011dffabb95c != 0x3ff0011dffabb950 (0x3ff0011dffabb950) (0x3ff0011dffabb950) +// 863370098.16819441318511962890625000000000000000000000000000000 ** 27: 0x7206b860614eb6df != 0x7206b860614eb6d9 (0x7206b860614eb6d9) (0x7206b860614eb6d9) +// -1.00011928123711690830077714053913950920104980468750000 ** 57: 0xbff01bf129d0ffab != 0xbff01bf129d0ffbf (0xbff01bf129d0ffbf) (0xbff01bf129d0ffbf) +// -1.14006037237328494704513559554470703005790710449218750 ** 30: 0x404983fd4d57c4aa != 0x404983fd4d57c4a0 (0x404983fd4d57c4a0) (0x404983fd4d57c4a0) +// -447.11057737163486081044538877904415130615234375000000000 ** 8: 0x4455a4e4be220fce != 0x4455a4e4be220fd0 (0x4455a4e4be220fd0) (0x4455a4e4be220fd0) +// -1.03656507831253685836259137431625276803970336914062500 ** 20: 0x4000681e0886d6db != 0x4000681e0886d6d9 (0x4000681e0886d6d9) (0x4000681e0886d6d9) +// -1.00000465330344945336094042431795969605445861816406250 ** 41: 0xbff000c81257efc1 != 0xbff000c81257efc6 (0xbff000c81257efc6) (0xbff000c81257efc6) +// -1.00002726631492944164847358479164540767669677734375000 ** 14: 0x3ff00190579a2f93 != 0x3ff00190579a2f90 (0x3ff00190579a2f90) (0x3ff00190579a2f90) +// 2512068.57641875604167580604553222656250000000000000000000000 ** 26: 0x627b50512391a46e != 0x627b50512391a46c (0x627b50512391a46c) (0x627b50512391a46c) +// 3309586784.85019683837890625000000000000000000000000000000000000 ** 30: 0x7b3a5b69a3a40717 != 0x7b3a5b69a3a40719 (0x7b3a5b69a3a40719) (0x7b3a5b69a3a40719) +// 1.40742719307547781149025922786677256226539611816406250 ** 19: 0x4084a6ad66b5f1ce != 0x4084a6ad66b5f1d1 (0x4084a6ad66b5f1d0) (0x4084a6ad66b5f1d1) +// 1.00035740860596344958821646287105977535247802734375000 ** 36: 0x3ff0350873b3189e != 0x3ff0350873b318a0 (0x3ff0350873b318a0) (0x3ff0350873b318a0) +testCases.push( + [7.86990183266223297664510027971118688583374023437500000 , 54], + [-1.00000018751738117828153917798772454261779785156250000 , 25], + [1.19908234423429393977755808009533211588859558105468750 , 58], + [1.00000649317438283780745678086532279849052429199218750 , 42], + [863370098.16819441318511962890625000000000000000000000000000000 , 27], + [-1.00011928123711690830077714053913950920104980468750000 , 57], + [-1.14006037237328494704513559554470703005790710449218750 , 30], + [-447.11057737163486081044538877904415130615234375000000000 , 8], + [-1.03656507831253685836259137431625276803970336914062500 , 20], + [-1.00000465330344945336094042431795969605445861816406250 , 41], + [-1.00002726631492944164847358479164540767669677734375000 , 14], + [2512068.57641875604167580604553222656250000000000000000000000 , 26], + [3309586784.85019683837890625000000000000000000000000000000000000 , 30], + [1.40742719307547781149025922786677256226539611816406250 , 19], + [1.00035740860596344958821646287105977535247802734375000 , 36], +); + +// Test program modified to only use small integer bases (< 20) and positive exponents. +// +// ```cpp +// f = static_cast<double>(x); +// f = std::fmod(f, 20); +// y &= 63; +// ``` +// +// 13.00000000000000000000000000000000000000000000000000000 ** 31: 0x471a3d23b248d522 != 0x471a3d23b248d520 (0x471a3d23b248d520) (0x471a3d23b248d520) +// 13.00000000000000000000000000000000000000000000000000000 ** 41: 0x496a51a4d0054bb2 != 0x496a51a4d0054bb1 (0x496a51a4d0054bb0) (0x496a51a4d0054bb1) +// 13.00000000000000000000000000000000000000000000000000000 ** 51: 0x4bba6635f3af40fa != 0x4bba6635f3af40f8 (0x4bba6635f3af40f8) (0x4bba6635f3af40f8) +// 13.00000000000000000000000000000000000000000000000000000 ** 58: 0x4d58af19e7576d60 != 0x4d58af19e7576d5e (0x4d58af19e7576d5e) (0x4d58af19e7576d5e) +// 13.00000000000000000000000000000000000000000000000000000 ** 63: 0x4e817b180a97c789 != 0x4e817b180a97c787 (0x4e817b180a97c787) (0x4e817b180a97c787) +// 11.00000000000000000000000000000000000000000000000000000 ** 63: 0x4d8ec9288a0088ce != 0x4d8ec9288a0088d0 (0x4d8ec9288a0088d0) (0x4d8ec9288a0088d0) +// 13.00000000000000000000000000000000000000000000000000000 ** 47: 0x4ace49afd4c20163 != 0x4ace49afd4c20161 (0x4ace49afd4c20161) (0x4ace49afd4c20161) +// 13.00000000000000000000000000000000000000000000000000000 ** 41: 0x496a51a4d0054bb2 != 0x496a51a4d0054bb1 (0x496a51a4d0054bb0) (0x496a51a4d0054bb1) +// 13.00000000000000000000000000000000000000000000000000000 ** 63: 0x4e817b180a97c789 != 0x4e817b180a97c787 (0x4e817b180a97c787) (0x4e817b180a97c787) +// 13.00000000000000000000000000000000000000000000000000000 ** 31: 0x471a3d23b248d522 != 0x471a3d23b248d520 (0x471a3d23b248d520) (0x471a3d23b248d520) +// 13.00000000000000000000000000000000000000000000000000000 ** 49: 0x4b43fea5137412eb != 0x4b43fea5137412e9 (0x4b43fea5137412e9) (0x4b43fea5137412e9) +// 13.00000000000000000000000000000000000000000000000000000 ** 58: 0x4d58af19e7576d60 != 0x4d58af19e7576d5e (0x4d58af19e7576d5e) (0x4d58af19e7576d5e) +// 13.00000000000000000000000000000000000000000000000000000 ** 31: 0x471a3d23b248d522 != 0x471a3d23b248d520 (0x471a3d23b248d520) (0x471a3d23b248d520) +// 11.00000000000000000000000000000000000000000000000000000 ** 63: 0x4d8ec9288a0088ce != 0x4d8ec9288a0088d0 (0x4d8ec9288a0088d0) (0x4d8ec9288a0088d0) +// 13.00000000000000000000000000000000000000000000000000000 ** 31: 0x471a3d23b248d522 != 0x471a3d23b248d520 (0x471a3d23b248d520) (0x471a3d23b248d520) +testCases.push( + [13.00000000000000000000000000000000000000000000000000000 , 31], + [13.00000000000000000000000000000000000000000000000000000 , 41], + [13.00000000000000000000000000000000000000000000000000000 , 51], + [13.00000000000000000000000000000000000000000000000000000 , 58], + [13.00000000000000000000000000000000000000000000000000000 , 63], + [11.00000000000000000000000000000000000000000000000000000 , 63], + [13.00000000000000000000000000000000000000000000000000000 , 47], + [13.00000000000000000000000000000000000000000000000000000 , 41], + [13.00000000000000000000000000000000000000000000000000000 , 63], + [13.00000000000000000000000000000000000000000000000000000 , 31], + [13.00000000000000000000000000000000000000000000000000000 , 49], + [13.00000000000000000000000000000000000000000000000000000 , 58], + [13.00000000000000000000000000000000000000000000000000000 , 31], + [11.00000000000000000000000000000000000000000000000000000 , 63], + [13.00000000000000000000000000000000000000000000000000000 , 31], +); + +// Test program modified to only use small integer bases (< 20) and negative exponents. +// +// ```cpp +// f = static_cast<double>(x); +// f = std::fmod(f, 20); +// y &= 63; +// y = -y; +// ``` +// +// 14.00000000000000000000000000000000000000000000000000000 ** -57: 0x325f938745f05e58 != 0x325f938745f05e5a (0x325f938745f05e5a) (0x325f938745f05e5a) +// 11.00000000000000000000000000000000000000000000000000000 ** -53: 0x34791bddc7b3025a != 0x34791bddc7b30259 (0x34791bddc7b30258) (0x34791bddc7b30259) +// 7.00000000000000000000000000000000000000000000000000000 ** -57: 0x35ef938745f05e58 != 0x35ef938745f05e5a (0x35ef938745f05e5a) (0x35ef938745f05e5a) +// 15.00000000000000000000000000000000000000000000000000000 ** -50: 0x33b933babb6d9cd8 != 0x33b933babb6d9cda (0x33b933babb6d9cda) (0x33b933babb6d9cda) +// 14.00000000000000000000000000000000000000000000000000000 ** -57: 0x325f938745f05e58 != 0x325f938745f05e5a (0x325f938745f05e5a) (0x325f938745f05e5a) +// 13.00000000000000000000000000000000000000000000000000000 ** -33: 0x384d8ee9f0edfd7c != 0x384d8ee9f0edfd7d (0x384d8ee9f0edfd7e) (0x384d8ee9f0edfd7d) +// 19.00000000000000000000000000000000000000000000000000000 ** -53: 0x31dd0994e8aaf4e0 != 0x31dd0994e8aaf4e1 (0x31dd0994e8aaf4e2) (0x31dd0994e8aaf4e1) +// 15.00000000000000000000000000000000000000000000000000000 ** -50: 0x33b933babb6d9cd8 != 0x33b933babb6d9cda (0x33b933babb6d9cda) (0x33b933babb6d9cda) +// 14.00000000000000000000000000000000000000000000000000000 ** -57: 0x325f938745f05e58 != 0x325f938745f05e5a (0x325f938745f05e5a) (0x325f938745f05e5a) +// 13.00000000000000000000000000000000000000000000000000000 ** -63: 0x315d4a0a2c8d4bd8 != 0x315d4a0a2c8d4bdb (0x315d4a0a2c8d4bdb) (0x315d4a0a2c8d4bdb) +// 11.00000000000000000000000000000000000000000000000000000 ** -53: 0x34791bddc7b3025a != 0x34791bddc7b30259 (0x34791bddc7b30258) (0x34791bddc7b30259) +// 15.00000000000000000000000000000000000000000000000000000 ** -50: 0x33b933babb6d9cd8 != 0x33b933babb6d9cda (0x33b933babb6d9cda) (0x33b933babb6d9cda) +// 13.00000000000000000000000000000000000000000000000000000 ** -53: 0x33ad60ed868e2926 != 0x33ad60ed868e2928 (0x33ad60ed868e2928) (0x33ad60ed868e2928) +// 19.00000000000000000000000000000000000000000000000000000 ** -53: 0x31dd0994e8aaf4e0 != 0x31dd0994e8aaf4e1 (0x31dd0994e8aaf4e2) (0x31dd0994e8aaf4e1) +// 13.00000000000000000000000000000000000000000000000000000 ** -33: 0x384d8ee9f0edfd7c != 0x384d8ee9f0edfd7d (0x384d8ee9f0edfd7e) (0x384d8ee9f0edfd7d) +testCases.push( + [14.00000000000000000000000000000000000000000000000000000 , -57], + [11.00000000000000000000000000000000000000000000000000000 , -53], + [7.00000000000000000000000000000000000000000000000000000 , -57], + [15.00000000000000000000000000000000000000000000000000000 , -50], + [14.00000000000000000000000000000000000000000000000000000 , -57], + [13.00000000000000000000000000000000000000000000000000000 , -33], + [19.00000000000000000000000000000000000000000000000000000 , -53], + [15.00000000000000000000000000000000000000000000000000000 , -50], + [14.00000000000000000000000000000000000000000000000000000 , -57], + [13.00000000000000000000000000000000000000000000000000000 , -63], + [11.00000000000000000000000000000000000000000000000000000 , -53], + [15.00000000000000000000000000000000000000000000000000000 , -50], + [13.00000000000000000000000000000000000000000000000000000 , -53], + [19.00000000000000000000000000000000000000000000000000000 , -53], + [13.00000000000000000000000000000000000000000000000000000 , -33], +); + +// std::pow is less precise on Windows. +const maxError = getBuildConfiguration("windows") ? 3 : 1; + +// Ensure the error is less-or-equal to |maxError| ULP when compared to fdlibm. +for (let [x, y] of testCases) { + let actual = Math.pow(x, y); + let expected = fdlibm.pow(x, y); + let error = errorInULP(actual, expected); + + assertEq(error <= maxError, true, + `${x} ** ${y}: ${actual} (${toBits(actual).toString(16)}) != ${expected} (${toBits(expected).toString(16)})`); +} + +// Test program modified to use 4 as the exponent: +// +// ```cpp +// y = 4; +// ``` +// +// -0.00000000000000000000000000000749666789562697097993956 ** 4: 0x27bfdbe3cf0b7e1d != 0x27bfdbe3cf0b7e1b (0x27bfdbe3cf0b7e1b) (0x27bfdbe3cf0b7e1b) +// 0.00000000000000000000000000000000000000000000000000000 ** 4: 0xd3e1e77bd0d8f5d != 0xd3e1e77bd0d8f5f (0xd3e1e77bd0d8f5f) (0xd3e1e77bd0d8f5f) +// -0.00000000000000000000000000023705601542216470968966009 ** 4: 0x28fe60d2f5131d02 != 0x28fe60d2f5131d04 (0x28fe60d2f5131d04) (0x28fe60d2f5131d04) +// 0.00000000000000000000000000000000000000000000000000441 ** 4: 0x161dad0fa681c66c != 0x161dad0fa681c66b (0x161dad0fa681c66a) (0x161dad0fa681c66b) +// 0.00000000000000537255761599995092558925668894011631095 ** 4: 0x3414eb4baea214b6 != 0x3414eb4baea214b5 (0x3414eb4baea214b4) (0x3414eb4baea214b5) +// 0.01225688384384779339164595057809492573142051696777344 ** 4: 0x3e583bd550871dfc != 0x3e583bd550871dfd (0x3e583bd550871dfe) (0x3e583bd550871dfd) +// -0.00000000000000000000000000000000000000000000000000000 ** 4: 0xa59292360f6d326 != 0xa59292360f6d324 (0xa59292360f6d324) (0xa59292360f6d324) +// -0.00000000000000000000000000000000000000000000000000000 ** 4: 0x109fb7a8459811ec != 0x109fb7a8459811ed (0x109fb7a8459811ee) (0x109fb7a8459811ed) +// -120834175976112453093144522854609799898808186321228136949237230085114691584.00000000000000000000000000000000000000000000000000000 ** 4: 0x7d74dcc37a2d7dc2 != 0x7d74dcc37a2d7dc3 (0x7d74dcc37a2d7dc4) (0x7d74dcc37a2d7dc3) +// -6676.83140968165753292851150035858154296875000000000000000 ** 4: 0x431c3e0ef48fe66a != 0x431c3e0ef48fe66c (0x431c3e0ef48fe66c) (0x431c3e0ef48fe66c) +// -0.00000000000000000000000000000000000000000000039753861 ** 4: 0x1a3a87f39f288766 != 0x1a3a87f39f288764 (0x1a3a87f39f288764) (0x1a3a87f39f288764) +// 129749516186492032220917661696.00000000000000000000000000000000000000000000000000000 ** 4: 0x581cc58a512bdd10 != 0x581cc58a512bdd12 (0x581cc58a512bdd12) (0x581cc58a512bdd12) +// -1888635225450734959219733085647207705818299180319259746124169216.00000000000000000000000000000000000000000000000000000 ** 4: 0x747bc423aba49de6 != 0x747bc423aba49de5 (0x747bc423aba49de4) (0x747bc423aba49de5) +// 7934926680560039158281691725824.00000000000000000000000000000000000000000000000000000 ** 4: 0x5997fceb5eed5c94 != 0x5997fceb5eed5c93 (0x5997fceb5eed5c92) (0x5997fceb5eed5c93) +// -0.00000000000000579868166379701264244398310517312073637 ** 4: 0x341c635a1a764ef2 != 0x341c635a1a764ef0 (0x341c635a1a764ef0) (0x341c635a1a764ef0) +// +// +// Test program modified to avoid bases with |abs(x) < 1| and large exponents. +// +// ```cpp +// // Skip over likely denormals. +// if (-1 < f && f < 0) { +// f -= 1; +// } else if (0 < f && f < 1) { +// f += 1; +// } +// +// f = std::fmod(f, 20); +// +// y = 4; +// ``` +// +// 4.73347349464893341064453125000000000000000000000000000 ** 4: 0x407f604c239c2323 != 0x407f604c239c2321 (0x407f604c239c2321) (0x407f604c239c2321) +// -12.35635152040049433708190917968750000000000000000000000 ** 4: 0x40d6c3c0652f0948 != 0x40d6c3c0652f0949 (0x40d6c3c0652f094a) (0x40d6c3c0652f0949) +// -1.50385549572482823954544528533006086945533752441406250 ** 4: 0x40147581145bc6e6 != 0x40147581145bc6e7 (0x40147581145bc6e8) (0x40147581145bc6e7) +// -8.93048901623114943504333496093750000000000000000000000 ** 4: 0x40b8d8a463c28bd6 != 0x40b8d8a463c28bd7 (0x40b8d8a463c28bd8) (0x40b8d8a463c28bd7) +// 19.02711385915608843788504600524902343750000000000000000 ** 4: 0x40ffffa7d5df2562 != 0x40ffffa7d5df2560 (0x40ffffa7d5df2560) (0x40ffffa7d5df2560) +// 17.83878016096969076897948980331420898437500000000000000 ** 4: 0x40f8b914a6acb498 != 0x40f8b914a6acb497 (0x40f8b914a6acb496) (0x40f8b914a6acb497) +// 12.90541613101959228515625000000000000000000000000000000 ** 4: 0x40db16b4c2dafa0a != 0x40db16b4c2dafa0c (0x40db16b4c2dafa0c) (0x40db16b4c2dafa0c) +// -18.34655402903445065021514892578125000000000000000000000 ** 4: 0x40fba90e5b7bbc6a != 0x40fba90e5b7bbc6b (0x40fba90e5b7bbc6c) (0x40fba90e5b7bbc6b) +// -13.28634420270100235939025878906250000000000000000000000 ** 4: 0x40de6e70b9ed821a != 0x40de6e70b9ed821c (0x40de6e70b9ed821c) (0x40de6e70b9ed821c) +// 18.52965961024165153503417968750000000000000000000000000 ** 4: 0x40fcc800b850b01a != 0x40fcc800b850b018 (0x40fcc800b850b018) (0x40fcc800b850b018) +// 13.32226210648514097556471824645996093750000000000000000 ** 4: 0x40dec3063a559350 != 0x40dec3063a55934e (0x40dec3063a55934e) (0x40dec3063a55934e) +// 1.09174693829848346027233674249146133661270141601562500 ** 4: 0x3ff6bafe5bbe7532 != 0x3ff6bafe5bbe7533 (0x3ff6bafe5bbe7534) (0x3ff6bafe5bbe7533) +// 9.35059530444141273619607090950012207031250000000000000 ** 4: 0x40bddca3dd9f5c8f != 0x40bddca3dd9f5c91 (0x40bddca3dd9f5c91) (0x40bddca3dd9f5c91) +// 17.59552449546754360198974609375000000000000000000000000 ** 4: 0x40f766db2706f434 != 0x40f766db2706f435 (0x40f766db2706f436) (0x40f766db2706f435) +// 17.94561576098203659057617187500000000000000000000000000 ** 4: 0x40f952110041965c != 0x40f952110041965a (0x40f952110041965a) (0x40f952110041965a) +const testCases4 = [ + [-0.00000000000000000000000000000749666789562697097993956 , 4], + [0.00000000000000000000000000000000000000000000000000000 , 4], + [-0.00000000000000000000000000023705601542216470968966009 , 4], + [0.00000000000000000000000000000000000000000000000000441 , 4], + [0.00000000000000537255761599995092558925668894011631095 , 4], + [0.01225688384384779339164595057809492573142051696777344 , 4], + [-0.00000000000000000000000000000000000000000000000000000 , 4], + [-0.00000000000000000000000000000000000000000000000000000 , 4], + [-120834175976112453093144522854609799898808186321228136949237230085114691584.00000000000000000000000000000000000000000000000000000 , 4], + [-6676.83140968165753292851150035858154296875000000000000000 , 4], + [-0.00000000000000000000000000000000000000000000039753861 , 4], + [129749516186492032220917661696.00000000000000000000000000000000000000000000000000000 , 4], + [-1888635225450734959219733085647207705818299180319259746124169216.00000000000000000000000000000000000000000000000000000 , 4], + [7934926680560039158281691725824.00000000000000000000000000000000000000000000000000000 , 4], + [-0.00000000000000579868166379701264244398310517312073637 , 4], + + [4.73347349464893341064453125000000000000000000000000000 , 4], + [-12.35635152040049433708190917968750000000000000000000000 , 4], + [-1.50385549572482823954544528533006086945533752441406250 , 4], + [-8.93048901623114943504333496093750000000000000000000000 , 4], + [19.02711385915608843788504600524902343750000000000000000 , 4], + [17.83878016096969076897948980331420898437500000000000000 , 4], + [12.90541613101959228515625000000000000000000000000000000 , 4], + [-18.34655402903445065021514892578125000000000000000000000 , 4], + [-13.28634420270100235939025878906250000000000000000000000 , 4], + [18.52965961024165153503417968750000000000000000000000000 , 4], + [13.32226210648514097556471824645996093750000000000000000 , 4], + [1.09174693829848346027233674249146133661270141601562500 , 4], + [9.35059530444141273619607090950012207031250000000000000 , 4], + [17.59552449546754360198974609375000000000000000000000000 , 4], + [17.94561576098203659057617187500000000000000000000000000 , 4], +]; + +// Ensure the error is less-or-equal to 2 ULP when compared to fdlibm. +// +// This can produce a larger error than std::pow, because we evaluate +// |x ** 4| as |(x * x) * (x * x)| to match Ion. +for (let [x, y] of testCases4) { + let actual = Math.pow(x, y); + let expected = fdlibm.pow(x, y); + let error = errorInULP(actual, expected); + assertEq(error <= 2, true, + `${x} ** ${y}: ${actual} (${toBits(actual).toString(16)}) != ${expected} (${toBits(expected).toString(16)})`); +} + +for (let [x, y] of testCases4) { + // Replace |y| with a constant to trigger Ion optimisations. + let actual = Math.pow(x, 4); + let expected = fdlibm.pow(x, y); + let error = errorInULP(actual, expected); + assertEq(error <= 2, true, + `${x} ** ${y}: ${actual} (${toBits(actual).toString(16)}) != ${expected} (${toBits(expected).toString(16)})`); +} + +// Test program modified to use 3 as the exponent: +// +// ```cpp +// y = 3; +// ``` +// +// 196194373276.42089843750000000000000000000000000000000000000000000 ** 3: 0x46f745720bc58e22 != 0x46f745720bc58e23 (0x46f745720bc58e24) (0x46f745720bc58e23) +// 17260025115986696435331651385474892363490876322742272.00000000000000000000000000000000000000000000000000000 ** 3: 0x6077f8040eb542fc != 0x6077f8040eb542fb (0x6077f8040eb542fa) (0x6077f8040eb542fb) +// -0.00000000000000000000000000000000000000000000000000000 ** 3: 0x9307c17ddf2c4af6 != 0x9307c17ddf2c4af7 (0x9307c17ddf2c4af8) (0x9307c17ddf2c4af7) +// 2359506498398344427475761591701240715936602989985583832867274752.00000000000000000000000000000000000000000000000000000 ** 3: 0x6767960b1076dc24 != 0x6767960b1076dc25 (0x6767960b1076dc26) (0x6767960b1076dc25) +// 22724457948673043906745552566513068013978508710758109286797554897659283949989408425377792.00000000000000000000000000000000000000000000000000000 ** 3: 0x76f74ab82115b372 != 0x76f74ab82115b373 (0x76f74ab82115b374) (0x76f74ab82115b373) +// -1024872849611580448634200763411882795753013248.00000000000000000000000000000000000000000000000000000 ** 3: 0xdbf7b2694dce1d6c != 0xdbf7b2694dce1d6b (0xdbf7b2694dce1d6a) (0xdbf7b2694dce1d6b) +// -918435268181356203923125447950336.00000000000000000000000000000000000000000000000000000 ** 3: 0xd476ab3173dbfcc0 != 0xd476ab3173dbfcbf (0xd476ab3173dbfcbe) (0xd476ab3173dbfcbf) +// 558545783776545344834655968246618719333738303286453207040.00000000000000000000000000000000000000000000000000000 ** 3: 0x634716045b3ee61c != 0x634716045b3ee61b (0x634716045b3ee61a) (0x634716045b3ee61b) +// 0.00000000000000000000000000000000000000000000000000000 ** 3: 0x1c6f3bddc90315c != 0x1c6f3bddc90315b (0x1c6f3bddc90315a) (0x1c6f3bddc90315b) +// -0.00000000000261062225071774409619236799548496917242058 ** 3: 0xb8b7a667f8b6344e != 0xb8b7a667f8b6344f (0xb8b7a667f8b63450) (0xb8b7a667f8b6344f) +// 0.00000000000000000000000000000000000000000000012475377 ** 3: 0x23571f25316bb01e != 0x23571f25316bb01f (0x23571f25316bb020) (0x23571f25316bb01f) +// -0.00000000000000000000000000000000000000000000000000000 ** 3: 0x93f6c04c12acc76c != 0x93f6c04c12acc76d (0x93f6c04c12acc76e) (0x93f6c04c12acc76d) +// 0.00000000000000000000000000000000000000000000000000000 ** 3: 0x676eb3aa0a63236 != 0x676eb3aa0a63237 (0x676eb3aa0a63238) (0x676eb3aa0a63237) +// 0.00000000000000000000000007454937961610833261396029146 ** 3: 0x3047fcbe59481112 != 0x3047fcbe59481111 (0x3047fcbe59481110) (0x3047fcbe59481111) +// 0.00000000000000000000000000000000000003326770580987513 ** 3: 0x2896aaec8bb845c8 != 0x2896aaec8bb845c9 (0x2896aaec8bb845ca) (0x2896aaec8bb845c9) +// +// +// Test program modified to avoid bases with |abs(x) < 1| and large exponents. +// +// ```cpp +// // Skip over likely denormals. +// if (-1 < f && f < 0) { +// f -= 1; +// } else if (0 < f && f < 1) { +// f += 1; +// } +// +// f = std::fmod(f, 20); +// +// y = 3; +// ``` +// +// -11.40858423709869384765625000000000000000000000000000000 ** 3: 0xc0973392c88cadcc != 0xc0973392c88cadcd (0xc0973392c88cadce) (0xc0973392c88cadcd) +// 11.42477834224700927734375000000000000000000000000000000 ** 3: 0x40974ce701d58518 != 0x40974ce701d58519 (0x40974ce701d5851a) (0x40974ce701d58519) +// -11.46123231985238533070514677092432975769042968750000000 ** 3: 0xc097862ed0211e58 != 0xc097862ed0211e59 (0xc097862ed0211e5a) (0xc097862ed0211e59) +// -11.40183842182159423828125000000000000000000000000000000 ** 3: 0xc097290b23fe8cdc != 0xc097290b23fe8cdd (0xc097290b23fe8cde) (0xc097290b23fe8cdd) +// 2.87109172078278795936512324260547757148742675781250000 ** 3: 0x4037aab95517cdd0 != 0x4037aab95517cdcf (0x4037aab95517cdce) (0x4037aab95517cdcf) +// -0.72109144181013107299804687500000000000000000000000000 ** 3: 0xbfd7ff25d4fd46bc != 0xbfd7ff25d4fd46bd (0xbfd7ff25d4fd46be) (0xbfd7ff25d4fd46bd) +// 5.70116788148880004882812500000000000000000000000000000 ** 3: 0x406729d1c53687b4 != 0x406729d1c53687b5 (0x406729d1c53687b6) (0x406729d1c53687b5) +// -11.32285048566092200417187996208667755126953125000000000 ** 3: 0xc096aeac14d25c0e != 0xc096aeac14d25c0f (0xc096aeac14d25c10) (0xc096aeac14d25c0f) +// 1.41961999237537384033203125000000000000000000000000000 ** 3: 0x4006e34ea8957732 != 0x4006e34ea8957733 (0x4006e34ea8957734) (0x4006e34ea8957733) +// -11.52091628707762538397219032049179077148437500000000000 ** 3: 0xc097e4c12ab5e96e != 0xc097e4c12ab5e96f (0xc097e4c12ab5e970) (0xc097e4c12ab5e96f) +// -5.73415940999984741210937500000000000000000000000000000 ** 3: 0xc067915c3febbeba != 0xc067915c3febbebb (0xc067915c3febbebc) (0xc067915c3febbebb) +// 1.41478560105390638312883311300538480281829833984375000 ** 3: 0x4006a7a69b402738 != 0x4006a7a69b402737 (0x4006a7a69b402736) (0x4006a7a69b402737) +// -2.88328036665916442871093750000000000000000000000000000 ** 3: 0xc037f8371e1d17ce != 0xc037f8371e1d17cf (0xc037f8371e1d17d0) (0xc037f8371e1d17cf) +// 1.42408178602072932328326260176254436373710632324218750 ** 3: 0x40071aba43b3bcea != 0x40071aba43b3bceb (0x40071aba43b3bcec) (0x40071aba43b3bceb) +// 11.48128501093015074729919433593750000000000000000000000 ** 3: 0x4097a5d8fdac3954 != 0x4097a5d8fdac3955 (0x4097a5d8fdac3956) (0x4097a5d8fdac3955) +const testCases3 = [ + [196194373276.42089843750000000000000000000000000000000000000000000 , 3], + [17260025115986696435331651385474892363490876322742272.00000000000000000000000000000000000000000000000000000 , 3], + [-0.00000000000000000000000000000000000000000000000000000 , 3], + [2359506498398344427475761591701240715936602989985583832867274752.00000000000000000000000000000000000000000000000000000 , 3], + [22724457948673043906745552566513068013978508710758109286797554897659283949989408425377792.00000000000000000000000000000000000000000000000000000 , 3], + [-1024872849611580448634200763411882795753013248.00000000000000000000000000000000000000000000000000000 , 3], + [-918435268181356203923125447950336.00000000000000000000000000000000000000000000000000000 , 3], + [558545783776545344834655968246618719333738303286453207040.00000000000000000000000000000000000000000000000000000 , 3], + [0.00000000000000000000000000000000000000000000000000000 , 3], + [-0.00000000000261062225071774409619236799548496917242058 , 3], + [0.00000000000000000000000000000000000000000000012475377 , 3], + [-0.00000000000000000000000000000000000000000000000000000 , 3], + [0.00000000000000000000000000000000000000000000000000000 , 3], + [0.00000000000000000000000007454937961610833261396029146 , 3], + [0.00000000000000000000000000000000000003326770580987513 , 3], + + [-11.40858423709869384765625000000000000000000000000000000 , 3], + [11.42477834224700927734375000000000000000000000000000000 , 3], + [-11.46123231985238533070514677092432975769042968750000000 , 3], + [-11.40183842182159423828125000000000000000000000000000000 , 3], + [2.87109172078278795936512324260547757148742675781250000 , 3], + [-0.72109144181013107299804687500000000000000000000000000 , 3], + [5.70116788148880004882812500000000000000000000000000000 , 3], + [-11.32285048566092200417187996208667755126953125000000000 , 3], + [1.41961999237537384033203125000000000000000000000000000 , 3], + [-11.52091628707762538397219032049179077148437500000000000 , 3], + [-5.73415940999984741210937500000000000000000000000000000 , 3], + [1.41478560105390638312883311300538480281829833984375000 , 3], + [-2.88328036665916442871093750000000000000000000000000000 , 3], + [1.42408178602072932328326260176254436373710632324218750 , 3], + [11.48128501093015074729919433593750000000000000000000000 , 3], +]; + +// Ensure the error is less-or-equal to 2 ULP when compared to fdlibm. +// +// This can produce a larger error than std::pow, because we evaluate +// |x ** 3| as |(x * x) * x| to match Ion. +for (let [x, y] of testCases3) { + let actual = Math.pow(x, y); + let expected = fdlibm.pow(x, y); + let error = errorInULP(actual, expected); + assertEq(error <= 2, true, + `${x} ** ${y}: ${actual} (${toBits(actual).toString(16)}) != ${expected} (${toBits(expected).toString(16)})`); +} + +for (let [x, y] of testCases3) { + // Replace |y| with a constant to trigger Ion optimisations. + let actual = Math.pow(x, 3); + let expected = fdlibm.pow(x, y); + let error = errorInULP(actual, expected); + assertEq(error <= 2, true, + `${x} ** ${y}: ${actual} (${toBits(actual).toString(16)}) != ${expected} (${toBits(expected).toString(16)})`); +} + +if (typeof reportCompare === "function") + reportCompare(true, true); |