1N/A#!/usr/bin/perl -w
1N/A
1N/Ause strict;
1N/Ause Test;
1N/A
1N/ABEGIN
1N/A {
1N/A $| = 1;
1N/A chdir 't' if -d 't';
1N/A unshift @INC, '../lib'; # for running manually
1N/A plan tests => 308;
1N/A }
1N/A
1N/Ause Math::BigInt::Calc;
1N/A
1N/Amy ($BASE_LEN, $AND_BITS, $XOR_BITS, $OR_BITS, $BASE_LEN_SMALL, $MAX_VAL) =
1N/A Math::BigInt::Calc->_base_len();
1N/A
1N/Aprint "# BASE_LEN = $BASE_LEN\n";
1N/Aprint "# MAX_VAL = $MAX_VAL\n";
1N/Aprint "# AND_BITS = $AND_BITS\n";
1N/Aprint "# XOR_BITS = $XOR_BITS\n";
1N/Aprint "# IOR_BITS = $OR_BITS\n";
1N/A
1N/A# testing of Math::BigInt::Calc
1N/A
1N/Amy $C = 'Math::BigInt::Calc'; # pass classname to sub's
1N/A
1N/A# _new and _str
1N/Amy $x = $C->_new("123"); my $y = $C->_new("321");
1N/Aok (ref($x),'ARRAY'); ok ($C->_str($x),123); ok ($C->_str($y),321);
1N/A
1N/A###############################################################################
1N/A# _add, _sub, _mul, _div
1N/Aok ($C->_str($C->_add($x,$y)),444);
1N/Aok ($C->_str($C->_sub($x,$y)),123);
1N/Aok ($C->_str($C->_mul($x,$y)),39483);
1N/Aok ($C->_str($C->_div($x,$y)),123);
1N/A
1N/A###############################################################################
1N/A# check that mul/div doesn't change $y
1N/A# and returns the same reference, not something new
1N/Aok ($C->_str($C->_mul($x,$y)),39483);
1N/Aok ($C->_str($x),39483); ok ($C->_str($y),321);
1N/A
1N/Aok ($C->_str($C->_div($x,$y)),123);
1N/Aok ($C->_str($x),123); ok ($C->_str($y),321);
1N/A
1N/A$x = $C->_new("39483");
1N/Amy ($x1,$r1) = $C->_div($x,$y);
1N/Aok ("$x1","$x");
1N/A$C->_inc($x1);
1N/Aok ("$x1","$x");
1N/Aok ($C->_str($r1),'0');
1N/A
1N/A$x = $C->_new("39483"); # reset
1N/A
1N/A###############################################################################
1N/Amy $z = $C->_new("2");
1N/Aok ($C->_str($C->_add($x,$z)),39485);
1N/Amy ($re,$rr) = $C->_div($x,$y);
1N/A
1N/Aok ($C->_str($re),123); ok ($C->_str($rr),2);
1N/A
1N/A# is_zero, _is_one, _one, _zero
1N/Aok ($C->_is_zero($x)||0,0);
1N/Aok ($C->_is_one($x)||0,0);
1N/A
1N/Aok ($C->_str($C->_zero()),"0");
1N/Aok ($C->_str($C->_one()),"1");
1N/A
1N/A# _two() and _ten()
1N/Aok ($C->_str($C->_two()),"2");
1N/Aok ($C->_str($C->_ten()),"10");
1N/Aok ($C->_is_ten($C->_two()),0);
1N/Aok ($C->_is_two($C->_two()),1);
1N/Aok ($C->_is_ten($C->_ten()),1);
1N/Aok ($C->_is_two($C->_ten()),0);
1N/A
1N/Aok ($C->_is_one($C->_one()),1);
1N/Aok ($C->_is_one($C->_two()),0);
1N/Aok ($C->_is_one($C->_ten()),0);
1N/A
1N/Aok ($C->_is_one($C->_zero()) || 0,0);
1N/A
1N/Aok ($C->_is_zero($C->_zero()),1);
1N/A
1N/Aok ($C->_is_zero($C->_one()) || 0,0);
1N/A
1N/A# is_odd, is_even
1N/Aok ($C->_is_odd($C->_one()),1); ok ($C->_is_odd($C->_zero())||0,0);
1N/Aok ($C->_is_even($C->_one()) || 0,0); ok ($C->_is_even($C->_zero()),1);
1N/A
1N/A# _len
1N/A$x = $C->_new("1"); ok ($C->_len($x),1);
1N/A$x = $C->_new("12"); ok ($C->_len($x),2);
1N/A$x = $C->_new("123"); ok ($C->_len($x),3);
1N/A$x = $C->_new("1234"); ok ($C->_len($x),4);
1N/A$x = $C->_new("12345"); ok ($C->_len($x),5);
1N/A$x = $C->_new("123456"); ok ($C->_len($x),6);
1N/A$x = $C->_new("1234567"); ok ($C->_len($x),7);
1N/A$x = $C->_new("12345678"); ok ($C->_len($x),8);
1N/A$x = $C->_new("123456789"); ok ($C->_len($x),9);
1N/A
1N/A$x = $C->_new("8"); ok ($C->_len($x),1);
1N/A$x = $C->_new("21"); ok ($C->_len($x),2);
1N/A$x = $C->_new("321"); ok ($C->_len($x),3);
1N/A$x = $C->_new("4321"); ok ($C->_len($x),4);
1N/A$x = $C->_new("54321"); ok ($C->_len($x),5);
1N/A$x = $C->_new("654321"); ok ($C->_len($x),6);
1N/A$x = $C->_new("7654321"); ok ($C->_len($x),7);
1N/A$x = $C->_new("87654321"); ok ($C->_len($x),8);
1N/A$x = $C->_new("987654321"); ok ($C->_len($x),9);
1N/A
1N/Afor (my $i = 1; $i < 9; $i++)
1N/A {
1N/A my $a = "$i" . '0' x ($i-1);
1N/A $x = $C->_new($a);
1N/A print "# Tried len '$a'\n" unless ok ($C->_len($x),$i);
1N/A }
1N/A
1N/A# _digit
1N/A$x = $C->_new("123456789");
1N/Aok ($C->_digit($x,0),9);
1N/Aok ($C->_digit($x,1),8);
1N/Aok ($C->_digit($x,2),7);
1N/Aok ($C->_digit($x,-1),1);
1N/Aok ($C->_digit($x,-2),2);
1N/Aok ($C->_digit($x,-3),3);
1N/A
1N/A# _copy
1N/Aforeach (qw/ 1 12 123 1234 12345 123456 1234567 12345678 123456789/)
1N/A {
1N/A $x = $C->_new("$_");
1N/A ok ($C->_str($C->_copy($x)),"$_");
1N/A ok ($C->_str($x),"$_"); # did _copy destroy original x?
1N/A }
1N/A
1N/A# _zeros
1N/A$x = $C->_new("1256000000"); ok ($C->_zeros($x),6);
1N/A$x = $C->_new("152"); ok ($C->_zeros($x),0);
1N/A$x = $C->_new("123000"); ok ($C->_zeros($x),3);
1N/A$x = $C->_new("0"); ok ($C->_zeros($x),0);
1N/A
1N/A# _lsft, _rsft
1N/A$x = $C->_new("10"); $y = $C->_new("3");
1N/Aok ($C->_str($C->_lsft($x,$y,10)),10000);
1N/A$x = $C->_new("20"); $y = $C->_new("3");
1N/Aok ($C->_str($C->_lsft($x,$y,10)),20000);
1N/A
1N/A$x = $C->_new("128"); $y = $C->_new("4");
1N/Aok ($C->_str($C->_lsft($x,$y,2)), 128 << 4);
1N/A
1N/A$x = $C->_new("1000"); $y = $C->_new("3");
1N/Aok ($C->_str($C->_rsft($x,$y,10)),1);
1N/A$x = $C->_new("20000"); $y = $C->_new("3");
1N/Aok ($C->_str($C->_rsft($x,$y,10)),20);
1N/A$x = $C->_new("256"); $y = $C->_new("4");
1N/Aok ($C->_str($C->_rsft($x,$y,2)),256 >> 4);
1N/A
1N/A$x = $C->_new("6411906467305339182857313397200584952398");
1N/A$y = $C->_new("45");
1N/Aok ($C->_str($C->_rsft($x,$y,10)),0);
1N/A
1N/A# _acmp
1N/A$x = $C->_new("123456789");
1N/A$y = $C->_new("987654321");
1N/Aok ($C->_acmp($x,$y),-1);
1N/Aok ($C->_acmp($y,$x),1);
1N/Aok ($C->_acmp($x,$x),0);
1N/Aok ($C->_acmp($y,$y),0);
1N/A$x = $C->_new("12");
1N/A$y = $C->_new("12");
1N/Aok ($C->_acmp($x,$y),0);
1N/A$x = $C->_new("21");
1N/Aok ($C->_acmp($x,$y),1);
1N/Aok ($C->_acmp($y,$x),-1);
1N/A$x = $C->_new("123456789");
1N/A$y = $C->_new("1987654321");
1N/Aok ($C->_acmp($x,$y),-1);
1N/Aok ($C->_acmp($y,$x),+1);
1N/A
1N/A$x = $C->_new("1234567890123456789");
1N/A$y = $C->_new("987654321012345678");
1N/Aok ($C->_acmp($x,$y),1);
1N/Aok ($C->_acmp($y,$x),-1);
1N/Aok ($C->_acmp($x,$x),0);
1N/Aok ($C->_acmp($y,$y),0);
1N/A
1N/A$x = $C->_new("1234");
1N/A$y = $C->_new("987654321012345678");
1N/Aok ($C->_acmp($x,$y),-1);
1N/Aok ($C->_acmp($y,$x),1);
1N/Aok ($C->_acmp($x,$x),0);
1N/Aok ($C->_acmp($y,$y),0);
1N/A
1N/A# _modinv
1N/A$x = $C->_new("8");
1N/A$y = $C->_new("5033");
1N/Amy ($xmod,$sign) = $C->_modinv($x,$y);
1N/Aok ($C->_str($xmod),'629'); # -629 % 5033 == 4404
1N/Aok ($sign, '-');
1N/A
1N/A# _div
1N/A$x = $C->_new("3333"); $y = $C->_new("1111");
1N/Aok ($C->_str(scalar $C->_div($x,$y)),3);
1N/A$x = $C->_new("33333"); $y = $C->_new("1111"); ($x,$y) = $C->_div($x,$y);
1N/Aok ($C->_str($x),30); ok ($C->_str($y),3);
1N/A$x = $C->_new("123"); $y = $C->_new("1111");
1N/A($x,$y) = $C->_div($x,$y); ok ($C->_str($x),0); ok ($C->_str($y),123);
1N/A
1N/A# _num
1N/Aforeach (qw/1 12 123 1234 12345 1234567 12345678 123456789 1234567890/)
1N/A {
1N/A $x = $C->_new("$_");
1N/A ok (ref($x)||'','ARRAY'); ok ($C->_str($x),"$_");
1N/A $x = $C->_num($x); ok (ref($x)||'',''); ok ($x,$_);
1N/A }
1N/A
1N/A# _sqrt
1N/A$x = $C->_new("144"); ok ($C->_str($C->_sqrt($x)),'12');
1N/A$x = $C->_new("144000000000000"); ok ($C->_str($C->_sqrt($x)),'12000000');
1N/A
1N/A# _root
1N/A$x = $C->_new("81"); my $n = $C->_new("3"); # 4*4*4 = 64, 5*5*5 = 125
1N/Aok ($C->_str($C->_root($x,$n)),'4'); # 4.xx => 4.0
1N/A$x = $C->_new("81"); $n = $C->_new("4"); # 3*3*3*3 == 81
1N/Aok ($C->_str($C->_root($x,$n)),'3');
1N/A
1N/A# _pow (and _root)
1N/A$x = $C->_new("0"); $n = $C->_new("3"); # 0 ** y => 0
1N/Aok ($C->_str($C->_pow($x,$n)), 0);
1N/A$x = $C->_new("3"); $n = $C->_new("0"); # x ** 0 => 1
1N/Aok ($C->_str($C->_pow($x,$n)), 1);
1N/A$x = $C->_new("1"); $n = $C->_new("3"); # 1 ** y => 1
1N/Aok ($C->_str($C->_pow($x,$n)), 1);
1N/A$x = $C->_new("5"); $n = $C->_new("1"); # x ** 1 => x
1N/Aok ($C->_str($C->_pow($x,$n)), 5);
1N/A
1N/A$x = $C->_new("81"); $n = $C->_new("3"); # 81 ** 3 == 531441
1N/Aok ($C->_str($C->_pow($x,$n)),81 ** 3);
1N/A
1N/Aok ($C->_str($C->_root($x,$n)),81);
1N/A
1N/A$x = $C->_new("81");
1N/Aok ($C->_str($C->_pow($x,$n)),81 ** 3);
1N/Aok ($C->_str($C->_pow($x,$n)),'150094635296999121'); # 531441 ** 3 ==
1N/A
1N/Aok ($C->_str($C->_root($x,$n)),'531441');
1N/Aok ($C->_str($C->_root($x,$n)),'81');
1N/A
1N/A$x = $C->_new("81"); $n = $C->_new("14");
1N/Aok ($C->_str($C->_pow($x,$n)),'523347633027360537213511521');
1N/Aok ($C->_str($C->_root($x,$n)),'81');
1N/A
1N/A$x = $C->_new("523347633027360537213511520");
1N/Aok ($C->_str($C->_root($x,$n)),'80');
1N/A
1N/A$x = $C->_new("523347633027360537213511522");
1N/Aok ($C->_str($C->_root($x,$n)),'81');
1N/A
1N/Amy $res = [ qw/ 9 31 99 316 999 3162 9999/ ];
1N/A
1N/A# 99 ** 2 = 9801, 999 ** 2 = 998001 etc
1N/Afor my $i (2 .. 9)
1N/A {
1N/A $x = '9' x $i; $x = $C->_new($x);
1N/A $n = $C->_new("2");
1N/A my $rc = '9' x ($i-1). '8' . '0' x ($i-1) . '1';
1N/A print "# _pow( ", '9' x $i, ", 2) \n" unless
1N/A ok ($C->_str($C->_pow($x,$n)),$rc);
1N/A
1N/A if ($i <= 7)
1N/A {
1N/A $x = '9' x $i; $x = $C->_new($x);
1N/A $n = '9' x $i; $n = $C->_new($n);
1N/A print "# _root( ", '9' x $i, ", ", 9 x $i, ") \n" unless
1N/A ok ($C->_str($C->_root($x,$n)),'1');
1N/A
1N/A $x = '9' x $i; $x = $C->_new($x);
1N/A $n = $C->_new("2");
1N/A print "# _root( ", '9' x $i, ", ", 9 x $i, ") \n" unless
1N/A ok ($C->_str($C->_root($x,$n)), $res->[$i-2]);
1N/A }
1N/A }
1N/A
1N/A##############################################################################
1N/A# _fac
1N/A$x = $C->_new("0"); ok ($C->_str($C->_fac($x)),'1');
1N/A$x = $C->_new("1"); ok ($C->_str($C->_fac($x)),'1');
1N/A$x = $C->_new("2"); ok ($C->_str($C->_fac($x)),'2');
1N/A$x = $C->_new("3"); ok ($C->_str($C->_fac($x)),'6');
1N/A$x = $C->_new("4"); ok ($C->_str($C->_fac($x)),'24');
1N/A$x = $C->_new("5"); ok ($C->_str($C->_fac($x)),'120');
1N/A$x = $C->_new("10"); ok ($C->_str($C->_fac($x)),'3628800');
1N/A$x = $C->_new("11"); ok ($C->_str($C->_fac($x)),'39916800');
1N/A$x = $C->_new("12"); ok ($C->_str($C->_fac($x)),'479001600');
1N/A$x = $C->_new("13"); ok ($C->_str($C->_fac($x)),'6227020800');
1N/A
1N/A# test that _fac modifes $x in place for small arguments
1N/A$x = $C->_new("3"); $C->_fac($x); ok ($C->_str($x),'6');
1N/A$x = $C->_new("13"); $C->_fac($x); ok ($C->_str($x),'6227020800');
1N/A
1N/A##############################################################################
1N/A# _inc and _dec
1N/Aforeach (qw/1 11 121 1231 12341 1234561 12345671 123456781 1234567891/)
1N/A {
1N/A $x = $C->_new("$_"); $C->_inc($x);
1N/A print "# \$x = ",$C->_str($x),"\n"
1N/A unless ok ($C->_str($x),substr($_,0,length($_)-1) . '2');
1N/A $C->_dec($x); ok ($C->_str($x),$_);
1N/A }
1N/Aforeach (qw/19 119 1219 12319 1234519 12345619 123456719 1234567819/)
1N/A {
1N/A $x = $C->_new("$_"); $C->_inc($x);
1N/A print "# \$x = ",$C->_str($x),"\n"
1N/A unless ok ($C->_str($x),substr($_,0,length($_)-2) . '20');
1N/A $C->_dec($x); ok ($C->_str($x),$_);
1N/A }
1N/Aforeach (qw/999 9999 99999 9999999 99999999 999999999 9999999999 99999999999/)
1N/A {
1N/A $x = $C->_new("$_"); $C->_inc($x);
1N/A print "# \$x = ",$C->_str($x),"\n"
1N/A unless ok ($C->_str($x), '1' . '0' x (length($_)));
1N/A $C->_dec($x); ok ($C->_str($x),$_);
1N/A }
1N/A
1N/A$x = $C->_new("1000"); $C->_inc($x); ok ($C->_str($x),'1001');
1N/A$C->_dec($x); ok ($C->_str($x),'1000');
1N/A
1N/Amy $BL;
1N/A{
1N/A no strict 'refs';
1N/A $BL = &{"$C"."::_base_len"}();
1N/A}
1N/A
1N/A$x = '1' . '0' x $BL;
1N/A$z = '1' . '0' x ($BL-1); $z .= '1';
1N/A$x = $C->_new($x); $C->_inc($x); ok ($C->_str($x),$z);
1N/A
1N/A$x = '1' . '0' x $BL; $z = '9' x $BL;
1N/A$x = $C->_new($x); $C->_dec($x); ok ($C->_str($x),$z);
1N/A
1N/A# should not happen:
1N/A# $x = $C->_new("-2"); $y = $C->_new("4"); ok ($C->_acmp($x,$y),-1);
1N/A
1N/A###############################################################################
1N/A# _mod
1N/A$x = $C->_new("1000"); $y = $C->_new("3");
1N/Aok ($C->_str(scalar $C->_mod($x,$y)),1);
1N/A$x = $C->_new("1000"); $y = $C->_new("2");
1N/Aok ($C->_str(scalar $C->_mod($x,$y)),0);
1N/A
1N/A# _and, _or, _xor
1N/A$x = $C->_new("5"); $y = $C->_new("2");
1N/Aok ($C->_str(scalar $C->_xor($x,$y)),7);
1N/A$x = $C->_new("5"); $y = $C->_new("2");
1N/Aok ($C->_str(scalar $C->_or($x,$y)),7);
1N/A$x = $C->_new("5"); $y = $C->_new("3");
1N/Aok ($C->_str(scalar $C->_and($x,$y)),1);
1N/A
1N/A# _from_hex, _from_bin
1N/Aok ($C->_str( $C->_from_hex("0xFf")),255);
1N/Aok ($C->_str( $C->_from_bin("0b10101011")),160+11);
1N/A
1N/A# _as_hex, _as_bin
1N/Aok ($C->_str( $C->_from_hex( $C->_as_hex( $C->_new("128")))), 128);
1N/Aok ($C->_str( $C->_from_bin( $C->_as_bin( $C->_new("128")))), 128);
1N/Aok ($C->_str( $C->_from_hex( $C->_as_hex( $C->_new("0")))), 0);
1N/Aok ($C->_str( $C->_from_bin( $C->_as_bin( $C->_new("0")))), 0);
1N/Aok ($C->_as_hex( $C->_new("0")), '0x0');
1N/Aok ($C->_as_bin( $C->_new("0")), '0b0');
1N/Aok ($C->_as_hex( $C->_new("12")), '0xc');
1N/Aok ($C->_as_bin( $C->_new("12")), '0b1100');
1N/A
1N/A# _check
1N/A$x = $C->_new("123456789");
1N/Aok ($C->_check($x),0);
1N/Aok ($C->_check(123),'123 is not a reference');
1N/A
1N/A###############################################################################
1N/A# __strip_zeros
1N/A
1N/A{
1N/A no strict 'refs';
1N/A # correct empty arrays
1N/A $x = &{$C."::__strip_zeros"}([]); ok (@$x,1); ok ($x->[0],0);
1N/A # don't strip single elements
1N/A $x = &{$C."::__strip_zeros"}([0]); ok (@$x,1); ok ($x->[0],0);
1N/A $x = &{$C."::__strip_zeros"}([1]); ok (@$x,1); ok ($x->[0],1);
1N/A # don't strip non-zero elements
1N/A $x = &{$C."::__strip_zeros"}([0,1]);
1N/A ok (@$x,2); ok ($x->[0],0); ok ($x->[1],1);
1N/A $x = &{$C."::__strip_zeros"}([0,1,2]);
1N/A ok (@$x,3); ok ($x->[0],0); ok ($x->[1],1); ok ($x->[2],2);
1N/A
1N/A # but strip leading zeros
1N/A $x = &{$C."::__strip_zeros"}([0,1,2,0]);
1N/A ok (@$x,3); ok ($x->[0],0); ok ($x->[1],1); ok ($x->[2],2);
1N/A
1N/A $x = &{$C."::__strip_zeros"}([0,1,2,0,0]);
1N/A ok (@$x,3); ok ($x->[0],0); ok ($x->[1],1); ok ($x->[2],2);
1N/A
1N/A $x = &{$C."::__strip_zeros"}([0,1,2,0,0,0]);
1N/A ok (@$x,3); ok ($x->[0],0); ok ($x->[1],1); ok ($x->[2],2);
1N/A
1N/A # collapse multiple zeros
1N/A $x = &{$C."::__strip_zeros"}([0,0,0,0]);
1N/A ok (@$x,1); ok ($x->[0],0);
1N/A}
1N/A
1N/A# done
1N/A
1N/A1;
1N/A