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