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