Raymond Hettinger | 7c85fa4 | 2004-07-01 11:01:35 +0000 | [diff] [blame] | 1 | ------------------------------------------------------------------------ |
| 2 | -- multiply.decTest -- decimal multiplication -- |
Benjamin Peterson | f17ff4e | 2008-07-31 16:32:12 +0000 | [diff] [blame] | 3 | -- Copyright (c) IBM Corporation, 1981, 2008. All rights reserved. -- |
Raymond Hettinger | 7c85fa4 | 2004-07-01 11:01:35 +0000 | [diff] [blame] | 4 | ------------------------------------------------------------------------ |
| 5 | -- Please see the document "General Decimal Arithmetic Testcases" -- |
| 6 | -- at http://www2.hursley.ibm.com/decimal for the description of -- |
| 7 | -- these testcases. -- |
| 8 | -- -- |
| 9 | -- These testcases are experimental ('beta' versions), and they -- |
| 10 | -- may contain errors. They are offered on an as-is basis. In -- |
| 11 | -- particular, achieving the same results as the tests here is not -- |
| 12 | -- a guarantee that an implementation complies with any Standard -- |
| 13 | -- or specification. The tests are not exhaustive. -- |
| 14 | -- -- |
| 15 | -- Please send comments, suggestions, and corrections to the author: -- |
| 16 | -- Mike Cowlishaw, IBM Fellow -- |
| 17 | -- IBM UK, PO Box 31, Birmingham Road, Warwick CV34 5JL, UK -- |
| 18 | -- mfc@uk.ibm.com -- |
| 19 | ------------------------------------------------------------------------ |
Mark Dickinson | 8a54653 | 2009-10-08 16:30:38 +0000 | [diff] [blame] | 20 | version: 2.59 |
Raymond Hettinger | 7c85fa4 | 2004-07-01 11:01:35 +0000 | [diff] [blame] | 21 | |
| 22 | extended: 1 |
| 23 | precision: 9 |
| 24 | rounding: half_up |
| 25 | maxExponent: 384 |
| 26 | minexponent: -383 |
| 27 | |
| 28 | -- sanity checks (as base, above) |
| 29 | mulx000 multiply 2 2 -> 4 |
| 30 | mulx001 multiply 2 3 -> 6 |
| 31 | mulx002 multiply 5 1 -> 5 |
| 32 | mulx003 multiply 5 2 -> 10 |
| 33 | mulx004 multiply 1.20 2 -> 2.40 |
| 34 | mulx005 multiply 1.20 0 -> 0.00 |
| 35 | mulx006 multiply 1.20 -2 -> -2.40 |
| 36 | mulx007 multiply -1.20 2 -> -2.40 |
| 37 | mulx008 multiply -1.20 0 -> -0.00 |
| 38 | mulx009 multiply -1.20 -2 -> 2.40 |
| 39 | mulx010 multiply 5.09 7.1 -> 36.139 |
| 40 | mulx011 multiply 2.5 4 -> 10.0 |
| 41 | mulx012 multiply 2.50 4 -> 10.00 |
| 42 | mulx013 multiply 1.23456789 1.00000000 -> 1.23456789 Rounded |
| 43 | mulx014 multiply 9.999999999 9.999999999 -> 100.000000 Inexact Rounded |
| 44 | mulx015 multiply 2.50 4 -> 10.00 |
| 45 | precision: 6 |
| 46 | mulx016 multiply 2.50 4 -> 10.00 |
Thomas Wouters | 1b7f891 | 2007-09-19 03:06:30 +0000 | [diff] [blame] | 47 | mulx017 multiply 9.999999999 9.999999999 -> 100.000 Inexact Rounded |
| 48 | mulx018 multiply 9.999999999 -9.999999999 -> -100.000 Inexact Rounded |
| 49 | mulx019 multiply -9.999999999 9.999999999 -> -100.000 Inexact Rounded |
| 50 | mulx020 multiply -9.999999999 -9.999999999 -> 100.000 Inexact Rounded |
Raymond Hettinger | 7c85fa4 | 2004-07-01 11:01:35 +0000 | [diff] [blame] | 51 | |
| 52 | -- 1999.12.21: next one is a edge case if intermediate longs are used |
| 53 | precision: 15 |
Thomas Wouters | 1b7f891 | 2007-09-19 03:06:30 +0000 | [diff] [blame] | 54 | mulx059 multiply 999999999999 9765625 -> 9.76562499999023E+18 Inexact Rounded |
Raymond Hettinger | 7c85fa4 | 2004-07-01 11:01:35 +0000 | [diff] [blame] | 55 | precision: 30 |
| 56 | mulx160 multiply 999999999999 9765625 -> 9765624999990234375 |
| 57 | precision: 9 |
| 58 | ----- |
| 59 | |
| 60 | -- zeros, etc. |
Thomas Wouters | 1b7f891 | 2007-09-19 03:06:30 +0000 | [diff] [blame] | 61 | mulx021 multiply 0 0 -> 0 |
| 62 | mulx022 multiply 0 -0 -> -0 |
| 63 | mulx023 multiply -0 0 -> -0 |
| 64 | mulx024 multiply -0 -0 -> 0 |
| 65 | mulx025 multiply -0.0 -0.0 -> 0.00 |
| 66 | mulx026 multiply -0.0 -0.0 -> 0.00 |
| 67 | mulx027 multiply -0.0 -0.0 -> 0.00 |
| 68 | mulx028 multiply -0.0 -0.0 -> 0.00 |
Raymond Hettinger | 7c85fa4 | 2004-07-01 11:01:35 +0000 | [diff] [blame] | 69 | mulx030 multiply 5.00 1E-3 -> 0.00500 |
| 70 | mulx031 multiply 00.00 0.000 -> 0.00000 |
| 71 | mulx032 multiply 00.00 0E-3 -> 0.00000 -- rhs is 0 |
| 72 | mulx033 multiply 0E-3 00.00 -> 0.00000 -- lhs is 0 |
| 73 | mulx034 multiply -5.00 1E-3 -> -0.00500 |
| 74 | mulx035 multiply -00.00 0.000 -> -0.00000 |
| 75 | mulx036 multiply -00.00 0E-3 -> -0.00000 -- rhs is 0 |
| 76 | mulx037 multiply -0E-3 00.00 -> -0.00000 -- lhs is 0 |
| 77 | mulx038 multiply 5.00 -1E-3 -> -0.00500 |
| 78 | mulx039 multiply 00.00 -0.000 -> -0.00000 |
| 79 | mulx040 multiply 00.00 -0E-3 -> -0.00000 -- rhs is 0 |
| 80 | mulx041 multiply 0E-3 -00.00 -> -0.00000 -- lhs is 0 |
| 81 | mulx042 multiply -5.00 -1E-3 -> 0.00500 |
| 82 | mulx043 multiply -00.00 -0.000 -> 0.00000 |
| 83 | mulx044 multiply -00.00 -0E-3 -> 0.00000 -- rhs is 0 |
| 84 | mulx045 multiply -0E-3 -00.00 -> 0.00000 -- lhs is 0 |
| 85 | |
| 86 | -- examples from decarith |
| 87 | mulx050 multiply 1.20 3 -> 3.60 |
| 88 | mulx051 multiply 7 3 -> 21 |
| 89 | mulx052 multiply 0.9 0.8 -> 0.72 |
| 90 | mulx053 multiply 0.9 -0 -> -0.0 |
| 91 | mulx054 multiply 654321 654321 -> 4.28135971E+11 Inexact Rounded |
| 92 | |
| 93 | mulx060 multiply 123.45 1e7 -> 1.2345E+9 |
| 94 | mulx061 multiply 123.45 1e8 -> 1.2345E+10 |
| 95 | mulx062 multiply 123.45 1e+9 -> 1.2345E+11 |
| 96 | mulx063 multiply 123.45 1e10 -> 1.2345E+12 |
| 97 | mulx064 multiply 123.45 1e11 -> 1.2345E+13 |
| 98 | mulx065 multiply 123.45 1e12 -> 1.2345E+14 |
| 99 | mulx066 multiply 123.45 1e13 -> 1.2345E+15 |
| 100 | |
| 101 | |
| 102 | -- test some intermediate lengths |
| 103 | precision: 9 |
Thomas Wouters | 1b7f891 | 2007-09-19 03:06:30 +0000 | [diff] [blame] | 104 | mulx080 multiply 0.1 123456789 -> 12345678.9 |
| 105 | mulx081 multiply 0.1 1234567891 -> 123456789 Inexact Rounded |
| 106 | mulx082 multiply 0.1 12345678912 -> 1.23456789E+9 Inexact Rounded |
| 107 | mulx083 multiply 0.1 12345678912345 -> 1.23456789E+12 Inexact Rounded |
| 108 | mulx084 multiply 0.1 123456789 -> 12345678.9 |
Raymond Hettinger | 7c85fa4 | 2004-07-01 11:01:35 +0000 | [diff] [blame] | 109 | precision: 8 |
Thomas Wouters | 1b7f891 | 2007-09-19 03:06:30 +0000 | [diff] [blame] | 110 | mulx085 multiply 0.1 12345678912 -> 1.2345679E+9 Inexact Rounded |
| 111 | mulx086 multiply 0.1 12345678912345 -> 1.2345679E+12 Inexact Rounded |
Raymond Hettinger | 7c85fa4 | 2004-07-01 11:01:35 +0000 | [diff] [blame] | 112 | precision: 7 |
Thomas Wouters | 1b7f891 | 2007-09-19 03:06:30 +0000 | [diff] [blame] | 113 | mulx087 multiply 0.1 12345678912 -> 1.234568E+9 Inexact Rounded |
| 114 | mulx088 multiply 0.1 12345678912345 -> 1.234568E+12 Inexact Rounded |
Raymond Hettinger | 7c85fa4 | 2004-07-01 11:01:35 +0000 | [diff] [blame] | 115 | |
| 116 | precision: 9 |
| 117 | mulx090 multiply 123456789 0.1 -> 12345678.9 |
| 118 | mulx091 multiply 1234567891 0.1 -> 123456789 Inexact Rounded |
| 119 | mulx092 multiply 12345678912 0.1 -> 1.23456789E+9 Inexact Rounded |
| 120 | mulx093 multiply 12345678912345 0.1 -> 1.23456789E+12 Inexact Rounded |
| 121 | mulx094 multiply 123456789 0.1 -> 12345678.9 |
| 122 | precision: 8 |
| 123 | mulx095 multiply 12345678912 0.1 -> 1.2345679E+9 Inexact Rounded |
| 124 | mulx096 multiply 12345678912345 0.1 -> 1.2345679E+12 Inexact Rounded |
| 125 | precision: 7 |
| 126 | mulx097 multiply 12345678912 0.1 -> 1.234568E+9 Inexact Rounded |
| 127 | mulx098 multiply 12345678912345 0.1 -> 1.234568E+12 Inexact Rounded |
| 128 | |
| 129 | -- test some more edge cases and carries |
| 130 | maxexponent: 9999 |
| 131 | minexponent: -9999 |
| 132 | precision: 33 |
| 133 | mulx101 multiply 9 9 -> 81 |
| 134 | mulx102 multiply 9 90 -> 810 |
| 135 | mulx103 multiply 9 900 -> 8100 |
| 136 | mulx104 multiply 9 9000 -> 81000 |
| 137 | mulx105 multiply 9 90000 -> 810000 |
| 138 | mulx106 multiply 9 900000 -> 8100000 |
| 139 | mulx107 multiply 9 9000000 -> 81000000 |
| 140 | mulx108 multiply 9 90000000 -> 810000000 |
| 141 | mulx109 multiply 9 900000000 -> 8100000000 |
| 142 | mulx110 multiply 9 9000000000 -> 81000000000 |
| 143 | mulx111 multiply 9 90000000000 -> 810000000000 |
| 144 | mulx112 multiply 9 900000000000 -> 8100000000000 |
| 145 | mulx113 multiply 9 9000000000000 -> 81000000000000 |
| 146 | mulx114 multiply 9 90000000000000 -> 810000000000000 |
| 147 | mulx115 multiply 9 900000000000000 -> 8100000000000000 |
| 148 | mulx116 multiply 9 9000000000000000 -> 81000000000000000 |
| 149 | mulx117 multiply 9 90000000000000000 -> 810000000000000000 |
| 150 | mulx118 multiply 9 900000000000000000 -> 8100000000000000000 |
| 151 | mulx119 multiply 9 9000000000000000000 -> 81000000000000000000 |
| 152 | mulx120 multiply 9 90000000000000000000 -> 810000000000000000000 |
| 153 | mulx121 multiply 9 900000000000000000000 -> 8100000000000000000000 |
| 154 | mulx122 multiply 9 9000000000000000000000 -> 81000000000000000000000 |
| 155 | mulx123 multiply 9 90000000000000000000000 -> 810000000000000000000000 |
| 156 | -- test some more edge cases without carries |
| 157 | mulx131 multiply 3 3 -> 9 |
| 158 | mulx132 multiply 3 30 -> 90 |
| 159 | mulx133 multiply 3 300 -> 900 |
| 160 | mulx134 multiply 3 3000 -> 9000 |
| 161 | mulx135 multiply 3 30000 -> 90000 |
| 162 | mulx136 multiply 3 300000 -> 900000 |
| 163 | mulx137 multiply 3 3000000 -> 9000000 |
| 164 | mulx138 multiply 3 30000000 -> 90000000 |
| 165 | mulx139 multiply 3 300000000 -> 900000000 |
| 166 | mulx140 multiply 3 3000000000 -> 9000000000 |
| 167 | mulx141 multiply 3 30000000000 -> 90000000000 |
| 168 | mulx142 multiply 3 300000000000 -> 900000000000 |
| 169 | mulx143 multiply 3 3000000000000 -> 9000000000000 |
| 170 | mulx144 multiply 3 30000000000000 -> 90000000000000 |
| 171 | mulx145 multiply 3 300000000000000 -> 900000000000000 |
| 172 | mulx146 multiply 3 3000000000000000 -> 9000000000000000 |
| 173 | mulx147 multiply 3 30000000000000000 -> 90000000000000000 |
| 174 | mulx148 multiply 3 300000000000000000 -> 900000000000000000 |
| 175 | mulx149 multiply 3 3000000000000000000 -> 9000000000000000000 |
| 176 | mulx150 multiply 3 30000000000000000000 -> 90000000000000000000 |
| 177 | mulx151 multiply 3 300000000000000000000 -> 900000000000000000000 |
| 178 | mulx152 multiply 3 3000000000000000000000 -> 9000000000000000000000 |
| 179 | mulx153 multiply 3 30000000000000000000000 -> 90000000000000000000000 |
| 180 | |
| 181 | maxexponent: 999999999 |
| 182 | minexponent: -999999999 |
| 183 | precision: 9 |
| 184 | -- test some cases that are close to exponent overflow/underflow |
| 185 | mulx170 multiply 1 9e999999999 -> 9E+999999999 |
| 186 | mulx171 multiply 1 9.9e999999999 -> 9.9E+999999999 |
| 187 | mulx172 multiply 1 9.99e999999999 -> 9.99E+999999999 |
| 188 | mulx173 multiply 9e999999999 1 -> 9E+999999999 |
| 189 | mulx174 multiply 9.9e999999999 1 -> 9.9E+999999999 |
| 190 | mulx176 multiply 9.99e999999999 1 -> 9.99E+999999999 |
| 191 | mulx177 multiply 1 9.99999999e999999999 -> 9.99999999E+999999999 |
| 192 | mulx178 multiply 9.99999999e999999999 1 -> 9.99999999E+999999999 |
| 193 | |
| 194 | mulx180 multiply 0.1 9e-999999998 -> 9E-999999999 |
| 195 | mulx181 multiply 0.1 99e-999999998 -> 9.9E-999999998 |
| 196 | mulx182 multiply 0.1 999e-999999998 -> 9.99E-999999997 |
| 197 | |
| 198 | mulx183 multiply 0.1 9e-999999998 -> 9E-999999999 |
| 199 | mulx184 multiply 0.1 99e-999999998 -> 9.9E-999999998 |
| 200 | mulx185 multiply 0.1 999e-999999998 -> 9.99E-999999997 |
| 201 | mulx186 multiply 0.1 999e-999999997 -> 9.99E-999999996 |
| 202 | mulx187 multiply 0.1 9999e-999999997 -> 9.999E-999999995 |
| 203 | mulx188 multiply 0.1 99999e-999999997 -> 9.9999E-999999994 |
| 204 | |
| 205 | mulx190 multiply 1 9e-999999998 -> 9E-999999998 |
| 206 | mulx191 multiply 1 99e-999999998 -> 9.9E-999999997 |
| 207 | mulx192 multiply 1 999e-999999998 -> 9.99E-999999996 |
| 208 | mulx193 multiply 9e-999999998 1 -> 9E-999999998 |
| 209 | mulx194 multiply 99e-999999998 1 -> 9.9E-999999997 |
| 210 | mulx195 multiply 999e-999999998 1 -> 9.99E-999999996 |
| 211 | |
| 212 | mulx196 multiply 1e-599999999 1e-400000000 -> 1E-999999999 |
| 213 | mulx197 multiply 1e-600000000 1e-399999999 -> 1E-999999999 |
| 214 | mulx198 multiply 1.2e-599999999 1.2e-400000000 -> 1.44E-999999999 |
| 215 | mulx199 multiply 1.2e-600000000 1.2e-399999999 -> 1.44E-999999999 |
| 216 | |
| 217 | mulx201 multiply 1e599999999 1e400000000 -> 1E+999999999 |
| 218 | mulx202 multiply 1e600000000 1e399999999 -> 1E+999999999 |
| 219 | mulx203 multiply 1.2e599999999 1.2e400000000 -> 1.44E+999999999 |
| 220 | mulx204 multiply 1.2e600000000 1.2e399999999 -> 1.44E+999999999 |
| 221 | |
| 222 | -- long operand triangle |
| 223 | precision: 33 |
| 224 | mulx246 multiply 30269.587755640502150977251770554 4.8046009735990873395936309640543 -> 145433.290801193369671916511992830 Inexact Rounded |
| 225 | precision: 32 |
| 226 | mulx247 multiply 30269.587755640502150977251770554 4.8046009735990873395936309640543 -> 145433.29080119336967191651199283 Inexact Rounded |
| 227 | precision: 31 |
| 228 | mulx248 multiply 30269.587755640502150977251770554 4.8046009735990873395936309640543 -> 145433.2908011933696719165119928 Inexact Rounded |
| 229 | precision: 30 |
| 230 | mulx249 multiply 30269.587755640502150977251770554 4.8046009735990873395936309640543 -> 145433.290801193369671916511993 Inexact Rounded |
| 231 | precision: 29 |
| 232 | mulx250 multiply 30269.587755640502150977251770554 4.8046009735990873395936309640543 -> 145433.29080119336967191651199 Inexact Rounded |
| 233 | precision: 28 |
| 234 | mulx251 multiply 30269.587755640502150977251770554 4.8046009735990873395936309640543 -> 145433.2908011933696719165120 Inexact Rounded |
| 235 | precision: 27 |
| 236 | mulx252 multiply 30269.587755640502150977251770554 4.8046009735990873395936309640543 -> 145433.290801193369671916512 Inexact Rounded |
| 237 | precision: 26 |
| 238 | mulx253 multiply 30269.587755640502150977251770554 4.8046009735990873395936309640543 -> 145433.29080119336967191651 Inexact Rounded |
| 239 | precision: 25 |
| 240 | mulx254 multiply 30269.587755640502150977251770554 4.8046009735990873395936309640543 -> 145433.2908011933696719165 Inexact Rounded |
| 241 | precision: 24 |
| 242 | mulx255 multiply 30269.587755640502150977251770554 4.8046009735990873395936309640543 -> 145433.290801193369671917 Inexact Rounded |
| 243 | precision: 23 |
| 244 | mulx256 multiply 30269.587755640502150977251770554 4.8046009735990873395936309640543 -> 145433.29080119336967192 Inexact Rounded |
| 245 | precision: 22 |
| 246 | mulx257 multiply 30269.587755640502150977251770554 4.8046009735990873395936309640543 -> 145433.2908011933696719 Inexact Rounded |
| 247 | precision: 21 |
| 248 | mulx258 multiply 30269.587755640502150977251770554 4.8046009735990873395936309640543 -> 145433.290801193369672 Inexact Rounded |
| 249 | precision: 20 |
| 250 | mulx259 multiply 30269.587755640502150977251770554 4.8046009735990873395936309640543 -> 145433.29080119336967 Inexact Rounded |
| 251 | precision: 19 |
| 252 | mulx260 multiply 30269.587755640502150977251770554 4.8046009735990873395936309640543 -> 145433.2908011933697 Inexact Rounded |
| 253 | precision: 18 |
| 254 | mulx261 multiply 30269.587755640502150977251770554 4.8046009735990873395936309640543 -> 145433.290801193370 Inexact Rounded |
| 255 | precision: 17 |
| 256 | mulx262 multiply 30269.587755640502150977251770554 4.8046009735990873395936309640543 -> 145433.29080119337 Inexact Rounded |
| 257 | precision: 16 |
| 258 | mulx263 multiply 30269.587755640502150977251770554 4.8046009735990873395936309640543 -> 145433.2908011934 Inexact Rounded |
| 259 | precision: 15 |
| 260 | mulx264 multiply 30269.587755640502150977251770554 4.8046009735990873395936309640543 -> 145433.290801193 Inexact Rounded |
| 261 | precision: 14 |
| 262 | mulx265 multiply 30269.587755640502150977251770554 4.8046009735990873395936309640543 -> 145433.29080119 Inexact Rounded |
| 263 | precision: 13 |
| 264 | mulx266 multiply 30269.587755640502150977251770554 4.8046009735990873395936309640543 -> 145433.2908012 Inexact Rounded |
| 265 | precision: 12 |
| 266 | mulx267 multiply 30269.587755640502150977251770554 4.8046009735990873395936309640543 -> 145433.290801 Inexact Rounded |
| 267 | precision: 11 |
| 268 | mulx268 multiply 30269.587755640502150977251770554 4.8046009735990873395936309640543 -> 145433.29080 Inexact Rounded |
| 269 | precision: 10 |
| 270 | mulx269 multiply 30269.587755640502150977251770554 4.8046009735990873395936309640543 -> 145433.2908 Inexact Rounded |
| 271 | precision: 9 |
| 272 | mulx270 multiply 30269.587755640502150977251770554 4.8046009735990873395936309640543 -> 145433.291 Inexact Rounded |
| 273 | precision: 8 |
| 274 | mulx271 multiply 30269.587755640502150977251770554 4.8046009735990873395936309640543 -> 145433.29 Inexact Rounded |
| 275 | precision: 7 |
| 276 | mulx272 multiply 30269.587755640502150977251770554 4.8046009735990873395936309640543 -> 145433.3 Inexact Rounded |
| 277 | precision: 6 |
| 278 | mulx273 multiply 30269.587755640502150977251770554 4.8046009735990873395936309640543 -> 145433 Inexact Rounded |
| 279 | precision: 5 |
| 280 | mulx274 multiply 30269.587755640502150977251770554 4.8046009735990873395936309640543 -> 1.4543E+5 Inexact Rounded |
| 281 | precision: 4 |
| 282 | mulx275 multiply 30269.587755640502150977251770554 4.8046009735990873395936309640543 -> 1.454E+5 Inexact Rounded |
| 283 | precision: 3 |
| 284 | mulx276 multiply 30269.587755640502150977251770554 4.8046009735990873395936309640543 -> 1.45E+5 Inexact Rounded |
| 285 | precision: 2 |
| 286 | mulx277 multiply 30269.587755640502150977251770554 4.8046009735990873395936309640543 -> 1.5E+5 Inexact Rounded |
| 287 | precision: 1 |
| 288 | mulx278 multiply 30269.587755640502150977251770554 4.8046009735990873395936309640543 -> 1E+5 Inexact Rounded |
| 289 | |
Thomas Wouters | 1b7f891 | 2007-09-19 03:06:30 +0000 | [diff] [blame] | 290 | -- test some edge cases with exact rounding |
| 291 | maxexponent: 9999 |
| 292 | minexponent: -9999 |
| 293 | precision: 9 |
| 294 | mulx301 multiply 9 9 -> 81 |
| 295 | mulx302 multiply 9 90 -> 810 |
| 296 | mulx303 multiply 9 900 -> 8100 |
| 297 | mulx304 multiply 9 9000 -> 81000 |
| 298 | mulx305 multiply 9 90000 -> 810000 |
| 299 | mulx306 multiply 9 900000 -> 8100000 |
| 300 | mulx307 multiply 9 9000000 -> 81000000 |
| 301 | mulx308 multiply 9 90000000 -> 810000000 |
| 302 | mulx309 multiply 9 900000000 -> 8.10000000E+9 Rounded |
| 303 | mulx310 multiply 9 9000000000 -> 8.10000000E+10 Rounded |
| 304 | mulx311 multiply 9 90000000000 -> 8.10000000E+11 Rounded |
| 305 | mulx312 multiply 9 900000000000 -> 8.10000000E+12 Rounded |
| 306 | mulx313 multiply 9 9000000000000 -> 8.10000000E+13 Rounded |
| 307 | mulx314 multiply 9 90000000000000 -> 8.10000000E+14 Rounded |
| 308 | mulx315 multiply 9 900000000000000 -> 8.10000000E+15 Rounded |
| 309 | mulx316 multiply 9 9000000000000000 -> 8.10000000E+16 Rounded |
| 310 | mulx317 multiply 9 90000000000000000 -> 8.10000000E+17 Rounded |
| 311 | mulx318 multiply 9 900000000000000000 -> 8.10000000E+18 Rounded |
| 312 | mulx319 multiply 9 9000000000000000000 -> 8.10000000E+19 Rounded |
| 313 | mulx320 multiply 9 90000000000000000000 -> 8.10000000E+20 Rounded |
| 314 | mulx321 multiply 9 900000000000000000000 -> 8.10000000E+21 Rounded |
| 315 | mulx322 multiply 9 9000000000000000000000 -> 8.10000000E+22 Rounded |
| 316 | mulx323 multiply 9 90000000000000000000000 -> 8.10000000E+23 Rounded |
| 317 | |
| 318 | -- fastpath breakers |
| 319 | precision: 29 |
| 320 | mulx330 multiply 1.491824697641270317824852952837224 1.105170918075647624811707826490246514675628614562883537345747603 -> 1.6487212707001281468486507878 Inexact Rounded |
| 321 | precision: 55 |
| 322 | mulx331 multiply 0.8958341352965282506768545828765117803873717284891040428 0.8958341352965282506768545828765117803873717284891040428 -> 0.8025187979624784829842553829934069955890983696752228299 Inexact Rounded |
| 323 | |
| 324 | |
Raymond Hettinger | 7c85fa4 | 2004-07-01 11:01:35 +0000 | [diff] [blame] | 325 | -- tryzeros cases |
| 326 | precision: 7 |
| 327 | rounding: half_up |
| 328 | maxExponent: 92 |
| 329 | minexponent: -92 |
| 330 | mulx504 multiply 0E-60 1000E-60 -> 0E-98 Clamped |
| 331 | mulx505 multiply 100E+60 0E+60 -> 0E+92 Clamped |
| 332 | |
| 333 | -- mixed with zeros |
| 334 | maxexponent: 999999999 |
| 335 | minexponent: -999999999 |
| 336 | precision: 9 |
| 337 | mulx541 multiply 0 -1 -> -0 |
| 338 | mulx542 multiply -0 -1 -> 0 |
| 339 | mulx543 multiply 0 1 -> 0 |
| 340 | mulx544 multiply -0 1 -> -0 |
| 341 | mulx545 multiply -1 0 -> -0 |
| 342 | mulx546 multiply -1 -0 -> 0 |
| 343 | mulx547 multiply 1 0 -> 0 |
| 344 | mulx548 multiply 1 -0 -> -0 |
| 345 | |
| 346 | mulx551 multiply 0.0 -1 -> -0.0 |
| 347 | mulx552 multiply -0.0 -1 -> 0.0 |
| 348 | mulx553 multiply 0.0 1 -> 0.0 |
| 349 | mulx554 multiply -0.0 1 -> -0.0 |
| 350 | mulx555 multiply -1.0 0 -> -0.0 |
| 351 | mulx556 multiply -1.0 -0 -> 0.0 |
| 352 | mulx557 multiply 1.0 0 -> 0.0 |
| 353 | mulx558 multiply 1.0 -0 -> -0.0 |
| 354 | |
| 355 | mulx561 multiply 0 -1.0 -> -0.0 |
| 356 | mulx562 multiply -0 -1.0 -> 0.0 |
| 357 | mulx563 multiply 0 1.0 -> 0.0 |
| 358 | mulx564 multiply -0 1.0 -> -0.0 |
| 359 | mulx565 multiply -1 0.0 -> -0.0 |
| 360 | mulx566 multiply -1 -0.0 -> 0.0 |
| 361 | mulx567 multiply 1 0.0 -> 0.0 |
| 362 | mulx568 multiply 1 -0.0 -> -0.0 |
| 363 | |
| 364 | mulx571 multiply 0.0 -1.0 -> -0.00 |
| 365 | mulx572 multiply -0.0 -1.0 -> 0.00 |
| 366 | mulx573 multiply 0.0 1.0 -> 0.00 |
| 367 | mulx574 multiply -0.0 1.0 -> -0.00 |
| 368 | mulx575 multiply -1.0 0.0 -> -0.00 |
| 369 | mulx576 multiply -1.0 -0.0 -> 0.00 |
| 370 | mulx577 multiply 1.0 0.0 -> 0.00 |
| 371 | mulx578 multiply 1.0 -0.0 -> -0.00 |
| 372 | |
| 373 | |
| 374 | -- Specials |
| 375 | mulx580 multiply Inf -Inf -> -Infinity |
| 376 | mulx581 multiply Inf -1000 -> -Infinity |
| 377 | mulx582 multiply Inf -1 -> -Infinity |
| 378 | mulx583 multiply Inf -0 -> NaN Invalid_operation |
| 379 | mulx584 multiply Inf 0 -> NaN Invalid_operation |
| 380 | mulx585 multiply Inf 1 -> Infinity |
| 381 | mulx586 multiply Inf 1000 -> Infinity |
| 382 | mulx587 multiply Inf Inf -> Infinity |
| 383 | mulx588 multiply -1000 Inf -> -Infinity |
| 384 | mulx589 multiply -Inf Inf -> -Infinity |
| 385 | mulx590 multiply -1 Inf -> -Infinity |
| 386 | mulx591 multiply -0 Inf -> NaN Invalid_operation |
| 387 | mulx592 multiply 0 Inf -> NaN Invalid_operation |
| 388 | mulx593 multiply 1 Inf -> Infinity |
| 389 | mulx594 multiply 1000 Inf -> Infinity |
| 390 | mulx595 multiply Inf Inf -> Infinity |
| 391 | |
| 392 | mulx600 multiply -Inf -Inf -> Infinity |
| 393 | mulx601 multiply -Inf -1000 -> Infinity |
| 394 | mulx602 multiply -Inf -1 -> Infinity |
| 395 | mulx603 multiply -Inf -0 -> NaN Invalid_operation |
| 396 | mulx604 multiply -Inf 0 -> NaN Invalid_operation |
| 397 | mulx605 multiply -Inf 1 -> -Infinity |
| 398 | mulx606 multiply -Inf 1000 -> -Infinity |
| 399 | mulx607 multiply -Inf Inf -> -Infinity |
| 400 | mulx608 multiply -1000 Inf -> -Infinity |
| 401 | mulx609 multiply -Inf -Inf -> Infinity |
| 402 | mulx610 multiply -1 -Inf -> Infinity |
| 403 | mulx611 multiply -0 -Inf -> NaN Invalid_operation |
| 404 | mulx612 multiply 0 -Inf -> NaN Invalid_operation |
| 405 | mulx613 multiply 1 -Inf -> -Infinity |
| 406 | mulx614 multiply 1000 -Inf -> -Infinity |
| 407 | mulx615 multiply Inf -Inf -> -Infinity |
| 408 | |
| 409 | mulx621 multiply NaN -Inf -> NaN |
| 410 | mulx622 multiply NaN -1000 -> NaN |
| 411 | mulx623 multiply NaN -1 -> NaN |
| 412 | mulx624 multiply NaN -0 -> NaN |
| 413 | mulx625 multiply NaN 0 -> NaN |
| 414 | mulx626 multiply NaN 1 -> NaN |
| 415 | mulx627 multiply NaN 1000 -> NaN |
| 416 | mulx628 multiply NaN Inf -> NaN |
| 417 | mulx629 multiply NaN NaN -> NaN |
| 418 | mulx630 multiply -Inf NaN -> NaN |
| 419 | mulx631 multiply -1000 NaN -> NaN |
| 420 | mulx632 multiply -1 NaN -> NaN |
| 421 | mulx633 multiply -0 NaN -> NaN |
| 422 | mulx634 multiply 0 NaN -> NaN |
| 423 | mulx635 multiply 1 NaN -> NaN |
| 424 | mulx636 multiply 1000 NaN -> NaN |
| 425 | mulx637 multiply Inf NaN -> NaN |
| 426 | |
| 427 | mulx641 multiply sNaN -Inf -> NaN Invalid_operation |
| 428 | mulx642 multiply sNaN -1000 -> NaN Invalid_operation |
| 429 | mulx643 multiply sNaN -1 -> NaN Invalid_operation |
| 430 | mulx644 multiply sNaN -0 -> NaN Invalid_operation |
| 431 | mulx645 multiply sNaN 0 -> NaN Invalid_operation |
| 432 | mulx646 multiply sNaN 1 -> NaN Invalid_operation |
| 433 | mulx647 multiply sNaN 1000 -> NaN Invalid_operation |
| 434 | mulx648 multiply sNaN NaN -> NaN Invalid_operation |
| 435 | mulx649 multiply sNaN sNaN -> NaN Invalid_operation |
| 436 | mulx650 multiply NaN sNaN -> NaN Invalid_operation |
| 437 | mulx651 multiply -Inf sNaN -> NaN Invalid_operation |
| 438 | mulx652 multiply -1000 sNaN -> NaN Invalid_operation |
| 439 | mulx653 multiply -1 sNaN -> NaN Invalid_operation |
| 440 | mulx654 multiply -0 sNaN -> NaN Invalid_operation |
| 441 | mulx655 multiply 0 sNaN -> NaN Invalid_operation |
| 442 | mulx656 multiply 1 sNaN -> NaN Invalid_operation |
| 443 | mulx657 multiply 1000 sNaN -> NaN Invalid_operation |
| 444 | mulx658 multiply Inf sNaN -> NaN Invalid_operation |
| 445 | mulx659 multiply NaN sNaN -> NaN Invalid_operation |
| 446 | |
| 447 | -- propagating NaNs |
| 448 | mulx661 multiply NaN9 -Inf -> NaN9 |
| 449 | mulx662 multiply NaN8 999 -> NaN8 |
| 450 | mulx663 multiply NaN71 Inf -> NaN71 |
| 451 | mulx664 multiply NaN6 NaN5 -> NaN6 |
| 452 | mulx665 multiply -Inf NaN4 -> NaN4 |
| 453 | mulx666 multiply -999 NaN33 -> NaN33 |
| 454 | mulx667 multiply Inf NaN2 -> NaN2 |
| 455 | |
| 456 | mulx671 multiply sNaN99 -Inf -> NaN99 Invalid_operation |
| 457 | mulx672 multiply sNaN98 -11 -> NaN98 Invalid_operation |
| 458 | mulx673 multiply sNaN97 NaN -> NaN97 Invalid_operation |
| 459 | mulx674 multiply sNaN16 sNaN94 -> NaN16 Invalid_operation |
| 460 | mulx675 multiply NaN95 sNaN93 -> NaN93 Invalid_operation |
| 461 | mulx676 multiply -Inf sNaN92 -> NaN92 Invalid_operation |
| 462 | mulx677 multiply 088 sNaN91 -> NaN91 Invalid_operation |
| 463 | mulx678 multiply Inf sNaN90 -> NaN90 Invalid_operation |
| 464 | mulx679 multiply NaN sNaN89 -> NaN89 Invalid_operation |
| 465 | |
| 466 | mulx681 multiply -NaN9 -Inf -> -NaN9 |
| 467 | mulx682 multiply -NaN8 999 -> -NaN8 |
| 468 | mulx683 multiply -NaN71 Inf -> -NaN71 |
| 469 | mulx684 multiply -NaN6 -NaN5 -> -NaN6 |
| 470 | mulx685 multiply -Inf -NaN4 -> -NaN4 |
| 471 | mulx686 multiply -999 -NaN33 -> -NaN33 |
| 472 | mulx687 multiply Inf -NaN2 -> -NaN2 |
| 473 | |
| 474 | mulx691 multiply -sNaN99 -Inf -> -NaN99 Invalid_operation |
| 475 | mulx692 multiply -sNaN98 -11 -> -NaN98 Invalid_operation |
| 476 | mulx693 multiply -sNaN97 NaN -> -NaN97 Invalid_operation |
| 477 | mulx694 multiply -sNaN16 -sNaN94 -> -NaN16 Invalid_operation |
| 478 | mulx695 multiply -NaN95 -sNaN93 -> -NaN93 Invalid_operation |
| 479 | mulx696 multiply -Inf -sNaN92 -> -NaN92 Invalid_operation |
| 480 | mulx697 multiply 088 -sNaN91 -> -NaN91 Invalid_operation |
| 481 | mulx698 multiply Inf -sNaN90 -> -NaN90 Invalid_operation |
| 482 | mulx699 multiply -NaN -sNaN89 -> -NaN89 Invalid_operation |
| 483 | |
| 484 | mulx701 multiply -NaN -Inf -> -NaN |
| 485 | mulx702 multiply -NaN 999 -> -NaN |
| 486 | mulx703 multiply -NaN Inf -> -NaN |
| 487 | mulx704 multiply -NaN -NaN -> -NaN |
| 488 | mulx705 multiply -Inf -NaN0 -> -NaN |
| 489 | mulx706 multiply -999 -NaN -> -NaN |
| 490 | mulx707 multiply Inf -NaN -> -NaN |
| 491 | |
| 492 | mulx711 multiply -sNaN -Inf -> -NaN Invalid_operation |
| 493 | mulx712 multiply -sNaN -11 -> -NaN Invalid_operation |
| 494 | mulx713 multiply -sNaN00 NaN -> -NaN Invalid_operation |
| 495 | mulx714 multiply -sNaN -sNaN -> -NaN Invalid_operation |
| 496 | mulx715 multiply -NaN -sNaN -> -NaN Invalid_operation |
| 497 | mulx716 multiply -Inf -sNaN -> -NaN Invalid_operation |
| 498 | mulx717 multiply 088 -sNaN -> -NaN Invalid_operation |
| 499 | mulx718 multiply Inf -sNaN -> -NaN Invalid_operation |
| 500 | mulx719 multiply -NaN -sNaN -> -NaN Invalid_operation |
| 501 | |
| 502 | -- overflow and underflow tests .. note subnormal results |
| 503 | maxexponent: 999999999 |
| 504 | minexponent: -999999999 |
| 505 | mulx730 multiply +1.23456789012345E-0 9E+999999999 -> Infinity Inexact Overflow Rounded |
| 506 | mulx731 multiply 9E+999999999 +1.23456789012345E-0 -> Infinity Inexact Overflow Rounded |
| 507 | mulx732 multiply +0.100 9E-999999999 -> 9.00E-1000000000 Subnormal |
| 508 | mulx733 multiply 9E-999999999 +0.100 -> 9.00E-1000000000 Subnormal |
| 509 | mulx735 multiply -1.23456789012345E-0 9E+999999999 -> -Infinity Inexact Overflow Rounded |
| 510 | mulx736 multiply 9E+999999999 -1.23456789012345E-0 -> -Infinity Inexact Overflow Rounded |
| 511 | mulx737 multiply -0.100 9E-999999999 -> -9.00E-1000000000 Subnormal |
| 512 | mulx738 multiply 9E-999999999 -0.100 -> -9.00E-1000000000 Subnormal |
| 513 | |
| 514 | mulx739 multiply 1e-599999999 1e-400000001 -> 1E-1000000000 Subnormal |
| 515 | mulx740 multiply 1e-599999999 1e-400000000 -> 1E-999999999 |
| 516 | mulx741 multiply 1e-600000000 1e-400000000 -> 1E-1000000000 Subnormal |
| 517 | mulx742 multiply 9e-999999998 0.01 -> 9E-1000000000 Subnormal |
| 518 | mulx743 multiply 9e-999999998 0.1 -> 9E-999999999 |
| 519 | mulx744 multiply 0.01 9e-999999998 -> 9E-1000000000 Subnormal |
| 520 | mulx745 multiply 1e599999999 1e400000001 -> Infinity Overflow Inexact Rounded |
| 521 | mulx746 multiply 1e599999999 1e400000000 -> 1E+999999999 |
| 522 | mulx747 multiply 1e600000000 1e400000000 -> Infinity Overflow Inexact Rounded |
| 523 | mulx748 multiply 9e999999998 100 -> Infinity Overflow Inexact Rounded |
| 524 | mulx749 multiply 9e999999998 10 -> 9.0E+999999999 |
| 525 | mulx750 multiply 100 9e999999998 -> Infinity Overflow Inexact Rounded |
| 526 | -- signs |
| 527 | mulx751 multiply 1e+777777777 1e+411111111 -> Infinity Overflow Inexact Rounded |
| 528 | mulx752 multiply 1e+777777777 -1e+411111111 -> -Infinity Overflow Inexact Rounded |
| 529 | mulx753 multiply -1e+777777777 1e+411111111 -> -Infinity Overflow Inexact Rounded |
| 530 | mulx754 multiply -1e+777777777 -1e+411111111 -> Infinity Overflow Inexact Rounded |
Thomas Wouters | 1b7f891 | 2007-09-19 03:06:30 +0000 | [diff] [blame] | 531 | mulx755 multiply 1e-777777777 1e-411111111 -> 0E-1000000007 Underflow Subnormal Inexact Rounded Clamped |
| 532 | mulx756 multiply 1e-777777777 -1e-411111111 -> -0E-1000000007 Underflow Subnormal Inexact Rounded Clamped |
| 533 | mulx757 multiply -1e-777777777 1e-411111111 -> -0E-1000000007 Underflow Subnormal Inexact Rounded Clamped |
| 534 | mulx758 multiply -1e-777777777 -1e-411111111 -> 0E-1000000007 Underflow Subnormal Inexact Rounded Clamped |
Raymond Hettinger | 7c85fa4 | 2004-07-01 11:01:35 +0000 | [diff] [blame] | 535 | |
| 536 | -- 'subnormal' boundary (all hard underflow or overflow in base arithemtic) |
| 537 | precision: 9 |
| 538 | mulx760 multiply 1e-600000000 1e-400000001 -> 1E-1000000001 Subnormal |
| 539 | mulx761 multiply 1e-600000000 1e-400000002 -> 1E-1000000002 Subnormal |
| 540 | mulx762 multiply 1e-600000000 1e-400000003 -> 1E-1000000003 Subnormal |
| 541 | mulx763 multiply 1e-600000000 1e-400000004 -> 1E-1000000004 Subnormal |
| 542 | mulx764 multiply 1e-600000000 1e-400000005 -> 1E-1000000005 Subnormal |
| 543 | mulx765 multiply 1e-600000000 1e-400000006 -> 1E-1000000006 Subnormal |
| 544 | mulx766 multiply 1e-600000000 1e-400000007 -> 1E-1000000007 Subnormal |
Thomas Wouters | 1b7f891 | 2007-09-19 03:06:30 +0000 | [diff] [blame] | 545 | mulx767 multiply 1e-600000000 1e-400000008 -> 0E-1000000007 Underflow Subnormal Inexact Rounded Clamped |
| 546 | mulx768 multiply 1e-600000000 1e-400000009 -> 0E-1000000007 Underflow Subnormal Inexact Rounded Clamped |
| 547 | mulx769 multiply 1e-600000000 1e-400000010 -> 0E-1000000007 Underflow Subnormal Inexact Rounded Clamped |
Raymond Hettinger | 7c85fa4 | 2004-07-01 11:01:35 +0000 | [diff] [blame] | 548 | -- [no equivalent of 'subnormal' for overflow] |
| 549 | mulx770 multiply 1e+600000000 1e+400000001 -> Infinity Overflow Inexact Rounded |
| 550 | mulx771 multiply 1e+600000000 1e+400000002 -> Infinity Overflow Inexact Rounded |
| 551 | mulx772 multiply 1e+600000000 1e+400000003 -> Infinity Overflow Inexact Rounded |
| 552 | mulx773 multiply 1e+600000000 1e+400000004 -> Infinity Overflow Inexact Rounded |
| 553 | mulx774 multiply 1e+600000000 1e+400000005 -> Infinity Overflow Inexact Rounded |
| 554 | mulx775 multiply 1e+600000000 1e+400000006 -> Infinity Overflow Inexact Rounded |
| 555 | mulx776 multiply 1e+600000000 1e+400000007 -> Infinity Overflow Inexact Rounded |
| 556 | mulx777 multiply 1e+600000000 1e+400000008 -> Infinity Overflow Inexact Rounded |
| 557 | mulx778 multiply 1e+600000000 1e+400000009 -> Infinity Overflow Inexact Rounded |
| 558 | mulx779 multiply 1e+600000000 1e+400000010 -> Infinity Overflow Inexact Rounded |
| 559 | |
| 560 | -- 'subnormal' test edge condition at higher precisions |
| 561 | precision: 99 |
| 562 | mulx780 multiply 1e-600000000 1e-400000007 -> 1E-1000000007 Subnormal |
| 563 | mulx781 multiply 1e-600000000 1e-400000008 -> 1E-1000000008 Subnormal |
| 564 | mulx782 multiply 1e-600000000 1e-400000097 -> 1E-1000000097 Subnormal |
Thomas Wouters | 1b7f891 | 2007-09-19 03:06:30 +0000 | [diff] [blame] | 565 | mulx783 multiply 1e-600000000 1e-400000098 -> 0E-1000000097 Underflow Subnormal Inexact Rounded Clamped |
Raymond Hettinger | 7c85fa4 | 2004-07-01 11:01:35 +0000 | [diff] [blame] | 566 | precision: 999 |
| 567 | mulx784 multiply 1e-600000000 1e-400000997 -> 1E-1000000997 Subnormal |
Thomas Wouters | 1b7f891 | 2007-09-19 03:06:30 +0000 | [diff] [blame] | 568 | mulx785 multiply 1e-600000000 1e-400000998 -> 0E-1000000997 Underflow Subnormal Inexact Rounded Clamped |
Raymond Hettinger | 7c85fa4 | 2004-07-01 11:01:35 +0000 | [diff] [blame] | 569 | |
| 570 | -- following testcases [through mulx800] not yet run against code |
| 571 | precision: 9999 |
| 572 | mulx786 multiply 1e-600000000 1e-400009997 -> 1E-1000009997 Subnormal |
Thomas Wouters | 1b7f891 | 2007-09-19 03:06:30 +0000 | [diff] [blame] | 573 | mulx787 multiply 1e-600000000 1e-400009998 -> 0E-1000009997 Underflow Subnormal Inexact Rounded Clamped |
Raymond Hettinger | 7c85fa4 | 2004-07-01 11:01:35 +0000 | [diff] [blame] | 574 | precision: 99999 |
| 575 | mulx788 multiply 1e-600000000 1e-400099997 -> 1E-1000099997 Subnormal |
Thomas Wouters | 1b7f891 | 2007-09-19 03:06:30 +0000 | [diff] [blame] | 576 | mulx789 multiply 1e-600000000 1e-400099998 -> 0E-1000099997 Underflow Subnormal Inexact Rounded Clamped |
Raymond Hettinger | 7c85fa4 | 2004-07-01 11:01:35 +0000 | [diff] [blame] | 577 | precision: 999999 |
| 578 | mulx790 multiply 1e-600000000 1e-400999997 -> 1E-1000999997 Subnormal |
Thomas Wouters | 1b7f891 | 2007-09-19 03:06:30 +0000 | [diff] [blame] | 579 | mulx791 multiply 1e-600000000 1e-400999998 -> 0E-1000999997 Underflow Subnormal Inexact Rounded Clamped |
Raymond Hettinger | 7c85fa4 | 2004-07-01 11:01:35 +0000 | [diff] [blame] | 580 | precision: 9999999 |
| 581 | mulx792 multiply 1e-600000000 1e-409999997 -> 1E-1009999997 Subnormal |
Thomas Wouters | 1b7f891 | 2007-09-19 03:06:30 +0000 | [diff] [blame] | 582 | mulx793 multiply 1e-600000000 1e-409999998 -> 0E-1009999997 Underflow Subnormal Inexact Rounded Clamped |
Raymond Hettinger | 7c85fa4 | 2004-07-01 11:01:35 +0000 | [diff] [blame] | 583 | precision: 99999999 |
| 584 | mulx794 multiply 1e-600000000 1e-499999997 -> 1E-1099999997 Subnormal |
Thomas Wouters | 1b7f891 | 2007-09-19 03:06:30 +0000 | [diff] [blame] | 585 | mulx795 multiply 1e-600000000 1e-499999998 -> 0E-1099999997 Underflow Subnormal Inexact Rounded Clamped |
Raymond Hettinger | 7c85fa4 | 2004-07-01 11:01:35 +0000 | [diff] [blame] | 586 | precision: 999999999 |
| 587 | mulx796 multiply 1e-999999999 1e-999999997 -> 1E-1999999996 Subnormal |
| 588 | mulx797 multiply 1e-999999999 1e-999999998 -> 1E-1999999997 Subnormal |
Thomas Wouters | 1b7f891 | 2007-09-19 03:06:30 +0000 | [diff] [blame] | 589 | mulx798 multiply 1e-999999999 1e-999999999 -> 0E-1999999997 Underflow Subnormal Inexact Rounded Clamped |
Raymond Hettinger | 7c85fa4 | 2004-07-01 11:01:35 +0000 | [diff] [blame] | 590 | mulx799 multiply 1e-600000000 1e-400000007 -> 1E-1000000007 Subnormal |
| 591 | mulx800 multiply 1e-600000000 1e-400000008 -> 1E-1000000008 Subnormal |
| 592 | |
| 593 | -- test subnormals rounding |
| 594 | precision: 5 |
| 595 | maxExponent: 999 |
| 596 | minexponent: -999 |
| 597 | rounding: half_even |
| 598 | |
| 599 | mulx801 multiply 1.0000E-999 1 -> 1.0000E-999 |
| 600 | mulx802 multiply 1.000E-999 1e-1 -> 1.000E-1000 Subnormal |
| 601 | mulx803 multiply 1.00E-999 1e-2 -> 1.00E-1001 Subnormal |
| 602 | mulx804 multiply 1.0E-999 1e-3 -> 1.0E-1002 Subnormal |
| 603 | mulx805 multiply 1.0E-999 1e-4 -> 1E-1003 Subnormal Rounded |
| 604 | mulx806 multiply 1.3E-999 1e-4 -> 1E-1003 Underflow Subnormal Inexact Rounded |
| 605 | mulx807 multiply 1.5E-999 1e-4 -> 2E-1003 Underflow Subnormal Inexact Rounded |
| 606 | mulx808 multiply 1.7E-999 1e-4 -> 2E-1003 Underflow Subnormal Inexact Rounded |
| 607 | mulx809 multiply 2.3E-999 1e-4 -> 2E-1003 Underflow Subnormal Inexact Rounded |
| 608 | mulx810 multiply 2.5E-999 1e-4 -> 2E-1003 Underflow Subnormal Inexact Rounded |
| 609 | mulx811 multiply 2.7E-999 1e-4 -> 3E-1003 Underflow Subnormal Inexact Rounded |
| 610 | mulx812 multiply 1.49E-999 1e-4 -> 1E-1003 Underflow Subnormal Inexact Rounded |
| 611 | mulx813 multiply 1.50E-999 1e-4 -> 2E-1003 Underflow Subnormal Inexact Rounded |
| 612 | mulx814 multiply 1.51E-999 1e-4 -> 2E-1003 Underflow Subnormal Inexact Rounded |
| 613 | mulx815 multiply 2.49E-999 1e-4 -> 2E-1003 Underflow Subnormal Inexact Rounded |
| 614 | mulx816 multiply 2.50E-999 1e-4 -> 2E-1003 Underflow Subnormal Inexact Rounded |
| 615 | mulx817 multiply 2.51E-999 1e-4 -> 3E-1003 Underflow Subnormal Inexact Rounded |
| 616 | |
| 617 | mulx818 multiply 1E-999 1e-4 -> 1E-1003 Subnormal |
Thomas Wouters | 1b7f891 | 2007-09-19 03:06:30 +0000 | [diff] [blame] | 618 | mulx819 multiply 3E-999 1e-5 -> 0E-1003 Underflow Subnormal Inexact Rounded Clamped |
| 619 | mulx820 multiply 5E-999 1e-5 -> 0E-1003 Underflow Subnormal Inexact Rounded Clamped |
Raymond Hettinger | 7c85fa4 | 2004-07-01 11:01:35 +0000 | [diff] [blame] | 620 | mulx821 multiply 7E-999 1e-5 -> 1E-1003 Underflow Subnormal Inexact Rounded |
| 621 | mulx822 multiply 9E-999 1e-5 -> 1E-1003 Underflow Subnormal Inexact Rounded |
| 622 | mulx823 multiply 9.9E-999 1e-5 -> 1E-1003 Underflow Subnormal Inexact Rounded |
| 623 | |
| 624 | mulx824 multiply 1E-999 -1e-4 -> -1E-1003 Subnormal |
Thomas Wouters | 1b7f891 | 2007-09-19 03:06:30 +0000 | [diff] [blame] | 625 | mulx825 multiply 3E-999 -1e-5 -> -0E-1003 Underflow Subnormal Inexact Rounded Clamped |
| 626 | mulx826 multiply -5E-999 1e-5 -> -0E-1003 Underflow Subnormal Inexact Rounded Clamped |
Raymond Hettinger | 7c85fa4 | 2004-07-01 11:01:35 +0000 | [diff] [blame] | 627 | mulx827 multiply 7E-999 -1e-5 -> -1E-1003 Underflow Subnormal Inexact Rounded |
| 628 | mulx828 multiply -9E-999 1e-5 -> -1E-1003 Underflow Subnormal Inexact Rounded |
| 629 | mulx829 multiply 9.9E-999 -1e-5 -> -1E-1003 Underflow Subnormal Inexact Rounded |
Thomas Wouters | 1b7f891 | 2007-09-19 03:06:30 +0000 | [diff] [blame] | 630 | mulx830 multiply 3.0E-999 -1e-5 -> -0E-1003 Underflow Subnormal Inexact Rounded Clamped |
Raymond Hettinger | 7c85fa4 | 2004-07-01 11:01:35 +0000 | [diff] [blame] | 631 | |
| 632 | mulx831 multiply 1.0E-501 1e-501 -> 1.0E-1002 Subnormal |
| 633 | mulx832 multiply 2.0E-501 2e-501 -> 4.0E-1002 Subnormal |
| 634 | mulx833 multiply 4.0E-501 4e-501 -> 1.60E-1001 Subnormal |
| 635 | mulx834 multiply 10.0E-501 10e-501 -> 1.000E-1000 Subnormal |
| 636 | mulx835 multiply 30.0E-501 30e-501 -> 9.000E-1000 Subnormal |
| 637 | mulx836 multiply 40.0E-501 40e-501 -> 1.6000E-999 |
| 638 | |
| 639 | -- squares |
Thomas Wouters | 1b7f891 | 2007-09-19 03:06:30 +0000 | [diff] [blame] | 640 | mulx840 multiply 1E-502 1e-502 -> 0E-1003 Underflow Subnormal Inexact Rounded Clamped |
Raymond Hettinger | 7c85fa4 | 2004-07-01 11:01:35 +0000 | [diff] [blame] | 641 | mulx841 multiply 1E-501 1e-501 -> 1E-1002 Subnormal |
| 642 | mulx842 multiply 2E-501 2e-501 -> 4E-1002 Subnormal |
| 643 | mulx843 multiply 4E-501 4e-501 -> 1.6E-1001 Subnormal |
| 644 | mulx844 multiply 10E-501 10e-501 -> 1.00E-1000 Subnormal |
| 645 | mulx845 multiply 30E-501 30e-501 -> 9.00E-1000 Subnormal |
| 646 | mulx846 multiply 40E-501 40e-501 -> 1.600E-999 |
| 647 | |
| 648 | -- cubes |
Thomas Wouters | 1b7f891 | 2007-09-19 03:06:30 +0000 | [diff] [blame] | 649 | mulx850 multiply 1E-670 1e-335 -> 0E-1003 Underflow Subnormal Inexact Rounded Clamped |
Raymond Hettinger | 7c85fa4 | 2004-07-01 11:01:35 +0000 | [diff] [blame] | 650 | mulx851 multiply 1E-668 1e-334 -> 1E-1002 Subnormal |
| 651 | mulx852 multiply 4E-668 2e-334 -> 8E-1002 Subnormal |
| 652 | mulx853 multiply 9E-668 3e-334 -> 2.7E-1001 Subnormal |
| 653 | mulx854 multiply 16E-668 4e-334 -> 6.4E-1001 Subnormal |
| 654 | mulx855 multiply 25E-668 5e-334 -> 1.25E-1000 Subnormal |
| 655 | mulx856 multiply 10E-668 100e-334 -> 1.000E-999 |
| 656 | |
Thomas Wouters | 1b7f891 | 2007-09-19 03:06:30 +0000 | [diff] [blame] | 657 | -- test derived from result of 0.099 ** 999 at 15 digits with unlimited exponent |
Raymond Hettinger | 7c85fa4 | 2004-07-01 11:01:35 +0000 | [diff] [blame] | 658 | precision: 19 |
| 659 | mulx860 multiply 6636851557994578716E-520 6636851557994578716E-520 -> 4.40477986028551E-1003 Underflow Subnormal Inexact Rounded |
| 660 | |
| 661 | -- Long operand overflow may be a different path |
| 662 | precision: 3 |
| 663 | maxExponent: 999999999 |
| 664 | minexponent: -999999999 |
| 665 | mulx870 multiply 1 9.999E+999999999 -> Infinity Inexact Overflow Rounded |
| 666 | mulx871 multiply 1 -9.999E+999999999 -> -Infinity Inexact Overflow Rounded |
| 667 | mulx872 multiply 9.999E+999999999 1 -> Infinity Inexact Overflow Rounded |
| 668 | mulx873 multiply -9.999E+999999999 1 -> -Infinity Inexact Overflow Rounded |
| 669 | |
| 670 | -- check for double-rounded subnormals |
| 671 | precision: 5 |
| 672 | maxexponent: 79 |
| 673 | minexponent: -79 |
Thomas Wouters | 1b7f891 | 2007-09-19 03:06:30 +0000 | [diff] [blame] | 674 | mulx881 multiply 1.2347E-40 1.2347E-40 -> 1.524E-80 Inexact Rounded Subnormal Underflow |
| 675 | mulx882 multiply 1.234E-40 1.234E-40 -> 1.523E-80 Inexact Rounded Subnormal Underflow |
| 676 | mulx883 multiply 1.23E-40 1.23E-40 -> 1.513E-80 Inexact Rounded Subnormal Underflow |
| 677 | mulx884 multiply 1.2E-40 1.2E-40 -> 1.44E-80 Subnormal |
| 678 | mulx885 multiply 1.2E-40 1.2E-41 -> 1.44E-81 Subnormal |
| 679 | mulx886 multiply 1.2E-40 1.2E-42 -> 1.4E-82 Subnormal Inexact Rounded Underflow |
| 680 | mulx887 multiply 1.2E-40 1.3E-42 -> 1.6E-82 Subnormal Inexact Rounded Underflow |
| 681 | mulx888 multiply 1.3E-40 1.3E-42 -> 1.7E-82 Subnormal Inexact Rounded Underflow |
| 682 | mulx889 multiply 1.3E-40 1.3E-43 -> 2E-83 Subnormal Inexact Rounded Underflow |
| 683 | mulx890 multiply 1.3E-41 1.3E-43 -> 0E-83 Clamped Subnormal Inexact Rounded Underflow |
Raymond Hettinger | 7c85fa4 | 2004-07-01 11:01:35 +0000 | [diff] [blame] | 684 | |
| 685 | mulx891 multiply 1.2345E-39 1.234E-40 -> 1.5234E-79 Inexact Rounded |
| 686 | mulx892 multiply 1.23456E-39 1.234E-40 -> 1.5234E-79 Inexact Rounded |
| 687 | mulx893 multiply 1.2345E-40 1.234E-40 -> 1.523E-80 Inexact Rounded Subnormal Underflow |
| 688 | mulx894 multiply 1.23456E-40 1.234E-40 -> 1.523E-80 Inexact Rounded Subnormal Underflow |
| 689 | mulx895 multiply 1.2345E-41 1.234E-40 -> 1.52E-81 Inexact Rounded Subnormal Underflow |
| 690 | mulx896 multiply 1.23456E-41 1.234E-40 -> 1.52E-81 Inexact Rounded Subnormal Underflow |
| 691 | |
Thomas Wouters | 1b7f891 | 2007-09-19 03:06:30 +0000 | [diff] [blame] | 692 | -- Now explore the case where we get a normal result with Underflow |
| 693 | precision: 16 |
| 694 | rounding: half_up |
| 695 | maxExponent: 384 |
| 696 | minExponent: -383 |
| 697 | |
| 698 | mulx900 multiply 0.3000000000E-191 0.3000000000E-191 -> 9.00000000000000E-384 Subnormal Rounded |
| 699 | mulx901 multiply 0.3000000001E-191 0.3000000001E-191 -> 9.00000000600000E-384 Underflow Inexact Subnormal Rounded |
| 700 | mulx902 multiply 9.999999999999999E-383 0.0999999999999 -> 9.99999999999000E-384 Underflow Inexact Subnormal Rounded |
| 701 | mulx903 multiply 9.999999999999999E-383 0.09999999999999 -> 9.99999999999900E-384 Underflow Inexact Subnormal Rounded |
| 702 | mulx904 multiply 9.999999999999999E-383 0.099999999999999 -> 9.99999999999990E-384 Underflow Inexact Subnormal Rounded |
| 703 | mulx905 multiply 9.999999999999999E-383 0.0999999999999999 -> 9.99999999999999E-384 Underflow Inexact Subnormal Rounded |
| 704 | -- prove operands are exact |
| 705 | mulx906 multiply 9.999999999999999E-383 1 -> 9.999999999999999E-383 |
| 706 | mulx907 multiply 1 0.09999999999999999 -> 0.09999999999999999 |
| 707 | -- the next rounds to Nmin |
| 708 | mulx908 multiply 9.999999999999999E-383 0.09999999999999999 -> 1.000000000000000E-383 Underflow Inexact Subnormal Rounded |
| 709 | mulx909 multiply 9.999999999999999E-383 0.099999999999999999 -> 1.000000000000000E-383 Underflow Inexact Subnormal Rounded |
| 710 | mulx910 multiply 9.999999999999999E-383 0.0999999999999999999 -> 1.000000000000000E-383 Underflow Inexact Subnormal Rounded |
| 711 | mulx911 multiply 9.999999999999999E-383 0.09999999999999999999 -> 1.000000000000000E-383 Underflow Inexact Subnormal Rounded |
| 712 | |
| 713 | |
| 714 | -- Examples from SQL proposal (Krishna Kulkarni) |
| 715 | precision: 34 |
| 716 | rounding: half_up |
| 717 | maxExponent: 6144 |
| 718 | minExponent: -6143 |
| 719 | mulx1001 multiply 130E-2 120E-2 -> 1.5600 |
| 720 | mulx1002 multiply 130E-2 12E-1 -> 1.560 |
| 721 | mulx1003 multiply 130E-2 1E0 -> 1.30 |
| 722 | mulx1004 multiply 1E2 1E4 -> 1E+6 |
| 723 | |
| 724 | -- payload decapitate |
| 725 | precision: 5 |
| 726 | mulx1010 multiply 11 -sNaN1234567890 -> -NaN67890 Invalid_operation |
| 727 | |
Raymond Hettinger | 7c85fa4 | 2004-07-01 11:01:35 +0000 | [diff] [blame] | 728 | -- Null tests |
Thomas Wouters | 1b7f891 | 2007-09-19 03:06:30 +0000 | [diff] [blame] | 729 | mulx990 multiply 10 # -> NaN Invalid_operation |
| 730 | mulx991 multiply # 10 -> NaN Invalid_operation |
Raymond Hettinger | 7c85fa4 | 2004-07-01 11:01:35 +0000 | [diff] [blame] | 731 | |