Commit | Line | Data |
56b9c951 |
1 | #!/usr/bin/perl -w |
2 | |
3 | use Test; |
4 | use strict; |
5 | |
6 | my $count; |
7 | |
8 | BEGIN |
9 | { |
10 | $| = 1; |
11 | unshift @INC, '../lib'; # for running manually |
12 | my $location = $0; $location =~ s/mbi_rand.t//; |
13 | unshift @INC, $location; # to locate the testing files |
14 | chdir 't' if -d 't'; |
15 | $count = 500; |
16 | plan tests => $count*2; |
17 | } |
18 | |
19 | use Math::BigInt; |
20 | my $c = 'Math::BigInt'; |
21 | |
22 | my $length = 200; |
23 | |
24 | # If you get a failure here, please re-run the test with the printed seed |
25 | # value as input: perl t/mbi_rand.t seed |
26 | |
27 | my $seed = int(rand(65537)); print "# seed: $seed\n"; srand($seed); |
28 | |
29 | my ($A,$B,$ADB,$AMB,$la,$lb); |
30 | for (my $i = 0; $i < $count; $i++) |
31 | { |
32 | # length of A and B |
33 | $la = int(rand($length)+1); $lb = int(rand($length)+1); |
34 | $A = ''; $B = ''; |
35 | # we create the numbers from "patterns", e.g. get a random number and a |
36 | # random count and string them together. This means things like |
37 | # "100000999999999999911122222222" are much more likely. If we just strung |
38 | # together digits, we would end up with "1272398823211223" etc. |
39 | while (length($A) < $la) { $A .= int(rand(100)) x int(rand(16)); } |
40 | while (length($B) < $lb) { $B .= int(rand(100)) x int(rand(16)); } |
41 | $A = $c->new($A); $B = $c->new($B); |
8f675a64 |
42 | # print "# A $A\n# B $B\n"; |
56b9c951 |
43 | if ($A->is_zero() || $B->is_zero()) |
44 | { |
45 | ok (1,1); ok (1,1); next; |
46 | } |
47 | # check that int(A/B)*B + A % B == A holds for all inputs |
48 | # $X = ($A/$B)*$B + 2 * ($A % $B) - ($A % $B); |
49 | ($ADB,$AMB) = $A->copy()->bdiv($B); |
8f675a64 |
50 | print "# ". join(' ',Math::BigInt::Calc->_base_len()),"\n" |
51 | unless ok ($ADB*$B+2*$AMB-$AMB,$A); |
56b9c951 |
52 | # swap 'em and try this, too |
53 | # $X = ($B/$A)*$A + $B % $A; |
54 | ($ADB,$AMB) = $B->copy()->bdiv($A); |
8f675a64 |
55 | print "# ". join(' ',Math::BigInt::Calc->_base_len()),"\n" |
56 | unless ok ($ADB*$A+2*$AMB-$AMB,$B); |
56b9c951 |
57 | } |
58 | |