1 #include this file into another test for subclass testing...
3 ok ($class->config()->{lib},$CL);
10 $_ =~ s/#.*$//; # remove comments
11 $_ =~ s/\s+$//; # trailing spaces
12 next if /^$/; # skip empty lines & comments
19 $setup = $_; $setup =~ s/\$/\$${class}::/g; # round_mode, div_scale
20 #print "\$setup== $setup\n";
27 @args = split(/:/,$1,99);
31 @args = split(/:/,$_,99); $ans = pop(@args);
33 $try = "\$x = $class->new('$args[0]');";
37 } elsif ($f eq "finf") {
38 $try .= "\$x->finf('$args[1]');";
39 } elsif ($f eq "is_inf") {
40 $try .= "\$x->is_inf('$args[1]');";
41 } elsif ($f eq "fone") {
42 $try .= "\$x->bone('$args[1]');";
43 } elsif ($f eq "fstr") {
44 $try .= "\$x->accuracy($args[1]); \$x->precision($args[2]);";
45 $try .= '$x->fstr();';
46 } elsif ($f eq "parts") {
47 # ->bstr() to see if an object is returned
48 $try .= '($a,$b) = $x->parts(); $a = $a->bstr(); $b = $b->bstr();';
50 } elsif ($f eq "exponent") {
51 # ->bstr() to see if an object is returned
52 $try .= '$x->exponent()->bstr();';
53 } elsif ($f eq "mantissa") {
54 # ->bstr() to see if an object is returned
55 $try .= '$x->mantissa()->bstr();';
56 } elsif ($f =~ /^(numify|length|as_number|as_hex|as_bin)$/) {
58 # some unary ops (test the fxxx form, since that is done by AUTOLOAD)
59 } elsif ($f =~ /^f(nan|sstr|neg|floor|ceil|abs)$/) {
60 $try .= "\$x->f$1();";
61 # some is_xxx test function
62 } elsif ($f =~ /^is_(zero|one|negative|positive|odd|even|nan|int)$/) {
64 } elsif ($f eq "finc") {
66 } elsif ($f eq "fdec") {
68 }elsif ($f eq "fround") {
69 $try .= "$setup; \$x->fround($args[1]);";
70 } elsif ($f eq "ffround") {
71 $try .= "$setup; \$x->ffround($args[1]);";
72 } elsif ($f eq "fsqrt") {
73 $try .= "$setup; \$x->fsqrt();";
74 } elsif ($f eq "ffac") {
75 $try .= "$setup; \$x->ffac();";
76 } elsif ($f eq "flog") {
79 $try .= "\$y = $class->new($args[1]);";
80 $try .= "$setup; \$x->flog(\$y);";
84 $try .= "$setup; \$x->flog();";
89 $try .= "\$y = $class->new(\"$args[1]\");";
92 } elsif ($f eq "facmp") {
93 $try .= '$x->facmp($y);';
94 } elsif ($f eq "fpow") {
96 } elsif ($f eq "froot") {
97 $try .= "$setup; \$x->froot(\$y);";
98 } elsif ($f eq "fadd") {
100 } elsif ($f eq "fsub") {
102 } elsif ($f eq "fmul") {
104 } elsif ($f eq "fdiv") {
105 $try .= "$setup; \$x / \$y;";
106 } elsif ($f eq "fdiv-list") {
107 $try .= "$setup; join(',',\$x->fdiv(\$y));";
108 } elsif ($f eq "frsft") {
110 } elsif ($f eq "flsft") {
112 } elsif ($f eq "fmod") {
114 } else { warn "Unknown op '$f'"; }
116 # print "# Trying: '$try'\n";
118 if ($ans =~ m|^/(.*)$|)
127 print "# '$try' expected: /$pat/ got: '$ans1'\n" if !ok(1,0);
138 print "# Tried: '$try'\n" if !ok ($ans1, $ans);
139 if (ref($ans1) eq "$class")
141 # float numbers are normalized (for now), so mantissa shouldn't have
143 #print $ans1->_trailing_zeros(),"\n";
144 print "# Has trailing zeros after '$try'\n"
145 if !ok ($ans1->{_m}->_trailing_zeros(), 0);
148 } # end pattern or string
152 # check whether $class->new( Math::BigInt->new()) destroys it
153 # ($y == 12 in this case)
154 $x = Math::BigInt->new(1200); $y = $class->new($x);
155 ok ($y,1200); ok ($x,1200);
157 ###############################################################################
158 # Really huge, big, ultra-mega-biggy-monster exponents
159 # Technically, the exponents should not be limited (they are BigInts), but
160 # practically there are a few places were they are limited to a Perl scalar.
161 # This is sometimes for speed, sometimes because otherwise the number wouldn't
162 # fit into your memory (just think of 1e123456789012345678901234567890 + 1!)
163 # anyway. We don't test everything here, but let's make sure it just basically
166 my $monster = '1e1234567890123456789012345678901234567890';
169 ok ($class->new($monster)->bsstr(),
170 '1e+1234567890123456789012345678901234567890');
172 ok ($class->new($monster) > 0,1);
175 ok ($class->new($monster)->bsub( $monster),0);
176 ok ($class->new($monster)->bmul(2)->bsstr(),
177 '2e+1234567890123456789012345678901234567890');
179 ###############################################################################
182 $x = $class->new(2); $x->fzero(); ok_undef ($x->{_a}); ok_undef ($x->{_p});
183 $x = $class->new(2); $x->finf(); ok_undef ($x->{_a}); ok_undef ($x->{_p});
184 $x = $class->new(2); $x->fone(); ok_undef ($x->{_a}); ok_undef ($x->{_p});
185 $x = $class->new(2); $x->fnan(); ok_undef ($x->{_a}); ok_undef ($x->{_p});
187 ###############################################################################
188 # bone/binf etc as plain calls (Lite failed them)
190 ok ($class->fzero(),0);
191 ok ($class->fone(),1);
192 ok ($class->fone('+'),1);
193 ok ($class->fone('-'),-1);
194 ok ($class->fnan(),'NaN');
195 ok ($class->finf(),'inf');
196 ok ($class->finf('+'),'inf');
197 ok ($class->finf('-'),'-inf');
198 ok ($class->finf('-inf'),'-inf');
200 $class->accuracy(undef); $class->precision(undef); # reset
202 ###############################################################################
203 # bug in bsstr()/numify() showed up in after-rounding in bdiv()
205 $x = $class->new('0.008'); $y = $class->new(2);
209 ###############################################################################
210 # fsqrt() with set global A/P or A/P enabled on $x, also a test whether fsqrt()
211 # correctly modifies $x
214 $x = $class->new(12); $class->precision(-2); $x->fsqrt(); ok ($x,'3.46');
216 $class->precision(undef);
217 $x = $class->new(12); $class->precision(0); $x->fsqrt(); ok ($x,'3');
219 $class->precision(-3); $x = $class->new(12); $x->fsqrt(); ok ($x,'3.464');
224 ${${class}.'::accuracy'} = 4; $x = $class->new(12);
225 $x->fsqrt(3); ok ($x,'NaN');
226 # supplied arg overrides set global
227 $class->precision(undef); $x = $class->new(12); $x->fsqrt(3); ok ($x,'3.46');
228 $class->accuracy(undef); $class->precision(undef); # reset for further tests
231 #############################################################################
232 # can we call objectify (broken until v1.52)
237 '@args' . " = $class" . "::objectify(2,$class,4,5);".'join(" ",@args);';
239 ok ($ans,"$class 4 5");
242 #############################################################################
243 # is_one('-') (broken until v1.64)
245 ok ($class->new(-1)->is_one(),0);
246 ok ($class->new(-1)->is_one('-'),1);
250 ###############################################################################
251 # Perl 5.005 does not like ok ($x,undef)
257 ok (1,1) and return if !defined $x;
267 # base > 0, base != 1
271 # log(1) is always 1, regardless of $base
275 # this is too slow for the testsuite
276 #2:0.6931471805599453094172321214581765680755
277 #2.718281828:0.9999999998311266953289851340574956564911
279 #2.718281828:0.99999999983112669533
284 #1000:0:6.90775527898214
285 #100:0:4.60517018598809
287 #3.1415:0:1.14470039286086
289 #12345:0:9.42100640177928
290 #0.001:0:-6.90775527898214
291 # reset for further tests
333 # test for bug in brsft() not handling cases that return 0
378 # uses bsstr() so 5 => 5e+0 to be compatible w/ Perls output
415 1234.567:9::1234.56700
416 1234.567::-6:1234.567000
418 0.001234:6::0.00123400
419 0.001234::-8:0.00123400
438 000000_0000000_00000:0
449 -123456789:-123456789
457 -.0000000004:-0.0000000004
472 -3e111:-3000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000
473 -4e-1111:-0.0000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000004
480 123.456:2:15241.383936
483 128:-2:0.00006103515625
490 # 2 ** 0.5 == sqrt(2)
491 # 1.41..7 and not 1.4170 since fallback (bsqrt(9) is '3', not 3.0...0)
492 2:0.5:1.41421356237309504880168872420969807857
493 #2:0.2:1.148698354997035006798626946777927589444
494 #6:1.5:14.6969384566990685891837044482353483518
496 #62.5:12.5:26447206647554886213592.3959144
505 +123456789:-123456789
507 +123.456789:-123.456789
508 -123456.789:123456.789
518 +123.456789:123.456789
519 -123456.789:123456.789
521 $round_mode = "trunc"
526 +10123456789:5:10123000000
527 -10123456789:5:-10123000000
528 +10123456789.123:5:10123000000
529 -10123456789.123:5:-10123000000
530 +10123456789:9:10123456700
531 -10123456789:9:-10123456700
532 +101234500:6:101234000
533 -101234500:6:-101234000
535 +20123456789:5:20123000000
536 -20123456789:5:-20123000000
537 +20123456789.123:5:20123000000
538 -20123456789.123:5:-20123000000
539 +20123456789:9:20123456800
540 -20123456789:9:-20123456800
541 +201234500:6:201234000
542 -201234500:6:-201234000
544 +30123456789:5:30123000000
545 -30123456789:5:-30123000000
546 +30123456789.123:5:30123000000
547 -30123456789.123:5:-30123000000
548 +30123456789:9:30123456800
549 -30123456789:9:-30123456800
550 +301234500:6:301235000
551 -301234500:6:-301234000
553 +40123456789:5:40123000000
554 -40123456789:5:-40123000000
555 +40123456789.123:5:40123000000
556 -40123456789.123:5:-40123000000
557 +40123456789:9:40123456800
558 -40123456789:9:-40123456800
559 +401234500:6:401234000
560 -401234500:6:-401235000
562 +50123456789:5:50123000000
563 -50123456789:5:-50123000000
564 +50123456789.123:5:50123000000
565 -50123456789.123:5:-50123000000
566 +50123456789:9:50123456800
567 -50123456789:9:-50123456800
568 +501234500:6:501235000
569 -501234500:6:-501235000
571 +60123456789:5:60123000000
572 -60123456789:5:-60123000000
573 +60123456789:9:60123456800
574 -60123456789:9:-60123456800
575 +601234500:6:601234000
576 -601234500:6:-601234000
577 +60123456789.0123:5:60123000000
578 -60123456789.0123:5:-60123000000
580 $round_mode = "trunc"
601 -0.0061234567890:-1:0.0
609 -0.0065:-3:/-0\.006|-6e-03
610 -0.0065:-4:/-0\.006(?:5|49{5}\d+)|-6\.5e-03
611 -0.0065:-5:/-0\.006(?:5|49{5}\d+)|-6\.5e-03
617 +2.23:-1:/2.2(?:0{5}\d+)?
618 -2.23:-1:/-2.2(?:0{5}\d+)?
619 +2.27:-1:/2.(?:3|29{5}\d+)
620 -2.27:-1:/-2.(?:3|29{5}\d+)
621 +2.25:-1:/2.2(?:0{5}\d+)?
622 -2.25:-1:/-2.2(?:0{5}\d+)?
623 +2.35:-1:/2.(?:3|29{5}\d+)
624 -2.35:-1:/-2.(?:3|29{5}\d+)
626 -0.0065:-2:/-0\.01|-1e-02
627 -0.0065:-3:/-0\.006|-6e-03
628 -0.0065:-4:/-0\.006(?:5|49{5}\d+)|-6\.5e-03
629 -0.0065:-5:/-0\.006(?:5|49{5}\d+)|-6\.5e-03
635 +3.23:-1:/3.2(?:0{5}\d+)?
636 -3.23:-1:/-3.2(?:0{5}\d+)?
637 +3.27:-1:/3.(?:3|29{5}\d+)
638 -3.27:-1:/-3.(?:3|29{5}\d+)
639 +3.25:-1:/3.(?:3|29{5}\d+)
640 -3.25:-1:/-3.2(?:0{5}\d+)?
641 +3.35:-1:/3.(?:4|39{5}\d+)
642 -3.35:-1:/-3.(?:3|29{5}\d+)
644 -0.0065:-2:/-0\.01|-1e-02
645 -0.0065:-3:/-0\.006|-6e-03
646 -0.0065:-4:/-0\.006(?:5|49{5}\d+)|-6\.5e-03
647 -0.0065:-5:/-0\.006(?:5|49{5}\d+)|-6\.5e-03
653 +4.23:-1:/4.2(?:0{5}\d+)?
654 -4.23:-1:/-4.2(?:0{5}\d+)?
655 +4.27:-1:/4.(?:3|29{5}\d+)
656 -4.27:-1:/-4.(?:3|29{5}\d+)
657 +4.25:-1:/4.2(?:0{5}\d+)?
658 -4.25:-1:/-4.(?:3|29{5}\d+)
659 +4.35:-1:/4.(?:3|29{5}\d+)
660 -4.35:-1:/-4.(?:4|39{5}\d+)
662 -0.0065:-2:/-0\.01|-1e-02
663 -0.0065:-3:/-0\.007|-7e-03
664 -0.0065:-4:/-0\.006(?:5|49{5}\d+)|-6\.5e-03
665 -0.0065:-5:/-0\.006(?:5|49{5}\d+)|-6\.5e-03
671 +5.23:-1:/5.2(?:0{5}\d+)?
672 -5.23:-1:/-5.2(?:0{5}\d+)?
673 +5.27:-1:/5.(?:3|29{5}\d+)
674 -5.27:-1:/-5.(?:3|29{5}\d+)
675 +5.25:-1:/5.(?:3|29{5}\d+)
676 -5.25:-1:/-5.(?:3|29{5}\d+)
677 +5.35:-1:/5.(?:3|29{5}\d+)
678 -5.35:-1:/-5.(?:3|29{5}\d+)
680 -0.0065:-2:/-0\.01|-1e-02
681 -0.0065:-3:/-0\.007|-7e-03
682 -0.0065:-4:/-0\.006(?:5|49{5}\d+)|-6\.5e-03
683 -0.0065:-5:/-0\.006(?:5|49{5}\d+)|-6\.5e-03
689 +6.23:-1:/6.2(?:0{5}\d+)?
690 -6.23:-1:/-6.2(?:0{5}\d+)?
691 +6.27:-1:/6.(?:3|29{5}\d+)
692 -6.27:-1:/-6.(?:3|29{5}\d+)
693 +6.25:-1:/6.(?:2(?:0{5}\d+)?|29{5}\d+)
694 -6.25:-1:/-6.(?:2(?:0{5}\d+)?|29{5}\d+)
695 +6.35:-1:/6.(?:4|39{5}\d+|29{8}\d+)
696 -6.35:-1:/-6.(?:4|39{5}\d+|29{8}\d+)
698 -0.0065:-2:/-0\.01|-1e-02
699 -0.0065:-3:/-0\.006|-7e-03
700 -0.0065:-4:/-0\.006(?:5|49{5}\d+)|-6\.5e-03
701 -0.0065:-5:/-0\.006(?:5|49{5}\d+)|-6\.5e-03
708 0.01234567:-5:0.01235
709 0.01234567:-6:0.012346
710 0.01234567:-7:0.0123457
711 0.01234567:-8:0.01234567
712 0.01234567:-9:0.012345670
713 0.01234567:-12:0.012345670000
915 +99999999:+1:100000000
916 +999999999:+1:1000000000
917 +9999999999:+1:10000000000
918 +99999999999:+1:100000000000
926 +100000000:-1:99999999
927 +1000000000:-1:999999999
928 +10000000000:-1:9999999999
929 +123456789:+987654321:1111111110
930 -123456789:+987654321:864197532
931 -123456789:-987654321:-1111111110
932 +123456789:-987654321:-864197532
933 0.001234:0.0001234:0.0013574
962 +99999999:+1:99999998
963 +999999999:+1:999999998
964 +9999999999:+1:9999999998
965 +99999999999:+1:99999999998
972 +10000000:-1:10000001
973 +100000000:-1:100000001
974 +1000000000:-1:1000000001
975 +10000000000:-1:10000000001
976 +123456789:+987654321:-864197532
977 -123456789:+987654321:-1111111110
978 -123456789:-987654321:864197532
979 +123456789:-987654321:1111111110
1005 +123456789123456789:+0:0
1006 +0:+123456789123456789:0
1016 +10101:+10101:102030201
1017 +1001001:+1001001:1002003002001
1018 +100010001:+100010001:10002000300020001
1019 +10000100001:+10000100001:100002000030000200001
1020 +11111111111:+9:99999999999
1021 +22222222222:+9:199999999998
1022 +33333333333:+9:299999999997
1023 +44444444444:+9:399999999996
1024 +55555555555:+9:499999999995
1025 +66666666666:+9:599999999994
1026 +77777777777:+9:699999999993
1027 +88888888888:+9:799999999992
1028 +99999999999:+9:899999999991
1037 $div_scale = 40; $round_mode = 'even'
1063 +999999999999:+9:111111111111
1064 +999999999999:+99:10101010101
1065 +999999999999:+999:1001001001
1066 +999999999999:+9999:100010001
1067 +999999999999999:+99999:10000100001
1068 +1000000000:+9:111111111.1111111111111111111111111111111
1069 +2000000000:+9:222222222.2222222222222222222222222222222
1070 +3000000000:+9:333333333.3333333333333333333333333333333
1071 +4000000000:+9:444444444.4444444444444444444444444444444
1072 +5000000000:+9:555555555.5555555555555555555555555555556
1073 +6000000000:+9:666666666.6666666666666666666666666666667
1074 +7000000000:+9:777777777.7777777777777777777777777777778
1075 +8000000000:+9:888888888.8888888888888888888888888888889
1076 +9000000000:+9:1000000000
1077 +35500000:+113:314159.2920353982300884955752212389380531
1078 +71000000:+226:314159.2920353982300884955752212389380531
1079 +106500000:+339:314159.2920353982300884955752212389380531
1080 +1000000000:+3:333333333.3333333333333333333333333333333
1081 2:25.024996000799840031993601279744051189762:0.07992009269196593320152084692285869265447
1084 +1000000000:+9:111111111.11111111111
1085 +2000000000:+9:222222222.22222222222
1086 +3000000000:+9:333333333.33333333333
1087 +4000000000:+9:444444444.44444444444
1088 +5000000000:+9:555555555.55555555556
1089 +6000000000:+9:666666666.66666666667
1090 +7000000000:+9:777777777.77777777778
1091 +8000000000:+9:888888888.88888888889
1092 +9000000000:+9:1000000000
1097 1:504:0.001984126984126984127
1098 2:1.987654321:1.0062111801179738436
1099 123456789.123456789123456789123456789:1:123456789.12345678912
1100 # the next two cases are the "old" behaviour, but are now (>v0.01) different
1101 #+35500000:+113:314159.292035398230088
1102 #+71000000:+226:314159.292035398230088
1103 +35500000:+113:314159.29203539823009
1104 +71000000:+226:314159.29203539823009
1105 +106500000:+339:314159.29203539823009
1106 +1000000000:+3:333333333.33333333333
1108 # round to accuracy 1 after bdiv
1110 123456789.1234:1:100000000
1111 # reset scale for further tests
1118 # inf handling, see table in doc
1137 # exceptions to reminder rule
1176 999999999999999:99999:0
1190 152403346:12345:4321
1192 # now some floating point tests
1215 # -$x ** (1/2) => -$y, but not in froot()
1219 2:2:1.41421356237309504880168872420969807857
1225 123.456:2:11.11107555549866648462149404118219234119
1226 15241.38393:2:123.4559999756998444766131352122991626468
1228 12:2:3.464101615137754587054892683011744733886
1269 2:1.41421356237309504880168872420969807857
1274 123.456:11.11107555549866648462149404118219234119
1275 15241.38393:123.4559999756998444766131352122991626468
1277 # sqrt(1.44) = 1.2, sqrt(e10) = e5 => 12e4
1279 2e10:141421.356237309504880168872420969807857
1281 # proved to be an endless loop under 7-9
1282 12:3.464101615137754587054892683011744733886
1301 # it must be exactly /^[+-]inf$/
1397 12345678901234567890:20