1 #include this file into another test for subclass testing...
3 ok ($class->config()->{lib},$CL);
8 $_ =~ s/#.*$//; # remove comments
9 $_ =~ s/\s+$//; # trailing spaces
10 next if /^$/; # skip empty lines & comments
17 $setup = $_; $setup =~ s/\$/\$${class}::/g; # round_mode, div_scale
18 #print "\$setup== $setup\n";
25 @args = split(/:/,$1,99);
29 @args = split(/:/,$_,99); $ans = pop(@args);
31 $try = "\$x = new $class \"$args[0]\";";
35 } elsif ($f eq "finf") {
36 $try .= "\$x->finf('$args[1]');";
37 } elsif ($f eq "is_inf") {
38 $try .= "\$x->is_inf('$args[1]');";
39 } elsif ($f eq "fone") {
40 $try .= "\$x->bone('$args[1]');";
41 } elsif ($f eq "fstr") {
42 $try .= "\$x->accuracy($args[1]); \$x->precision($args[2]);";
43 $try .= '$x->fstr();';
44 } elsif ($f eq "parts") {
45 # ->bstr() to see if an object is returned
46 $try .= '($a,$b) = $x->parts(); $a = $a->bstr(); $b = $b->bstr();';
48 } elsif ($f eq "exponent") {
49 # ->bstr() to see if an object is returned
50 $try .= '$x->exponent()->bstr();';
51 } elsif ($f eq "mantissa") {
52 # ->bstr() to see if an object is returned
53 $try .= '$x->mantissa()->bstr();';
54 } elsif ($f eq "numify") {
55 $try .= "\$x->numify();";
56 } elsif ($f eq "length") {
57 $try .= "\$x->length();";
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 "as_number") {
65 $try .= '$x->as_number();';
66 } elsif ($f eq "finc") {
68 } elsif ($f eq "fdec") {
70 }elsif ($f eq "fround") {
71 $try .= "$setup; \$x->fround($args[1]);";
72 } elsif ($f eq "ffround") {
73 $try .= "$setup; \$x->ffround($args[1]);";
74 } elsif ($f eq "fsqrt") {
75 $try .= "$setup; \$x->fsqrt();";
76 } elsif ($f eq "flog") {
77 $try .= "$setup; \$x->flog();";
78 } elsif ($f eq "ffac") {
79 $try .= "$setup; \$x->ffac();";
83 $try .= "\$y = new $class \"$args[1]\";";
86 } elsif ($f eq "facmp") {
87 $try .= '$x->facmp($y);';
88 } elsif ($f eq "fpow") {
90 } elsif ($f eq "fadd") {
92 } elsif ($f eq "fsub") {
94 } elsif ($f eq "fmul") {
96 } elsif ($f eq "fdiv") {
97 $try .= "$setup; \$x / \$y;";
98 } elsif ($f eq "fdiv-list") {
99 $try .= "$setup; join(',',\$x->fdiv(\$y));";
100 } elsif ($f eq "frsft") {
102 } elsif ($f eq "flsft") {
104 } elsif ($f eq "fmod") {
106 } else { warn "Unknown op '$f'"; }
108 # print "# Trying: '$try'\n";
110 if ($ans =~ m|^/(.*)$|)
119 print "# '$try' expected: /$pat/ got: '$ans1'\n" if !ok(1,0);
130 print "# Tried: '$try'\n" if !ok ($ans1, $ans);
131 if (ref($ans1) eq "$class")
133 # float numbers are normalized (for now), so mantissa shouldn't have
135 #print $ans1->_trailing_zeros(),"\n";
136 print "# Has trailing zeros after '$try'\n"
137 if !ok ($ans1->{_m}->_trailing_zeros(), 0);
140 } # end pattern or string
144 # check whether $class->new( Math::BigInt->new()) destroys it
145 # ($y == 12 in this case)
146 $x = Math::BigInt->new(1200); $y = $class->new($x);
147 ok ($y,1200); ok ($x,1200);
149 ###############################################################################
150 # Really huge, big, ultra-mega-biggy-monster exponents
151 # Technically, the exponents should not be limited (they are BigInts), but
152 # practically there are a few places were they are limited to a Perl scalar.
153 # This is sometimes for speed, sometimes because otherwise the number wouldn't
154 # fit into your memory (just think of 1e123456789012345678901234567890 + 1!)
155 # anyway. We don't test everything here, but let's make sure it just basically
158 my $monster = '1e1234567890123456789012345678901234567890';
161 ok ($class->new($monster)->bsstr(),
162 '1e+1234567890123456789012345678901234567890');
164 ok ($class->new($monster) > 0,1);
167 ok ($class->new($monster)->bsub( $monster),0);
168 ok ($class->new($monster)->bmul(2)->bsstr(),
169 '2e+1234567890123456789012345678901234567890');
171 ###############################################################################
174 $x = $class->new(2); $x->fzero(); ok_undef ($x->{_a}); ok_undef ($x->{_p});
175 $x = $class->new(2); $x->finf(); ok_undef ($x->{_a}); ok_undef ($x->{_p});
176 $x = $class->new(2); $x->fone(); ok_undef ($x->{_a}); ok_undef ($x->{_p});
177 $x = $class->new(2); $x->fnan(); ok_undef ($x->{_a}); ok_undef ($x->{_p});
179 ###############################################################################
180 # bone/binf etc as plain calls (Lite failed them)
182 ok ($class->fzero(),0);
183 ok ($class->fone(),1);
184 ok ($class->fone('+'),1);
185 ok ($class->fone('-'),-1);
186 ok ($class->fnan(),'NaN');
187 ok ($class->finf(),'inf');
188 ok ($class->finf('+'),'inf');
189 ok ($class->finf('-'),'-inf');
190 ok ($class->finf('-inf'),'-inf');
192 ###############################################################################
193 # fsqrt() with set global A/P or A/P enabled on $x, also a test whether fsqrt()
194 # correctly modifies $x
196 $class->accuracy(undef); $class->precision(undef); # reset
198 $x = $class->new(12); $class->precision(-2); $x->fsqrt(); ok ($x,'3.46');
200 $class->precision(undef);
201 $x = $class->new(12); $class->precision(0); $x->fsqrt(); ok ($x,'3');
203 $class->precision(-3); $x = $class->new(12); $x->fsqrt(); ok ($x,'3.464');
206 ${${class}.'::accuracy'} = 4; $x = $class->new(12); $x->fsqrt(3); ok ($x,'NaN');
207 # supplied arg overrides set global
208 $class->precision(undef); $x = $class->new(12); $x->fsqrt(3); ok ($x,'3.46');
210 $class->accuracy(undef); $class->precision(undef); # reset for further tests
212 ###############################################################################
213 # can we call objectify (broken until v1.52)
215 $try = '@args' . " = $class" . "::objectify(2,$class,4,5);".'join(" ",@args);';
217 ok ($ans,"$class 4 5");
219 ###############################################################################
220 # test whether an opp calls objectify properly or not (or at least does what
221 # it should do given non-objects, w/ or w/o objectify())
223 ok ($class->new(123)->badd(123),246);
224 ok ($class->badd(123,321),444);
225 ok ($class->badd(123,$class->new(321)),444);
227 ok ($class->new(123)->bsub(122),1);
228 ok ($class->bsub(321,123),198);
229 ok ($class->bsub(321,$class->new(123)),198);
231 ok ($class->new(123)->bmul(123),15129);
232 ok ($class->bmul(123,123),15129);
233 ok ($class->bmul(123,$class->new(123)),15129);
235 ok ($class->new(15129)->bdiv(123),123);
236 ok ($class->bdiv(15129,123),123);
237 ok ($class->bdiv(15129,$class->new(123)),123);
239 ok ($class->new(15131)->bmod(123),2);
240 ok ($class->bmod(15131,123),2);
241 ok ($class->bmod(15131,$class->new(123)),2);
243 ok ($class->new(2)->bpow(16),65536);
244 ok ($class->bpow(2,16),65536);
245 ok ($class->bpow(2,$class->new(16)),65536);
247 ok ($class->new(2**15)->brsft(1),2**14);
248 ok ($class->brsft(2**15,1),2**14);
249 ok ($class->brsft(2**15,$class->new(1)),2**14);
251 ok ($class->new(2**13)->blsft(1),2**14);
252 ok ($class->blsft(2**13,1),2**14);
253 ok ($class->blsft(2**13,$class->new(1)),2**14);
257 ###############################################################################
258 # Perl 5.005 does not like ok ($x,undef)
264 ok (1,1) and return if !defined $x;
275 # this is too slow for the testsuite
276 #2:0.6931471805599453094172321214581765680755
277 #2.718281828:0.9999999998311266953289851340574956564911
279 #2.718281828:0.99999999983112669533
280 # too slow, too (or hangs?)
284 #1000:0:6.90775527898214
285 #100:0:4.60517018598809
287 #3.1415:0:1.14470039286086
288 #12345:0:9.42100640177928
289 #0.001:0:-6.90775527898214
290 # reset for further tests
366 1234.567:9::1234.56700
367 1234.567::-6:1234.567000
369 0.001234:6::0.00123400
370 0.001234::-8:0.00123400
389 000000_0000000_00000:0
400 -123456789:-123456789
408 -.0000000004:-0.0000000004
423 -3e111:-3000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000
424 -4e-1111:-0.0000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000004
431 123.456:2:15241.383936
434 128:-2:0.00006103515625
441 # 2 ** 0.5 == sqrt(2)
442 # 1.41..7 and not 1.4170 since fallback (bsqrt(9) is '3', not 3.0...0)
443 2:0.5:1.41421356237309504880168872420969807857
444 #2:0.2:1.148698354997035006798626946777927589444
445 #6:1.5:14.6969384566990685891837044482353483518
447 #62.5:12.5:26447206647554886213592.3959144
456 +123456789:-123456789
458 +123.456789:-123.456789
459 -123456.789:123456.789
469 +123.456789:123.456789
470 -123456.789:123456.789
472 $round_mode = "trunc"
477 +10123456789:5:10123000000
478 -10123456789:5:-10123000000
479 +10123456789.123:5:10123000000
480 -10123456789.123:5:-10123000000
481 +10123456789:9:10123456700
482 -10123456789:9:-10123456700
483 +101234500:6:101234000
484 -101234500:6:-101234000
486 +20123456789:5:20123000000
487 -20123456789:5:-20123000000
488 +20123456789.123:5:20123000000
489 -20123456789.123:5:-20123000000
490 +20123456789:9:20123456800
491 -20123456789:9:-20123456800
492 +201234500:6:201234000
493 -201234500:6:-201234000
495 +30123456789:5:30123000000
496 -30123456789:5:-30123000000
497 +30123456789.123:5:30123000000
498 -30123456789.123:5:-30123000000
499 +30123456789:9:30123456800
500 -30123456789:9:-30123456800
501 +301234500:6:301235000
502 -301234500:6:-301234000
504 +40123456789:5:40123000000
505 -40123456789:5:-40123000000
506 +40123456789.123:5:40123000000
507 -40123456789.123:5:-40123000000
508 +40123456789:9:40123456800
509 -40123456789:9:-40123456800
510 +401234500:6:401234000
511 -401234500:6:-401235000
513 +50123456789:5:50123000000
514 -50123456789:5:-50123000000
515 +50123456789.123:5:50123000000
516 -50123456789.123:5:-50123000000
517 +50123456789:9:50123456800
518 -50123456789:9:-50123456800
519 +501234500:6:501235000
520 -501234500:6:-501235000
522 +60123456789:5:60123000000
523 -60123456789:5:-60123000000
524 +60123456789:9:60123456800
525 -60123456789:9:-60123456800
526 +601234500:6:601234000
527 -601234500:6:-601234000
528 +60123456789.0123:5:60123000000
529 -60123456789.0123:5:-60123000000
531 $round_mode = "trunc"
552 -0.0061234567890:-1:0.0
560 -0.0065:-3:/-0\.006|-6e-03
561 -0.0065:-4:/-0\.006(?:5|49{5}\d+)|-6\.5e-03
562 -0.0065:-5:/-0\.006(?:5|49{5}\d+)|-6\.5e-03
568 +2.23:-1:/2.2(?:0{5}\d+)?
569 -2.23:-1:/-2.2(?:0{5}\d+)?
570 +2.27:-1:/2.(?:3|29{5}\d+)
571 -2.27:-1:/-2.(?:3|29{5}\d+)
572 +2.25:-1:/2.2(?:0{5}\d+)?
573 -2.25:-1:/-2.2(?:0{5}\d+)?
574 +2.35:-1:/2.(?:3|29{5}\d+)
575 -2.35:-1:/-2.(?:3|29{5}\d+)
577 -0.0065:-2:/-0\.01|-1e-02
578 -0.0065:-3:/-0\.006|-6e-03
579 -0.0065:-4:/-0\.006(?:5|49{5}\d+)|-6\.5e-03
580 -0.0065:-5:/-0\.006(?:5|49{5}\d+)|-6\.5e-03
586 +3.23:-1:/3.2(?:0{5}\d+)?
587 -3.23:-1:/-3.2(?:0{5}\d+)?
588 +3.27:-1:/3.(?:3|29{5}\d+)
589 -3.27:-1:/-3.(?:3|29{5}\d+)
590 +3.25:-1:/3.(?:3|29{5}\d+)
591 -3.25:-1:/-3.2(?:0{5}\d+)?
592 +3.35:-1:/3.(?:4|39{5}\d+)
593 -3.35:-1:/-3.(?:3|29{5}\d+)
595 -0.0065:-2:/-0\.01|-1e-02
596 -0.0065:-3:/-0\.006|-6e-03
597 -0.0065:-4:/-0\.006(?:5|49{5}\d+)|-6\.5e-03
598 -0.0065:-5:/-0\.006(?:5|49{5}\d+)|-6\.5e-03
604 +4.23:-1:/4.2(?:0{5}\d+)?
605 -4.23:-1:/-4.2(?:0{5}\d+)?
606 +4.27:-1:/4.(?:3|29{5}\d+)
607 -4.27:-1:/-4.(?:3|29{5}\d+)
608 +4.25:-1:/4.2(?:0{5}\d+)?
609 -4.25:-1:/-4.(?:3|29{5}\d+)
610 +4.35:-1:/4.(?:3|29{5}\d+)
611 -4.35:-1:/-4.(?:4|39{5}\d+)
613 -0.0065:-2:/-0\.01|-1e-02
614 -0.0065:-3:/-0\.007|-7e-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 +5.23:-1:/5.2(?:0{5}\d+)?
623 -5.23:-1:/-5.2(?:0{5}\d+)?
624 +5.27:-1:/5.(?:3|29{5}\d+)
625 -5.27:-1:/-5.(?:3|29{5}\d+)
626 +5.25:-1:/5.(?:3|29{5}\d+)
627 -5.25:-1:/-5.(?:3|29{5}\d+)
628 +5.35:-1:/5.(?:3|29{5}\d+)
629 -5.35:-1:/-5.(?:3|29{5}\d+)
631 -0.0065:-2:/-0\.01|-1e-02
632 -0.0065:-3:/-0\.007|-7e-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 +6.23:-1:/6.2(?:0{5}\d+)?
641 -6.23:-1:/-6.2(?:0{5}\d+)?
642 +6.27:-1:/6.(?:3|29{5}\d+)
643 -6.27:-1:/-6.(?:3|29{5}\d+)
644 +6.25:-1:/6.(?:2(?:0{5}\d+)?|29{5}\d+)
645 -6.25:-1:/-6.(?:2(?:0{5}\d+)?|29{5}\d+)
646 +6.35:-1:/6.(?:4|39{5}\d+|29{8}\d+)
647 -6.35:-1:/-6.(?:4|39{5}\d+|29{8}\d+)
649 -0.0065:-2:/-0\.01|-1e-02
650 -0.0065:-3:/-0\.006|-7e-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
659 0.01234567:-5:0.01235
660 0.01234567:-6:0.012346
661 0.01234567:-7:0.0123457
662 0.01234567:-8:0.01234567
663 0.01234567:-9:0.012345670
664 0.01234567:-12:0.012345670000
866 +99999999:+1:100000000
867 +999999999:+1:1000000000
868 +9999999999:+1:10000000000
869 +99999999999:+1:100000000000
877 +100000000:-1:99999999
878 +1000000000:-1:999999999
879 +10000000000:-1:9999999999
880 +123456789:+987654321:1111111110
881 -123456789:+987654321:864197532
882 -123456789:-987654321:-1111111110
883 +123456789:-987654321:-864197532
884 0.001234:0.0001234:0.0013574
913 +99999999:+1:99999998
914 +999999999:+1:999999998
915 +9999999999:+1:9999999998
916 +99999999999:+1:99999999998
923 +10000000:-1:10000001
924 +100000000:-1:100000001
925 +1000000000:-1:1000000001
926 +10000000000:-1:10000000001
927 +123456789:+987654321:-864197532
928 -123456789:+987654321:-1111111110
929 -123456789:-987654321:864197532
930 +123456789:-987654321:1111111110
956 +123456789123456789:+0:0
957 +0:+123456789123456789:0
967 +10101:+10101:102030201
968 +1001001:+1001001:1002003002001
969 +100010001:+100010001:10002000300020001
970 +10000100001:+10000100001:100002000030000200001
971 +11111111111:+9:99999999999
972 +22222222222:+9:199999999998
973 +33333333333:+9:299999999997
974 +44444444444:+9:399999999996
975 +55555555555:+9:499999999995
976 +66666666666:+9:599999999994
977 +77777777777:+9:699999999993
978 +88888888888:+9:799999999992
979 +99999999999:+9:899999999991
988 $div_scale = 40; $round_mode = 'even'
1014 +999999999999:+9:111111111111
1015 +999999999999:+99:10101010101
1016 +999999999999:+999:1001001001
1017 +999999999999:+9999:100010001
1018 +999999999999999:+99999:10000100001
1019 +1000000000:+9:111111111.1111111111111111111111111111111
1020 +2000000000:+9:222222222.2222222222222222222222222222222
1021 +3000000000:+9:333333333.3333333333333333333333333333333
1022 +4000000000:+9:444444444.4444444444444444444444444444444
1023 +5000000000:+9:555555555.5555555555555555555555555555556
1024 +6000000000:+9:666666666.6666666666666666666666666666667
1025 +7000000000:+9:777777777.7777777777777777777777777777778
1026 +8000000000:+9:888888888.8888888888888888888888888888889
1027 +9000000000:+9:1000000000
1028 +35500000:+113:314159.2920353982300884955752212389380531
1029 +71000000:+226:314159.2920353982300884955752212389380531
1030 +106500000:+339:314159.2920353982300884955752212389380531
1031 +1000000000:+3:333333333.3333333333333333333333333333333
1032 2:25.024996000799840031993601279744051189762:0.07992009269196593320152084692285869265447
1035 +1000000000:+9:111111111.11111111111
1036 +2000000000:+9:222222222.22222222222
1037 +3000000000:+9:333333333.33333333333
1038 +4000000000:+9:444444444.44444444444
1039 +5000000000:+9:555555555.55555555556
1040 +6000000000:+9:666666666.66666666667
1041 +7000000000:+9:777777777.77777777778
1042 +8000000000:+9:888888888.88888888889
1043 +9000000000:+9:1000000000
1048 1:504:0.001984126984126984127
1049 2:1.987654321:1.0062111801179738436
1050 123456789.123456789123456789123456789:1:123456789.12345678912
1051 # the next two cases are the "old" behaviour, but are now (>v0.01) different
1052 #+35500000:+113:314159.292035398230088
1053 #+71000000:+226:314159.292035398230088
1054 +35500000:+113:314159.29203539823009
1055 +71000000:+226:314159.29203539823009
1056 +106500000:+339:314159.29203539823009
1057 +1000000000:+3:333333333.33333333333
1059 # round to accuracy 1 after bdiv
1061 123456789.1234:1:100000000
1062 # reset scale for further tests
1069 # inf handling, see table in doc
1088 # exceptions to reminder rule
1127 999999999999999:99999:0
1141 152403346:12345:4321
1143 # now some floating point tests
1171 2:1.41421356237309504880168872420969807857
1176 123.456:11.11107555549866648462149404118219234119
1177 15241.38393:123.4559999756998444766131352122991626468
1179 # sqrt(1.44) = 1.2, sqrt(e10) = e5 => 12e4
1181 2e10:141421.356237309504880168872420969807857
1183 # proved to be an endless loop under 7-9
1184 12:3.464101615137754587054892683011744733886
1201 # it must be exactly /^[+-]inf$/
1297 12345678901234567890:20