Gecko [ www.blackthorn-inc.com ]

Name	Size	Permission	Date
Math	[ DIR ]	drwxr-xr-x	2015-02-14 16:56
_e_math.t	2.22 KB	-r--r--r--	2014-12-27 11:48
alias.inc	401 B	-r--r--r--	2014-12-27 11:48
bare_mbf.t	321 B	-r--r--r--	2014-12-27 11:48
bare_mbi.t	447 B	-r--r--r--	2014-12-27 11:48
bare_mif.t	502 B	-r--r--r--	2014-12-27 11:48
big_pi_e.t	689 B	-r--r--r--	2014-12-27 11:48
bigfltpm.inc	35.83 KB	-r--r--r--	2014-12-27 11:49
bigfltpm.t	846 B	-r--r--r--	2014-12-27 11:48
bigintc.t	14.24 KB	-r--r--r--	2014-12-27 11:48
bigintpm.inc	55.94 KB	-r--r--r--	2014-12-27 11:48
bigintpm.t	1.1 KB	-r--r--r--	2014-12-27 11:48
bigints.t	2.66 KB	-r--r--r--	2014-12-27 11:48
biglog.t	6.69 KB	-r--r--r--	2014-12-27 11:49
bigroot.t	1.41 KB	-r--r--r--	2014-12-27 11:48
calling.t	2.81 KB	-r--r--r--	2014-12-27 11:48
config.t	3.36 KB	-r--r--r--	2014-12-27 11:48
const_mbf.t	301 B	-r--r--r--	2014-12-27 11:48
constant.t	907 B	-r--r--r--	2014-12-27 11:48
downgrade.t	1.49 KB	-r--r--r--	2014-12-27 11:48
inf_nan.t	5.5 KB	-r--r--r--	2014-12-27 11:48
isa.t	886 B	-r--r--r--	2014-12-27 11:48
lib_load.t	881 B	-r--r--r--	2014-12-27 11:48
mbf_ali.t	189 B	-r--r--r--	2014-12-27 11:48
mbi_ali.t	185 B	-r--r--r--	2014-12-27 11:48
mbi_rand.t	2.81 KB	-r--r--r--	2014-12-27 11:48
mbimbf.inc	27.61 KB	-r--r--r--	2014-12-27 11:48
mbimbf.t	2.23 KB	-r--r--r--	2014-12-27 11:48
nan_cmp.t	805 B	-r--r--r--	2014-12-27 11:48
new_overloaded.t	739 B	-r--r--r--	2014-12-27 11:48
req_mbf0.t	228 B	-r--r--r--	2014-12-27 11:48
req_mbf1.t	208 B	-r--r--r--	2014-12-27 11:48
req_mbfa.t	213 B	-r--r--r--	2014-12-27 11:48
req_mbfi.t	212 B	-r--r--r--	2014-12-27 11:48
req_mbfn.t	213 B	-r--r--r--	2014-12-27 11:48
req_mbfw.t	632 B	-r--r--r--	2014-12-27 11:48
require.t	200 B	-r--r--r--	2014-12-27 11:48
round.t	2.7 KB	-r--r--r--	2014-12-27 11:48
rt-16221.t	1.52 KB	-r--r--r--	2014-12-27 11:48
sub_ali.t	234 B	-r--r--r--	2014-12-27 11:48
sub_mbf.t	1.03 KB	-r--r--r--	2014-12-27 11:48
sub_mbi.t	1 KB	-r--r--r--	2014-12-27 11:48
sub_mif.t	356 B	-r--r--r--	2014-12-27 11:48
trap.t	2.35 KB	-r--r--r--	2014-12-27 11:48
upgrade.inc	26.83 KB	-r--r--r--	2014-12-27 11:49
upgrade.t	475 B	-r--r--r--	2014-12-27 11:48
upgrade2.t	372 B	-r--r--r--	2014-12-27 11:48
upgradef.t	1.22 KB	-r--r--r--	2014-12-27 11:48
use.t	406 B	-r--r--r--	2014-12-27 11:48
use_lib1.t	346 B	-r--r--r--	2014-12-27 11:48
use_lib2.t	364 B	-r--r--r--	2014-12-27 11:48
use_lib3.t	364 B	-r--r--r--	2014-12-27 11:48
use_lib4.t	441 B	-r--r--r--	2014-12-27 11:48
use_mbfw.t	915 B	-r--r--r--	2014-12-27 11:48
with_sub.t	479 B	-r--r--r--	2014-12-27 11:48

Rename

#!/usr/bin/perl -w

use strict;

my $count;

BEGIN {
  $count = 128;
}

use Test::More tests => $count*4;

use Math::BigInt;
my $c = 'Math::BigInt';

my $length = 128;

# If you get a failure here, please re-run the test with the printed seed
# value as input "perl t/mbi_rand.t seed" and send me the output

my $seed = ($#ARGV == 0) ? $ARGV[0] : int(rand(1165537));
print "# seed: $seed\n"; srand($seed);

print "# lib: ", Math::BigInt->config()->{lib},"\n";
if (Math::BigInt->config()->{lib} =~ /::Calc/)
  {
  print "# base len: ", scalar Math::BigInt::Calc->_base_len(),"\n";
  }

my ($A,$B,$As,$Bs,$ADB,$AMB,$la,$lb);
my $two = Math::BigInt->new(2);
for (my $i = 0; $i < $count; $i++)
  {
  # length of A and B
  $la = int(rand($length)+1); $lb = int(rand($length)+1);
  $As = ''; $Bs = '';

# we create the numbers from "patterns", e.g. get a random number and a
  # random count and string them together. This means things like
  # "100000999999999999911122222222" are much more likely. If we just strung
  # together digits, we would end up with "1272398823211223" etc. It also means
  # that we get more frequently equal numbers or other special cases.
  while (length($As) < $la) { $As .= int(rand(100)) x int(rand(16)); }
  while (length($Bs) < $lb) { $Bs .= int(rand(100)) x int(rand(16)); }

$As =~ s/^0+//; $Bs =~ s/^0+//; 
  $As = $As || '0'; $Bs = $Bs || '0';
#  print "# As $As\n# Bs $Bs\n";
  $A = $c->new($As); $B = $c->new($Bs);
  print "# A $A\n# B $B\n";
  if ($A->is_zero() || $B->is_zero())
    {
    for (1..4) { is (1,1, 'skipped this test'); } next;
    }

# check that int(A/B)*B + A % B == A holds for all inputs

# $X = ($A/$B)*$B + 2 * ($A % $B) - ($A % $B);

($ADB,$AMB) = $A->copy()->bdiv($B);
  print "# ($A / $B, $A % $B ) = $ADB $AMB\n";

print "# seed $seed, ". join(' ',Math::BigInt::Calc->_base_len()),"\n".
        "# tried $ADB * $B + $two*$AMB - $AMB\n"
   unless is ($ADB*$B+$two*$AMB-$AMB,$As, "ADB * B + 2 * AMB - AMB == A");
  if (is ($ADB*$B/$B,$ADB, "ADB * B / B == ADB"))
    {
    print "# seed: $seed, \$ADB * \$B / \$B = ", $ADB * $B / $B, " != $ADB (\$B=$B)\n";
    if (Math::BigInt->config()->{lib} =~ /::Calc/)
      {
      print "# ADB->[-1]: ", $ADB->{value}->[-1], " B->[-1]: ", $B->{value}->[-1],"\n";
      }
    }
  # swap 'em and try this, too
  # $X = ($B/$A)*$A + $B % $A;
  ($ADB,$AMB) = $B->copy()->bdiv($A);
  # print "check: $ADB $AMB";
  print "# seed $seed, ". join(' ',Math::BigInt::Calc->_base_len()),"\n".
        "# tried $ADB * $A + $two*$AMB - $AMB\n"
   unless is ($ADB*$A+$two*$AMB-$AMB,$Bs, "ADB * A + 2 * AMB - AMB == B");
  print "# +$two * $AMB = ",$ADB * $A + $two * $AMB,"\n";
  print "# -$AMB = ",$ADB * $A + $two * $AMB - $AMB,"\n";
  print "# seed $seed, \$ADB * \$A / \$A = ", $ADB * $A / $A, " != $ADB (\$A=$A)\n"
   unless is ($ADB*$A/$A,$ADB, "ADB * A/A == ADB");
  }