t/lib/bigfltpm.t

   1 #!./perl
   2
   3 BEGIN {
   4     chdir 't' if -d 't';
   5     @INC = '../lib';
   6 }
   7
   8 use Math::BigFloat;
   9
  10 $test = 0;
  11 $| = 1;
  12 print "1..358\n";
  13 while (<DATA>) {
  14        chop;
  15        if (s/^&//) {
  16                $f = $_;
  17         } elsif (/^\$.*/) {
  18                 eval "$_;";
  19        } else {
  20                ++$test;
  21                if (m|^(.*?):(/.+)$|) {
  22                    $ans = $2;
  23                    @args = split(/:/,$1,99);
  24                }
  25                else {
  26                    @args = split(/:/,$_,99);
  27                    $ans = pop(@args);
  28                }
  29                $try = "\$x = new Math::BigFloat \"$args[0]\";";
  30                if ($f eq "fnorm"){
  31                    $try .= "\$x+0;";
  32                } elsif ($f eq "fneg") {
  33                    $try .= "-\$x;";
  34                } elsif ($f eq "fabs") {
  35                    $try .= "abs \$x;";
  36                } elsif ($f eq "fround") {
  37                    $try .= "0+\$x->fround($args[1]);";
  38                } elsif ($f eq "ffround") {
  39                    $try .= "0+\$x->ffround($args[1]);";
  40                } elsif ($f eq "fsqrt") {
  41                    $try .= "0+\$x->fsqrt;";
  42                } else {
  43                    $try .= "\$y = new Math::BigFloat \"$args[1]\";";
  44                    if ($f eq fcmp){
  45                        $try .= "\$x <=> \$y;";
  46                    }elsif ($f eq fadd){
  47                        $try .= "\$x + \$y;";
  48                    }elsif ($f eq fsub){
  49                        $try .= "\$x - \$y;";
  50                    }elsif ($f eq fmul){
  51                        $try .= "\$x * \$y;";
  52                    }elsif ($f eq fdiv){
  53                        $try .= "\$x / \$y;";
  54                    } else { warn "Unknown op"; }
  55                }
  56                #print ">>>",$try,"<<<\n";
  57                $ans1 = eval $try;
  58                if ($ans =~ m|^/(.*)$|) {
  59                    my $pat = $1;
  60                    if ($ans1 =~ /$pat/) {
  61                        print "ok $test\n";
  62                    }
  63                    else {
  64                        print "not ok $test\n";
  65                        print "# '$try' expected: /$pat/ got: '$ans1'\n";
  66                    }
  67                }
  68                elsif ("$ans1" eq $ans) { #bug!
  69                        print "ok $test\n";
  70                } else {
  71                        print "not ok $test\n";
  72                        print "# '$try' expected: '$ans' got: '$ans1'\n";
  73                }
  74        }
  75 }
  76 __END__
  77 &fnorm
  78 abc:NaN.
  79    1 a:NaN.
  80 1bcd2:NaN.
  81 11111b:NaN.
  82 +1z:NaN.
  83 -1z:NaN.
  84 0:0.
  85 +0:0.
  86 +00:0.
  87 +0 0 0:0.
  88 000000  0000000   00000:0.
  89 -0:0.
  90 -0000:0.
  91 +1:1.
  92 +01:1.
  93 +001:1.
  94 +00000100000:100000.
  95 123456789:123456789.
  96 -1:-1.
  97 -01:-1.
  98 -001:-1.
  99 -123456789:-123456789.
 100 -00000100000:-100000.
 101 123.456a:NaN.
 102 123.456:123.456
 103 0.01:.01
 104 .002:.002
 105 -0.0003:-.0003
 106 -.0000000004:-.0000000004
 107 123456E2:12345600.
 108 123456E-2:1234.56
 109 -123456E2:-12345600.
 110 -123456E-2:-1234.56
 111 1e1:10.
 112 2e-11:.00000000002
 113 -3e111:-3000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000.
 114 -4e-1111:-.0000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000004
 115 &fneg
 116 abd:NaN.
 117 +0:0.
 118 +1:-1.
 119 -1:1.
 120 +123456789:-123456789.
 121 -123456789:123456789.
 122 +123.456789:-123.456789
 123 -123456.789:123456.789
 124 &fabs
 125 abc:NaN.
 126 +0:0.
 127 +1:1.
 128 -1:1.
 129 +123456789:123456789.
 130 -123456789:123456789.
 131 +123.456789:123.456789
 132 -123456.789:123456.789
 133 &fround
 134 $Math::BigFloat::rnd_mode = 'trunc'
 135 +10123456789:5:10123000000
 136 -10123456789:5:-10123000000
 137 +10123456789:9:10123456700
 138 -10123456789:9:-10123456700
 139 +101234500:6:101234000
 140 -101234500:6:-101234000
 141 $Math::BigFloat::rnd_mode = 'zero'
 142 +20123456789:5:20123000000
 143 -20123456789:5:-20123000000
 144 +20123456789:9:20123456800
 145 -20123456789:9:-20123456800
 146 +201234500:6:201234000
 147 -201234500:6:-201234000
 148 $Math::BigFloat::rnd_mode = '+inf'
 149 +30123456789:5:30123000000
 150 -30123456789:5:-30123000000
 151 +30123456789:9:30123456800
 152 -30123456789:9:-30123456800
 153 +301234500:6:301235000
 154 -301234500:6:-301234000
 155 $Math::BigFloat::rnd_mode = '-inf'
 156 +40123456789:5:40123000000
 157 -40123456789:5:-40123000000
 158 +40123456789:9:40123456800
 159 -40123456789:9:-40123456800
 160 +401234500:6:401234000
 161 -401234500:6:-401235000
 162 $Math::BigFloat::rnd_mode = 'odd'
 163 +50123456789:5:50123000000
 164 -50123456789:5:-50123000000
 165 +50123456789:9:50123456800
 166 -50123456789:9:-50123456800
 167 +501234500:6:501235000
 168 -501234500:6:-501235000
 169 $Math::BigFloat::rnd_mode = 'even'
 170 +60123456789:5:60123000000
 171 -60123456789:5:-60123000000
 172 +60123456789:9:60123456800
 173 -60123456789:9:-60123456800
 174 +601234500:6:601234000
 175 -601234500:6:-601234000
 176 &ffround
 177 $Math::BigFloat::rnd_mode = 'trunc'
 178 +1.23:-1:1.2
 179 -1.23:-1:-1.2
 180 +1.27:-1:1.2
 181 -1.27:-1:-1.2
 182 +1.25:-1:1.2
 183 -1.25:-1:-1.2
 184 +1.35:-1:1.3
 185 -1.35:-1:-1.3
 186 -0.006:-1:0
 187 -0.006:-2:0
 188 -0.0065:-3:/-0\.006|-6e-03
 189 -0.0065:-4:/-0\.0065|-6\.5e-03
 190 -0.0065:-5:/-0\.0065|-6\.5e-03
 191 $Math::BigFloat::rnd_mode = 'zero'
 192 +2.23:-1:2.2
 193 -2.23:-1:-2.2
 194 +2.27:-1:2.3
 195 -2.27:-1:-2.3
 196 +2.25:-1:2.2
 197 -2.25:-1:-2.2
 198 +2.35:-1:2.3
 199 -2.35:-1:-2.3
 200 -0.0065:-1:0
 201 -0.0065:-2:/-0\.01|-1e-02
 202 -0.0065:-3:/-0\.006|-6e-03
 203 -0.0065:-4:/-0\.0065|-6\.5e-03
 204 -0.0065:-5:/-0\.0065|-6\.5e-03
 205 $Math::BigFloat::rnd_mode = '+inf'
 206 +3.23:-1:3.2
 207 -3.23:-1:-3.2
 208 +3.27:-1:3.3
 209 -3.27:-1:-3.3
 210 +3.25:-1:3.3
 211 -3.25:-1:-3.2
 212 +3.35:-1:3.4
 213 -3.35:-1:-3.3
 214 -0.0065:-1:0
 215 -0.0065:-2:/-0\.01|-1e-02
 216 -0.0065:-3:/-0\.006|-6e-03
 217 -0.0065:-4:/-0\.0065|-6\.5e-03
 218 -0.0065:-5:/-0\.0065|-6\.5e-03
 219 $Math::BigFloat::rnd_mode = '-inf'
 220 +4.23:-1:4.2
 221 -4.23:-1:-4.2
 222 +4.27:-1:4.3
 223 -4.27:-1:-4.3
 224 +4.25:-1:4.2
 225 -4.25:-1:-4.3
 226 +4.35:-1:4.3
 227 -4.35:-1:-4.4
 228 -0.0065:-1:0
 229 -0.0065:-2:/-0\.01|-1e-02
 230 -0.0065:-3:/-0\.007|-7e-03
 231 -0.0065:-4:/-0\.0065|-6\.5e-03
 232 -0.0065:-5:/-0\.0065|-6\.5e-03
 233 $Math::BigFloat::rnd_mode = 'odd'
 234 +5.23:-1:5.2
 235 -5.23:-1:-5.2
 236 +5.27:-1:5.3
 237 -5.27:-1:-5.3
 238 +5.25:-1:5.3
 239 -5.25:-1:-5.3
 240 +5.35:-1:5.3
 241 -5.35:-1:-5.3
 242 -0.0065:-1:0
 243 -0.0065:-2:/-0\.01|-1e-02
 244 -0.0065:-3:/-0\.007|-7e-03
 245 -0.0065:-4:/-0\.0065|-6\.5e-03
 246 -0.0065:-5:/-0\.0065|-6\.5e-03
 247 $Math::BigFloat::rnd_mode = 'even'
 248 +6.23:-1:6.2
 249 -6.23:-1:-6.2
 250 +6.27:-1:6.3
 251 -6.27:-1:-6.3
 252 +6.25:-1:6.2
 253 -6.25:-1:-6.2
 254 +6.35:-1:6.4
 255 -6.35:-1:-6.4
 256 -0.0065:-1:0
 257 -0.0065:-2:/-0\.01|-1e-02
 258 -0.0065:-3:/-0\.006|-6e-03
 259 -0.0065:-4:/-0\.0065|-6\.5e-03
 260 -0.0065:-5:/-0\.0065|-6\.5e-03
 261 &fcmp
 262 abc:abc:
 263 abc:+0:
 264 +0:abc:
 265 +0:+0:0
 266 -1:+0:-1
 267 +0:-1:1
 268 +1:+0:1
 269 +0:+1:-1
 270 -1:+1:-1
 271 +1:-1:1
 272 -1:-1:0
 273 +1:+1:0
 274 +123:+123:0
 275 +123:+12:1
 276 +12:+123:-1
 277 -123:-123:0
 278 -123:-12:-1
 279 -12:-123:1
 280 +123:+124:-1
 281 +124:+123:1
 282 -123:-124:1
 283 -124:-123:-1
 284 &fadd
 285 abc:abc:NaN.
 286 abc:+0:NaN.
 287 +0:abc:NaN.
 288 +0:+0:0.
 289 +1:+0:1.
 290 +0:+1:1.
 291 +1:+1:2.
 292 -1:+0:-1.
 293 +0:-1:-1.
 294 -1:-1:-2.
 295 -1:+1:0.
 296 +1:-1:0.
 297 +9:+1:10.
 298 +99:+1:100.
 299 +999:+1:1000.
 300 +9999:+1:10000.
 301 +99999:+1:100000.
 302 +999999:+1:1000000.
 303 +9999999:+1:10000000.
 304 +99999999:+1:100000000.
 305 +999999999:+1:1000000000.
 306 +9999999999:+1:10000000000.
 307 +99999999999:+1:100000000000.
 308 +10:-1:9.
 309 +100:-1:99.
 310 +1000:-1:999.
 311 +10000:-1:9999.
 312 +100000:-1:99999.
 313 +1000000:-1:999999.
 314 +10000000:-1:9999999.
 315 +100000000:-1:99999999.
 316 +1000000000:-1:999999999.
 317 +10000000000:-1:9999999999.
 318 +123456789:+987654321:1111111110.
 319 -123456789:+987654321:864197532.
 320 -123456789:-987654321:-1111111110.
 321 +123456789:-987654321:-864197532.
 322 &fsub
 323 abc:abc:NaN.
 324 abc:+0:NaN.
 325 +0:abc:NaN.
 326 +0:+0:0.
 327 +1:+0:1.
 328 +0:+1:-1.
 329 +1:+1:0.
 330 -1:+0:-1.
 331 +0:-1:1.
 332 -1:-1:0.
 333 -1:+1:-2.
 334 +1:-1:2.
 335 +9:+1:8.
 336 +99:+1:98.
 337 +999:+1:998.
 338 +9999:+1:9998.
 339 +99999:+1:99998.
 340 +999999:+1:999998.
 341 +9999999:+1:9999998.
 342 +99999999:+1:99999998.
 343 +999999999:+1:999999998.
 344 +9999999999:+1:9999999998.
 345 +99999999999:+1:99999999998.
 346 +10:-1:11.
 347 +100:-1:101.
 348 +1000:-1:1001.
 349 +10000:-1:10001.
 350 +100000:-1:100001.
 351 +1000000:-1:1000001.
 352 +10000000:-1:10000001.
 353 +100000000:-1:100000001.
 354 +1000000000:-1:1000000001.
 355 +10000000000:-1:10000000001.
 356 +123456789:+987654321:-864197532.
 357 -123456789:+987654321:-1111111110.
 358 -123456789:-987654321:864197532.
 359 +123456789:-987654321:1111111110.
 360 &fmul
 361 abc:abc:NaN.
 362 abc:+0:NaN.
 363 +0:abc:NaN.
 364 +0:+0:0.
 365 +0:+1:0.
 366 +1:+0:0.
 367 +0:-1:0.
 368 -1:+0:0.
 369 +123456789123456789:+0:0.
 370 +0:+123456789123456789:0.
 371 -1:-1:1.
 372 -1:+1:-1.
 373 +1:-1:-1.
 374 +1:+1:1.
 375 +2:+3:6.
 376 -2:+3:-6.
 377 +2:-3:-6.
 378 -2:-3:6.
 379 +111:+111:12321.
 380 +10101:+10101:102030201.
 381 +1001001:+1001001:1002003002001.
 382 +100010001:+100010001:10002000300020001.
 383 +10000100001:+10000100001:100002000030000200001.
 384 +11111111111:+9:99999999999.
 385 +22222222222:+9:199999999998.
 386 +33333333333:+9:299999999997.
 387 +44444444444:+9:399999999996.
 388 +55555555555:+9:499999999995.
 389 +66666666666:+9:599999999994.
 390 +77777777777:+9:699999999993.
 391 +88888888888:+9:799999999992.
 392 +99999999999:+9:899999999991.
 393 &fdiv
 394 abc:abc:NaN.
 395 abc:+1:abc:NaN.
 396 +1:abc:NaN.
 397 +0:+0:NaN.
 398 +0:+1:0.
 399 +1:+0:NaN.
 400 +0:-1:0.
 401 -1:+0:NaN.
 402 +1:+1:1.
 403 -1:-1:1.
 404 +1:-1:-1.
 405 -1:+1:-1.
 406 +1:+2:.5
 407 +2:+1:2.
 408 +10:+5:2.
 409 +100:+4:25.
 410 +1000:+8:125.
 411 +10000:+16:625.
 412 +10000:-16:-625.
 413 +999999999999:+9:111111111111.
 414 +999999999999:+99:10101010101.
 415 +999999999999:+999:1001001001.
 416 +999999999999:+9999:100010001.
 417 +999999999999999:+99999:10000100001.
 418 +1000000000:+9:111111111.1111111111111111111111111111111
 419 +2000000000:+9:222222222.2222222222222222222222222222222
 420 +3000000000:+9:333333333.3333333333333333333333333333333
 421 +4000000000:+9:444444444.4444444444444444444444444444444
 422 +5000000000:+9:555555555.5555555555555555555555555555556
 423 +6000000000:+9:666666666.6666666666666666666666666666667
 424 +7000000000:+9:777777777.7777777777777777777777777777778
 425 +8000000000:+9:888888888.8888888888888888888888888888889
 426 +9000000000:+9:1000000000.
 427 +35500000:+113:314159.2920353982300884955752212389380531
 428 +71000000:+226:314159.2920353982300884955752212389380531
 429 +106500000:+339:314159.2920353982300884955752212389380531
 430 +1000000000:+3:333333333.3333333333333333333333333333333
 431 $Math::BigFloat::div_scale = 20
 432 +1000000000:+9:111111111.11111111111
 433 +2000000000:+9:222222222.22222222222
 434 +3000000000:+9:333333333.33333333333
 435 +4000000000:+9:444444444.44444444444
 436 +5000000000:+9:555555555.55555555556
 437 +6000000000:+9:666666666.66666666667
 438 +7000000000:+9:777777777.77777777778
 439 +8000000000:+9:888888888.88888888889
 440 +9000000000:+9:1000000000.
 441 +35500000:+113:314159.292035398230088
 442 +71000000:+226:314159.292035398230088
 443 +106500000:+339:314159.29203539823009
 444 +1000000000:+3:333333333.33333333333
 445 $Math::BigFloat::div_scale = 40
 446 &fsqrt
 447 +0:0
 448 -1:/^(?i:0|\?|NaNQ?)$
 449 -2:/^(?i:0|\?|NaNQ?)$
 450 -16:/^(?i:0|\?|NaNQ?)$
 451 -123.456:/^(?i:0|\?|NaNQ?)$
 452 +1:1.
 453 +1.44:1.2
 454 +2:1.41421356237309504880168872420969807857
 455 +4:2.
 456 +16:4.
 457 +100:10.
 458 +123.456:11.11107555549866648462149404118219234119
 459 +15241.383936:123.456