| Thomas Wouters | 1b7f891 | 2007-09-19 03:06:30 +0000 | [diff] [blame] | 1 | ------------------------------------------------------------------------
 | 
 | 2 | -- logb.decTest -- return integral adjusted exponent as per 754r      --
 | 
| Mark Dickinson | 8a54653 | 2009-10-08 16:30:38 +0000 | [diff] [blame] | 3 | -- Copyright (c) IBM Corporation, 2005, 2009.  All rights reserved.   --
 | 
| Thomas Wouters | 1b7f891 | 2007-09-19 03:06:30 +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
 | 
| Thomas Wouters | 1b7f891 | 2007-09-19 03:06:30 +0000 | [diff] [blame] | 21 | 
 | 
 | 22 | -- This emphasises the testing of notable cases, as they will often
 | 
 | 23 | -- have unusual paths (especially the 10**n results).
 | 
 | 24 | 
 | 
 | 25 | extended:    1
 | 
 | 26 | rounding:    half_even
 | 
 | 27 | maxExponent: 999
 | 
 | 28 | minexponent: -999
 | 
 | 29 | 
 | 
 | 30 | -- basics & examples
 | 
 | 31 | precision:   9
 | 
 | 32 | logbx001 logb  0                 -> -Infinity  Division_by_zero
 | 
 | 33 | logbx002 logb  1E-999            -> -999
 | 
 | 34 | logbx003 logb  9E-999            -> -999
 | 
 | 35 | logbx004 logb  0.001             -> -3
 | 
 | 36 | logbx005 logb  0.03              -> -2
 | 
 | 37 | logbx006 logb  1                 ->  0
 | 
 | 38 | logbx007 logb  2                 ->  0
 | 
 | 39 | logbx008 logb  2.5               ->  0
 | 
 | 40 | logbx009 logb  2.50              ->  0
 | 
 | 41 | logbx010 logb  10                ->  1
 | 
 | 42 | logbx011 logb  70                ->  1
 | 
 | 43 | logbx012 logb  100               ->  2
 | 
 | 44 | logbx013 logb  250               ->  2
 | 
 | 45 | logbx014 logb +Infinity          ->  Infinity
 | 
 | 46 | 
 | 
 | 47 | -- negatives are treated as positives
 | 
 | 48 | logbx021 logb -0                 -> -Infinity  Division_by_zero
 | 
 | 49 | logbx022 logb -1E-999            -> -999
 | 
 | 50 | logbx023 logb -9E-999            -> -999
 | 
 | 51 | logbx024 logb -0.001             -> -3
 | 
 | 52 | logbx025 logb -1                 ->  0
 | 
 | 53 | logbx026 logb -2                 ->  0
 | 
 | 54 | logbx027 logb -10                ->  1
 | 
 | 55 | logbx028 logb -70                ->  1
 | 
 | 56 | logbx029 logb -100               ->  2
 | 
 | 57 | logbx030 logb -100000000         ->  8
 | 
 | 58 | logbx031 logb -Infinity          ->  Infinity
 | 
 | 59 | 
 | 
 | 60 | -- zeros
 | 
 | 61 | logbx111 logb          0   -> -Infinity  Division_by_zero
 | 
 | 62 | logbx112 logb         -0   -> -Infinity  Division_by_zero
 | 
 | 63 | logbx113 logb       0E+4   -> -Infinity  Division_by_zero
 | 
 | 64 | logbx114 logb      -0E+4   -> -Infinity  Division_by_zero
 | 
 | 65 | logbx115 logb     0.0000   -> -Infinity  Division_by_zero
 | 
 | 66 | logbx116 logb    -0.0000   -> -Infinity  Division_by_zero
 | 
 | 67 | logbx117 logb      0E-141  -> -Infinity  Division_by_zero
 | 
 | 68 | logbx118 logb     -0E-141  -> -Infinity  Division_by_zero
 | 
 | 69 | 
 | 
 | 70 | -- full coefficients, alternating bits
 | 
 | 71 | logbx121 logb   268268268        -> 8
 | 
 | 72 | logbx122 logb  -268268268        -> 8
 | 
 | 73 | logbx123 logb   134134134        -> 8
 | 
 | 74 | logbx124 logb  -134134134        -> 8
 | 
 | 75 | 
 | 
 | 76 | -- Nmax, Nmin, Ntiny
 | 
 | 77 | logbx131 logb  9.99999999E+999   -> 999
 | 
 | 78 | logbx132 logb  1E-999            -> -999
 | 
 | 79 | logbx133 logb  1.00000000E-999   -> -999
 | 
 | 80 | logbx134 logb  1E-1007           -> -1007
 | 
 | 81 | 
 | 
 | 82 | logbx135 logb  -1E-1007          -> -1007
 | 
 | 83 | logbx136 logb  -1.00000000E-999  -> -999
 | 
 | 84 | logbx137 logb  -1E-999           -> -999
 | 
 | 85 | logbx138 logb  -9.99999999E+999  ->  999
 | 
 | 86 | 
 | 
 | 87 | -- ones
 | 
 | 88 | logbx0061 logb  1                 ->   0
 | 
 | 89 | logbx0062 logb  1.0               ->   0
 | 
 | 90 | logbx0063 logb  1.000000000000000 ->   0
 | 
 | 91 | logbx0064 logb  1.000000000000000000 ->   0
 | 
 | 92 | 
 | 
 | 93 | -- notable cases -- exact powers of 10
 | 
 | 94 | logbx1100 logb 1             -> 0
 | 
 | 95 | logbx1101 logb 10            -> 1
 | 
 | 96 | logbx1102 logb 100           -> 2
 | 
 | 97 | logbx1103 logb 1000          -> 3
 | 
 | 98 | logbx1104 logb 10000         -> 4
 | 
 | 99 | logbx1105 logb 100000        -> 5
 | 
 | 100 | logbx1106 logb 1000000       -> 6
 | 
 | 101 | logbx1107 logb 10000000      -> 7
 | 
 | 102 | logbx1108 logb 100000000     -> 8
 | 
 | 103 | logbx1109 logb 1000000000    -> 9
 | 
 | 104 | logbx1110 logb 10000000000   -> 10
 | 
 | 105 | logbx1111 logb 100000000000  -> 11
 | 
 | 106 | logbx1112 logb 1000000000000 -> 12
 | 
 | 107 | logbx1113 logb 0.00000000001 -> -11
 | 
 | 108 | logbx1114 logb 0.0000000001 -> -10
 | 
 | 109 | logbx1115 logb 0.000000001 -> -9
 | 
 | 110 | logbx1116 logb 0.00000001 -> -8
 | 
 | 111 | logbx1117 logb 0.0000001 -> -7
 | 
 | 112 | logbx1118 logb 0.000001 -> -6
 | 
 | 113 | logbx1119 logb 0.00001 -> -5
 | 
 | 114 | logbx1120 logb 0.0001 -> -4
 | 
 | 115 | logbx1121 logb 0.001 -> -3
 | 
 | 116 | logbx1122 logb 0.01 -> -2
 | 
 | 117 | logbx1123 logb 0.1 -> -1
 | 
 | 118 | logbx1124 logb 1E-99  -> -99
 | 
 | 119 | logbx1125 logb 1E-100 -> -100
 | 
 | 120 | logbx1126 logb 1E-383 -> -383
 | 
 | 121 | logbx1127 logb 1E-999 -> -999
 | 
 | 122 | 
 | 
 | 123 | -- suggestions from Ilan Nehama
 | 
 | 124 | logbx1400 logb 10E-3    -> -2
 | 
 | 125 | logbx1401 logb 10E-2    -> -1
 | 
 | 126 | logbx1402 logb 100E-2   ->  0
 | 
 | 127 | logbx1403 logb 1000E-2  ->  1
 | 
 | 128 | logbx1404 logb 10000E-2 ->  2
 | 
 | 129 | logbx1405 logb 10E-1    ->  0
 | 
 | 130 | logbx1406 logb 100E-1   ->  1
 | 
 | 131 | logbx1407 logb 1000E-1  ->  2
 | 
 | 132 | logbx1408 logb 10000E-1 ->  3
 | 
 | 133 | logbx1409 logb 10E0     ->  1
 | 
 | 134 | logbx1410 logb 100E0    ->  2
 | 
 | 135 | logbx1411 logb 1000E0   ->  3
 | 
 | 136 | logbx1412 logb 10000E0  ->  4
 | 
 | 137 | logbx1413 logb 10E1     ->  2
 | 
 | 138 | logbx1414 logb 100E1    ->  3
 | 
 | 139 | logbx1415 logb 1000E1   ->  4
 | 
 | 140 | logbx1416 logb 10000E1  ->  5
 | 
 | 141 | logbx1417 logb 10E2     ->  3
 | 
 | 142 | logbx1418 logb 100E2    ->  4
 | 
 | 143 | logbx1419 logb 1000E2   ->  5
 | 
 | 144 | logbx1420 logb 10000E2  ->  6
 | 
 | 145 | 
 | 
| Mark Dickinson | 8a54653 | 2009-10-08 16:30:38 +0000 | [diff] [blame] | 146 | -- inexacts
 | 
 | 147 | precision: 2
 | 
 | 148 | logbx1500 logb 10000E2       ->  6
 | 
 | 149 | logbx1501 logb 1E+99         ->  99
 | 
 | 150 | logbx1502 logb 1E-99         -> -99
 | 
 | 151 | logbx1503 logb 1E+100        ->  1.0E+2  Rounded
 | 
 | 152 | logbx1504 logb 1E+999        ->  1.0E+3  Inexact Rounded
 | 
 | 153 | logbx1505 logb 1E-100        -> -1.0E+2  Rounded
 | 
 | 154 | logbx1506 logb 1E-999        -> -1.0E+3  Inexact Rounded
 | 
 | 155 | logbx1507 logb 1E-1111       -> -1.1E+3  Inexact Rounded
 | 
 | 156 | logbx1508 logb 1E-3333       -> -3.3E+3  Inexact Rounded
 | 
 | 157 | logbx1509 logb 1E-6666       -> -6.7E+3  Inexact Rounded
 | 
 | 158 | logbx1510 logb 1E+999999999  ->  1.0E+9  Inexact Rounded
 | 
 | 159 | logbx1511 logb 1E-999999999  -> -1.0E+9  Inexact Rounded
 | 
 | 160 | precision: 1
 | 
 | 161 | logbx1517 logb 1E-1111       -> -1E+3    Inexact Rounded
 | 
 | 162 | logbx1518 logb 1E-3333       -> -3E+3    Inexact Rounded
 | 
 | 163 | logbx1519 logb 1E-6666       -> -7E+3    Inexact Rounded
 | 
 | 164 | precision: 8
 | 
 | 165 | logbx1520 logb 1E+999999999  ->  1.0000000E+9 Inexact Rounded
 | 
 | 166 | logbx1521 logb 1E-999999999  -> -1.0000000E+9 Inexact Rounded
 | 
 | 167 | precision: 9
 | 
 | 168 | logbx1523 logb 1E+999999999  ->  999999999
 | 
 | 169 | logbx1524 logb 1E-999999999  -> -999999999
 | 
 | 170 | 
 | 
| Thomas Wouters | 1b7f891 | 2007-09-19 03:06:30 +0000 | [diff] [blame] | 171 | -- special values
 | 
| Mark Dickinson | 8a54653 | 2009-10-08 16:30:38 +0000 | [diff] [blame] | 172 | precision: 9
 | 
| Thomas Wouters | 1b7f891 | 2007-09-19 03:06:30 +0000 | [diff] [blame] | 173 | logbx820  logb   Infinity ->   Infinity
 | 
 | 174 | logbx821  logb  -Infinity ->   Infinity
 | 
 | 175 | logbx822  logb   0        ->  -Infinity Division_by_zero
 | 
 | 176 | logbx823  logb   NaN      ->   NaN
 | 
 | 177 | logbx824  logb   sNaN     ->   NaN     Invalid_operation
 | 
 | 178 | -- propagating NaNs
 | 
 | 179 | logbx825  logb   sNaN123  ->   NaN123  Invalid_operation
 | 
 | 180 | logbx826  logb   -sNaN321 ->  -NaN321  Invalid_operation
 | 
 | 181 | logbx827  logb   NaN456   ->   NaN456
 | 
 | 182 | logbx828  logb   -NaN654  ->  -NaN654
 | 
 | 183 | logbx829  logb   NaN1     ->   NaN1
 | 
 | 184 | 
 | 
 | 185 | -- Null test
 | 
 | 186 | logbx900  logb #   -> NaN Invalid_operation
 | 
 | 187 | 
 | 
 | 188 | 
 |