1#!/usr/bin/perl -w 2 3use strict; 4use Test; 5 6BEGIN 7 { 8 $| = 1; 9 chdir 't' if -d 't' && !$ENV{PERL_CORE}; 10 unshift @INC, '../lib'; # for running manually 11 unshift @INC, '../blib/arch'; # for running manually 12 plan tests => 359; 13 } 14 15use Math::BigInt::FastCalc; 16 17my ($BASE_LEN, $AND_BITS, $XOR_BITS, $OR_BITS, $BASE_LEN_SMALL, $MAX_VAL) = 18 Math::BigInt::FastCalc->_base_len(); 19 20print "# BASE_LEN = $BASE_LEN\n"; 21print "# MAX_VAL = $MAX_VAL\n"; 22print "# AND_BITS = $AND_BITS\n"; 23print "# XOR_BITS = $XOR_BITS\n"; 24print "# IOR_BITS = $OR_BITS\n"; 25 26# testing of Math::BigInt::FastCalc 27 28my $C = 'Math::BigInt::FastCalc'; # pass classname to sub's 29 30# _new and _str 31my $x = $C->_new("123"); my $y = $C->_new("321"); 32ok (ref($x),'ARRAY'); ok ($C->_str($x),123); ok ($C->_str($y),321); 33 34############################################################################### 35# _add, _sub, _mul, _div 36ok ($C->_str($C->_add($x,$y)),444); 37ok ($C->_str($C->_sub($x,$y)),123); 38ok ($C->_str($C->_mul($x,$y)),39483); 39ok ($C->_str($C->_div($x,$y)),123); 40 41############################################################################### 42# check that mul/div doesn't change $y 43# and returns the same reference, not something new 44ok ($C->_str($C->_mul($x,$y)),39483); 45ok ($C->_str($x),39483); ok ($C->_str($y),321); 46 47ok ($C->_str($C->_div($x,$y)),123); 48ok ($C->_str($x),123); ok ($C->_str($y),321); 49 50$x = $C->_new("39483"); 51my ($x1,$r1) = $C->_div($x,$y); 52ok ("$x1","$x"); 53$C->_inc($x1); 54ok ("$x1","$x"); 55ok ($C->_str($r1),'0'); 56 57$x = $C->_new("39483"); # reset 58 59############################################################################### 60my $z = $C->_new("2"); 61ok ($C->_str($C->_add($x,$z)),39485); 62my ($re,$rr) = $C->_div($x,$y); 63 64ok ($C->_str($re),123); ok ($C->_str($rr),2); 65 66# is_zero, _is_one, _one, _zero 67ok ($C->_is_zero($x)||0,0); 68ok ($C->_is_one($x)||0,0); 69 70ok ($C->_str($C->_zero()),"0"); 71ok ($C->_str($C->_one()),"1"); 72 73# _two() and _ten() 74ok ($C->_str($C->_two()),"2"); 75ok ($C->_str($C->_ten()),"10"); 76ok ($C->_is_ten($C->_two())||0,0); 77ok ($C->_is_two($C->_two()),1); 78ok ($C->_is_ten($C->_ten()),1); 79ok ($C->_is_two($C->_ten())||0,0); 80 81ok ($C->_is_one($C->_one()),1); 82ok ($C->_is_one($C->_two()) || 0,0); 83ok ($C->_is_one($C->_ten()) || 0,0); 84 85ok ($C->_is_one($C->_zero()) || 0,0); 86 87ok ($C->_is_zero($C->_zero()),1); 88 89ok ($C->_is_zero($C->_one()) || 0,0); 90 91# is_odd, is_even 92ok ($C->_is_odd($C->_one()),1); ok ($C->_is_odd($C->_zero())||0,0); 93ok ($C->_is_even($C->_one()) || 0,0); ok ($C->_is_even($C->_zero()),1); 94 95# _len 96for my $method (qw/_alen _len/) 97 { 98 $x = $C->_new("1"); ok ($C->$method($x),1); 99 $x = $C->_new("12"); ok ($C->$method($x),2); 100 $x = $C->_new("123"); ok ($C->$method($x),3); 101 $x = $C->_new("1234"); ok ($C->$method($x),4); 102 $x = $C->_new("12345"); ok ($C->$method($x),5); 103 $x = $C->_new("123456"); ok ($C->$method($x),6); 104 $x = $C->_new("1234567"); ok ($C->$method($x),7); 105 $x = $C->_new("12345678"); ok ($C->$method($x),8); 106 $x = $C->_new("123456789"); ok ($C->$method($x),9); 107 108 $x = $C->_new("8"); ok ($C->$method($x),1); 109 $x = $C->_new("21"); ok ($C->$method($x),2); 110 $x = $C->_new("321"); ok ($C->$method($x),3); 111 $x = $C->_new("4321"); ok ($C->$method($x),4); 112 $x = $C->_new("54321"); ok ($C->$method($x),5); 113 $x = $C->_new("654321"); ok ($C->$method($x),6); 114 $x = $C->_new("7654321"); ok ($C->$method($x),7); 115 $x = $C->_new("87654321"); ok ($C->$method($x),8); 116 $x = $C->_new("987654321"); ok ($C->$method($x),9); 117 118 $x = $C->_new("0"); ok ($C->$method($x),1); 119 $x = $C->_new("20"); ok ($C->$method($x),2); 120 $x = $C->_new("320"); ok ($C->$method($x),3); 121 $x = $C->_new("4320"); ok ($C->$method($x),4); 122 $x = $C->_new("54320"); ok ($C->$method($x),5); 123 $x = $C->_new("654320"); ok ($C->$method($x),6); 124 $x = $C->_new("7654320"); ok ($C->$method($x),7); 125 $x = $C->_new("87654320"); ok ($C->$method($x),8); 126 $x = $C->_new("987654320"); ok ($C->$method($x),9); 127 128 for (my $i = 1; $i < 9; $i++) 129 { 130 my $a = "$i" . '0' x ($i-1); 131 $x = $C->_new($a); 132 print "# Tried len '$a'\n" unless ok ($C->_len($x),$i); 133 } 134 } 135 136# _digit 137$x = $C->_new("123456789"); 138ok ($C->_digit($x,0),9); 139ok ($C->_digit($x,1),8); 140ok ($C->_digit($x,2),7); 141ok ($C->_digit($x,-1),1); 142ok ($C->_digit($x,-2),2); 143ok ($C->_digit($x,-3),3); 144 145# _copy 146foreach (qw/ 1 12 123 1234 12345 123456 1234567 12345678 123456789/) 147 { 148 $x = $C->_new("$_"); 149 ok ($C->_str($C->_copy($x)),"$_"); 150 ok ($C->_str($x),"$_"); # did _copy destroy original x? 151 } 152 153# _zeros 154$x = $C->_new("1256000000"); ok ($C->_zeros($x),6); 155$x = $C->_new("152"); ok ($C->_zeros($x),0); 156$x = $C->_new("123000"); ok ($C->_zeros($x),3); 157$x = $C->_new("0"); ok ($C->_zeros($x),0); 158 159# _lsft, _rsft 160$x = $C->_new("10"); $y = $C->_new("3"); 161ok ($C->_str($C->_lsft($x,$y,10)),10000); 162$x = $C->_new("20"); $y = $C->_new("3"); 163ok ($C->_str($C->_lsft($x,$y,10)),20000); 164 165$x = $C->_new("128"); $y = $C->_new("4"); 166ok ($C->_str($C->_lsft($x,$y,2)), 128 << 4); 167 168$x = $C->_new("1000"); $y = $C->_new("3"); 169ok ($C->_str($C->_rsft($x,$y,10)),1); 170$x = $C->_new("20000"); $y = $C->_new("3"); 171ok ($C->_str($C->_rsft($x,$y,10)),20); 172$x = $C->_new("256"); $y = $C->_new("4"); 173ok ($C->_str($C->_rsft($x,$y,2)),256 >> 4); 174 175$x = $C->_new("6411906467305339182857313397200584952398"); 176$y = $C->_new("45"); 177ok ($C->_str($C->_rsft($x,$y,10)),0); 178 179# _acmp 180$x = $C->_new("123456789"); 181$y = $C->_new("987654321"); 182ok ($C->_acmp($x,$y),-1); 183ok ($C->_acmp($y,$x),1); 184ok ($C->_acmp($x,$x),0); 185ok ($C->_acmp($y,$y),0); 186$x = $C->_new("12"); 187$y = $C->_new("12"); 188ok ($C->_acmp($x,$y),0); 189$x = $C->_new("21"); 190ok ($C->_acmp($x,$y),1); 191ok ($C->_acmp($y,$x),-1); 192$x = $C->_new("123456789"); 193$y = $C->_new("1987654321"); 194ok ($C->_acmp($x,$y),-1); 195ok ($C->_acmp($y,$x),+1); 196 197$x = $C->_new("1234567890123456789"); 198$y = $C->_new("987654321012345678"); 199ok ($C->_acmp($x,$y),1); 200ok ($C->_acmp($y,$x),-1); 201ok ($C->_acmp($x,$x),0); 202ok ($C->_acmp($y,$y),0); 203 204$x = $C->_new("1234"); 205$y = $C->_new("987654321012345678"); 206ok ($C->_acmp($x,$y),-1); 207ok ($C->_acmp($y,$x),1); 208ok ($C->_acmp($x,$x),0); 209ok ($C->_acmp($y,$y),0); 210 211# _modinv 212$x = $C->_new("8"); 213$y = $C->_new("5033"); 214my ($xmod,$sign) = $C->_modinv($x,$y); 215ok ($C->_str($xmod),'629'); # -629 % 5033 == 4404 216ok ($sign, '-'); 217 218# _div 219$x = $C->_new("3333"); $y = $C->_new("1111"); 220ok ($C->_str(scalar $C->_div($x,$y)),3); 221$x = $C->_new("33333"); $y = $C->_new("1111"); ($x,$y) = $C->_div($x,$y); 222ok ($C->_str($x),30); ok ($C->_str($y),3); 223$x = $C->_new("123"); $y = $C->_new("1111"); 224($x,$y) = $C->_div($x,$y); ok ($C->_str($x),0); ok ($C->_str($y),123); 225 226# _num 227foreach (qw/1 12 123 1234 12345 1234567 12345678 123456789 1234567890/) 228 { 229 $x = $C->_new("$_"); 230 ok (ref($x)||'','ARRAY'); ok ($C->_str($x),"$_"); 231 $x = $C->_num($x); ok (ref($x)||'',''); ok ($x,$_); 232 } 233 234# _sqrt 235$x = $C->_new("144"); ok ($C->_str($C->_sqrt($x)),'12'); 236$x = $C->_new("144000000000000"); ok ($C->_str($C->_sqrt($x)),'12000000'); 237 238# _root 239$x = $C->_new("81"); my $n = $C->_new("3"); # 4*4*4 = 64, 5*5*5 = 125 240ok ($C->_str($C->_root($x,$n)),'4'); # 4.xx => 4.0 241$x = $C->_new("81"); $n = $C->_new("4"); # 3*3*3*3 == 81 242ok ($C->_str($C->_root($x,$n)),'3'); 243 244# _pow (and _root) 245$x = $C->_new("0"); $n = $C->_new("3"); # 0 ** y => 0 246ok ($C->_str($C->_pow($x,$n)), 0); 247$x = $C->_new("3"); $n = $C->_new("0"); # x ** 0 => 1 248ok ($C->_str($C->_pow($x,$n)), 1); 249$x = $C->_new("1"); $n = $C->_new("3"); # 1 ** y => 1 250ok ($C->_str($C->_pow($x,$n)), 1); 251$x = $C->_new("5"); $n = $C->_new("1"); # x ** 1 => x 252ok ($C->_str($C->_pow($x,$n)), 5); 253 254$x = $C->_new("81"); $n = $C->_new("3"); # 81 ** 3 == 531441 255ok ($C->_str($C->_pow($x,$n)),81 ** 3); 256 257ok ($C->_str($C->_root($x,$n)),81); 258 259$x = $C->_new("81"); 260ok ($C->_str($C->_pow($x,$n)),81 ** 3); 261ok ($C->_str($C->_pow($x,$n)),'150094635296999121'); # 531441 ** 3 == 262 263ok ($C->_str($C->_root($x,$n)),'531441'); 264ok ($C->_str($C->_root($x,$n)),'81'); 265 266$x = $C->_new("81"); $n = $C->_new("14"); 267ok ($C->_str($C->_pow($x,$n)),'523347633027360537213511521'); 268ok ($C->_str($C->_root($x,$n)),'81'); 269 270$x = $C->_new("523347633027360537213511520"); 271ok ($C->_str($C->_root($x,$n)),'80'); 272 273$x = $C->_new("523347633027360537213511522"); 274ok ($C->_str($C->_root($x,$n)),'81'); 275 276my $res = [ qw/ 9 31 99 316 999 3162 9999/ ]; 277 278# 99 ** 2 = 9801, 999 ** 2 = 998001 etc 279for my $i (2 .. 9) 280 { 281 $x = '9' x $i; $x = $C->_new($x); 282 $n = $C->_new("2"); 283 my $rc = '9' x ($i-1). '8' . '0' x ($i-1) . '1'; 284 print "# _pow( ", '9' x $i, ", 2) \n" unless 285 ok ($C->_str($C->_pow($x,$n)),$rc); 286 287 if ($i <= 7) 288 { 289 $x = '9' x $i; $x = $C->_new($x); 290 $n = '9' x $i; $n = $C->_new($n); 291 print "# _root( ", '9' x $i, ", ", 9 x $i, ") \n" unless 292 ok ($C->_str($C->_root($x,$n)),'1'); 293 294 $x = '9' x $i; $x = $C->_new($x); 295 $n = $C->_new("2"); 296 print "# _root( ", '9' x $i, ", ", 9 x $i, ") \n" unless 297 ok ($C->_str($C->_root($x,$n)), $res->[$i-2]); 298 } 299 } 300 301############################################################################## 302# _fac 303$x = $C->_new("0"); ok ($C->_str($C->_fac($x)),'1'); 304$x = $C->_new("1"); ok ($C->_str($C->_fac($x)),'1'); 305$x = $C->_new("2"); ok ($C->_str($C->_fac($x)),'2'); 306$x = $C->_new("3"); ok ($C->_str($C->_fac($x)),'6'); 307$x = $C->_new("4"); ok ($C->_str($C->_fac($x)),'24'); 308$x = $C->_new("5"); ok ($C->_str($C->_fac($x)),'120'); 309$x = $C->_new("10"); ok ($C->_str($C->_fac($x)),'3628800'); 310$x = $C->_new("11"); ok ($C->_str($C->_fac($x)),'39916800'); 311$x = $C->_new("12"); ok ($C->_str($C->_fac($x)),'479001600'); 312$x = $C->_new("13"); ok ($C->_str($C->_fac($x)),'6227020800'); 313 314# test that _fac modifes $x in place for small arguments 315$x = $C->_new("3"); $C->_fac($x); ok ($C->_str($x),'6'); 316$x = $C->_new("13"); $C->_fac($x); ok ($C->_str($x),'6227020800'); 317 318############################################################################## 319# _inc and _dec 320foreach (qw/1 11 121 1231 12341 1234561 12345671 123456781 1234567891/) 321 { 322 $x = $C->_new("$_"); $C->_inc($x); 323 print "# \$x = ",$C->_str($x),"\n" 324 unless ok ($C->_str($x),substr($_,0,length($_)-1) . '2'); 325 $C->_dec($x); ok ($C->_str($x),$_); 326 } 327foreach (qw/19 119 1219 12319 1234519 12345619 123456719 1234567819/) 328 { 329 $x = $C->_new("$_"); $C->_inc($x); 330 print "# \$x = ",$C->_str($x),"\n" 331 unless ok ($C->_str($x),substr($_,0,length($_)-2) . '20'); 332 $C->_dec($x); ok ($C->_str($x),$_); 333 } 334foreach (qw/999 9999 99999 9999999 99999999 999999999 9999999999 99999999999/) 335 { 336 $x = $C->_new("$_"); $C->_inc($x); 337 print "# \$x = ",$C->_str($x),"\n" 338 unless ok ($C->_str($x), '1' . '0' x (length($_))); 339 $C->_dec($x); ok ($C->_str($x),$_); 340 } 341 342$x = $C->_new("1000"); $C->_inc($x); ok ($C->_str($x),'1001'); 343$C->_dec($x); ok ($C->_str($x),'1000'); 344 345my $BL; 346{ 347 no strict 'refs'; 348 $BL = &{"$C"."::_base_len"}(); 349} 350 351$x = '1' . '0' x $BL; 352$z = '1' . '0' x ($BL-1); $z .= '1'; 353$x = $C->_new($x); $C->_inc($x); ok ($C->_str($x),$z); 354 355$x = '1' . '0' x $BL; $z = '9' x $BL; 356$x = $C->_new($x); $C->_dec($x); ok ($C->_str($x),$z); 357 358# should not happen: 359# $x = $C->_new("-2"); $y = $C->_new("4"); ok ($C->_acmp($x,$y),-1); 360 361############################################################################### 362# _mod 363$x = $C->_new("1000"); $y = $C->_new("3"); 364ok ($C->_str(scalar $C->_mod($x,$y)),1); 365$x = $C->_new("1000"); $y = $C->_new("2"); 366ok ($C->_str(scalar $C->_mod($x,$y)),0); 367 368# _and, _or, _xor 369$x = $C->_new("5"); $y = $C->_new("2"); 370ok ($C->_str(scalar $C->_xor($x,$y)),7); 371$x = $C->_new("5"); $y = $C->_new("2"); 372ok ($C->_str(scalar $C->_or($x,$y)),7); 373$x = $C->_new("5"); $y = $C->_new("3"); 374ok ($C->_str(scalar $C->_and($x,$y)),1); 375 376# _from_hex, _from_bin, _from_oct 377ok ($C->_str( $C->_from_hex("0xFf")),255); 378ok ($C->_str( $C->_from_bin("0b10101011")),160+11); 379ok ($C->_str( $C->_from_oct("0100")), 8*8); 380ok ($C->_str( $C->_from_oct("01000")), 8*8*8); 381ok ($C->_str( $C->_from_oct("010001")), 8*8*8*8+1); 382ok ($C->_str( $C->_from_oct("010007")), 8*8*8*8+7); 383 384# _as_hex, _as_bin, as_oct 385ok ($C->_str( $C->_from_hex( $C->_as_hex( $C->_new("128")))), 128); 386ok ($C->_str( $C->_from_bin( $C->_as_bin( $C->_new("128")))), 128); 387ok ($C->_str( $C->_from_oct( $C->_as_oct( $C->_new("128")))), 128); 388 389ok ($C->_str( $C->_from_oct( $C->_as_oct( $C->_new("123456")))), 123456); 390ok ($C->_str( $C->_from_oct( $C->_as_oct( $C->_new("123456789")))), "123456789"); 391ok ($C->_str( $C->_from_oct( $C->_as_oct( $C->_new("1234567890123")))), "1234567890123"); 392 393# _1ex 394ok ($C->_str($C->_1ex(0)), "1"); 395ok ($C->_str($C->_1ex(1)), "10"); 396ok ($C->_str($C->_1ex(2)), "100"); 397ok ($C->_str($C->_1ex(12)), "1000000000000"); 398ok ($C->_str($C->_1ex(16)), "10000000000000000"); 399 400# _check 401$x = $C->_new("123456789"); 402ok ($C->_check($x),0); 403ok ($C->_check(123),'123 is not a reference'); 404 405############################################################################### 406# __strip_zeros 407 408{ 409 no strict 'refs'; 410 # correct empty arrays 411 $x = &{$C."::__strip_zeros"}([]); ok (@$x,1); ok ($x->[0],0); 412 # don't strip single elements 413 $x = &{$C."::__strip_zeros"}([0]); ok (@$x,1); ok ($x->[0],0); 414 $x = &{$C."::__strip_zeros"}([1]); ok (@$x,1); ok ($x->[0],1); 415 # don't strip non-zero elements 416 $x = &{$C."::__strip_zeros"}([0,1]); 417 ok (@$x,2); ok ($x->[0],0); ok ($x->[1],1); 418 $x = &{$C."::__strip_zeros"}([0,1,2]); 419 ok (@$x,3); ok ($x->[0],0); ok ($x->[1],1); ok ($x->[2],2); 420 421 # but strip leading zeros 422 $x = &{$C."::__strip_zeros"}([0,1,2,0]); 423 ok (@$x,3); ok ($x->[0],0); ok ($x->[1],1); ok ($x->[2],2); 424 425 $x = &{$C."::__strip_zeros"}([0,1,2,0,0]); 426 ok (@$x,3); ok ($x->[0],0); ok ($x->[1],1); ok ($x->[2],2); 427 428 $x = &{$C."::__strip_zeros"}([0,1,2,0,0,0]); 429 ok (@$x,3); ok ($x->[0],0); ok ($x->[1],1); ok ($x->[2],2); 430 431 # collapse multiple zeros 432 $x = &{$C."::__strip_zeros"}([0,0,0,0]); 433 ok (@$x,1); ok ($x->[0],0); 434} 435 436# done 437 4381; 439 440