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 ref($ans) eq 'HASH' && exists $ans->{_m} && !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
167 #my $monster = '1e1234567890123456789012345678901234567890';
170 #ok ($class->new($monster)->bsstr(),
171 # '1e+1234567890123456789012345678901234567890');
173 #ok ($class->new($monster) > 0,1);
176 #ok ($class->new($monster)->bsub( $monster),0);
177 #ok ($class->new($monster)->bmul(2)->bsstr(),
178 # '2e+1234567890123456789012345678901234567890');
180 ###############################################################################
183 $x = $class->new(2); $x->fzero(); ok_undef ($x->{_a}); ok_undef ($x->{_p});
184 $x = $class->new(2); $x->finf(); ok_undef ($x->{_a}); ok_undef ($x->{_p});
185 $x = $class->new(2); $x->fone(); ok_undef ($x->{_a}); ok_undef ($x->{_p});
186 $x = $class->new(2); $x->fnan(); ok_undef ($x->{_a}); ok_undef ($x->{_p});
188 ###############################################################################
189 # bone/binf etc as plain calls (Lite failed them)
191 ok ($class->fzero(),0);
192 ok ($class->fone(),1);
193 ok ($class->fone('+'),1);
194 ok ($class->fone('-'),-1);
195 ok ($class->fnan(),'NaN');
196 ok ($class->finf(),'inf');
197 ok ($class->finf('+'),'inf');
198 ok ($class->finf('-'),'-inf');
199 ok ($class->finf('-inf'),'-inf');
201 $class->accuracy(undef); $class->precision(undef); # reset
203 ###############################################################################
204 # bug in bsstr()/numify() showed up in after-rounding in bdiv()
206 $x = $class->new('0.008'); $y = $class->new(2);
210 ###############################################################################
211 # fsqrt() with set global A/P or A/P enabled on $x, also a test whether fsqrt()
212 # correctly modifies $x
215 $x = $class->new(12); $class->precision(-2); $x->fsqrt(); ok ($x,'3.46');
217 $class->precision(undef);
218 $x = $class->new(12); $class->precision(0); $x->fsqrt(); ok ($x,'3');
220 $class->precision(-3); $x = $class->new(12); $x->fsqrt(); ok ($x,'3.464');
225 ${${class}.'::accuracy'} = 4; $x = $class->new(12);
226 $x->fsqrt(3); ok ($x,'NaN');
227 # supplied arg overrides set global
228 $class->precision(undef); $x = $class->new(12); $x->fsqrt(3); ok ($x,'3.46');
229 $class->accuracy(undef); $class->precision(undef); # reset for further tests
232 #############################################################################
233 # can we call objectify (broken until v1.52)
238 '@args' . " = $class" . "::objectify(2,$class,4,5);".'join(" ",@args);';
240 ok ($ans,"$class 4 5");
243 #############################################################################
244 # is_one('-') (broken until v1.64)
246 ok ($class->new(-1)->is_one(),0);
247 ok ($class->new(-1)->is_one('-'),1);
251 ###############################################################################
252 # Perl 5.005 does not like ok ($x,undef)
258 ok (1,1) and return if !defined $x;
268 # base > 0, base != 1
272 # log(1) is always 1, regardless of $base
276 # this is too slow for the testsuite
277 #2:0.6931471805599453094172321214581765680755
278 #2.718281828:0.9999999998311266953289851340574956564911
280 #2.718281828:0.99999999983112669533
285 #1000:0:6.90775527898214
286 #100:0:4.60517018598809
288 #3.1415:0:1.14470039286086
290 #12345:0:9.42100640177928
291 #0.001:0:-6.90775527898214
292 # reset for further tests
334 # test for bug in brsft() not handling cases that return 0
379 # uses bsstr() so 5 => 5e+0 to be compatible w/ Perls output
416 1234.567:9::1234.56700
417 1234.567::-6:1234.567000
419 0.001234:6::0.00123400
420 0.001234::-8:0.00123400
443 000000_0000000_00000:0
454 -123456789:-123456789
462 -.0000000004:-0.0000000004
477 -3e111:-3000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000
478 -4e-1111:-0.0000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000004
485 123.456:2:15241.383936
488 128:-2:0.00006103515625
495 # 2 ** 0.5 == sqrt(2)
496 # 1.41..7 and not 1.4170 since fallback (bsqrt(9) is '3', not 3.0...0)
497 2:0.5:1.41421356237309504880168872420969807857
498 #2:0.2:1.148698354997035006798626946777927589444
499 #6:1.5:14.6969384566990685891837044482353483518
501 #62.5:12.5:26447206647554886213592.3959144
510 +123456789:-123456789
512 +123.456789:-123.456789
513 -123456.789:123456.789
523 +123.456789:123.456789
524 -123456.789:123456.789
526 $round_mode = "trunc"
531 +10123456789:5:10123000000
532 -10123456789:5:-10123000000
533 +10123456789.123:5:10123000000
534 -10123456789.123:5:-10123000000
535 +10123456789:9:10123456700
536 -10123456789:9:-10123456700
537 +101234500:6:101234000
538 -101234500:6:-101234000
540 +20123456789:5:20123000000
541 -20123456789:5:-20123000000
542 +20123456789.123:5:20123000000
543 -20123456789.123:5:-20123000000
544 +20123456789:9:20123456800
545 -20123456789:9:-20123456800
546 +201234500:6:201234000
547 -201234500:6:-201234000
549 +30123456789:5:30123000000
550 -30123456789:5:-30123000000
551 +30123456789.123:5:30123000000
552 -30123456789.123:5:-30123000000
553 +30123456789:9:30123456800
554 -30123456789:9:-30123456800
555 +301234500:6:301235000
556 -301234500:6:-301234000
558 +40123456789:5:40123000000
559 -40123456789:5:-40123000000
560 +40123456789.123:5:40123000000
561 -40123456789.123:5:-40123000000
562 +40123456789:9:40123456800
563 -40123456789:9:-40123456800
564 +401234500:6:401234000
565 -401234500:6:-401235000
567 +50123456789:5:50123000000
568 -50123456789:5:-50123000000
569 +50123456789.123:5:50123000000
570 -50123456789.123:5:-50123000000
571 +50123456789:9:50123456800
572 -50123456789:9:-50123456800
573 +501234500:6:501235000
574 -501234500:6:-501235000
576 +60123456789:5:60123000000
577 -60123456789:5:-60123000000
578 +60123456789:9:60123456800
579 -60123456789:9:-60123456800
580 +601234500:6:601234000
581 -601234500:6:-601234000
582 +60123456789.0123:5:60123000000
583 -60123456789.0123:5:-60123000000
585 $round_mode = "trunc"
606 -0.0061234567890:-1:0.0
614 -0.0065:-3:/-0\.006|-6e-03
615 -0.0065:-4:/-0\.006(?:5|49{5}\d+)|-6\.5e-03
616 -0.0065:-5:/-0\.006(?:5|49{5}\d+)|-6\.5e-03
622 +2.23:-1:/2.2(?:0{5}\d+)?
623 -2.23:-1:/-2.2(?:0{5}\d+)?
624 +2.27:-1:/2.(?:3|29{5}\d+)
625 -2.27:-1:/-2.(?:3|29{5}\d+)
626 +2.25:-1:/2.2(?:0{5}\d+)?
627 -2.25:-1:/-2.2(?:0{5}\d+)?
628 +2.35:-1:/2.(?:3|29{5}\d+)
629 -2.35:-1:/-2.(?:3|29{5}\d+)
631 -0.0065:-2:/-0\.01|-1e-02
632 -0.0065:-3:/-0\.006|-6e-03
633 -0.0065:-4:/-0\.006(?:5|49{5}\d+)|-6\.5e-03
634 -0.0065:-5:/-0\.006(?:5|49{5}\d+)|-6\.5e-03
640 +3.23:-1:/3.2(?:0{5}\d+)?
641 -3.23:-1:/-3.2(?:0{5}\d+)?
642 +3.27:-1:/3.(?:3|29{5}\d+)
643 -3.27:-1:/-3.(?:3|29{5}\d+)
644 +3.25:-1:/3.(?:3|29{5}\d+)
645 -3.25:-1:/-3.2(?:0{5}\d+)?
646 +3.35:-1:/3.(?:4|39{5}\d+)
647 -3.35:-1:/-3.(?:3|29{5}\d+)
649 -0.0065:-2:/-0\.01|-1e-02
650 -0.0065:-3:/-0\.006|-6e-03
651 -0.0065:-4:/-0\.006(?:5|49{5}\d+)|-6\.5e-03
652 -0.0065:-5:/-0\.006(?:5|49{5}\d+)|-6\.5e-03
658 +4.23:-1:/4.2(?:0{5}\d+)?
659 -4.23:-1:/-4.2(?:0{5}\d+)?
660 +4.27:-1:/4.(?:3|29{5}\d+)
661 -4.27:-1:/-4.(?:3|29{5}\d+)
662 +4.25:-1:/4.2(?:0{5}\d+)?
663 -4.25:-1:/-4.(?:3|29{5}\d+)
664 +4.35:-1:/4.(?:3|29{5}\d+)
665 -4.35:-1:/-4.(?:4|39{5}\d+)
667 -0.0065:-2:/-0\.01|-1e-02
668 -0.0065:-3:/-0\.007|-7e-03
669 -0.0065:-4:/-0\.006(?:5|49{5}\d+)|-6\.5e-03
670 -0.0065:-5:/-0\.006(?:5|49{5}\d+)|-6\.5e-03
676 +5.23:-1:/5.2(?:0{5}\d+)?
677 -5.23:-1:/-5.2(?:0{5}\d+)?
678 +5.27:-1:/5.(?:3|29{5}\d+)
679 -5.27:-1:/-5.(?:3|29{5}\d+)
680 +5.25:-1:/5.(?:3|29{5}\d+)
681 -5.25:-1:/-5.(?:3|29{5}\d+)
682 +5.35:-1:/5.(?:3|29{5}\d+)
683 -5.35:-1:/-5.(?:3|29{5}\d+)
685 -0.0065:-2:/-0\.01|-1e-02
686 -0.0065:-3:/-0\.007|-7e-03
687 -0.0065:-4:/-0\.006(?:5|49{5}\d+)|-6\.5e-03
688 -0.0065:-5:/-0\.006(?:5|49{5}\d+)|-6\.5e-03
694 +6.23:-1:/6.2(?:0{5}\d+)?
695 -6.23:-1:/-6.2(?:0{5}\d+)?
696 +6.27:-1:/6.(?:3|29{5}\d+)
697 -6.27:-1:/-6.(?:3|29{5}\d+)
698 +6.25:-1:/6.(?:2(?:0{5}\d+)?|29{5}\d+)
699 -6.25:-1:/-6.(?:2(?:0{5}\d+)?|29{5}\d+)
700 +6.35:-1:/6.(?:4|39{5}\d+|29{8}\d+)
701 -6.35:-1:/-6.(?:4|39{5}\d+|29{8}\d+)
703 -0.0065:-2:/-0\.01|-1e-02
704 -0.0065:-3:/-0\.006|-7e-03
705 -0.0065:-4:/-0\.006(?:5|49{5}\d+)|-6\.5e-03
706 -0.0065:-5:/-0\.006(?:5|49{5}\d+)|-6\.5e-03
713 0.01234567:-5:0.01235
714 0.01234567:-6:0.012346
715 0.01234567:-7:0.0123457
716 0.01234567:-8:0.01234567
717 0.01234567:-9:0.012345670
718 0.01234567:-12:0.012345670000
920 +99999999:+1:100000000
921 +999999999:+1:1000000000
922 +9999999999:+1:10000000000
923 +99999999999:+1:100000000000
931 +100000000:-1:99999999
932 +1000000000:-1:999999999
933 +10000000000:-1:9999999999
934 +123456789:+987654321:1111111110
935 -123456789:+987654321:864197532
936 -123456789:-987654321:-1111111110
937 +123456789:-987654321:-864197532
938 0.001234:0.0001234:0.0013574
967 +99999999:+1:99999998
968 +999999999:+1:999999998
969 +9999999999:+1:9999999998
970 +99999999999:+1:99999999998
977 +10000000:-1:10000001
978 +100000000:-1:100000001
979 +1000000000:-1:1000000001
980 +10000000000:-1:10000000001
981 +123456789:+987654321:-864197532
982 -123456789:+987654321:-1111111110
983 -123456789:-987654321:864197532
984 +123456789:-987654321:1111111110
1010 +123456789123456789:+0:0
1011 +0:+123456789123456789:0
1021 +10101:+10101:102030201
1022 +1001001:+1001001:1002003002001
1023 +100010001:+100010001:10002000300020001
1024 +10000100001:+10000100001:100002000030000200001
1025 +11111111111:+9:99999999999
1026 +22222222222:+9:199999999998
1027 +33333333333:+9:299999999997
1028 +44444444444:+9:399999999996
1029 +55555555555:+9:499999999995
1030 +66666666666:+9:599999999994
1031 +77777777777:+9:699999999993
1032 +88888888888:+9:799999999992
1033 +99999999999:+9:899999999991
1042 $div_scale = 40; $round_mode = 'even'
1068 +999999999999:+9:111111111111
1069 +999999999999:+99:10101010101
1070 +999999999999:+999:1001001001
1071 +999999999999:+9999:100010001
1072 +999999999999999:+99999:10000100001
1073 +1000000000:+9:111111111.1111111111111111111111111111111
1074 +2000000000:+9:222222222.2222222222222222222222222222222
1075 +3000000000:+9:333333333.3333333333333333333333333333333
1076 +4000000000:+9:444444444.4444444444444444444444444444444
1077 +5000000000:+9:555555555.5555555555555555555555555555556
1078 +6000000000:+9:666666666.6666666666666666666666666666667
1079 +7000000000:+9:777777777.7777777777777777777777777777778
1080 +8000000000:+9:888888888.8888888888888888888888888888889
1081 +9000000000:+9:1000000000
1082 +35500000:+113:314159.2920353982300884955752212389380531
1083 +71000000:+226:314159.2920353982300884955752212389380531
1084 +106500000:+339:314159.2920353982300884955752212389380531
1085 +1000000000:+3:333333333.3333333333333333333333333333333
1086 2:25.024996000799840031993601279744051189762:0.07992009269196593320152084692285869265447
1089 +1000000000:+9:111111111.11111111111
1090 +2000000000:+9:222222222.22222222222
1091 +3000000000:+9:333333333.33333333333
1092 +4000000000:+9:444444444.44444444444
1093 +5000000000:+9:555555555.55555555556
1094 +6000000000:+9:666666666.66666666667
1095 +7000000000:+9:777777777.77777777778
1096 +8000000000:+9:888888888.88888888889
1097 +9000000000:+9:1000000000
1102 1:504:0.001984126984126984127
1103 2:1.987654321:1.0062111801179738436
1104 123456789.123456789123456789123456789:1:123456789.12345678912
1105 # the next two cases are the "old" behaviour, but are now (>v0.01) different
1106 #+35500000:+113:314159.292035398230088
1107 #+71000000:+226:314159.292035398230088
1108 +35500000:+113:314159.29203539823009
1109 +71000000:+226:314159.29203539823009
1110 +106500000:+339:314159.29203539823009
1111 +1000000000:+3:333333333.33333333333
1113 # round to accuracy 1 after bdiv
1115 123456789.1234:1:100000000
1116 # reset scale for further tests
1123 # inf handling, see table in doc
1142 # exceptions to reminder rule
1181 999999999999999:99999:0
1195 152403346:12345:4321
1197 # now some floating point tests
1222 # -$x ** (1/2) => -$y, but not in froot()
1226 2:2:1.41421356237309504880168872420969807857
1232 123.456:2:11.11107555549866648462149404118219234119
1233 15241.38393:2:123.4559999756998444766131352122991626468
1235 12:2:3.464101615137754587054892683011744733886
1266 # see t/bigroot() for more tests
1277 2:1.41421356237309504880168872420969807857
1282 123.456:11.11107555549866648462149404118219234119
1283 15241.38393:123.4559999756998444766131352122991626468
1285 # sqrt(1.44) = 1.2, sqrt(e10) = e5 => 12e4
1287 2e10:141421.356237309504880168872420969807857
1289 # proved to be an endless loop under 7-9
1290 12:3.464101615137754587054892683011744733886
1309 # it must be exactly /^[+-]inf$/
1405 12345678901234567890:20