diff options
Diffstat (limited to 'Lib/test/decimaltestdata/power.decTest')
| -rw-r--r-- | Lib/test/decimaltestdata/power.decTest | 651 |
1 files changed, 651 insertions, 0 deletions
diff --git a/Lib/test/decimaltestdata/power.decTest b/Lib/test/decimaltestdata/power.decTest new file mode 100644 index 0000000..d7357e8 --- /dev/null +++ b/Lib/test/decimaltestdata/power.decTest @@ -0,0 +1,651 @@ +---------------------------------------------------------------------- +-- power.decTest -- decimal exponentiation -- +-- Copyright (c) IBM Corporation, 1981, 2003. All rights reserved. -- +------------------------------------------------------------------------ +-- Please see the document "General Decimal Arithmetic Testcases" -- +-- at http://www2.hursley.ibm.com/decimal for the description of -- +-- these testcases. -- +-- -- +-- These testcases are experimental ('beta' versions), and they -- +-- may contain errors. They are offered on an as-is basis. In -- +-- particular, achieving the same results as the tests here is not -- +-- a guarantee that an implementation complies with any Standard -- +-- or specification. The tests are not exhaustive. -- +-- -- +-- Please send comments, suggestions, and corrections to the author: -- +-- Mike Cowlishaw, IBM Fellow -- +-- IBM UK, PO Box 31, Birmingham Road, Warwick CV34 5JL, UK -- +-- mfc@uk.ibm.com -- +------------------------------------------------------------------------ +version: 2.38 + +-- This set of testcases tests raising numbers to an integer power only. +-- If arbitrary powers were supported, 1 ulp differences would be +-- permitted. + +extended: 1 +precision: 9 +rounding: half_up +maxExponent: 999 +minexponent: -999 + +-- base checks. Note 0**0 is an error. +powx001 power '0' '0' -> NaN Invalid_operation +powx002 power '0' '1' -> '0' +powx003 power '0' '2' -> '0' +powx004 power '1' '0' -> '1' +powx005 power '1' '1' -> '1' +powx006 power '1' '2' -> '1' + +powx010 power '2' '0' -> '1' +powx011 power '2' '1' -> '2' +powx012 power '2' '2' -> '4' +powx013 power '2' '3' -> '8' +powx014 power '2' '4' -> '16' +powx015 power '2' '5' -> '32' +powx016 power '2' '6' -> '64' +powx017 power '2' '7' -> '128' +powx018 power '2' '8' -> '256' +powx019 power '2' '9' -> '512' +powx020 power '2' '10' -> '1024' +powx021 power '2' '11' -> '2048' +powx022 power '2' '12' -> '4096' +powx023 power '2' '15' -> '32768' +powx024 power '2' '16' -> '65536' +powx025 power '2' '31' -> '2.14748365E+9' Inexact Rounded +-- NB 0 not stripped in next +powx026 power '2' '32' -> '4.29496730E+9' Inexact Rounded +precision: 10 +powx027 power '2' '31' -> '2147483648' +powx028 power '2' '32' -> '4294967296' +precision: 9 + +powx030 power '3' '2' -> 9 +powx031 power '4' '2' -> 16 +powx032 power '5' '2' -> 25 +powx033 power '6' '2' -> 36 +powx034 power '7' '2' -> 49 +powx035 power '8' '2' -> 64 +powx036 power '9' '2' -> 81 +powx037 power '10' '2' -> 100 +powx038 power '11' '2' -> 121 +powx039 power '12' '2' -> 144 + +powx040 power '3' '3' -> 27 +powx041 power '4' '3' -> 64 +powx042 power '5' '3' -> 125 +powx043 power '6' '3' -> 216 +powx044 power '7' '3' -> 343 + +powx050 power '10' '0' -> 1 +powx051 power '10' '1' -> 10 +powx052 power '10' '2' -> 100 +powx053 power '10' '3' -> 1000 +powx054 power '10' '4' -> 10000 +powx055 power '10' '5' -> 100000 +powx056 power '10' '6' -> 1000000 +powx057 power '10' '7' -> 10000000 +powx058 power '10' '8' -> 100000000 +powx059 power '10' '9' -> 1.00000000E+9 Rounded +powx060 power '10' '22' -> 1.00000000E+22 Rounded +powx061 power '10' '77' -> 1.00000000E+77 Rounded +powx062 power '10' '99' -> 1.00000000E+99 Rounded + +maxexponent: 999999999 +minexponent: -999999999 +powx063 power '10' '999999999' -> '1.00000000E+999999999' Rounded +powx064 power '10' '999999998' -> '1.00000000E+999999998' Rounded +powx065 power '10' '999999997' -> '1.00000000E+999999997' Rounded +powx066 power '10' '333333333' -> '1.00000000E+333333333' Rounded + +powx070 power '0.3' '0' -> '1' +powx071 power '0.3' '1' -> '0.3' +powx072 power '0.3' '1.00' -> '0.3' +powx073 power '0.3' '2.00' -> '0.09' +powx074 power '0.3' '2.000000000' -> '0.09' +powx075 power '6.0' '1' -> '6.0' -- NB zeros not stripped +powx076 power '6.0' '2' -> '36.00' -- .. +powx077 power '-3' '2' -> '9' -- from NetRexx book +powx078 power '4' '3' -> '64' -- .. (sort of) + +powx080 power 0.1 0 -> 1 +powx081 power 0.1 1 -> 0.1 +powx082 power 0.1 2 -> 0.01 +powx083 power 0.1 3 -> 0.001 +powx084 power 0.1 4 -> 0.0001 +powx085 power 0.1 5 -> 0.00001 +powx086 power 0.1 6 -> 0.000001 +powx087 power 0.1 7 -> 1E-7 +powx088 power 0.1 8 -> 1E-8 +powx089 power 0.1 9 -> 1E-9 + +powx090 power 101 2 -> 10201 +powx091 power 101 3 -> 1030301 +powx092 power 101 4 -> 104060401 +powx093 power 101 5 -> 1.05101005E+10 Inexact Rounded +powx094 power 101 6 -> 1.06152015E+12 Inexact Rounded +powx095 power 101 7 -> 1.07213535E+14 Inexact Rounded + +-- negative powers +powx101 power '2' '-1' -> 0.5 +powx102 power '2' '-2' -> 0.25 +powx103 power '2' '-4' -> 0.0625 +powx104 power '2' '-8' -> 0.00390625 +powx105 power '2' '-16' -> 0.0000152587891 Inexact Rounded +powx106 power '2' '-32' -> 2.32830644E-10 Inexact Rounded +powx108 power '2' '-64' -> 5.42101086E-20 Inexact Rounded +powx110 power '10' '-8' -> 1E-8 +powx111 power '10' '-7' -> 1E-7 +powx112 power '10' '-6' -> 0.000001 +powx113 power '10' '-5' -> 0.00001 +powx114 power '10' '-4' -> 0.0001 +powx115 power '10' '-3' -> 0.001 +powx116 power '10' '-2' -> 0.01 +powx117 power '10' '-1' -> 0.1 + +powx118 power '10' '-333333333' -> 1E-333333333 +powx119 power '10' '-999999998' -> 1E-999999998 +powx120 power '10' '-999999999' -> 1E-999999999 +powx121 power '10' '-77' -> '1E-77' +powx122 power '10' '-22' -> '1E-22' + +powx123 power '2' '-1' -> '0.5' +powx124 power '2' '-2' -> '0.25' +powx125 power '2' '-4' -> '0.0625' +powx126 power '0' '-1' -> Infinity Division_by_zero +powx127 power '0' '-2' -> Infinity Division_by_zero +powx128 power -0 '-1' -> -Infinity Division_by_zero +powx129 power -0 '-2' -> Infinity Division_by_zero + +-- out-of-range edge cases +powx181 power '7' '999999998' -> 2.10892313E+845098038 Inexact Rounded +powx182 power '7' '999999999' -> 1.47624619E+845098039 Inexact Rounded +powx183 power '7' '1000000000' -> NaN Invalid_operation +powx184 power '7' '1000000001' -> NaN Invalid_operation +powx185 power '7' '10000000000' -> NaN Invalid_operation +powx186 power '7' '-1000000001' -> NaN Invalid_operation +powx187 power '7' '-1000000000' -> NaN Invalid_operation +powx189 power '7' '-999999999' -> 6.77393787E-845098040 Inexact Rounded +powx190 power '7' '-999999998' -> 4.74175651E-845098039 Inexact Rounded + +-- some baddies [more below] +powx191 power '2' '2.000001' -> NaN Invalid_operation +powx192 power '2' '2.00000000' -> 4 +powx193 power '2' '2.000000001' -> NaN Invalid_operation +powx194 power '2' '2.0000000001' -> NaN Invalid_operation + +-- "0.5" tests from original Rexx diagnostics [loop unrolled] +powx200 power 0.5 0 -> 1 +powx201 power 0.5 1 -> 0.5 +powx202 power 0.5 2 -> 0.25 +powx203 power 0.5 3 -> 0.125 +powx204 power 0.5 4 -> 0.0625 +powx205 power 0.5 5 -> 0.03125 +powx206 power 0.5 6 -> 0.015625 +powx207 power 0.5 7 -> 0.0078125 +powx208 power 0.5 8 -> 0.00390625 +powx209 power 0.5 9 -> 0.001953125 +powx210 power 0.5 10 -> 0.0009765625 + +-- A (rare) case where the last digit is not within 0.5 ULP +precision: 9 +powx215 power "-21971575.0E+31454441" "-7" -> "-4.04549503E-220181139" Inexact Rounded +precision: 20 +powx216 power "-21971575.0E+31454441" "-7" -> "-4.0454950249324891788E-220181139" Inexact Rounded + +-- The Vienna case. Checks both setup and 1/acc working precision +-- Modified 1998.12.14 as RHS no longer rounded before use (must fit) +-- Modified 1990.02.04 as LHS is now rounded (instead of truncated to guard) +-- '123456789E+10' -- lhs .. rounded to 1.23E+18 +-- '-1.23000e+2' -- rhs .. [was: -1.23455e+2, rounds to -123] +-- Modified 2002.10.06 -- finally, no input rounding +-- With input rounding, result would be 8.74E-2226 +precision: 3 +powx219 power '123456789E+10' '-1.23000e+2' -> '5.54E-2226' Inexact Rounded + +-- whole number checks +precision: 9 +powx221 power 1 1234 -> 1 +precision: 4 +powx222 power 1 1234 -> 1 +precision: 3 +powx223 power 1 1234 -> 1 +powx224 power 1 12.34e+2 -> 1 +powx225 power 1 12.3 -> NaN Invalid_operation +powx226 power 1 12.0 -> 1 +powx227 power 1 1.01 -> NaN Invalid_operation +powx228 power 2 1.00 -> 2 +powx229 power 2 2.00 -> 4 +precision: 9 +powx230 power 1 1.0001 -> NaN Invalid_operation +powx231 power 1 1.0000001 -> NaN Invalid_operation +powx232 power 1 1.0000000001 -> NaN Invalid_operation +powx233 power 1 1.0000000000001 -> NaN Invalid_operation +precision: 5 +powx234 power 1 1.0001 -> NaN Invalid_operation +powx235 power 1 1.0000001 -> NaN Invalid_operation +powx236 power 1 1.0000000001 -> NaN Invalid_operation +powx237 power 1 1.0000000000001 -> NaN Invalid_operation +powx238 power 1 1.0000000000001 -> NaN Invalid_operation + +maxexponent: 999999999 +minexponent: -999999999 +powx239 power 1 5.67E-987654321 -> NaN Invalid_operation + +powx240 power 1 100000000 -> 1 +powx241 power 1 999999998 -> 1 +powx242 power 1 999999999 -> 1 +powx243 power 1 1000000000 -> NaN Invalid_operation +powx244 power 1 9999999999 -> NaN Invalid_operation + +-- Checks for 'Too much precision needed' +-- For x^12, digits+elength+1 = digits+3 +precision: 999999999 +powx249 add 1 1 -> 2 -- check basic operation at this precision +powx250 power 2 12 -> Infinity Overflow +precision: 999999998 +powx251 power 2 12 -> Infinity Overflow +precision: 999999997 +powx252 power 2 12 -> Infinity Overflow +precision: 999999996 +powx253 power 2 12 -> 4096 +precision: 999999995 +powx254 power 2 12 -> 4096 + +-- zeros +maxexponent: +96 +minexponent: -95 +precision: 7 +powx260 power 0E-34 3 -> 0E-101 Clamped +powx261 power 0E-33 3 -> 0E-99 +powx262 power 0E-32 3 -> 0E-96 +powx263 power 0E-30 3 -> 0E-90 +powx264 power 0E-10 3 -> 0E-30 +powx265 power 0E-1 3 -> 0.000 +powx266 power 0E+0 3 -> 0 +powx267 power 0 3 -> 0 +powx268 power 0E+1 3 -> 0E+3 +powx269 power 0E+10 3 -> 0E+30 +powx270 power 0E+30 3 -> 0E+90 +powx271 power 0E+32 3 -> 0E+96 +powx272 power 0E+33 3 -> 0E+96 Clamped + +-- overflow and underflow tests +maxexponent: 999999999 +minexponent: -999999999 +precision: 9 +powx280 power 9 999999999 -> 3.05550054E+954242508 Inexact Rounded +powx281 power 10 999999999 -> 1.00000000E+999999999 Rounded +powx282 power 10.0001 999999999 -> Infinity Overflow Inexact Rounded +powx283 power 10.1 999999999 -> Infinity Overflow Inexact Rounded +powx284 power 11 999999999 -> Infinity Overflow Inexact Rounded +powx285 power 12 999999999 -> Infinity Overflow Inexact Rounded +powx286 power 999 999999999 -> Infinity Overflow Inexact Rounded +powx287 power 999999 999999999 -> Infinity Overflow Inexact Rounded +powx288 power 999999999 999999999 -> Infinity Overflow Inexact Rounded +powx289 power 9.9E999999999 999999999 -> Infinity Overflow Inexact Rounded + +powx290 power 0.5 999999999 -> 4.33559594E-301029996 Inexact Rounded +powx291 power 0.1 999999999 -> 1E-999999999 -- unrounded +powx292 power 0.09 999999999 -> 0E-1000000007 Underflow Subnormal Inexact Rounded Clamped +powx293 power 0.05 999999999 -> 0E-1000000007 Underflow Subnormal Inexact Rounded Clamped +powx294 power 0.01 999999999 -> 0E-1000000007 Underflow Subnormal Inexact Rounded Clamped +powx295 power 0.0001 999999999 -> 0E-1000000007 Underflow Subnormal Inexact Rounded Clamped +powx297 power 0.0000001 999999999 -> 0E-1000000007 Underflow Subnormal Inexact Rounded Clamped +powx298 power 0.0000000001 999999999 -> 0E-1000000007 Underflow Subnormal Inexact Rounded Clamped +powx299 power 1E-999999999 999999999 -> 0E-1000000007 Underflow Subnormal Inexact Rounded Clamped + +powx310 power -9 999999999 -> -3.05550054E+954242508 Inexact Rounded +powx311 power -10 999999999 -> -1.00000000E+999999999 Rounded +powx312 power -10.0001 999999999 -> -Infinity Overflow Inexact Rounded +powx313 power -10.1 999999999 -> -Infinity Overflow Inexact Rounded +powx314 power -11 999999999 -> -Infinity Overflow Inexact Rounded +powx315 power -12 999999999 -> -Infinity Overflow Inexact Rounded +powx316 power -999 999999999 -> -Infinity Overflow Inexact Rounded +powx317 power -999999 999999999 -> -Infinity Overflow Inexact Rounded +powx318 power -999999999 999999999 -> -Infinity Overflow Inexact Rounded +powx319 power -9.9E999999999 999999999 -> -Infinity Overflow Inexact Rounded + +powx320 power -0.5 999999999 -> -4.33559594E-301029996 Inexact Rounded +powx321 power -0.1 999999999 -> -1E-999999999 +powx322 power -0.09 999999999 -> -0E-1000000007 Underflow Subnormal Inexact Rounded Clamped +powx323 power -0.05 999999999 -> -0E-1000000007 Underflow Subnormal Inexact Rounded Clamped +powx324 power -0.01 999999999 -> -0E-1000000007 Underflow Subnormal Inexact Rounded Clamped +powx325 power -0.0001 999999999 -> -0E-1000000007 Underflow Subnormal Inexact Rounded Clamped +powx327 power -0.0000001 999999999 -> -0E-1000000007 Underflow Subnormal Inexact Rounded Clamped +powx328 power -0.0000000001 999999999 -> -0E-1000000007 Underflow Subnormal Inexact Rounded Clamped +powx329 power -1E-999999999 999999999 -> -0E-1000000007 Underflow Subnormal Inexact Rounded Clamped + +-- note no trim of next result +powx330 power -9 999999998 -> 3.39500060E+954242507 Inexact Rounded +powx331 power -10 999999998 -> 1.00000000E+999999998 Rounded +powx332 power -10.0001 999999998 -> Infinity Overflow Inexact Rounded +powx333 power -10.1 999999998 -> Infinity Overflow Inexact Rounded +powx334 power -11 999999998 -> Infinity Overflow Inexact Rounded +powx335 power -12 999999998 -> Infinity Overflow Inexact Rounded +powx336 power -999 999999998 -> Infinity Overflow Inexact Rounded +powx337 power -999999 999999998 -> Infinity Overflow Inexact Rounded +powx338 power -999999999 999999998 -> Infinity Overflow Inexact Rounded +powx339 power -9.9E999999999 999999998 -> Infinity Overflow Inexact Rounded + +powx340 power -0.5 999999998 -> 8.67119187E-301029996 Inexact Rounded +powx341 power -0.1 999999998 -> 1E-999999998 -- NB exact unrounded +powx342 power -0.09 999999998 -> 0E-1000000007 Underflow Subnormal Inexact Rounded Clamped +powx343 power -0.05 999999998 -> 0E-1000000007 Underflow Subnormal Inexact Rounded Clamped +powx344 power -0.01 999999998 -> 0E-1000000007 Underflow Subnormal Inexact Rounded Clamped +powx345 power -0.0001 999999998 -> 0E-1000000007 Underflow Subnormal Inexact Rounded Clamped +powx347 power -0.0000001 999999998 -> 0E-1000000007 Underflow Subnormal Inexact Rounded Clamped +powx348 power -0.0000000001 999999998 -> 0E-1000000007 Underflow Subnormal Inexact Rounded Clamped +powx349 power -1E-999999999 999999998 -> 0E-1000000007 Underflow Subnormal Inexact Rounded Clamped + +-- some subnormals +precision: 9 +-- [precision is 9, so smallest exponent is -1000000007 +powx350 power 1e-1 500000000 -> 1E-500000000 +powx351 power 1e-2 999999999 -> 0E-1000000007 Underflow Subnormal Inexact Rounded Clamped +powx352 power 1e-2 500000000 -> 1E-1000000000 Subnormal +powx353 power 1e-2 500000001 -> 1E-1000000002 Subnormal +powx354 power 1e-2 500000002 -> 1E-1000000004 Subnormal +powx355 power 1e-2 500000003 -> 1E-1000000006 Subnormal +powx356 power 1e-2 500000004 -> 0E-1000000007 Underflow Subnormal Inexact Rounded + +powx360 power 0.010001 500000000 -> 4.34941988E-999978287 Inexact Rounded +powx361 power 0.010000001 500000000 -> 5.18469257E-999999979 Inexact Rounded +powx362 power 0.010000001 500000001 -> 5.18469309E-999999981 Inexact Rounded +powx363 power 0.0100000009 500000000 -> 3.49342003E-999999981 Inexact Rounded +powx364 power 0.0100000001 500000000 -> 1.48413155E-999999998 Inexact Rounded +powx365 power 0.01 500000000 -> 1E-1000000000 Subnormal +powx366 power 0.0099999999 500000000 -> 6.7379E-1000000003 Underflow Subnormal Inexact Rounded +powx367 power 0.0099999998 500000000 -> 4.54E-1000000005 Underflow Subnormal Inexact Rounded +powx368 power 0.0099999997 500000000 -> 3E-1000000007 Underflow Subnormal Inexact Rounded +powx369 power 0.0099999996 500000000 -> 0E-1000000007 Underflow Subnormal Inexact Rounded +powx370 power 0.009 500000000 -> 0E-1000000007 Underflow Subnormal Inexact Rounded Clamped + +-- 1/subnormal -> overflow +powx371 power 1e-1 -500000000 -> 1E+500000000 +powx372 power 1e-2 -999999999 -> Infinity Overflow Inexact Rounded +powx373 power 1e-2 -500000000 -> Infinity Overflow Inexact Rounded +powx374 power 1e-2 -500000001 -> Infinity Overflow Inexact Rounded +powx375 power 1e-2 -500000002 -> Infinity Overflow Inexact Rounded +powx376 power 1e-2 -500000003 -> Infinity Overflow Inexact Rounded +powx377 power 1e-2 -500000004 -> Infinity Overflow Inexact Rounded + +powx381 power 0.010001 -500000000 -> 2.29915719E+999978286 Inexact Rounded +powx382 power 0.010000001 -500000000 -> 1.92875467E+999999978 Inexact Rounded +powx383 power 0.010000001 -500000001 -> 1.92875448E+999999980 Inexact Rounded +powx384 power 0.0100000009 -500000000 -> 2.86252438E+999999980 Inexact Rounded +powx385 power 0.0100000001 -500000000 -> 6.73794717E+999999997 Inexact Rounded +powx386 power 0.01 -500000000 -> Infinity Overflow Inexact Rounded +powx387 power 0.009999 -500000000 -> Infinity Overflow Inexact Rounded + +-- negative power giving subnormal +powx388 power 100.000001 -500000000 -> 6.7379E-1000000003 Underflow Subnormal Inexact Rounded + +-- some more edge cases +precision: 15 +maxExponent: 999 +minexponent: -999 +powx391 power 0.1 999 -> 1E-999 +powx392 power 0.099 999 -> 4.360732062E-1004 Underflow Subnormal Inexact Rounded +powx393 power 0.098 999 -> 1.71731E-1008 Underflow Subnormal Inexact Rounded +powx394 power 0.097 999 -> 6E-1013 Underflow Subnormal Inexact Rounded +powx395 power 0.096 999 -> 0E-1013 Underflow Subnormal Inexact Rounded +powx396 power 0.01 999 -> 0E-1013 Underflow Subnormal Inexact Rounded Clamped + +-- multiply tests are here to aid checking and test for consistent handling +-- of underflow +precision: 5 +maxexponent: 999 +minexponent: -999 + +-- squares +mulx400 multiply 1E-502 1e-502 -> 0E-1003 Subnormal Inexact Underflow Rounded +mulx401 multiply 1E-501 1e-501 -> 1E-1002 Subnormal +mulx402 multiply 2E-501 2e-501 -> 4E-1002 Subnormal +mulx403 multiply 4E-501 4e-501 -> 1.6E-1001 Subnormal +mulx404 multiply 10E-501 10e-501 -> 1.00E-1000 Subnormal +mulx405 multiply 30E-501 30e-501 -> 9.00E-1000 Subnormal +mulx406 multiply 40E-501 40e-501 -> 1.600E-999 + +powx400 power 1E-502 2 -> 0E-1003 Underflow Subnormal Inexact Rounded +powx401 power 1E-501 2 -> 1E-1002 Subnormal +powx402 power 2E-501 2 -> 4E-1002 Subnormal +powx403 power 4E-501 2 -> 1.6E-1001 Subnormal +powx404 power 10E-501 2 -> 1.00E-1000 Subnormal +powx405 power 30E-501 2 -> 9.00E-1000 Subnormal +powx406 power 40E-501 2 -> 1.600E-999 + +-- cubes +mulx410 multiply 1E-670 1e-335 -> 0E-1003 Underflow Subnormal Inexact Rounded +mulx411 multiply 1E-668 1e-334 -> 1E-1002 Subnormal +mulx412 multiply 4E-668 2e-334 -> 8E-1002 Subnormal +mulx413 multiply 9E-668 3e-334 -> 2.7E-1001 Subnormal +mulx414 multiply 16E-668 4e-334 -> 6.4E-1001 Subnormal +mulx415 multiply 25E-668 5e-334 -> 1.25E-1000 Subnormal +mulx416 multiply 10E-668 100e-334 -> 1.000E-999 + +powx410 power 1E-335 3 -> 0E-1003 Underflow Subnormal Inexact Rounded +powx411 power 1E-334 3 -> 1E-1002 Subnormal +powx412 power 2E-334 3 -> 8E-1002 Subnormal +powx413 power 3E-334 3 -> 2.7E-1001 Subnormal +powx414 power 4E-334 3 -> 6.4E-1001 Subnormal +powx415 power 5E-334 3 -> 1.25E-1000 Subnormal +powx416 power 10E-334 3 -> 1.000E-999 + +-- negative powers, testing subnormals +precision: 5 +maxExponent: 999 +minexponent: -999 +powx421 power 2.5E-501 -2 -> Infinity Overflow Inexact Rounded +powx422 power 2.5E-500 -2 -> 1.6E+999 + +powx423 power 2.5E+499 -2 -> 1.6E-999 +powx424 power 2.5E+500 -2 -> 1.6E-1001 Subnormal +powx425 power 2.5E+501 -2 -> 2E-1003 Underflow Subnormal Inexact Rounded +powx426 power 2.5E+502 -2 -> 0E-1003 Underflow Subnormal Inexact Rounded + +powx427 power 0.25E+499 -2 -> 1.6E-997 +powx428 power 0.25E+500 -2 -> 1.6E-999 +powx429 power 0.25E+501 -2 -> 1.6E-1001 Subnormal +powx430 power 0.25E+502 -2 -> 2E-1003 Underflow Subnormal Inexact Rounded +powx431 power 0.25E+503 -2 -> 0E-1003 Underflow Subnormal Inexact Rounded + +powx432 power 0.04E+499 -2 -> 6.25E-996 +powx433 power 0.04E+500 -2 -> 6.25E-998 +powx434 power 0.04E+501 -2 -> 6.25E-1000 Subnormal +powx435 power 0.04E+502 -2 -> 6.3E-1002 Underflow Subnormal Inexact Rounded +powx436 power 0.04E+503 -2 -> 1E-1003 Underflow Subnormal Inexact Rounded +powx437 power 0.04E+504 -2 -> 0E-1003 Underflow Subnormal Inexact Rounded + +powx441 power 0.04E+334 -3 -> 1.5625E-998 +powx442 power 0.04E+335 -3 -> 1.56E-1001 Underflow Subnormal Inexact Rounded +powx443 power 0.04E+336 -3 -> 0E-1003 Underflow Subnormal Inexact Rounded +powx444 power 0.25E+333 -3 -> 6.4E-998 +powx445 power 0.25E+334 -3 -> 6.4E-1001 Subnormal +powx446 power 0.25E+335 -3 -> 1E-1003 Underflow Subnormal Inexact Rounded +powx447 power 0.25E+336 -3 -> 0E-1003 Underflow Subnormal Inexact Rounded Clamped +-- check sign for cubes and a few squares +powx448 power -0.04E+334 -3 -> -1.5625E-998 +powx449 power -0.04E+335 -3 -> -1.56E-1001 Underflow Subnormal Inexact Rounded +powx450 power -0.04E+336 -3 -> -0E-1003 Underflow Subnormal Inexact Rounded +powx451 power -0.25E+333 -3 -> -6.4E-998 +powx452 power -0.25E+334 -3 -> -6.4E-1001 Subnormal +powx453 power -0.25E+335 -3 -> -1E-1003 Underflow Subnormal Inexact Rounded +powx454 power -0.25E+336 -3 -> -0E-1003 Underflow Subnormal Inexact Rounded Clamped +powx455 power -0.04E+499 -2 -> 6.25E-996 +powx456 power -0.04E+500 -2 -> 6.25E-998 +powx457 power -0.04E+501 -2 -> 6.25E-1000 Subnormal +powx458 power -0.04E+502 -2 -> 6.3E-1002 Underflow Subnormal Inexact Rounded + +-- test -0s +precision: 9 +powx560 power 0 0 -> NaN Invalid_operation +powx561 power 0 -0 -> NaN Invalid_operation +powx562 power -0 0 -> NaN Invalid_operation +powx563 power -0 -0 -> NaN Invalid_operation +powx564 power 1 0 -> 1 +powx565 power 1 -0 -> 1 +powx566 power -1 0 -> 1 +powx567 power -1 -0 -> 1 +powx568 power 0 1 -> 0 +powx569 power 0 -1 -> Infinity Division_by_zero +powx570 power -0 1 -> -0 +powx571 power -0 -1 -> -Infinity Division_by_zero +powx572 power 0 2 -> 0 +powx573 power 0 -2 -> Infinity Division_by_zero +powx574 power -0 2 -> 0 +powx575 power -0 -2 -> Infinity Division_by_zero +powx576 power 0 3 -> 0 +powx577 power 0 -3 -> Infinity Division_by_zero +powx578 power -0 3 -> -0 +powx579 power -0 -3 -> -Infinity Division_by_zero + +-- Specials +powx580 power Inf -Inf -> NaN Invalid_operation +powx581 power Inf -1000 -> 0 +powx582 power Inf -1 -> 0 +powx583 power Inf -0 -> 1 +powx584 power Inf 0 -> 1 +powx585 power Inf 1 -> Infinity +powx586 power Inf 1000 -> Infinity +powx587 power Inf Inf -> NaN Invalid_operation +powx588 power -1000 Inf -> NaN Invalid_operation +powx589 power -Inf Inf -> NaN Invalid_operation +powx590 power -1 Inf -> NaN Invalid_operation +powx591 power -0 Inf -> NaN Invalid_operation +powx592 power 0 Inf -> NaN Invalid_operation +powx593 power 1 Inf -> NaN Invalid_operation +powx594 power 1000 Inf -> NaN Invalid_operation +powx595 power Inf Inf -> NaN Invalid_operation + +powx600 power -Inf -Inf -> NaN Invalid_operation +powx601 power -Inf -1000 -> 0 +powx602 power -Inf -1 -> -0 +powx603 power -Inf -0 -> 1 +powx604 power -Inf 0 -> 1 +powx605 power -Inf 1 -> -Infinity +powx606 power -Inf 1000 -> Infinity +powx607 power -Inf Inf -> NaN Invalid_operation +powx608 power -1000 Inf -> NaN Invalid_operation +powx609 power -Inf -Inf -> NaN Invalid_operation +powx610 power -1 -Inf -> NaN Invalid_operation +powx611 power -0 -Inf -> NaN Invalid_operation +powx612 power 0 -Inf -> NaN Invalid_operation +powx613 power 1 -Inf -> NaN Invalid_operation +powx614 power 1000 -Inf -> NaN Invalid_operation +powx615 power Inf -Inf -> NaN Invalid_operation + +powx621 power NaN -Inf -> NaN Invalid_operation +powx622 power NaN -1000 -> NaN +powx623 power NaN -1 -> NaN +powx624 power NaN -0 -> NaN +powx625 power NaN 0 -> NaN +powx626 power NaN 1 -> NaN +powx627 power NaN 1000 -> NaN +powx628 power NaN Inf -> NaN Invalid_operation +powx629 power NaN NaN -> NaN +powx630 power -Inf NaN -> NaN +powx631 power -1000 NaN -> NaN +powx632 power -1 NaN -> NaN +powx633 power -0 NaN -> NaN +powx634 power 0 NaN -> NaN +powx635 power 1 NaN -> NaN +powx636 power 1000 NaN -> NaN +powx637 power Inf NaN -> NaN + +powx641 power sNaN -Inf -> NaN Invalid_operation +powx642 power sNaN -1000 -> NaN Invalid_operation +powx643 power sNaN -1 -> NaN Invalid_operation +powx644 power sNaN -0 -> NaN Invalid_operation +powx645 power sNaN 0 -> NaN Invalid_operation +powx646 power sNaN 1 -> NaN Invalid_operation +powx647 power sNaN 1000 -> NaN Invalid_operation +powx648 power sNaN NaN -> NaN Invalid_operation +powx649 power sNaN sNaN -> NaN Invalid_operation +powx650 power NaN sNaN -> NaN Invalid_operation +powx651 power -Inf sNaN -> NaN Invalid_operation +powx652 power -1000 sNaN -> NaN Invalid_operation +powx653 power -1 sNaN -> NaN Invalid_operation +powx654 power -0 sNaN -> NaN Invalid_operation +powx655 power 0 sNaN -> NaN Invalid_operation +powx656 power 1 sNaN -> NaN Invalid_operation +powx657 power 1000 sNaN -> NaN Invalid_operation +powx658 power Inf sNaN -> NaN Invalid_operation +powx659 power NaN sNaN -> NaN Invalid_operation + +-- NaN propagation +powx660 power NaN3 sNaN7 -> NaN7 Invalid_operation +powx661 power sNaN8 NaN6 -> NaN8 Invalid_operation +powx662 power 1 sNaN7 -> NaN7 Invalid_operation +powx663 power sNaN8 1 -> NaN8 Invalid_operation +powx664 power Inf sNaN7 -> NaN7 Invalid_operation +powx665 power sNaN8 Inf -> NaN Invalid_operation +powx666 power Inf NaN9 -> NaN9 +powx667 power NaN6 Inf -> NaN Invalid_operation +powx668 power 1 NaN5 -> NaN5 +powx669 power NaN2 1 -> NaN2 +powx670 power NaN2 Nan4 -> NaN2 +powx671 power NaN Nan4 -> NaN +powx672 power NaN345 Nan -> NaN345 +powx673 power Inf -sNaN7 -> -NaN7 Invalid_operation +powx674 power -sNaN8 Inf -> NaN Invalid_operation +powx675 power Inf -NaN9 -> -NaN9 +powx676 power -NaN6 Inf -> NaN Invalid_operation +powx677 power -NaN2 -Nan4 -> -NaN2 + +-- Examples from extended specification +powx690 power Inf -2 -> 0 +powx691 power Inf -1 -> 0 +powx692 power Inf 0 -> 1 +powx693 power Inf 1 -> Infinity +powx694 power Inf 2 -> Infinity +powx695 power -Inf -2 -> 0 +powx696 power -Inf -1 -> -0 +powx697 power -Inf 0 -> 1 +powx698 power -Inf 1 -> -Infinity +powx699 power -Inf 2 -> Infinity +powx700 power 0 0 -> NaN Invalid_operation + +-- long operand and RHS range checks +maxexponent: 999 +minexponent: -999 +precision: 9 +powx701 power 12345678000 1 -> 1.23456780E+10 Rounded +powx702 power 1234567800 1 -> 1.23456780E+9 Rounded +powx703 power 1234567890 1 -> 1.23456789E+9 Rounded +powx704 power 1234567891 1 -> 1.23456789E+9 Inexact Rounded +powx705 power 12345678901 1 -> 1.23456789E+10 Inexact Rounded +powx706 power 1234567896 1 -> 1.23456790E+9 Inexact Rounded +powx707 power 1 12345678000 -> NaN Invalid_operation +powx708 power 1 1234567800 -> NaN Invalid_operation +powx709 power 1 1234567890 -> NaN Invalid_operation +powx710 power 1 11234567891 -> NaN Invalid_operation +powx711 power 1 12345678901 -> NaN Invalid_operation +powx712 power 1 1234567896 -> NaN Invalid_operation +powx713 power 1 -1234567896 -> NaN Invalid_operation +powx714 power 1 1000000000 -> NaN Invalid_operation +powx715 power 1 -1000000000 -> NaN Invalid_operation + +precision: 15 +-- still checking +powx741 power 12345678000 1 -> 12345678000 +powx742 power 1234567800 1 -> 1234567800 +powx743 power 1234567890 1 -> 1234567890 +powx744 power 1234567891 1 -> 1234567891 +powx745 power 12345678901 1 -> 12345678901 +powx746 power 1234567896 1 -> 1234567896 +powx747 power 1 12345678000 -> NaN Invalid_operation +powx748 power 1 -1234567896 -> NaN Invalid_operation +powx749 power 1 1000000000 -> NaN Invalid_operation +powx740 power 1 -1000000000 -> NaN Invalid_operation + +-- check for double-rounded subnormals +precision: 5 +maxexponent: 79 +minexponent: -79 +powx750 power 1.2347E-40 2 -> 1.524E-80 Inexact Rounded Subnormal Underflow + +-- Null tests +powx900 power 1 # -> NaN Invalid_operation +powx901 power # 1 -> NaN Invalid_operation + |