1------------------------------------------------------------------------
2-- subtract.decTest -- decimal subtraction                            --
3-- Copyright (c) Mike Cowlishaw,  1981, 2010.  All rights reserved.   --
4-- Parts copyright (c) IBM Corporation, 1981, 2008.                   --
5------------------------------------------------------------------------
6-- Please see the document "General Decimal Arithmetic Testcases"     --
7-- at http://speleotrove.com/decimal for the description of           --
8-- these testcases.                                                   --
9--                                                                    --
10-- These testcases are experimental ('beta' versions), and they       --
11-- may contain errors.  They are offered on an as-is basis.  In       --
12-- particular, achieving the same results as the tests here is not    --
13-- a guarantee that an implementation complies with any Standard      --
14-- or specification.  The tests are not exhaustive.                   --
15--                                                                    --
16-- Please send comments, suggestions, and corrections to the author:  --
17--   Mike Cowlishaw, mfc@speleotrove.com                              --
18------------------------------------------------------------------------
19version: 2.62
20
21extended:    1
22precision:   9
23rounding:    half_up
24maxExponent: 384
25minexponent: -383
26
27-- [first group are 'quick confidence check']
28subx001 subtract  0   0  -> '0'
29subx002 subtract  1   1  -> '0'
30subx003 subtract  1   2  -> '-1'
31subx004 subtract  2   1  -> '1'
32subx005 subtract  2   2  -> '0'
33subx006 subtract  3   2  -> '1'
34subx007 subtract  2   3  -> '-1'
35
36subx011 subtract -0   0  -> '-0'
37subx012 subtract -1   1  -> '-2'
38subx013 subtract -1   2  -> '-3'
39subx014 subtract -2   1  -> '-3'
40subx015 subtract -2   2  -> '-4'
41subx016 subtract -3   2  -> '-5'
42subx017 subtract -2   3  -> '-5'
43
44subx021 subtract  0  -0  -> '0'
45subx022 subtract  1  -1  -> '2'
46subx023 subtract  1  -2  -> '3'
47subx024 subtract  2  -1  -> '3'
48subx025 subtract  2  -2  -> '4'
49subx026 subtract  3  -2  -> '5'
50subx027 subtract  2  -3  -> '5'
51
52subx030 subtract  11  1  -> 10
53subx031 subtract  10  1  ->  9
54subx032 subtract  9   1  ->  8
55subx033 subtract  1   1  ->  0
56subx034 subtract  0   1  -> -1
57subx035 subtract -1   1  -> -2
58subx036 subtract -9   1  -> -10
59subx037 subtract -10  1  -> -11
60subx038 subtract -11  1  -> -12
61
62subx040 subtract '5.75' '3.3'  -> '2.45'
63subx041 subtract '5'    '-3'   -> '8'
64subx042 subtract '-5'   '-3'   -> '-2'
65subx043 subtract '-7'   '2.5'  -> '-9.5'
66subx044 subtract '0.7'  '0.3'  -> '0.4'
67subx045 subtract '1.3'  '0.3'  -> '1.0'
68subx046 subtract '1.25' '1.25' -> '0.00'
69
70subx050 subtract '1.23456789'    '1.00000000' -> '0.23456789'
71subx051 subtract '1.23456789'    '1.00000089' -> '0.23456700'
72subx052 subtract '0.5555555559'    '0.0000000001' -> '0.555555556' Inexact Rounded
73subx053 subtract '0.5555555559'    '0.0000000005' -> '0.555555555' Inexact Rounded
74subx054 subtract '0.4444444444'    '0.1111111111' -> '0.333333333' Inexact Rounded
75subx055 subtract '1.0000000000'    '0.00000001' -> '0.999999990' Rounded
76subx056 subtract '0.4444444444999'    '0' -> '0.444444444' Inexact Rounded
77subx057 subtract '0.4444444445000'    '0' -> '0.444444445' Inexact Rounded
78
79subx060 subtract '70'    '10000e+9' -> '-1.00000000E+13' Inexact Rounded
80subx061 subtract '700'    '10000e+9' -> '-1.00000000E+13' Inexact Rounded
81subx062 subtract '7000'    '10000e+9' -> '-9.99999999E+12' Inexact Rounded
82subx063 subtract '70000'    '10000e+9' -> '-9.99999993E+12' Rounded
83subx064 subtract '700000'    '10000e+9' -> '-9.99999930E+12' Rounded
84  -- symmetry:
85subx065 subtract '10000e+9'    '70' -> '1.00000000E+13' Inexact Rounded
86subx066 subtract '10000e+9'    '700' -> '1.00000000E+13' Inexact Rounded
87subx067 subtract '10000e+9'    '7000' -> '9.99999999E+12' Inexact Rounded
88subx068 subtract '10000e+9'    '70000' -> '9.99999993E+12' Rounded
89subx069 subtract '10000e+9'    '700000' -> '9.99999930E+12' Rounded
90
91  -- change precision
92subx080 subtract '10000e+9'    '70000' -> '9.99999993E+12' Rounded
93precision: 6
94subx081 subtract '10000e+9'    '70000' -> '1.00000E+13' Inexact Rounded
95precision: 9
96
97  -- some of the next group are really constructor tests
98subx090 subtract '00.0'    '0.0'  -> '0.0'
99subx091 subtract '00.0'    '0.00' -> '0.00'
100subx092 subtract '0.00'    '00.0' -> '0.00'
101subx093 subtract '00.0'    '0.00' -> '0.00'
102subx094 subtract '0.00'    '00.0' -> '0.00'
103subx095 subtract '3'    '.3'   -> '2.7'
104subx096 subtract '3.'   '.3'   -> '2.7'
105subx097 subtract '3.0'  '.3'   -> '2.7'
106subx098 subtract '3.00' '.3'   -> '2.70'
107subx099 subtract '3'    '3'    -> '0'
108subx100 subtract '3'    '+3'   -> '0'
109subx101 subtract '3'    '-3'   -> '6'
110subx102 subtract '3'    '0.3'  -> '2.7'
111subx103 subtract '3.'   '0.3'  -> '2.7'
112subx104 subtract '3.0'  '0.3'  -> '2.7'
113subx105 subtract '3.00' '0.3'  -> '2.70'
114subx106 subtract '3'    '3.0'  -> '0.0'
115subx107 subtract '3'    '+3.0' -> '0.0'
116subx108 subtract '3'    '-3.0' -> '6.0'
117
118-- the above all from add; massaged and extended.  Now some new ones...
119-- [particularly important for comparisons]
120-- NB: -xE-8 below were non-exponents pre-ANSI X3-274, and -1E-7 or 0E-7
121-- with input rounding.
122subx120 subtract  '10.23456784'    '10.23456789'  -> '-5E-8'
123subx121 subtract  '10.23456785'    '10.23456789'  -> '-4E-8'
124subx122 subtract  '10.23456786'    '10.23456789'  -> '-3E-8'
125subx123 subtract  '10.23456787'    '10.23456789'  -> '-2E-8'
126subx124 subtract  '10.23456788'    '10.23456789'  -> '-1E-8'
127subx125 subtract  '10.23456789'    '10.23456789'  -> '0E-8'
128subx126 subtract  '10.23456790'    '10.23456789'  -> '1E-8'
129subx127 subtract  '10.23456791'    '10.23456789'  -> '2E-8'
130subx128 subtract  '10.23456792'    '10.23456789'  -> '3E-8'
131subx129 subtract  '10.23456793'    '10.23456789'  -> '4E-8'
132subx130 subtract  '10.23456794'    '10.23456789'  -> '5E-8'
133subx131 subtract  '10.23456781'    '10.23456786'  -> '-5E-8'
134subx132 subtract  '10.23456782'    '10.23456786'  -> '-4E-8'
135subx133 subtract  '10.23456783'    '10.23456786'  -> '-3E-8'
136subx134 subtract  '10.23456784'    '10.23456786'  -> '-2E-8'
137subx135 subtract  '10.23456785'    '10.23456786'  -> '-1E-8'
138subx136 subtract  '10.23456786'    '10.23456786'  -> '0E-8'
139subx137 subtract  '10.23456787'    '10.23456786'  -> '1E-8'
140subx138 subtract  '10.23456788'    '10.23456786'  -> '2E-8'
141subx139 subtract  '10.23456789'    '10.23456786'  -> '3E-8'
142subx140 subtract  '10.23456790'    '10.23456786'  -> '4E-8'
143subx141 subtract  '10.23456791'    '10.23456786'  -> '5E-8'
144subx142 subtract  '1'              '0.999999999'  -> '1E-9'
145subx143 subtract  '0.999999999'    '1'            -> '-1E-9'
146subx144 subtract  '-10.23456780'   '-10.23456786' -> '6E-8'
147subx145 subtract  '-10.23456790'   '-10.23456786' -> '-4E-8'
148subx146 subtract  '-10.23456791'   '-10.23456786' -> '-5E-8'
149
150precision: 3
151subx150 subtract '12345678900000' '9999999999999'  -> 2.35E+12 Inexact Rounded
152subx151 subtract '9999999999999'  '12345678900000' -> -2.35E+12 Inexact Rounded
153precision: 6
154subx152 subtract '12345678900000' '9999999999999'  -> 2.34568E+12 Inexact Rounded
155subx153 subtract '9999999999999'  '12345678900000' -> -2.34568E+12 Inexact Rounded
156precision: 9
157subx154 subtract '12345678900000' '9999999999999'  -> 2.34567890E+12 Inexact Rounded
158subx155 subtract '9999999999999'  '12345678900000' -> -2.34567890E+12 Inexact Rounded
159precision: 12
160subx156 subtract '12345678900000' '9999999999999'  -> 2.34567890000E+12 Inexact Rounded
161subx157 subtract '9999999999999'  '12345678900000' -> -2.34567890000E+12 Inexact Rounded
162precision: 15
163subx158 subtract '12345678900000' '9999999999999'  -> 2345678900001
164subx159 subtract '9999999999999'  '12345678900000' -> -2345678900001
165precision: 9
166
167-- additional scaled arithmetic tests [0.97 problem]
168subx160 subtract '0'     '.1'      -> '-0.1'
169subx161 subtract '00'    '.97983'  -> '-0.97983'
170subx162 subtract '0'     '.9'      -> '-0.9'
171subx163 subtract '0'     '0.102'   -> '-0.102'
172subx164 subtract '0'     '.4'      -> '-0.4'
173subx165 subtract '0'     '.307'    -> '-0.307'
174subx166 subtract '0'     '.43822'  -> '-0.43822'
175subx167 subtract '0'     '.911'    -> '-0.911'
176subx168 subtract '.0'    '.02'     -> '-0.02'
177subx169 subtract '00'    '.392'    -> '-0.392'
178subx170 subtract '0'     '.26'     -> '-0.26'
179subx171 subtract '0'     '0.51'    -> '-0.51'
180subx172 subtract '0'     '.2234'   -> '-0.2234'
181subx173 subtract '0'     '.2'      -> '-0.2'
182subx174 subtract '.0'    '.0008'   -> '-0.0008'
183-- 0. on left
184subx180 subtract '0.0'     '-.1'      -> '0.1'
185subx181 subtract '0.00'    '-.97983'  -> '0.97983'
186subx182 subtract '0.0'     '-.9'      -> '0.9'
187subx183 subtract '0.0'     '-0.102'   -> '0.102'
188subx184 subtract '0.0'     '-.4'      -> '0.4'
189subx185 subtract '0.0'     '-.307'    -> '0.307'
190subx186 subtract '0.0'     '-.43822'  -> '0.43822'
191subx187 subtract '0.0'     '-.911'    -> '0.911'
192subx188 subtract '0.0'     '-.02'     -> '0.02'
193subx189 subtract '0.00'    '-.392'    -> '0.392'
194subx190 subtract '0.0'     '-.26'     -> '0.26'
195subx191 subtract '0.0'     '-0.51'    -> '0.51'
196subx192 subtract '0.0'     '-.2234'   -> '0.2234'
197subx193 subtract '0.0'     '-.2'      -> '0.2'
198subx194 subtract '0.0'     '-.0008'   -> '0.0008'
199-- negatives of same
200subx200 subtract '0'     '-.1'      -> '0.1'
201subx201 subtract '00'    '-.97983'  -> '0.97983'
202subx202 subtract '0'     '-.9'      -> '0.9'
203subx203 subtract '0'     '-0.102'   -> '0.102'
204subx204 subtract '0'     '-.4'      -> '0.4'
205subx205 subtract '0'     '-.307'    -> '0.307'
206subx206 subtract '0'     '-.43822'  -> '0.43822'
207subx207 subtract '0'     '-.911'    -> '0.911'
208subx208 subtract '.0'    '-.02'     -> '0.02'
209subx209 subtract '00'    '-.392'    -> '0.392'
210subx210 subtract '0'     '-.26'     -> '0.26'
211subx211 subtract '0'     '-0.51'    -> '0.51'
212subx212 subtract '0'     '-.2234'   -> '0.2234'
213subx213 subtract '0'     '-.2'      -> '0.2'
214subx214 subtract '.0'    '-.0008'   -> '0.0008'
215
216-- more fixed, LHS swaps [really the same as testcases under add]
217subx220 subtract '-56267E-12' 0  -> '-5.6267E-8'
218subx221 subtract '-56267E-11' 0  -> '-5.6267E-7'
219subx222 subtract '-56267E-10' 0  -> '-0.0000056267'
220subx223 subtract '-56267E-9'  0  -> '-0.000056267'
221subx224 subtract '-56267E-8'  0  -> '-0.00056267'
222subx225 subtract '-56267E-7'  0  -> '-0.0056267'
223subx226 subtract '-56267E-6'  0  -> '-0.056267'
224subx227 subtract '-56267E-5'  0  -> '-0.56267'
225subx228 subtract '-56267E-2'  0  -> '-562.67'
226subx229 subtract '-56267E-1'  0  -> '-5626.7'
227subx230 subtract '-56267E-0'  0  -> '-56267'
228-- symmetry ...
229subx240 subtract 0 '-56267E-12'  -> '5.6267E-8'
230subx241 subtract 0 '-56267E-11'  -> '5.6267E-7'
231subx242 subtract 0 '-56267E-10'  -> '0.0000056267'
232subx243 subtract 0 '-56267E-9'   -> '0.000056267'
233subx244 subtract 0 '-56267E-8'   -> '0.00056267'
234subx245 subtract 0 '-56267E-7'   -> '0.0056267'
235subx246 subtract 0 '-56267E-6'   -> '0.056267'
236subx247 subtract 0 '-56267E-5'   -> '0.56267'
237subx248 subtract 0 '-56267E-2'   -> '562.67'
238subx249 subtract 0 '-56267E-1'   -> '5626.7'
239subx250 subtract 0 '-56267E-0'   -> '56267'
240
241-- now some more from the 'new' add
242precision: 9
243subx301 subtract '1.23456789'  '1.00000000' -> '0.23456789'
244subx302 subtract '1.23456789'  '1.00000011' -> '0.23456778'
245
246subx311 subtract '0.4444444444'  '0.5555555555' -> '-0.111111111' Inexact Rounded
247subx312 subtract '0.4444444440'  '0.5555555555' -> '-0.111111112' Inexact Rounded
248subx313 subtract '0.4444444444'  '0.5555555550' -> '-0.111111111' Inexact Rounded
249subx314 subtract '0.44444444449'    '0' -> '0.444444444' Inexact Rounded
250subx315 subtract '0.444444444499'   '0' -> '0.444444444' Inexact Rounded
251subx316 subtract '0.4444444444999'  '0' -> '0.444444444' Inexact Rounded
252subx317 subtract '0.4444444445000'  '0' -> '0.444444445' Inexact Rounded
253subx318 subtract '0.4444444445001'  '0' -> '0.444444445' Inexact Rounded
254subx319 subtract '0.444444444501'   '0' -> '0.444444445' Inexact Rounded
255subx320 subtract '0.44444444451'    '0' -> '0.444444445' Inexact Rounded
256
257-- some carrying effects
258subx321 subtract '0.9998'  '0.0000' -> '0.9998'
259subx322 subtract '0.9998'  '0.0001' -> '0.9997'
260subx323 subtract '0.9998'  '0.0002' -> '0.9996'
261subx324 subtract '0.9998'  '0.0003' -> '0.9995'
262subx325 subtract '0.9998'  '-0.0000' -> '0.9998'
263subx326 subtract '0.9998'  '-0.0001' -> '0.9999'
264subx327 subtract '0.9998'  '-0.0002' -> '1.0000'
265subx328 subtract '0.9998'  '-0.0003' -> '1.0001'
266
267subx330 subtract '70'  '10000e+9' -> '-1.00000000E+13' Inexact Rounded
268subx331 subtract '700'  '10000e+9' -> '-1.00000000E+13' Inexact Rounded
269subx332 subtract '7000'  '10000e+9' -> '-9.99999999E+12' Inexact Rounded
270subx333 subtract '70000'  '10000e+9' -> '-9.99999993E+12' Rounded
271subx334 subtract '700000'  '10000e+9' -> '-9.99999930E+12' Rounded
272subx335 subtract '7000000'  '10000e+9' -> '-9.99999300E+12' Rounded
273-- symmetry:
274subx340 subtract '10000e+9'  '70' -> '1.00000000E+13' Inexact Rounded
275subx341 subtract '10000e+9'  '700' -> '1.00000000E+13' Inexact Rounded
276subx342 subtract '10000e+9'  '7000' -> '9.99999999E+12' Inexact Rounded
277subx343 subtract '10000e+9'  '70000' -> '9.99999993E+12' Rounded
278subx344 subtract '10000e+9'  '700000' -> '9.99999930E+12' Rounded
279subx345 subtract '10000e+9'  '7000000' -> '9.99999300E+12' Rounded
280
281-- same, higher precision
282precision: 15
283subx346 subtract '10000e+9'  '7'   -> '9999999999993'
284subx347 subtract '10000e+9'  '70'   -> '9999999999930'
285subx348 subtract '10000e+9'  '700'   -> '9999999999300'
286subx349 subtract '10000e+9'  '7000'   -> '9999999993000'
287subx350 subtract '10000e+9'  '70000'   -> '9999999930000'
288subx351 subtract '10000e+9'  '700000'   -> '9999999300000'
289subx352 subtract '7' '10000e+9'   -> '-9999999999993'
290subx353 subtract '70' '10000e+9'   -> '-9999999999930'
291subx354 subtract '700' '10000e+9'   -> '-9999999999300'
292subx355 subtract '7000' '10000e+9'   -> '-9999999993000'
293subx356 subtract '70000' '10000e+9'   -> '-9999999930000'
294subx357 subtract '700000' '10000e+9'   -> '-9999999300000'
295
296-- zero preservation
297precision: 6
298subx360 subtract '10000e+9'  '70000' -> '1.00000E+13' Inexact Rounded
299subx361 subtract 1 '0.0001' -> '0.9999'
300subx362 subtract 1 '0.00001' -> '0.99999'
301subx363 subtract 1 '0.000001' -> '0.999999'
302subx364 subtract 1 '0.0000001' -> '1.00000' Inexact Rounded
303subx365 subtract 1 '0.00000001' -> '1.00000' Inexact Rounded
304
305-- some funny zeros [in case of bad signum]
306subx370 subtract 1  0  -> 1
307subx371 subtract 1 0.  -> 1
308subx372 subtract 1  .0 -> 1.0
309subx373 subtract 1 0.0 -> 1.0
310subx374 subtract  0  1 -> -1
311subx375 subtract 0.  1 -> -1
312subx376 subtract  .0 1 -> -1.0
313subx377 subtract 0.0 1 -> -1.0
314
315precision: 9
316
317-- leading 0 digit before round
318subx910 subtract -103519362 -51897955.3 -> -51621406.7
319subx911 subtract 159579.444 89827.5229 -> 69751.9211
320
321subx920 subtract 333.123456 33.1234566 -> 299.999999 Inexact Rounded
322subx921 subtract 333.123456 33.1234565 -> 300.000000 Inexact Rounded
323subx922 subtract 133.123456 33.1234565 ->  99.9999995
324subx923 subtract 133.123456 33.1234564 ->  99.9999996
325subx924 subtract 133.123456 33.1234540 -> 100.000002 Rounded
326subx925 subtract 133.123456 43.1234560 ->  90.0000000
327subx926 subtract 133.123456 43.1234561 ->  89.9999999
328subx927 subtract 133.123456 43.1234566 ->  89.9999994
329subx928 subtract 101.123456 91.1234566 ->   9.9999994
330subx929 subtract 101.123456 99.1234566 ->   1.9999994
331
332-- more of the same; probe for cluster boundary problems
333precision: 1
334subx930 subtract  11 2           -> 9
335precision: 2
336subx932 subtract 101 2           -> 99
337precision: 3
338subx934 subtract 101 2.1         -> 98.9
339subx935 subtract 101 92.01       ->  8.99
340precision: 4
341subx936 subtract 101 2.01        -> 98.99
342subx937 subtract 101 92.01       ->  8.99
343subx938 subtract 101 92.006      ->  8.994
344precision: 5
345subx939 subtract 101 2.001       -> 98.999
346subx940 subtract 101 92.001      ->  8.999
347subx941 subtract 101 92.0006     ->  8.9994
348precision: 6
349subx942 subtract 101 2.0001      -> 98.9999
350subx943 subtract 101 92.0001     ->  8.9999
351subx944 subtract 101 92.00006    ->  8.99994
352precision: 7
353subx945 subtract 101 2.00001     -> 98.99999
354subx946 subtract 101 92.00001    ->  8.99999
355subx947 subtract 101 92.000006   ->  8.999994
356precision: 8
357subx948 subtract 101 2.000001    -> 98.999999
358subx949 subtract 101 92.000001   ->  8.999999
359subx950 subtract 101 92.0000006  ->  8.9999994
360precision: 9
361subx951 subtract 101 2.0000001   -> 98.9999999
362subx952 subtract 101 92.0000001  ->  8.9999999
363subx953 subtract 101 92.00000006 ->  8.99999994
364
365precision: 9
366
367-- more LHS swaps [were fixed]
368subx390 subtract '-56267E-10'   0 ->  '-0.0000056267'
369subx391 subtract '-56267E-6'    0 ->  '-0.056267'
370subx392 subtract '-56267E-5'    0 ->  '-0.56267'
371subx393 subtract '-56267E-4'    0 ->  '-5.6267'
372subx394 subtract '-56267E-3'    0 ->  '-56.267'
373subx395 subtract '-56267E-2'    0 ->  '-562.67'
374subx396 subtract '-56267E-1'    0 ->  '-5626.7'
375subx397 subtract '-56267E-0'    0 ->  '-56267'
376subx398 subtract '-5E-10'       0 ->  '-5E-10'
377subx399 subtract '-5E-7'        0 ->  '-5E-7'
378subx400 subtract '-5E-6'        0 ->  '-0.000005'
379subx401 subtract '-5E-5'        0 ->  '-0.00005'
380subx402 subtract '-5E-4'        0 ->  '-0.0005'
381subx403 subtract '-5E-1'        0 ->  '-0.5'
382subx404 subtract '-5E0'         0 ->  '-5'
383subx405 subtract '-5E1'         0 ->  '-50'
384subx406 subtract '-5E5'         0 ->  '-500000'
385subx407 subtract '-5E8'         0 ->  '-500000000'
386subx408 subtract '-5E9'         0 ->  '-5.00000000E+9'   Rounded
387subx409 subtract '-5E10'        0 ->  '-5.00000000E+10'  Rounded
388subx410 subtract '-5E11'        0 ->  '-5.00000000E+11'  Rounded
389subx411 subtract '-5E100'       0 ->  '-5.00000000E+100' Rounded
390
391-- more RHS swaps [were fixed]
392subx420 subtract 0  '-56267E-10' ->  '0.0000056267'
393subx421 subtract 0  '-56267E-6'  ->  '0.056267'
394subx422 subtract 0  '-56267E-5'  ->  '0.56267'
395subx423 subtract 0  '-56267E-4'  ->  '5.6267'
396subx424 subtract 0  '-56267E-3'  ->  '56.267'
397subx425 subtract 0  '-56267E-2'  ->  '562.67'
398subx426 subtract 0  '-56267E-1'  ->  '5626.7'
399subx427 subtract 0  '-56267E-0'  ->  '56267'
400subx428 subtract 0  '-5E-10'     ->  '5E-10'
401subx429 subtract 0  '-5E-7'      ->  '5E-7'
402subx430 subtract 0  '-5E-6'      ->  '0.000005'
403subx431 subtract 0  '-5E-5'      ->  '0.00005'
404subx432 subtract 0  '-5E-4'      ->  '0.0005'
405subx433 subtract 0  '-5E-1'      ->  '0.5'
406subx434 subtract 0  '-5E0'       ->  '5'
407subx435 subtract 0  '-5E1'       ->  '50'
408subx436 subtract 0  '-5E5'       ->  '500000'
409subx437 subtract 0  '-5E8'       ->  '500000000'
410subx438 subtract 0  '-5E9'       ->  '5.00000000E+9'    Rounded
411subx439 subtract 0  '-5E10'      ->  '5.00000000E+10'   Rounded
412subx440 subtract 0  '-5E11'      ->  '5.00000000E+11'   Rounded
413subx441 subtract 0  '-5E100'     ->  '5.00000000E+100'  Rounded
414
415
416-- try borderline precision, with carries, etc.
417precision: 15
418subx461 subtract '1E+12' '1'       -> '999999999999'
419subx462 subtract '1E+12' '-1.11'   -> '1000000000001.11'
420subx463 subtract '1.11'  '-1E+12'  -> '1000000000001.11'
421subx464 subtract '-1'    '-1E+12'  -> '999999999999'
422subx465 subtract '7E+12' '1'       -> '6999999999999'
423subx466 subtract '7E+12' '-1.11'   -> '7000000000001.11'
424subx467 subtract '1.11'  '-7E+12'  -> '7000000000001.11'
425subx468 subtract '-1'    '-7E+12'  -> '6999999999999'
426
427--                 123456789012345       123456789012345      1 23456789012345
428subx470 subtract '0.444444444444444'  '-0.555555555555563' -> '1.00000000000001' Inexact Rounded
429subx471 subtract '0.444444444444444'  '-0.555555555555562' -> '1.00000000000001' Inexact Rounded
430subx472 subtract '0.444444444444444'  '-0.555555555555561' -> '1.00000000000001' Inexact Rounded
431subx473 subtract '0.444444444444444'  '-0.555555555555560' -> '1.00000000000000' Inexact Rounded
432subx474 subtract '0.444444444444444'  '-0.555555555555559' -> '1.00000000000000' Inexact Rounded
433subx475 subtract '0.444444444444444'  '-0.555555555555558' -> '1.00000000000000' Inexact Rounded
434subx476 subtract '0.444444444444444'  '-0.555555555555557' -> '1.00000000000000' Inexact Rounded
435subx477 subtract '0.444444444444444'  '-0.555555555555556' -> '1.00000000000000' Rounded
436subx478 subtract '0.444444444444444'  '-0.555555555555555' -> '0.999999999999999'
437subx479 subtract '0.444444444444444'  '-0.555555555555554' -> '0.999999999999998'
438subx480 subtract '0.444444444444444'  '-0.555555555555553' -> '0.999999999999997'
439subx481 subtract '0.444444444444444'  '-0.555555555555552' -> '0.999999999999996'
440subx482 subtract '0.444444444444444'  '-0.555555555555551' -> '0.999999999999995'
441subx483 subtract '0.444444444444444'  '-0.555555555555550' -> '0.999999999999994'
442
443-- and some more, including residue effects and different roundings
444precision: 9
445rounding: half_up
446subx500 subtract '123456789' 0             -> '123456789'
447subx501 subtract '123456789' 0.000000001   -> '123456789' Inexact Rounded
448subx502 subtract '123456789' 0.000001      -> '123456789' Inexact Rounded
449subx503 subtract '123456789' 0.1           -> '123456789' Inexact Rounded
450subx504 subtract '123456789' 0.4           -> '123456789' Inexact Rounded
451subx505 subtract '123456789' 0.49          -> '123456789' Inexact Rounded
452subx506 subtract '123456789' 0.499999      -> '123456789' Inexact Rounded
453subx507 subtract '123456789' 0.499999999   -> '123456789' Inexact Rounded
454subx508 subtract '123456789' 0.5           -> '123456789' Inexact Rounded
455subx509 subtract '123456789' 0.500000001   -> '123456788' Inexact Rounded
456subx510 subtract '123456789' 0.500001      -> '123456788' Inexact Rounded
457subx511 subtract '123456789' 0.51          -> '123456788' Inexact Rounded
458subx512 subtract '123456789' 0.6           -> '123456788' Inexact Rounded
459subx513 subtract '123456789' 0.9           -> '123456788' Inexact Rounded
460subx514 subtract '123456789' 0.99999       -> '123456788' Inexact Rounded
461subx515 subtract '123456789' 0.999999999   -> '123456788' Inexact Rounded
462subx516 subtract '123456789' 1             -> '123456788'
463subx517 subtract '123456789' 1.000000001   -> '123456788' Inexact Rounded
464subx518 subtract '123456789' 1.00001       -> '123456788' Inexact Rounded
465subx519 subtract '123456789' 1.1           -> '123456788' Inexact Rounded
466
467rounding: half_even
468subx520 subtract '123456789' 0             -> '123456789'
469subx521 subtract '123456789' 0.000000001   -> '123456789' Inexact Rounded
470subx522 subtract '123456789' 0.000001      -> '123456789' Inexact Rounded
471subx523 subtract '123456789' 0.1           -> '123456789' Inexact Rounded
472subx524 subtract '123456789' 0.4           -> '123456789' Inexact Rounded
473subx525 subtract '123456789' 0.49          -> '123456789' Inexact Rounded
474subx526 subtract '123456789' 0.499999      -> '123456789' Inexact Rounded
475subx527 subtract '123456789' 0.499999999   -> '123456789' Inexact Rounded
476subx528 subtract '123456789' 0.5           -> '123456788' Inexact Rounded
477subx529 subtract '123456789' 0.500000001   -> '123456788' Inexact Rounded
478subx530 subtract '123456789' 0.500001      -> '123456788' Inexact Rounded
479subx531 subtract '123456789' 0.51          -> '123456788' Inexact Rounded
480subx532 subtract '123456789' 0.6           -> '123456788' Inexact Rounded
481subx533 subtract '123456789' 0.9           -> '123456788' Inexact Rounded
482subx534 subtract '123456789' 0.99999       -> '123456788' Inexact Rounded
483subx535 subtract '123456789' 0.999999999   -> '123456788' Inexact Rounded
484subx536 subtract '123456789' 1             -> '123456788'
485subx537 subtract '123456789' 1.00000001    -> '123456788' Inexact Rounded
486subx538 subtract '123456789' 1.00001       -> '123456788' Inexact Rounded
487subx539 subtract '123456789' 1.1           -> '123456788' Inexact Rounded
488-- critical few with even bottom digit...
489subx540 subtract '123456788' 0.499999999   -> '123456788' Inexact Rounded
490subx541 subtract '123456788' 0.5           -> '123456788' Inexact Rounded
491subx542 subtract '123456788' 0.500000001   -> '123456787' Inexact Rounded
492
493rounding: down
494subx550 subtract '123456789' 0             -> '123456789'
495subx551 subtract '123456789' 0.000000001   -> '123456788' Inexact Rounded
496subx552 subtract '123456789' 0.000001      -> '123456788' Inexact Rounded
497subx553 subtract '123456789' 0.1           -> '123456788' Inexact Rounded
498subx554 subtract '123456789' 0.4           -> '123456788' Inexact Rounded
499subx555 subtract '123456789' 0.49          -> '123456788' Inexact Rounded
500subx556 subtract '123456789' 0.499999      -> '123456788' Inexact Rounded
501subx557 subtract '123456789' 0.499999999   -> '123456788' Inexact Rounded
502subx558 subtract '123456789' 0.5           -> '123456788' Inexact Rounded
503subx559 subtract '123456789' 0.500000001   -> '123456788' Inexact Rounded
504subx560 subtract '123456789' 0.500001      -> '123456788' Inexact Rounded
505subx561 subtract '123456789' 0.51          -> '123456788' Inexact Rounded
506subx562 subtract '123456789' 0.6           -> '123456788' Inexact Rounded
507subx563 subtract '123456789' 0.9           -> '123456788' Inexact Rounded
508subx564 subtract '123456789' 0.99999       -> '123456788' Inexact Rounded
509subx565 subtract '123456789' 0.999999999   -> '123456788' Inexact Rounded
510subx566 subtract '123456789' 1             -> '123456788'
511subx567 subtract '123456789' 1.00000001    -> '123456787' Inexact Rounded
512subx568 subtract '123456789' 1.00001       -> '123456787' Inexact Rounded
513subx569 subtract '123456789' 1.1           -> '123456787' Inexact Rounded
514
515-- symmetry...
516rounding: half_up
517subx600 subtract 0             '123456789' -> '-123456789'
518subx601 subtract 0.000000001   '123456789' -> '-123456789' Inexact Rounded
519subx602 subtract 0.000001      '123456789' -> '-123456789' Inexact Rounded
520subx603 subtract 0.1           '123456789' -> '-123456789' Inexact Rounded
521subx604 subtract 0.4           '123456789' -> '-123456789' Inexact Rounded
522subx605 subtract 0.49          '123456789' -> '-123456789' Inexact Rounded
523subx606 subtract 0.499999      '123456789' -> '-123456789' Inexact Rounded
524subx607 subtract 0.499999999   '123456789' -> '-123456789' Inexact Rounded
525subx608 subtract 0.5           '123456789' -> '-123456789' Inexact Rounded
526subx609 subtract 0.500000001   '123456789' -> '-123456788' Inexact Rounded
527subx610 subtract 0.500001      '123456789' -> '-123456788' Inexact Rounded
528subx611 subtract 0.51          '123456789' -> '-123456788' Inexact Rounded
529subx612 subtract 0.6           '123456789' -> '-123456788' Inexact Rounded
530subx613 subtract 0.9           '123456789' -> '-123456788' Inexact Rounded
531subx614 subtract 0.99999       '123456789' -> '-123456788' Inexact Rounded
532subx615 subtract 0.999999999   '123456789' -> '-123456788' Inexact Rounded
533subx616 subtract 1             '123456789' -> '-123456788'
534subx617 subtract 1.000000001   '123456789' -> '-123456788' Inexact Rounded
535subx618 subtract 1.00001       '123456789' -> '-123456788' Inexact Rounded
536subx619 subtract 1.1           '123456789' -> '-123456788' Inexact Rounded
537
538rounding: half_even
539subx620 subtract 0             '123456789' -> '-123456789'
540subx621 subtract 0.000000001   '123456789' -> '-123456789' Inexact Rounded
541subx622 subtract 0.000001      '123456789' -> '-123456789' Inexact Rounded
542subx623 subtract 0.1           '123456789' -> '-123456789' Inexact Rounded
543subx624 subtract 0.4           '123456789' -> '-123456789' Inexact Rounded
544subx625 subtract 0.49          '123456789' -> '-123456789' Inexact Rounded
545subx626 subtract 0.499999      '123456789' -> '-123456789' Inexact Rounded
546subx627 subtract 0.499999999   '123456789' -> '-123456789' Inexact Rounded
547subx628 subtract 0.5           '123456789' -> '-123456788' Inexact Rounded
548subx629 subtract 0.500000001   '123456789' -> '-123456788' Inexact Rounded
549subx630 subtract 0.500001      '123456789' -> '-123456788' Inexact Rounded
550subx631 subtract 0.51          '123456789' -> '-123456788' Inexact Rounded
551subx632 subtract 0.6           '123456789' -> '-123456788' Inexact Rounded
552subx633 subtract 0.9           '123456789' -> '-123456788' Inexact Rounded
553subx634 subtract 0.99999       '123456789' -> '-123456788' Inexact Rounded
554subx635 subtract 0.999999999   '123456789' -> '-123456788' Inexact Rounded
555subx636 subtract 1             '123456789' -> '-123456788'
556subx637 subtract 1.00000001    '123456789' -> '-123456788' Inexact Rounded
557subx638 subtract 1.00001       '123456789' -> '-123456788' Inexact Rounded
558subx639 subtract 1.1           '123456789' -> '-123456788' Inexact Rounded
559-- critical few with even bottom digit...
560subx640 subtract 0.499999999   '123456788' -> '-123456788' Inexact Rounded
561subx641 subtract 0.5           '123456788' -> '-123456788' Inexact Rounded
562subx642 subtract 0.500000001   '123456788' -> '-123456787' Inexact Rounded
563
564rounding: down
565subx650 subtract 0             '123456789' -> '-123456789'
566subx651 subtract 0.000000001   '123456789' -> '-123456788' Inexact Rounded
567subx652 subtract 0.000001      '123456789' -> '-123456788' Inexact Rounded
568subx653 subtract 0.1           '123456789' -> '-123456788' Inexact Rounded
569subx654 subtract 0.4           '123456789' -> '-123456788' Inexact Rounded
570subx655 subtract 0.49          '123456789' -> '-123456788' Inexact Rounded
571subx656 subtract 0.499999      '123456789' -> '-123456788' Inexact Rounded
572subx657 subtract 0.499999999   '123456789' -> '-123456788' Inexact Rounded
573subx658 subtract 0.5           '123456789' -> '-123456788' Inexact Rounded
574subx659 subtract 0.500000001   '123456789' -> '-123456788' Inexact Rounded
575subx660 subtract 0.500001      '123456789' -> '-123456788' Inexact Rounded
576subx661 subtract 0.51          '123456789' -> '-123456788' Inexact Rounded
577subx662 subtract 0.6           '123456789' -> '-123456788' Inexact Rounded
578subx663 subtract 0.9           '123456789' -> '-123456788' Inexact Rounded
579subx664 subtract 0.99999       '123456789' -> '-123456788' Inexact Rounded
580subx665 subtract 0.999999999   '123456789' -> '-123456788' Inexact Rounded
581subx666 subtract 1             '123456789' -> '-123456788'
582subx667 subtract 1.00000001    '123456789' -> '-123456787' Inexact Rounded
583subx668 subtract 1.00001       '123456789' -> '-123456787' Inexact Rounded
584subx669 subtract 1.1           '123456789' -> '-123456787' Inexact Rounded
585
586
587-- lots of leading zeros in intermediate result, and showing effects of
588-- input rounding would have affected the following
589precision: 9
590rounding: half_up
591subx670 subtract '123456789' '123456788.1' -> 0.9
592subx671 subtract '123456789' '123456788.9' -> 0.1
593subx672 subtract '123456789' '123456789.1' -> -0.1
594subx673 subtract '123456789' '123456789.5' -> -0.5
595subx674 subtract '123456789' '123456789.9' -> -0.9
596
597rounding: half_even
598subx680 subtract '123456789' '123456788.1' -> 0.9
599subx681 subtract '123456789' '123456788.9' -> 0.1
600subx682 subtract '123456789' '123456789.1' -> -0.1
601subx683 subtract '123456789' '123456789.5' -> -0.5
602subx684 subtract '123456789' '123456789.9' -> -0.9
603
604subx685 subtract '123456788' '123456787.1' -> 0.9
605subx686 subtract '123456788' '123456787.9' -> 0.1
606subx687 subtract '123456788' '123456788.1' -> -0.1
607subx688 subtract '123456788' '123456788.5' -> -0.5
608subx689 subtract '123456788' '123456788.9' -> -0.9
609
610rounding: down
611subx690 subtract '123456789' '123456788.1' -> 0.9
612subx691 subtract '123456789' '123456788.9' -> 0.1
613subx692 subtract '123456789' '123456789.1' -> -0.1
614subx693 subtract '123456789' '123456789.5' -> -0.5
615subx694 subtract '123456789' '123456789.9' -> -0.9
616
617-- input preparation tests
618rounding: half_up
619precision: 3
620
621subx700 subtract '12345678900000'  -9999999999999 ->  '2.23E+13' Inexact Rounded
622subx701 subtract  '9999999999999' -12345678900000 ->  '2.23E+13' Inexact Rounded
623subx702 subtract '12E+3'  '-3456' ->  '1.55E+4' Inexact Rounded
624subx703 subtract '12E+3'  '-3446' ->  '1.54E+4' Inexact Rounded
625subx704 subtract '12E+3'  '-3454' ->  '1.55E+4' Inexact Rounded
626subx705 subtract '12E+3'  '-3444' ->  '1.54E+4' Inexact Rounded
627
628subx706 subtract '3456'  '-12E+3' ->  '1.55E+4' Inexact Rounded
629subx707 subtract '3446'  '-12E+3' ->  '1.54E+4' Inexact Rounded
630subx708 subtract '3454'  '-12E+3' ->  '1.55E+4' Inexact Rounded
631subx709 subtract '3444'  '-12E+3' ->  '1.54E+4' Inexact Rounded
632
633-- overflow and underflow tests [subnormals now possible]
634maxexponent: 999999999
635minexponent: -999999999
636precision: 9
637rounding: down
638subx710 subtract 1E+999999999    -9E+999999999   -> 9.99999999E+999999999 Overflow Inexact Rounded
639subx711 subtract 9E+999999999    -1E+999999999   -> 9.99999999E+999999999 Overflow Inexact Rounded
640rounding: half_up
641subx712 subtract 1E+999999999    -9E+999999999   -> Infinity Overflow Inexact Rounded
642subx713 subtract 9E+999999999    -1E+999999999   -> Infinity Overflow Inexact Rounded
643subx714 subtract -1.1E-999999999 -1E-999999999   -> -1E-1000000000 Subnormal
644subx715 subtract 1E-999999999    +1.1e-999999999 -> -1E-1000000000 Subnormal
645subx716 subtract -1E+999999999   +9E+999999999   -> -Infinity Overflow Inexact Rounded
646subx717 subtract -9E+999999999   +1E+999999999   -> -Infinity Overflow Inexact Rounded
647subx718 subtract +1.1E-999999999 +1E-999999999   -> 1E-1000000000 Subnormal
648subx719 subtract -1E-999999999   -1.1e-999999999 -> 1E-1000000000 Subnormal
649
650precision: 3
651subx720 subtract 1  9.999E+999999999   -> -Infinity Inexact Overflow Rounded
652subx721 subtract 1 -9.999E+999999999   ->  Infinity Inexact Overflow Rounded
653subx722 subtract    9.999E+999999999 1 ->  Infinity Inexact Overflow Rounded
654subx723 subtract   -9.999E+999999999 1 -> -Infinity Inexact Overflow Rounded
655subx724 subtract 1  9.999E+999999999   -> -Infinity Inexact Overflow Rounded
656subx725 subtract 1 -9.999E+999999999   ->  Infinity Inexact Overflow Rounded
657subx726 subtract    9.999E+999999999 1 ->  Infinity Inexact Overflow Rounded
658subx727 subtract   -9.999E+999999999 1 -> -Infinity Inexact Overflow Rounded
659
660-- [more below]
661
662-- long operand checks
663maxexponent: 999
664minexponent: -999
665precision: 9
666sub731 subtract 12345678000 0 ->  1.23456780E+10 Rounded
667sub732 subtract 0 12345678000 -> -1.23456780E+10 Rounded
668sub733 subtract 1234567800  0 ->  1.23456780E+9 Rounded
669sub734 subtract 0 1234567800  -> -1.23456780E+9 Rounded
670sub735 subtract 1234567890  0 ->  1.23456789E+9 Rounded
671sub736 subtract 0 1234567890  -> -1.23456789E+9 Rounded
672sub737 subtract 1234567891  0 ->  1.23456789E+9 Inexact Rounded
673sub738 subtract 0 1234567891  -> -1.23456789E+9 Inexact Rounded
674sub739 subtract 12345678901 0 ->  1.23456789E+10 Inexact Rounded
675sub740 subtract 0 12345678901 -> -1.23456789E+10 Inexact Rounded
676sub741 subtract 1234567896  0 ->  1.23456790E+9 Inexact Rounded
677sub742 subtract 0 1234567896  -> -1.23456790E+9 Inexact Rounded
678
679precision: 15
680sub751 subtract 12345678000 0 ->  12345678000
681sub752 subtract 0 12345678000 -> -12345678000
682sub753 subtract 1234567800  0 ->  1234567800
683sub754 subtract 0 1234567800  -> -1234567800
684sub755 subtract 1234567890  0 ->  1234567890
685sub756 subtract 0 1234567890  -> -1234567890
686sub757 subtract 1234567891  0 ->  1234567891
687sub758 subtract 0 1234567891  -> -1234567891
688sub759 subtract 12345678901 0 ->  12345678901
689sub760 subtract 0 12345678901 -> -12345678901
690sub761 subtract 1234567896  0 ->  1234567896
691sub762 subtract 0 1234567896  -> -1234567896
692
693-- Specials
694subx780 subtract -Inf   Inf   -> -Infinity
695subx781 subtract -Inf   1000  -> -Infinity
696subx782 subtract -Inf   1     -> -Infinity
697subx783 subtract -Inf  -0     -> -Infinity
698subx784 subtract -Inf  -1     -> -Infinity
699subx785 subtract -Inf  -1000  -> -Infinity
700subx787 subtract -1000  Inf   -> -Infinity
701subx788 subtract -Inf   Inf   -> -Infinity
702subx789 subtract -1     Inf   -> -Infinity
703subx790 subtract  0     Inf   -> -Infinity
704subx791 subtract  1     Inf   -> -Infinity
705subx792 subtract  1000  Inf   -> -Infinity
706
707subx800 subtract  Inf   Inf   ->  NaN  Invalid_operation
708subx801 subtract  Inf   1000  ->  Infinity
709subx802 subtract  Inf   1     ->  Infinity
710subx803 subtract  Inf   0     ->  Infinity
711subx804 subtract  Inf  -0     ->  Infinity
712subx805 subtract  Inf  -1     ->  Infinity
713subx806 subtract  Inf  -1000  ->  Infinity
714subx807 subtract  Inf  -Inf   ->  Infinity
715subx808 subtract -1000 -Inf   ->  Infinity
716subx809 subtract -Inf  -Inf   ->  NaN  Invalid_operation
717subx810 subtract -1    -Inf   ->  Infinity
718subx811 subtract -0    -Inf   ->  Infinity
719subx812 subtract  0    -Inf   ->  Infinity
720subx813 subtract  1    -Inf   ->  Infinity
721subx814 subtract  1000 -Inf   ->  Infinity
722subx815 subtract  Inf  -Inf   ->  Infinity
723
724subx821 subtract  NaN   Inf   ->  NaN
725subx822 subtract -NaN   1000  -> -NaN
726subx823 subtract  NaN   1     ->  NaN
727subx824 subtract  NaN   0     ->  NaN
728subx825 subtract  NaN  -0     ->  NaN
729subx826 subtract  NaN  -1     ->  NaN
730subx827 subtract  NaN  -1000  ->  NaN
731subx828 subtract  NaN  -Inf   ->  NaN
732subx829 subtract -NaN   NaN   -> -NaN
733subx830 subtract -Inf   NaN   ->  NaN
734subx831 subtract -1000  NaN   ->  NaN
735subx832 subtract -1     NaN   ->  NaN
736subx833 subtract -0     NaN   ->  NaN
737subx834 subtract  0     NaN   ->  NaN
738subx835 subtract  1     NaN   ->  NaN
739subx836 subtract  1000 -NaN   -> -NaN
740subx837 subtract  Inf   NaN   ->  NaN
741
742subx841 subtract  sNaN  Inf   ->  NaN  Invalid_operation
743subx842 subtract -sNaN  1000  -> -NaN  Invalid_operation
744subx843 subtract  sNaN  1     ->  NaN  Invalid_operation
745subx844 subtract  sNaN  0     ->  NaN  Invalid_operation
746subx845 subtract  sNaN -0     ->  NaN  Invalid_operation
747subx846 subtract  sNaN -1     ->  NaN  Invalid_operation
748subx847 subtract  sNaN -1000  ->  NaN  Invalid_operation
749subx848 subtract  sNaN  NaN   ->  NaN  Invalid_operation
750subx849 subtract  sNaN sNaN   ->  NaN  Invalid_operation
751subx850 subtract  NaN  sNaN   ->  NaN  Invalid_operation
752subx851 subtract -Inf -sNaN   -> -NaN  Invalid_operation
753subx852 subtract -1000 sNaN   ->  NaN  Invalid_operation
754subx853 subtract -1    sNaN   ->  NaN  Invalid_operation
755subx854 subtract -0    sNaN   ->  NaN  Invalid_operation
756subx855 subtract  0    sNaN   ->  NaN  Invalid_operation
757subx856 subtract  1    sNaN   ->  NaN  Invalid_operation
758subx857 subtract  1000 sNaN   ->  NaN  Invalid_operation
759subx858 subtract  Inf  sNaN   ->  NaN  Invalid_operation
760subx859 subtract  NaN  sNaN   ->  NaN  Invalid_operation
761
762-- propagating NaNs
763subx861 subtract  NaN01   -Inf     ->  NaN1
764subx862 subtract -NaN02   -1000    -> -NaN2
765subx863 subtract  NaN03    1000    ->  NaN3
766subx864 subtract  NaN04    Inf     ->  NaN4
767subx865 subtract  NaN05    NaN61   ->  NaN5
768subx866 subtract -Inf     -NaN71   -> -NaN71
769subx867 subtract -1000     NaN81   ->  NaN81
770subx868 subtract  1000     NaN91   ->  NaN91
771subx869 subtract  Inf      NaN101  ->  NaN101
772subx871 subtract  sNaN011  -Inf    ->  NaN11  Invalid_operation
773subx872 subtract  sNaN012  -1000   ->  NaN12  Invalid_operation
774subx873 subtract -sNaN013   1000   -> -NaN13  Invalid_operation
775subx874 subtract  sNaN014   NaN171 ->  NaN14  Invalid_operation
776subx875 subtract  sNaN015  sNaN181 ->  NaN15  Invalid_operation
777subx876 subtract  NaN016   sNaN191 ->  NaN191 Invalid_operation
778subx877 subtract -Inf      sNaN201 ->  NaN201 Invalid_operation
779subx878 subtract -1000     sNaN211 ->  NaN211 Invalid_operation
780subx879 subtract  1000    -sNaN221 -> -NaN221 Invalid_operation
781subx880 subtract  Inf      sNaN231 ->  NaN231 Invalid_operation
782subx881 subtract  NaN025   sNaN241 ->  NaN241 Invalid_operation
783
784-- edge case spills
785subx901 subtract  2.E-3  1.002  -> -1.000
786subx902 subtract  2.0E-3  1.002  -> -1.0000
787subx903 subtract  2.00E-3  1.0020  -> -1.00000
788subx904 subtract  2.000E-3  1.00200  -> -1.000000
789subx905 subtract  2.0000E-3  1.002000  -> -1.0000000
790subx906 subtract  2.00000E-3  1.0020000  -> -1.00000000
791subx907 subtract  2.000000E-3  1.00200000  -> -1.000000000
792subx908 subtract  2.0000000E-3  1.002000000  -> -1.0000000000
793
794-- subnormals and underflows
795precision: 3
796maxexponent: 999
797minexponent: -999
798subx1010 subtract  0  1.00E-999       ->  -1.00E-999
799subx1011 subtract  0  0.1E-999        ->  -1E-1000   Subnormal
800subx1012 subtract  0  0.10E-999       ->  -1.0E-1000 Subnormal
801subx1013 subtract  0  0.100E-999      ->  -1.0E-1000 Subnormal Rounded
802subx1014 subtract  0  0.01E-999       ->  -1E-1001   Subnormal
803-- next is rounded to Emin
804subx1015 subtract  0  0.999E-999      ->  -1.00E-999 Inexact Rounded Subnormal Underflow
805subx1016 subtract  0  0.099E-999      ->  -1.0E-1000 Inexact Rounded Subnormal Underflow
806subx1017 subtract  0  0.009E-999      ->  -1E-1001   Inexact Rounded Subnormal Underflow
807subx1018 subtract  0  0.001E-999      ->  -0E-1001   Inexact Rounded Subnormal Underflow Clamped
808subx1019 subtract  0  0.0009E-999     ->  -0E-1001   Inexact Rounded Subnormal Underflow Clamped
809subx1020 subtract  0  0.0001E-999     ->  -0E-1001   Inexact Rounded Subnormal Underflow Clamped
810
811subx1030 subtract  0 -1.00E-999       ->   1.00E-999
812subx1031 subtract  0 -0.1E-999        ->   1E-1000   Subnormal
813subx1032 subtract  0 -0.10E-999       ->   1.0E-1000 Subnormal
814subx1033 subtract  0 -0.100E-999      ->   1.0E-1000 Subnormal Rounded
815subx1034 subtract  0 -0.01E-999       ->   1E-1001   Subnormal
816-- next is rounded to Emin
817subx1035 subtract  0 -0.999E-999      ->   1.00E-999 Inexact Rounded Subnormal Underflow
818subx1036 subtract  0 -0.099E-999      ->   1.0E-1000 Inexact Rounded Subnormal Underflow
819subx1037 subtract  0 -0.009E-999      ->   1E-1001   Inexact Rounded Subnormal Underflow
820subx1038 subtract  0 -0.001E-999      ->   0E-1001   Inexact Rounded Subnormal Underflow Clamped
821subx1039 subtract  0 -0.0009E-999     ->   0E-1001   Inexact Rounded Subnormal Underflow Clamped
822subx1040 subtract  0 -0.0001E-999     ->   0E-1001   Inexact Rounded Subnormal Underflow Clamped
823
824-- some non-zero subnormal subtracts
825-- subx1056 is a tricky case
826rounding: half_up
827subx1050 subtract  1.00E-999   0.1E-999  ->   9.0E-1000  Subnormal
828subx1051 subtract  0.1E-999    0.1E-999  ->   0E-1000
829subx1052 subtract  0.10E-999   0.1E-999  ->   0E-1001
830subx1053 subtract  0.100E-999  0.1E-999  ->   0E-1001    Clamped
831subx1054 subtract  0.01E-999   0.1E-999  ->   -9E-1001   Subnormal
832subx1055 subtract  0.999E-999  0.1E-999  ->   9.0E-1000  Inexact Rounded Subnormal Underflow
833subx1056 subtract  0.099E-999  0.1E-999  ->   -0E-1001   Inexact Rounded Subnormal Underflow Clamped
834subx1057 subtract  0.009E-999  0.1E-999  ->   -9E-1001   Inexact Rounded Subnormal Underflow
835subx1058 subtract  0.001E-999  0.1E-999  ->   -1.0E-1000 Inexact Rounded Subnormal Underflow
836subx1059 subtract  0.0009E-999 0.1E-999  ->   -1.0E-1000 Inexact Rounded Subnormal Underflow
837subx1060 subtract  0.0001E-999 0.1E-999  ->   -1.0E-1000 Inexact Rounded Subnormal Underflow
838
839
840-- check for double-rounded subnormals
841precision:   5
842maxexponent: 79
843minexponent: -79
844subx1101 subtract  0 1.52444E-80 -> -1.524E-80 Inexact Rounded Subnormal Underflow
845subx1102 subtract  0 1.52445E-80 -> -1.524E-80 Inexact Rounded Subnormal Underflow
846subx1103 subtract  0 1.52446E-80 -> -1.524E-80 Inexact Rounded Subnormal Underflow
847subx1104 subtract  1.52444E-80 0 ->  1.524E-80 Inexact Rounded Subnormal Underflow
848subx1105 subtract  1.52445E-80 0 ->  1.524E-80 Inexact Rounded Subnormal Underflow
849subx1106 subtract  1.52446E-80 0 ->  1.524E-80 Inexact Rounded Subnormal Underflow
850
851subx1111 subtract  1.2345678E-80  1.2345671E-80 ->  0E-83 Inexact Rounded Subnormal Underflow Clamped
852subx1112 subtract  1.2345678E-80  1.2345618E-80 ->  0E-83 Inexact Rounded Subnormal Underflow Clamped
853subx1113 subtract  1.2345678E-80  1.2345178E-80 ->  0E-83 Inexact Rounded Subnormal Underflow Clamped
854subx1114 subtract  1.2345678E-80  1.2341678E-80 ->  0E-83 Inexact Rounded Subnormal Underflow Clamped
855subx1115 subtract  1.2345678E-80  1.2315678E-80 ->  3E-83         Rounded Subnormal
856subx1116 subtract  1.2345678E-80  1.2145678E-80 ->  2.0E-82       Rounded Subnormal
857subx1117 subtract  1.2345678E-80  1.1345678E-80 ->  1.00E-81      Rounded Subnormal
858subx1118 subtract  1.2345678E-80  0.2345678E-80 ->  1.000E-80     Rounded Subnormal
859
860precision:   34
861rounding:    half_up
862maxExponent: 6144
863minExponent: -6143
864-- Examples from SQL proposal (Krishna Kulkarni)
865subx1125  subtract 130E-2  120E-2 -> 0.10
866subx1126  subtract 130E-2  12E-1  -> 0.10
867subx1127  subtract 130E-2  1E0    -> 0.30
868subx1128  subtract 1E2     1E4    -> -9.9E+3
869
870-- Null tests
871subx9990 subtract 10  # -> NaN Invalid_operation
872subx9991 subtract  # 10 -> NaN Invalid_operation
873