------------------------------------------------------------------------

-- subtract.decTest -- decimal subtraction                            --

-- Copyright (c) IBM Corporation, 1981, 2008.  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.58



extended:    1

precision:   9

rounding:    half_up

maxExponent: 384

minexponent: -383



-- [first group are 'quick confidence check']

subx001 subtract  0   0  -> '0'

subx002 subtract  1   1  -> '0'

subx003 subtract  1   2  -> '-1'

subx004 subtract  2   1  -> '1'

subx005 subtract  2   2  -> '0'

subx006 subtract  3   2  -> '1'

subx007 subtract  2   3  -> '-1'



subx011 subtract -0   0  -> '-0'

subx012 subtract -1   1  -> '-2'

subx013 subtract -1   2  -> '-3'

subx014 subtract -2   1  -> '-3'

subx015 subtract -2   2  -> '-4'

subx016 subtract -3   2  -> '-5'

subx017 subtract -2   3  -> '-5'



subx021 subtract  0  -0  -> '0'

subx022 subtract  1  -1  -> '2'

subx023 subtract  1  -2  -> '3'

subx024 subtract  2  -1  -> '3'

subx025 subtract  2  -2  -> '4'

subx026 subtract  3  -2  -> '5'

subx027 subtract  2  -3  -> '5'



subx030 subtract  11  1  -> 10

subx031 subtract  10  1  ->  9

subx032 subtract  9   1  ->  8

subx033 subtract  1   1  ->  0

subx034 subtract  0   1  -> -1

subx035 subtract -1   1  -> -2

subx036 subtract -9   1  -> -10

subx037 subtract -10  1  -> -11

subx038 subtract -11  1  -> -12



subx040 subtract '5.75' '3.3'  -> '2.45'

subx041 subtract '5'    '-3'   -> '8'

subx042 subtract '-5'   '-3'   -> '-2'

subx043 subtract '-7'   '2.5'  -> '-9.5'

subx044 subtract '0.7'  '0.3'  -> '0.4'

subx045 subtract '1.3'  '0.3'  -> '1.0'

subx046 subtract '1.25' '1.25' -> '0.00'



subx050 subtract '1.23456789'    '1.00000000' -> '0.23456789'

subx051 subtract '1.23456789'    '1.00000089' -> '0.23456700'

subx052 subtract '0.5555555559'    '0.0000000001' -> '0.555555556' Inexact Rounded

subx053 subtract '0.5555555559'    '0.0000000005' -> '0.555555555' Inexact Rounded

subx054 subtract '0.4444444444'    '0.1111111111' -> '0.333333333' Inexact Rounded

subx055 subtract '1.0000000000'    '0.00000001' -> '0.999999990' Rounded

subx056 subtract '0.4444444444999'    '0' -> '0.444444444' Inexact Rounded

subx057 subtract '0.4444444445000'    '0' -> '0.444444445' Inexact Rounded



subx060 subtract '70'    '10000e+9' -> '-1.00000000E+13' Inexact Rounded

subx061 subtract '700'    '10000e+9' -> '-1.00000000E+13' Inexact Rounded

subx062 subtract '7000'    '10000e+9' -> '-9.99999999E+12' Inexact Rounded

subx063 subtract '70000'    '10000e+9' -> '-9.99999993E+12' Rounded

subx064 subtract '700000'    '10000e+9' -> '-9.99999930E+12' Rounded

  -- symmetry:

subx065 subtract '10000e+9'    '70' -> '1.00000000E+13' Inexact Rounded

subx066 subtract '10000e+9'    '700' -> '1.00000000E+13' Inexact Rounded

subx067 subtract '10000e+9'    '7000' -> '9.99999999E+12' Inexact Rounded

subx068 subtract '10000e+9'    '70000' -> '9.99999993E+12' Rounded

subx069 subtract '10000e+9'    '700000' -> '9.99999930E+12' Rounded



  -- change precision

subx080 subtract '10000e+9'    '70000' -> '9.99999993E+12' Rounded

precision: 6

subx081 subtract '10000e+9'    '70000' -> '1.00000E+13' Inexact Rounded

precision: 9



  -- some of the next group are really constructor tests

subx090 subtract '00.0'    '0.0'  -> '0.0'

subx091 subtract '00.0'    '0.00' -> '0.00'

subx092 subtract '0.00'    '00.0' -> '0.00'

subx093 subtract '00.0'    '0.00' -> '0.00'

subx094 subtract '0.00'    '00.0' -> '0.00'

subx095 subtract '3'    '.3'   -> '2.7'

subx096 subtract '3.'   '.3'   -> '2.7'

subx097 subtract '3.0'  '.3'   -> '2.7'

subx098 subtract '3.00' '.3'   -> '2.70'

subx099 subtract '3'    '3'    -> '0'

subx100 subtract '3'    '+3'   -> '0'

subx101 subtract '3'    '-3'   -> '6'

subx102 subtract '3'    '0.3'  -> '2.7'

subx103 subtract '3.'   '0.3'  -> '2.7'

subx104 subtract '3.0'  '0.3'  -> '2.7'

subx105 subtract '3.00' '0.3'  -> '2.70'

subx106 subtract '3'    '3.0'  -> '0.0'

subx107 subtract '3'    '+3.0' -> '0.0'

subx108 subtract '3'    '-3.0' -> '6.0'



-- the above all from add; massaged and extended.  Now some new ones...

-- [particularly important for comparisons]

-- NB: -xE-8 below were non-exponents pre-ANSI X3-274, and -1E-7 or 0E-7

-- with input rounding.

subx120 subtract  '10.23456784'    '10.23456789'  -> '-5E-8'

subx121 subtract  '10.23456785'    '10.23456789'  -> '-4E-8'

subx122 subtract  '10.23456786'    '10.23456789'  -> '-3E-8'

subx123 subtract  '10.23456787'    '10.23456789'  -> '-2E-8'

subx124 subtract  '10.23456788'    '10.23456789'  -> '-1E-8'

subx125 subtract  '10.23456789'    '10.23456789'  -> '0E-8'

subx126 subtract  '10.23456790'    '10.23456789'  -> '1E-8'

subx127 subtract  '10.23456791'    '10.23456789'  -> '2E-8'

subx128 subtract  '10.23456792'    '10.23456789'  -> '3E-8'

subx129 subtract  '10.23456793'    '10.23456789'  -> '4E-8'

subx130 subtract  '10.23456794'    '10.23456789'  -> '5E-8'

subx131 subtract  '10.23456781'    '10.23456786'  -> '-5E-8'

subx132 subtract  '10.23456782'    '10.23456786'  -> '-4E-8'

subx133 subtract  '10.23456783'    '10.23456786'  -> '-3E-8'

subx134 subtract  '10.23456784'    '10.23456786'  -> '-2E-8'

subx135 subtract  '10.23456785'    '10.23456786'  -> '-1E-8'

subx136 subtract  '10.23456786'    '10.23456786'  -> '0E-8'

subx137 subtract  '10.23456787'    '10.23456786'  -> '1E-8'

subx138 subtract  '10.23456788'    '10.23456786'  -> '2E-8'

subx139 subtract  '10.23456789'    '10.23456786'  -> '3E-8'

subx140 subtract  '10.23456790'    '10.23456786'  -> '4E-8'

subx141 subtract  '10.23456791'    '10.23456786'  -> '5E-8'

subx142 subtract  '1'              '0.999999999'  -> '1E-9'

subx143 subtract  '0.999999999'    '1'            -> '-1E-9'

subx144 subtract  '-10.23456780'   '-10.23456786' -> '6E-8'

subx145 subtract  '-10.23456790'   '-10.23456786' -> '-4E-8'

subx146 subtract  '-10.23456791'   '-10.23456786' -> '-5E-8'



precision: 3

subx150 subtract '12345678900000' '9999999999999'  -> 2.35E+12 Inexact Rounded

subx151 subtract '9999999999999'  '12345678900000' -> -2.35E+12 Inexact Rounded

precision: 6

subx152 subtract '12345678900000' '9999999999999'  -> 2.34568E+12 Inexact Rounded

subx153 subtract '9999999999999'  '12345678900000' -> -2.34568E+12 Inexact Rounded

precision: 9

subx154 subtract '12345678900000' '9999999999999'  -> 2.34567890E+12 Inexact Rounded

subx155 subtract '9999999999999'  '12345678900000' -> -2.34567890E+12 Inexact Rounded

precision: 12

subx156 subtract '12345678900000' '9999999999999'  -> 2.34567890000E+12 Inexact Rounded

subx157 subtract '9999999999999'  '12345678900000' -> -2.34567890000E+12 Inexact Rounded

precision: 15

subx158 subtract '12345678900000' '9999999999999'  -> 2345678900001

subx159 subtract '9999999999999'  '12345678900000' -> -2345678900001

precision: 9



-- additional scaled arithmetic tests [0.97 problem]

subx160 subtract '0'     '.1'      -> '-0.1'

subx161 subtract '00'    '.97983'  -> '-0.97983'

subx162 subtract '0'     '.9'      -> '-0.9'

subx163 subtract '0'     '0.102'   -> '-0.102'

subx164 subtract '0'     '.4'      -> '-0.4'

subx165 subtract '0'     '.307'    -> '-0.307'

subx166 subtract '0'     '.43822'  -> '-0.43822'

subx167 subtract '0'     '.911'    -> '-0.911'

subx168 subtract '.0'    '.02'     -> '-0.02'

subx169 subtract '00'    '.392'    -> '-0.392'

subx170 subtract '0'     '.26'     -> '-0.26'

subx171 subtract '0'     '0.51'    -> '-0.51'

subx172 subtract '0'     '.2234'   -> '-0.2234'

subx173 subtract '0'     '.2'      -> '-0.2'

subx174 subtract '.0'    '.0008'   -> '-0.0008'

-- 0. on left

subx180 subtract '0.0'     '-.1'      -> '0.1'

subx181 subtract '0.00'    '-.97983'  -> '0.97983'

subx182 subtract '0.0'     '-.9'      -> '0.9'

subx183 subtract '0.0'     '-0.102'   -> '0.102'

subx184 subtract '0.0'     '-.4'      -> '0.4'

subx185 subtract '0.0'     '-.307'    -> '0.307'

subx186 subtract '0.0'     '-.43822'  -> '0.43822'

subx187 subtract '0.0'     '-.911'    -> '0.911'

subx188 subtract '0.0'     '-.02'     -> '0.02'

subx189 subtract '0.00'    '-.392'    -> '0.392'

subx190 subtract '0.0'     '-.26'     -> '0.26'

subx191 subtract '0.0'     '-0.51'    -> '0.51'

subx192 subtract '0.0'     '-.2234'   -> '0.2234'

subx193 subtract '0.0'     '-.2'      -> '0.2'

subx194 subtract '0.0'     '-.0008'   -> '0.0008'

-- negatives of same

subx200 subtract '0'     '-.1'      -> '0.1'

subx201 subtract '00'    '-.97983'  -> '0.97983'

subx202 subtract '0'     '-.9'      -> '0.9'

subx203 subtract '0'     '-0.102'   -> '0.102'

subx204 subtract '0'     '-.4'      -> '0.4'

subx205 subtract '0'     '-.307'    -> '0.307'

subx206 subtract '0'     '-.43822'  -> '0.43822'

subx207 subtract '0'     '-.911'    -> '0.911'

subx208 subtract '.0'    '-.02'     -> '0.02'

subx209 subtract '00'    '-.392'    -> '0.392'

subx210 subtract '0'     '-.26'     -> '0.26'

subx211 subtract '0'     '-0.51'    -> '0.51'

subx212 subtract '0'     '-.2234'   -> '0.2234'

subx213 subtract '0'     '-.2'      -> '0.2'

subx214 subtract '.0'    '-.0008'   -> '0.0008'



-- more fixed, LHS swaps [really the same as testcases under add]

subx220 subtract '-56267E-12' 0  -> '-5.6267E-8'

subx221 subtract '-56267E-11' 0  -> '-5.6267E-7'

subx222 subtract '-56267E-10' 0  -> '-0.0000056267'

subx223 subtract '-56267E-9'  0  -> '-0.000056267'

subx224 subtract '-56267E-8'  0  -> '-0.00056267'

subx225 subtract '-56267E-7'  0  -> '-0.0056267'

subx226 subtract '-56267E-6'  0  -> '-0.056267'

subx227 subtract '-56267E-5'  0  -> '-0.56267'

subx228 subtract '-56267E-2'  0  -> '-562.67'

subx229 subtract '-56267E-1'  0  -> '-5626.7'

subx230 subtract '-56267E-0'  0  -> '-56267'

-- symmetry ...

subx240 subtract 0 '-56267E-12'  -> '5.6267E-8'

subx241 subtract 0 '-56267E-11'  -> '5.6267E-7'

subx242 subtract 0 '-56267E-10'  -> '0.0000056267'

subx243 subtract 0 '-56267E-9'   -> '0.000056267'

subx244 subtract 0 '-56267E-8'   -> '0.00056267'

subx245 subtract 0 '-56267E-7'   -> '0.0056267'

subx246 subtract 0 '-56267E-6'   -> '0.056267'

subx247 subtract 0 '-56267E-5'   -> '0.56267'

subx248 subtract 0 '-56267E-2'   -> '562.67'

subx249 subtract 0 '-56267E-1'   -> '5626.7'

subx250 subtract 0 '-56267E-0'   -> '56267'



-- now some more from the 'new' add

precision: 9

subx301 subtract '1.23456789'  '1.00000000' -> '0.23456789'

subx302 subtract '1.23456789'  '1.00000011' -> '0.23456778'



subx311 subtract '0.4444444444'  '0.5555555555' -> '-0.111111111' Inexact Rounded

subx312 subtract '0.4444444440'  '0.5555555555' -> '-0.111111112' Inexact Rounded

subx313 subtract '0.4444444444'  '0.5555555550' -> '-0.111111111' Inexact Rounded

subx314 subtract '0.44444444449'    '0' -> '0.444444444' Inexact Rounded

subx315 subtract '0.444444444499'   '0' -> '0.444444444' Inexact Rounded

subx316 subtract '0.4444444444999'  '0' -> '0.444444444' Inexact Rounded

subx317 subtract '0.4444444445000'  '0' -> '0.444444445' Inexact Rounded

subx318 subtract '0.4444444445001'  '0' -> '0.444444445' Inexact Rounded

subx319 subtract '0.444444444501'   '0' -> '0.444444445' Inexact Rounded

subx320 subtract '0.44444444451'    '0' -> '0.444444445' Inexact Rounded



-- some carrying effects

subx321 subtract '0.9998'  '0.0000' -> '0.9998'

subx322 subtract '0.9998'  '0.0001' -> '0.9997'

subx323 subtract '0.9998'  '0.0002' -> '0.9996'

subx324 subtract '0.9998'  '0.0003' -> '0.9995'

subx325 subtract '0.9998'  '-0.0000' -> '0.9998'

subx326 subtract '0.9998'  '-0.0001' -> '0.9999'

subx327 subtract '0.9998'  '-0.0002' -> '1.0000'

subx328 subtract '0.9998'  '-0.0003' -> '1.0001'



subx330 subtract '70'  '10000e+9' -> '-1.00000000E+13' Inexact Rounded

subx331 subtract '700'  '10000e+9' -> '-1.00000000E+13' Inexact Rounded

subx332 subtract '7000'  '10000e+9' -> '-9.99999999E+12' Inexact Rounded

subx333 subtract '70000'  '10000e+9' -> '-9.99999993E+12' Rounded

subx334 subtract '700000'  '10000e+9' -> '-9.99999930E+12' Rounded

subx335 subtract '7000000'  '10000e+9' -> '-9.99999300E+12' Rounded

-- symmetry:

subx340 subtract '10000e+9'  '70' -> '1.00000000E+13' Inexact Rounded

subx341 subtract '10000e+9'  '700' -> '1.00000000E+13' Inexact Rounded

subx342 subtract '10000e+9'  '7000' -> '9.99999999E+12' Inexact Rounded

subx343 subtract '10000e+9'  '70000' -> '9.99999993E+12' Rounded

subx344 subtract '10000e+9'  '700000' -> '9.99999930E+12' Rounded

subx345 subtract '10000e+9'  '7000000' -> '9.99999300E+12' Rounded



-- same, higher precision

precision: 15

subx346 subtract '10000e+9'  '7'   -> '9999999999993'

subx347 subtract '10000e+9'  '70'   -> '9999999999930'

subx348 subtract '10000e+9'  '700'   -> '9999999999300'

subx349 subtract '10000e+9'  '7000'   -> '9999999993000'

subx350 subtract '10000e+9'  '70000'   -> '9999999930000'

subx351 subtract '10000e+9'  '700000'   -> '9999999300000'

subx352 subtract '7' '10000e+9'   -> '-9999999999993'

subx353 subtract '70' '10000e+9'   -> '-9999999999930'

subx354 subtract '700' '10000e+9'   -> '-9999999999300'

subx355 subtract '7000' '10000e+9'   -> '-9999999993000'

subx356 subtract '70000' '10000e+9'   -> '-9999999930000'

subx357 subtract '700000' '10000e+9'   -> '-9999999300000'



-- zero preservation

precision: 6

subx360 subtract '10000e+9'  '70000' -> '1.00000E+13' Inexact Rounded

subx361 subtract 1 '0.0001' -> '0.9999'

subx362 subtract 1 '0.00001' -> '0.99999'

subx363 subtract 1 '0.000001' -> '0.999999'

subx364 subtract 1 '0.0000001' -> '1.00000' Inexact Rounded

subx365 subtract 1 '0.00000001' -> '1.00000' Inexact Rounded



-- some funny zeros [in case of bad signum]

subx370 subtract 1  0  -> 1

subx371 subtract 1 0.  -> 1

subx372 subtract 1  .0 -> 1.0

subx373 subtract 1 0.0 -> 1.0

subx374 subtract  0  1 -> -1

subx375 subtract 0.  1 -> -1

subx376 subtract  .0 1 -> -1.0

subx377 subtract 0.0 1 -> -1.0



precision: 9



-- leading 0 digit before round

subx910 subtract -103519362 -51897955.3 -> -51621406.7

subx911 subtract 159579.444 89827.5229 -> 69751.9211



subx920 subtract 333.123456 33.1234566 -> 299.999999 Inexact Rounded

subx921 subtract 333.123456 33.1234565 -> 300.000000 Inexact Rounded

subx922 subtract 133.123456 33.1234565 ->  99.9999995

subx923 subtract 133.123456 33.1234564 ->  99.9999996

subx924 subtract 133.123456 33.1234540 -> 100.000002 Rounded

subx925 subtract 133.123456 43.1234560 ->  90.0000000

subx926 subtract 133.123456 43.1234561 ->  89.9999999

subx927 subtract 133.123456 43.1234566 ->  89.9999994

subx928 subtract 101.123456 91.1234566 ->   9.9999994

subx929 subtract 101.123456 99.1234566 ->   1.9999994



-- more of the same; probe for cluster boundary problems

precision: 1

subx930 subtract  11 2           -> 9

precision: 2

subx932 subtract 101 2           -> 99

precision: 3

subx934 subtract 101 2.1         -> 98.9

subx935 subtract 101 92.01       ->  8.99

precision: 4

subx936 subtract 101 2.01        -> 98.99

subx937 subtract 101 92.01       ->  8.99

subx938 subtract 101 92.006      ->  8.994

precision: 5

subx939 subtract 101 2.001       -> 98.999

subx940 subtract 101 92.001      ->  8.999

subx941 subtract 101 92.0006     ->  8.9994

precision: 6

subx942 subtract 101 2.0001      -> 98.9999

subx943 subtract 101 92.0001     ->  8.9999

subx944 subtract 101 92.00006    ->  8.99994

precision: 7

subx945 subtract 101 2.00001     -> 98.99999

subx946 subtract 101 92.00001    ->  8.99999

subx947 subtract 101 92.000006   ->  8.999994

precision: 8

subx948 subtract 101 2.000001    -> 98.999999

subx949 subtract 101 92.000001   ->  8.999999

subx950 subtract 101 92.0000006  ->  8.9999994

precision: 9

subx951 subtract 101 2.0000001   -> 98.9999999

subx952 subtract 101 92.0000001  ->  8.9999999

subx953 subtract 101 92.00000006 ->  8.99999994



precision: 9



-- more LHS swaps [were fixed]

subx390 subtract '-56267E-10'   0 ->  '-0.0000056267'

subx391 subtract '-56267E-6'    0 ->  '-0.056267'

subx392 subtract '-56267E-5'    0 ->  '-0.56267'

subx393 subtract '-56267E-4'    0 ->  '-5.6267'

subx394 subtract '-56267E-3'    0 ->  '-56.267'

subx395 subtract '-56267E-2'    0 ->  '-562.67'

subx396 subtract '-56267E-1'    0 ->  '-5626.7'

subx397 subtract '-56267E-0'    0 ->  '-56267'

subx398 subtract '-5E-10'       0 ->  '-5E-10'

subx399 subtract '-5E-7'        0 ->  '-5E-7'

subx400 subtract '-5E-6'        0 ->  '-0.000005'

subx401 subtract '-5E-5'        0 ->  '-0.00005'

subx402 subtract '-5E-4'        0 ->  '-0.0005'

subx403 subtract '-5E-1'        0 ->  '-0.5'

subx404 subtract '-5E0'         0 ->  '-5'

subx405 subtract '-5E1'         0 ->  '-50'

subx406 subtract '-5E5'         0 ->  '-500000'

subx407 subtract '-5E8'         0 ->  '-500000000'

subx408 subtract '-5E9'         0 ->  '-5.00000000E+9'   Rounded

subx409 subtract '-5E10'        0 ->  '-5.00000000E+10'  Rounded

subx410 subtract '-5E11'        0 ->  '-5.00000000E+11'  Rounded

subx411 subtract '-5E100'       0 ->  '-5.00000000E+100' Rounded



-- more RHS swaps [were fixed]

subx420 subtract 0  '-56267E-10' ->  '0.0000056267'

subx421 subtract 0  '-56267E-6'  ->  '0.056267'

subx422 subtract 0  '-56267E-5'  ->  '0.56267'

subx423 subtract 0  '-56267E-4'  ->  '5.6267'

subx424 subtract 0  '-56267E-3'  ->  '56.267'

subx425 subtract 0  '-56267E-2'  ->  '562.67'

subx426 subtract 0  '-56267E-1'  ->  '5626.7'

subx427 subtract 0  '-56267E-0'  ->  '56267'

subx428 subtract 0  '-5E-10'     ->  '5E-10'

subx429 subtract 0  '-5E-7'      ->  '5E-7'

subx430 subtract 0  '-5E-6'      ->  '0.000005'

subx431 subtract 0  '-5E-5'      ->  '0.00005'

subx432 subtract 0  '-5E-4'      ->  '0.0005'

subx433 subtract 0  '-5E-1'      ->  '0.5'

subx434 subtract 0  '-5E0'       ->  '5'

subx435 subtract 0  '-5E1'       ->  '50'

subx436 subtract 0  '-5E5'       ->  '500000'

subx437 subtract 0  '-5E8'       ->  '500000000'

subx438 subtract 0  '-5E9'       ->  '5.00000000E+9'    Rounded

subx439 subtract 0  '-5E10'      ->  '5.00000000E+10'   Rounded

subx440 subtract 0  '-5E11'      ->  '5.00000000E+11'   Rounded

subx441 subtract 0  '-5E100'     ->  '5.00000000E+100'  Rounded





-- try borderline precision, with carries, etc.

precision: 15

subx461 subtract '1E+12' '1'       -> '999999999999'

subx462 subtract '1E+12' '-1.11'   -> '1000000000001.11'

subx463 subtract '1.11'  '-1E+12'  -> '1000000000001.11'

subx464 subtract '-1'    '-1E+12'  -> '999999999999'

subx465 subtract '7E+12' '1'       -> '6999999999999'

subx466 subtract '7E+12' '-1.11'   -> '7000000000001.11'

subx467 subtract '1.11'  '-7E+12'  -> '7000000000001.11'

subx468 subtract '-1'    '-7E+12'  -> '6999999999999'



--                 123456789012345       123456789012345      1 23456789012345

subx470 subtract '0.444444444444444'  '-0.555555555555563' -> '1.00000000000001' Inexact Rounded

subx471 subtract '0.444444444444444'  '-0.555555555555562' -> '1.00000000000001' Inexact Rounded

subx472 subtract '0.444444444444444'  '-0.555555555555561' -> '1.00000000000001' Inexact Rounded

subx473 subtract '0.444444444444444'  '-0.555555555555560' -> '1.00000000000000' Inexact Rounded

subx474 subtract '0.444444444444444'  '-0.555555555555559' -> '1.00000000000000' Inexact Rounded

subx475 subtract '0.444444444444444'  '-0.555555555555558' -> '1.00000000000000' Inexact Rounded

subx476 subtract '0.444444444444444'  '-0.555555555555557' -> '1.00000000000000' Inexact Rounded

subx477 subtract '0.444444444444444'  '-0.555555555555556' -> '1.00000000000000' Rounded

subx478 subtract '0.444444444444444'  '-0.555555555555555' -> '0.999999999999999'

subx479 subtract '0.444444444444444'  '-0.555555555555554' -> '0.999999999999998'

subx480 subtract '0.444444444444444'  '-0.555555555555553' -> '0.999999999999997'

subx481 subtract '0.444444444444444'  '-0.555555555555552' -> '0.999999999999996'

subx482 subtract '0.444444444444444'  '-0.555555555555551' -> '0.999999999999995'

subx483 subtract '0.444444444444444'  '-0.555555555555550' -> '0.999999999999994'



-- and some more, including residue effects and different roundings

precision: 9

rounding: half_up

subx500 subtract '123456789' 0             -> '123456789'

subx501 subtract '123456789' 0.000000001   -> '123456789' Inexact Rounded

subx502 subtract '123456789' 0.000001      -> '123456789' Inexact Rounded

subx503 subtract '123456789' 0.1           -> '123456789' Inexact Rounded

subx504 subtract '123456789' 0.4           -> '123456789' Inexact Rounded

subx505 subtract '123456789' 0.49          -> '123456789' Inexact Rounded

subx506 subtract '123456789' 0.499999      -> '123456789' Inexact Rounded

subx507 subtract '123456789' 0.499999999   -> '123456789' Inexact Rounded

subx508 subtract '123456789' 0.5           -> '123456789' Inexact Rounded

subx509 subtract '123456789' 0.500000001   -> '123456788' Inexact Rounded

subx510 subtract '123456789' 0.500001      -> '123456788' Inexact Rounded

subx511 subtract '123456789' 0.51          -> '123456788' Inexact Rounded

subx512 subtract '123456789' 0.6           -> '123456788' Inexact Rounded

subx513 subtract '123456789' 0.9           -> '123456788' Inexact Rounded

subx514 subtract '123456789' 0.99999       -> '123456788' Inexact Rounded

subx515 subtract '123456789' 0.999999999   -> '123456788' Inexact Rounded

subx516 subtract '123456789' 1             -> '123456788'

subx517 subtract '123456789' 1.000000001   -> '123456788' Inexact Rounded

subx518 subtract '123456789' 1.00001       -> '123456788' Inexact Rounded

subx519 subtract '123456789' 1.1           -> '123456788' Inexact Rounded



rounding: half_even

subx520 subtract '123456789' 0             -> '123456789'

subx521 subtract '123456789' 0.000000001   -> '123456789' Inexact Rounded

subx522 subtract '123456789' 0.000001      -> '123456789' Inexact Rounded

subx523 subtract '123456789' 0.1           -> '123456789' Inexact Rounded

subx524 subtract '123456789' 0.4           -> '123456789' Inexact Rounded

subx525 subtract '123456789' 0.49          -> '123456789' Inexact Rounded

subx526 subtract '123456789' 0.499999      -> '123456789' Inexact Rounded

subx527 subtract '123456789' 0.499999999   -> '123456789' Inexact Rounded

subx528 subtract '123456789' 0.5           -> '123456788' Inexact Rounded

subx529 subtract '123456789' 0.500000001   -> '123456788' Inexact Rounded

subx530 subtract '123456789' 0.500001      -> '123456788' Inexact Rounded

subx531 subtract '123456789' 0.51          -> '123456788' Inexact Rounded

subx532 subtract '123456789' 0.6           -> '123456788' Inexact Rounded

subx533 subtract '123456789' 0.9           -> '123456788' Inexact Rounded

subx534 subtract '123456789' 0.99999       -> '123456788' Inexact Rounded

subx535 subtract '123456789' 0.999999999   -> '123456788' Inexact Rounded

subx536 subtract '123456789' 1             -> '123456788'

subx537 subtract '123456789' 1.00000001    -> '123456788' Inexact Rounded

subx538 subtract '123456789' 1.00001       -> '123456788' Inexact Rounded

subx539 subtract '123456789' 1.1           -> '123456788' Inexact Rounded

-- critical few with even bottom digit...

subx540 subtract '123456788' 0.499999999   -> '123456788' Inexact Rounded

subx541 subtract '123456788' 0.5           -> '123456788' Inexact Rounded

subx542 subtract '123456788' 0.500000001   -> '123456787' Inexact Rounded



rounding: down

subx550 subtract '123456789' 0             -> '123456789'

subx551 subtract '123456789' 0.000000001   -> '123456788' Inexact Rounded

subx552 subtract '123456789' 0.000001      -> '123456788' Inexact Rounded

subx553 subtract '123456789' 0.1           -> '123456788' Inexact Rounded

subx554 subtract '123456789' 0.4           -> '123456788' Inexact Rounded

subx555 subtract '123456789' 0.49          -> '123456788' Inexact Rounded

subx556 subtract '123456789' 0.499999      -> '123456788' Inexact Rounded

subx557 subtract '123456789' 0.499999999   -> '123456788' Inexact Rounded

subx558 subtract '123456789' 0.5           -> '123456788' Inexact Rounded

subx559 subtract '123456789' 0.500000001   -> '123456788' Inexact Rounded

subx560 subtract '123456789' 0.500001      -> '123456788' Inexact Rounded

subx561 subtract '123456789' 0.51          -> '123456788' Inexact Rounded

subx562 subtract '123456789' 0.6           -> '123456788' Inexact Rounded

subx563 subtract '123456789' 0.9           -> '123456788' Inexact Rounded

subx564 subtract '123456789' 0.99999       -> '123456788' Inexact Rounded

subx565 subtract '123456789' 0.999999999   -> '123456788' Inexact Rounded

subx566 subtract '123456789' 1             -> '123456788'

subx567 subtract '123456789' 1.00000001    -> '123456787' Inexact Rounded

subx568 subtract '123456789' 1.00001       -> '123456787' Inexact Rounded

subx569 subtract '123456789' 1.1           -> '123456787' Inexact Rounded



-- symmetry...

rounding: half_up

subx600 subtract 0             '123456789' -> '-123456789'

subx601 subtract 0.000000001   '123456789' -> '-123456789' Inexact Rounded

subx602 subtract 0.000001      '123456789' -> '-123456789' Inexact Rounded

subx603 subtract 0.1           '123456789' -> '-123456789' Inexact Rounded

subx604 subtract 0.4           '123456789' -> '-123456789' Inexact Rounded

subx605 subtract 0.49          '123456789' -> '-123456789' Inexact Rounded

subx606 subtract 0.499999      '123456789' -> '-123456789' Inexact Rounded

subx607 subtract 0.499999999   '123456789' -> '-123456789' Inexact Rounded

subx608 subtract 0.5           '123456789' -> '-123456789' Inexact Rounded

subx609 subtract 0.500000001   '123456789' -> '-123456788' Inexact Rounded

subx610 subtract 0.500001      '123456789' -> '-123456788' Inexact Rounded

subx611 subtract 0.51          '123456789' -> '-123456788' Inexact Rounded

subx612 subtract 0.6           '123456789' -> '-123456788' Inexact Rounded

subx613 subtract 0.9           '123456789' -> '-123456788' Inexact Rounded

subx614 subtract 0.99999       '123456789' -> '-123456788' Inexact Rounded

subx615 subtract 0.999999999   '123456789' -> '-123456788' Inexact Rounded

subx616 subtract 1             '123456789' -> '-123456788'

subx617 subtract 1.000000001   '123456789' -> '-123456788' Inexact Rounded

subx618 subtract 1.00001       '123456789' -> '-123456788' Inexact Rounded

subx619 subtract 1.1           '123456789' -> '-123456788' Inexact Rounded



rounding: half_even

subx620 subtract 0             '123456789' -> '-123456789'

subx621 subtract 0.000000001   '123456789' -> '-123456789' Inexact Rounded

subx622 subtract 0.000001      '123456789' -> '-123456789' Inexact Rounded

subx623 subtract 0.1           '123456789' -> '-123456789' Inexact Rounded

subx624 subtract 0.4           '123456789' -> '-123456789' Inexact Rounded

subx625 subtract 0.49          '123456789' -> '-123456789' Inexact Rounded

subx626 subtract 0.499999      '123456789' -> '-123456789' Inexact Rounded

subx627 subtract 0.499999999   '123456789' -> '-123456789' Inexact Rounded

subx628 subtract 0.5           '123456789' -> '-123456788' Inexact Rounded

subx629 subtract 0.500000001   '123456789' -> '-123456788' Inexact Rounded

subx630 subtract 0.500001      '123456789' -> '-123456788' Inexact Rounded

subx631 subtract 0.51          '123456789' -> '-123456788' Inexact Rounded

subx632 subtract 0.6           '123456789' -> '-123456788' Inexact Rounded

subx633 subtract 0.9           '123456789' -> '-123456788' Inexact Rounded

subx634 subtract 0.99999       '123456789' -> '-123456788' Inexact Rounded

subx635 subtract 0.999999999   '123456789' -> '-123456788' Inexact Rounded

subx636 subtract 1             '123456789' -> '-123456788'

subx637 subtract 1.00000001    '123456789' -> '-123456788' Inexact Rounded

subx638 subtract 1.00001       '123456789' -> '-123456788' Inexact Rounded

subx639 subtract 1.1           '123456789' -> '-123456788' Inexact Rounded

-- critical few with even bottom digit...

subx640 subtract 0.499999999   '123456788' -> '-123456788' Inexact Rounded

subx641 subtract 0.5           '123456788' -> '-123456788' Inexact Rounded

subx642 subtract 0.500000001   '123456788' -> '-123456787' Inexact Rounded



rounding: down

subx650 subtract 0             '123456789' -> '-123456789'

subx651 subtract 0.000000001   '123456789' -> '-123456788' Inexact Rounded

subx652 subtract 0.000001      '123456789' -> '-123456788' Inexact Rounded

subx653 subtract 0.1           '123456789' -> '-123456788' Inexact Rounded

subx654 subtract 0.4           '123456789' -> '-123456788' Inexact Rounded

subx655 subtract 0.49          '123456789' -> '-123456788' Inexact Rounded

subx656 subtract 0.499999      '123456789' -> '-123456788' Inexact Rounded

subx657 subtract 0.499999999   '123456789' -> '-123456788' Inexact Rounded

subx658 subtract 0.5           '123456789' -> '-123456788' Inexact Rounded

subx659 subtract 0.500000001   '123456789' -> '-123456788' Inexact Rounded

subx660 subtract 0.500001      '123456789' -> '-123456788' Inexact Rounded

subx661 subtract 0.51          '123456789' -> '-123456788' Inexact Rounded

subx662 subtract 0.6           '123456789' -> '-123456788' Inexact Rounded

subx663 subtract 0.9           '123456789' -> '-123456788' Inexact Rounded

subx664 subtract 0.99999       '123456789' -> '-123456788' Inexact Rounded

subx665 subtract 0.999999999   '123456789' -> '-123456788' Inexact Rounded

subx666 subtract 1             '123456789' -> '-123456788'

subx667 subtract 1.00000001    '123456789' -> '-123456787' Inexact Rounded

subx668 subtract 1.00001       '123456789' -> '-123456787' Inexact Rounded

subx669 subtract 1.1           '123456789' -> '-123456787' Inexact Rounded





-- lots of leading zeros in intermediate result, and showing effects of

-- input rounding would have affected the following

precision: 9

rounding: half_up

subx670 subtract '123456789' '123456788.1' -> 0.9

subx671 subtract '123456789' '123456788.9' -> 0.1

subx672 subtract '123456789' '123456789.1' -> -0.1

subx673 subtract '123456789' '123456789.5' -> -0.5

subx674 subtract '123456789' '123456789.9' -> -0.9



rounding: half_even

subx680 subtract '123456789' '123456788.1' -> 0.9

subx681 subtract '123456789' '123456788.9' -> 0.1

subx682 subtract '123456789' '123456789.1' -> -0.1

subx683 subtract '123456789' '123456789.5' -> -0.5

subx684 subtract '123456789' '123456789.9' -> -0.9



subx685 subtract '123456788' '123456787.1' -> 0.9

subx686 subtract '123456788' '123456787.9' -> 0.1

subx687 subtract '123456788' '123456788.1' -> -0.1

subx688 subtract '123456788' '123456788.5' -> -0.5

subx689 subtract '123456788' '123456788.9' -> -0.9



rounding: down

subx690 subtract '123456789' '123456788.1' -> 0.9

subx691 subtract '123456789' '123456788.9' -> 0.1

subx692 subtract '123456789' '123456789.1' -> -0.1

subx693 subtract '123456789' '123456789.5' -> -0.5

subx694 subtract '123456789' '123456789.9' -> -0.9



-- input preparation tests

rounding: half_up

precision: 3



subx700 subtract '12345678900000'  -9999999999999 ->  '2.23E+13' Inexact Rounded

subx701 subtract  '9999999999999' -12345678900000 ->  '2.23E+13' Inexact Rounded

subx702 subtract '12E+3'  '-3456' ->  '1.55E+4' Inexact Rounded

subx703 subtract '12E+3'  '-3446' ->  '1.54E+4' Inexact Rounded

subx704 subtract '12E+3'  '-3454' ->  '1.55E+4' Inexact Rounded

subx705 subtract '12E+3'  '-3444' ->  '1.54E+4' Inexact Rounded



subx706 subtract '3456'  '-12E+3' ->  '1.55E+4' Inexact Rounded

subx707 subtract '3446'  '-12E+3' ->  '1.54E+4' Inexact Rounded

subx708 subtract '3454'  '-12E+3' ->  '1.55E+4' Inexact Rounded

subx709 subtract '3444'  '-12E+3' ->  '1.54E+4' Inexact Rounded



-- overflow and underflow tests [subnormals now possible]

maxexponent: 999999999

minexponent: -999999999

precision: 9

rounding: down

subx710 subtract 1E+999999999    -9E+999999999   -> 9.99999999E+999999999 Overflow Inexact Rounded

subx711 subtract 9E+999999999    -1E+999999999   -> 9.99999999E+999999999 Overflow Inexact Rounded

rounding: half_up

subx712 subtract 1E+999999999    -9E+999999999   -> Infinity Overflow Inexact Rounded

subx713 subtract 9E+999999999    -1E+999999999   -> Infinity Overflow Inexact Rounded

subx714 subtract -1.1E-999999999 -1E-999999999   -> -1E-1000000000 Subnormal

subx715 subtract 1E-999999999    +1.1e-999999999 -> -1E-1000000000 Subnormal

subx716 subtract -1E+999999999   +9E+999999999   -> -Infinity Overflow Inexact Rounded

subx717 subtract -9E+999999999   +1E+999999999   -> -Infinity Overflow Inexact Rounded

subx718 subtract +1.1E-999999999 +1E-999999999   -> 1E-1000000000 Subnormal

subx719 subtract -1E-999999999   -1.1e-999999999 -> 1E-1000000000 Subnormal



precision: 3

subx720 subtract 1  9.999E+999999999   -> -Infinity Inexact Overflow Rounded

subx721 subtract 1 -9.999E+999999999   ->  Infinity Inexact Overflow Rounded

subx722 subtract    9.999E+999999999 1 ->  Infinity Inexact Overflow Rounded

subx723 subtract   -9.999E+999999999 1 -> -Infinity Inexact Overflow Rounded

subx724 subtract 1  9.999E+999999999   -> -Infinity Inexact Overflow Rounded

subx725 subtract 1 -9.999E+999999999   ->  Infinity Inexact Overflow Rounded

subx726 subtract    9.999E+999999999 1 ->  Infinity Inexact Overflow Rounded

subx727 subtract   -9.999E+999999999 1 -> -Infinity Inexact Overflow Rounded



-- [more below]



-- long operand checks

maxexponent: 999

minexponent: -999

precision: 9

sub731 subtract 12345678000 0 ->  1.23456780E+10 Rounded

sub732 subtract 0 12345678000 -> -1.23456780E+10 Rounded

sub733 subtract 1234567800  0 ->  1.23456780E+9 Rounded

sub734 subtract 0 1234567800  -> -1.23456780E+9 Rounded

sub735 subtract 1234567890  0 ->  1.23456789E+9 Rounded

sub736 subtract 0 1234567890  -> -1.23456789E+9 Rounded

sub737 subtract 1234567891  0 ->  1.23456789E+9 Inexact Rounded

sub738 subtract 0 1234567891  -> -1.23456789E+9 Inexact Rounded

sub739 subtract 12345678901 0 ->  1.23456789E+10 Inexact Rounded

sub740 subtract 0 12345678901 -> -1.23456789E+10 Inexact Rounded

sub741 subtract 1234567896  0 ->  1.23456790E+9 Inexact Rounded

sub742 subtract 0 1234567896  -> -1.23456790E+9 Inexact Rounded



precision: 15

sub751 subtract 12345678000 0 ->  12345678000

sub752 subtract 0 12345678000 -> -12345678000

sub753 subtract 1234567800  0 ->  1234567800

sub754 subtract 0 1234567800  -> -1234567800

sub755 subtract 1234567890  0 ->  1234567890

sub756 subtract 0 1234567890  -> -1234567890

sub757 subtract 1234567891  0 ->  1234567891

sub758 subtract 0 1234567891  -> -1234567891

sub759 subtract 12345678901 0 ->  12345678901

sub760 subtract 0 12345678901 -> -12345678901

sub761 subtract 1234567896  0 ->  1234567896

sub762 subtract 0 1234567896  -> -1234567896



-- Specials

subx780 subtract -Inf   Inf   -> -Infinity

subx781 subtract -Inf   1000  -> -Infinity

subx782 subtract -Inf   1     -> -Infinity

subx783 subtract -Inf  -0     -> -Infinity

subx784 subtract -Inf  -1     -> -Infinity

subx785 subtract -Inf  -1000  -> -Infinity

subx787 subtract -1000  Inf   -> -Infinity

subx788 subtract -Inf   Inf   -> -Infinity

subx789 subtract -1     Inf   -> -Infinity

subx790 subtract  0     Inf   -> -Infinity

subx791 subtract  1     Inf   -> -Infinity

subx792 subtract  1000  Inf   -> -Infinity



subx800 subtract  Inf   Inf   ->  NaN  Invalid_operation

subx801 subtract  Inf   1000  ->  Infinity

subx802 subtract  Inf   1     ->  Infinity

subx803 subtract  Inf   0     ->  Infinity

subx804 subtract  Inf  -0     ->  Infinity

subx805 subtract  Inf  -1     ->  Infinity

subx806 subtract  Inf  -1000  ->  Infinity

subx807 subtract  Inf  -Inf   ->  Infinity

subx808 subtract -1000 -Inf   ->  Infinity

subx809 subtract -Inf  -Inf   ->  NaN  Invalid_operation

subx810 subtract -1    -Inf   ->  Infinity

subx811 subtract -0    -Inf   ->  Infinity

subx812 subtract  0    -Inf   ->  Infinity

subx813 subtract  1    -Inf   ->  Infinity

subx814 subtract  1000 -Inf   ->  Infinity

subx815 subtract  Inf  -Inf   ->  Infinity



subx821 subtract  NaN   Inf   ->  NaN

subx822 subtract -NaN   1000  -> -NaN

subx823 subtract  NaN   1     ->  NaN

subx824 subtract  NaN   0     ->  NaN

subx825 subtract  NaN  -0     ->  NaN

subx826 subtract  NaN  -1     ->  NaN

subx827 subtract  NaN  -1000  ->  NaN

subx828 subtract  NaN  -Inf   ->  NaN

subx829 subtract -NaN   NaN   -> -NaN

subx830 subtract -Inf   NaN   ->  NaN

subx831 subtract -1000  NaN   ->  NaN

subx832 subtract -1     NaN   ->  NaN

subx833 subtract -0     NaN   ->  NaN

subx834 subtract  0     NaN   ->  NaN

subx835 subtract  1     NaN   ->  NaN

subx836 subtract  1000 -NaN   -> -NaN

subx837 subtract  Inf   NaN   ->  NaN



subx841 subtract  sNaN  Inf   ->  NaN  Invalid_operation

subx842 subtract -sNaN  1000  -> -NaN  Invalid_operation

subx843 subtract  sNaN  1     ->  NaN  Invalid_operation

subx844 subtract  sNaN  0     ->  NaN  Invalid_operation

subx845 subtract  sNaN -0     ->  NaN  Invalid_operation

subx846 subtract  sNaN -1     ->  NaN  Invalid_operation

subx847 subtract  sNaN -1000  ->  NaN  Invalid_operation

subx848 subtract  sNaN  NaN   ->  NaN  Invalid_operation

subx849 subtract  sNaN sNaN   ->  NaN  Invalid_operation

subx850 subtract  NaN  sNaN   ->  NaN  Invalid_operation

subx851 subtract -Inf -sNaN   -> -NaN  Invalid_operation

subx852 subtract -1000 sNaN   ->  NaN  Invalid_operation

subx853 subtract -1    sNaN   ->  NaN  Invalid_operation

subx854 subtract -0    sNaN   ->  NaN  Invalid_operation

subx855 subtract  0    sNaN   ->  NaN  Invalid_operation

subx856 subtract  1    sNaN   ->  NaN  Invalid_operation

subx857 subtract  1000 sNaN   ->  NaN  Invalid_operation

subx858 subtract  Inf  sNaN   ->  NaN  Invalid_operation

subx859 subtract  NaN  sNaN   ->  NaN  Invalid_operation



-- propagating NaNs

subx861 subtract  NaN01   -Inf     ->  NaN1

subx862 subtract -NaN02   -1000    -> -NaN2

subx863 subtract  NaN03    1000    ->  NaN3

subx864 subtract  NaN04    Inf     ->  NaN4

subx865 subtract  NaN05    NaN61   ->  NaN5

subx866 subtract -Inf     -NaN71   -> -NaN71

subx867 subtract -1000     NaN81   ->  NaN81

subx868 subtract  1000     NaN91   ->  NaN91

subx869 subtract  Inf      NaN101  ->  NaN101

subx871 subtract  sNaN011  -Inf    ->  NaN11  Invalid_operation

subx872 subtract  sNaN012  -1000   ->  NaN12  Invalid_operation

subx873 subtract -sNaN013   1000   -> -NaN13  Invalid_operation

subx874 subtract  sNaN014   NaN171 ->  NaN14  Invalid_operation

subx875 subtract  sNaN015  sNaN181 ->  NaN15  Invalid_operation

subx876 subtract  NaN016   sNaN191 ->  NaN191 Invalid_operation

subx877 subtract -Inf      sNaN201 ->  NaN201 Invalid_operation

subx878 subtract -1000     sNaN211 ->  NaN211 Invalid_operation

subx879 subtract  1000    -sNaN221 -> -NaN221 Invalid_operation

subx880 subtract  Inf      sNaN231 ->  NaN231 Invalid_operation

subx881 subtract  NaN025   sNaN241 ->  NaN241 Invalid_operation



-- edge case spills

subx901 subtract  2.E-3  1.002  -> -1.000

subx902 subtract  2.0E-3  1.002  -> -1.0000

subx903 subtract  2.00E-3  1.0020  -> -1.00000

subx904 subtract  2.000E-3  1.00200  -> -1.000000

subx905 subtract  2.0000E-3  1.002000  -> -1.0000000

subx906 subtract  2.00000E-3  1.0020000  -> -1.00000000

subx907 subtract  2.000000E-3  1.00200000  -> -1.000000000

subx908 subtract  2.0000000E-3  1.002000000  -> -1.0000000000



-- subnormals and underflows

precision: 3

maxexponent: 999

minexponent: -999

subx1010 subtract  0  1.00E-999       ->  -1.00E-999

subx1011 subtract  0  0.1E-999        ->  -1E-1000   Subnormal

subx1012 subtract  0  0.10E-999       ->  -1.0E-1000 Subnormal

subx1013 subtract  0  0.100E-999      ->  -1.0E-1000 Subnormal Rounded

subx1014 subtract  0  0.01E-999       ->  -1E-1001   Subnormal

-- next is rounded to Emin

subx1015 subtract  0  0.999E-999      ->  -1.00E-999 Inexact Rounded Subnormal Underflow

subx1016 subtract  0  0.099E-999      ->  -1.0E-1000 Inexact Rounded Subnormal Underflow

subx1017 subtract  0  0.009E-999      ->  -1E-1001   Inexact Rounded Subnormal Underflow

subx1018 subtract  0  0.001E-999      ->  -0E-1001   Inexact Rounded Subnormal Underflow Clamped

subx1019 subtract  0  0.0009E-999     ->  -0E-1001   Inexact Rounded Subnormal Underflow Clamped

subx1020 subtract  0  0.0001E-999     ->  -0E-1001   Inexact Rounded Subnormal Underflow Clamped



subx1030 subtract  0 -1.00E-999       ->   1.00E-999

subx1031 subtract  0 -0.1E-999        ->   1E-1000   Subnormal

subx1032 subtract  0 -0.10E-999       ->   1.0E-1000 Subnormal

subx1033 subtract  0 -0.100E-999      ->   1.0E-1000 Subnormal Rounded

subx1034 subtract  0 -0.01E-999       ->   1E-1001   Subnormal

-- next is rounded to Emin

subx1035 subtract  0 -0.999E-999      ->   1.00E-999 Inexact Rounded Subnormal Underflow

subx1036 subtract  0 -0.099E-999      ->   1.0E-1000 Inexact Rounded Subnormal Underflow

subx1037 subtract  0 -0.009E-999      ->   1E-1001   Inexact Rounded Subnormal Underflow

subx1038 subtract  0 -0.001E-999      ->   0E-1001   Inexact Rounded Subnormal Underflow Clamped

subx1039 subtract  0 -0.0009E-999     ->   0E-1001   Inexact Rounded Subnormal Underflow Clamped

subx1040 subtract  0 -0.0001E-999     ->   0E-1001   Inexact Rounded Subnormal Underflow Clamped



-- some non-zero subnormal subtracts

-- subx1056 is a tricky case

rounding: half_up

subx1050 subtract  1.00E-999   0.1E-999  ->   9.0E-1000  Subnormal

subx1051 subtract  0.1E-999    0.1E-999  ->   0E-1000

subx1052 subtract  0.10E-999   0.1E-999  ->   0E-1001

subx1053 subtract  0.100E-999  0.1E-999  ->   0E-1001    Clamped

subx1054 subtract  0.01E-999   0.1E-999  ->   -9E-1001   Subnormal

subx1055 subtract  0.999E-999  0.1E-999  ->   9.0E-1000  Inexact Rounded Subnormal Underflow

subx1056 subtract  0.099E-999  0.1E-999  ->   -0E-1001   Inexact Rounded Subnormal Underflow Clamped

subx1057 subtract  0.009E-999  0.1E-999  ->   -9E-1001   Inexact Rounded Subnormal Underflow

subx1058 subtract  0.001E-999  0.1E-999  ->   -1.0E-1000 Inexact Rounded Subnormal Underflow

subx1059 subtract  0.0009E-999 0.1E-999  ->   -1.0E-1000 Inexact Rounded Subnormal Underflow

subx1060 subtract  0.0001E-999 0.1E-999  ->   -1.0E-1000 Inexact Rounded Subnormal Underflow





-- check for double-rounded subnormals

precision:   5

maxexponent: 79

minexponent: -79

subx1101 subtract  0 1.52444E-80 -> -1.524E-80 Inexact Rounded Subnormal Underflow

subx1102 subtract  0 1.52445E-80 -> -1.524E-80 Inexact Rounded Subnormal Underflow

subx1103 subtract  0 1.52446E-80 -> -1.524E-80 Inexact Rounded Subnormal Underflow

subx1104 subtract  1.52444E-80 0 ->  1.524E-80 Inexact Rounded Subnormal Underflow

subx1105 subtract  1.52445E-80 0 ->  1.524E-80 Inexact Rounded Subnormal Underflow

subx1106 subtract  1.52446E-80 0 ->  1.524E-80 Inexact Rounded Subnormal Underflow



subx1111 subtract  1.2345678E-80  1.2345671E-80 ->  0E-83 Inexact Rounded Subnormal Underflow Clamped

subx1112 subtract  1.2345678E-80  1.2345618E-80 ->  0E-83 Inexact Rounded Subnormal Underflow Clamped

subx1113 subtract  1.2345678E-80  1.2345178E-80 ->  0E-83 Inexact Rounded Subnormal Underflow Clamped

subx1114 subtract  1.2345678E-80  1.2341678E-80 ->  0E-83 Inexact Rounded Subnormal Underflow Clamped

subx1115 subtract  1.2345678E-80  1.2315678E-80 ->  3E-83         Rounded Subnormal

subx1116 subtract  1.2345678E-80  1.2145678E-80 ->  2.0E-82       Rounded Subnormal

subx1117 subtract  1.2345678E-80  1.1345678E-80 ->  1.00E-81      Rounded Subnormal

subx1118 subtract  1.2345678E-80  0.2345678E-80 ->  1.000E-80     Rounded Subnormal



precision:   34

rounding:    half_up

maxExponent: 6144

minExponent: -6143

-- Examples from SQL proposal (Krishna Kulkarni)

subx1125  subtract 130E-2  120E-2 -> 0.10

subx1126  subtract 130E-2  12E-1  -> 0.10

subx1127  subtract 130E-2  1E0    -> 0.30

subx1128  subtract 1E2     1E4    -> -9.9E+3



-- Null tests

subx9990 subtract 10  # -> NaN Invalid_operation

subx9991 subtract  # 10 -> NaN Invalid_operation

