2013-07-07 10:01:40 -04:00
|
|
|
require 'test/unit'
|
|
|
|
require "-test-/bignum"
|
|
|
|
|
|
|
|
class TestBignum < Test::Unit::TestCase
|
|
|
|
class TestMul < Test::Unit::TestCase
|
|
|
|
|
|
|
|
SIZEOF_BDIGITS = Bignum::SIZEOF_BDIGITS
|
|
|
|
BITSPERDIG = Bignum::BITSPERDIG
|
2013-07-07 11:21:26 -04:00
|
|
|
BDIGMAX = (1 << BITSPERDIG) - 1
|
2013-07-07 10:01:40 -04:00
|
|
|
|
|
|
|
def test_mul_normal
|
|
|
|
x = (1 << BITSPERDIG) | 1
|
|
|
|
y = (1 << BITSPERDIG) | 1
|
|
|
|
z = (1 << (BITSPERDIG*2)) | (2 << BITSPERDIG) | 1
|
|
|
|
assert_equal(z, x.big_mul_normal(y))
|
|
|
|
end
|
|
|
|
|
|
|
|
def test_mul_normal_zero_in_x
|
|
|
|
x = (1 << (2*BITSPERDIG)) | 1
|
|
|
|
y = (1 << BITSPERDIG) | 1
|
|
|
|
z = (1 << (BITSPERDIG*3)) | (1 << (BITSPERDIG*2)) | (1 << BITSPERDIG) | 1
|
|
|
|
assert_equal(z, x.big_mul_normal(y))
|
|
|
|
end
|
|
|
|
|
|
|
|
def test_mul_normal_zero_in_y
|
|
|
|
x = (1 << BITSPERDIG) | 1
|
|
|
|
y = (1 << (2*BITSPERDIG)) | 1
|
|
|
|
z = (1 << (BITSPERDIG*3)) | (1 << (BITSPERDIG*2)) | (1 << BITSPERDIG) | 1
|
|
|
|
assert_equal(z, x.big_mul_normal(y))
|
|
|
|
end
|
|
|
|
|
|
|
|
def test_mul_normal_max_max
|
|
|
|
x = (1 << (2*BITSPERDIG)) - 1
|
|
|
|
y = (1 << (2*BITSPERDIG)) - 1
|
|
|
|
z = (1 << (4*BITSPERDIG)) - (1 << (2*BITSPERDIG+1)) + 1
|
|
|
|
assert_equal(z, x.big_mul_normal(y))
|
|
|
|
end
|
|
|
|
|
|
|
|
def test_mul_balance
|
|
|
|
x = (1 << BITSPERDIG) | 1
|
|
|
|
y = (1 << BITSPERDIG) | 1
|
|
|
|
z = (1 << (BITSPERDIG*2)) | (2 << BITSPERDIG) | 1
|
|
|
|
assert_equal(z, x.big_mul_balance(y))
|
|
|
|
end
|
|
|
|
|
2013-07-08 07:56:55 -04:00
|
|
|
def test_mul_balance_2x16
|
|
|
|
x = (1 << Bignum::BITSPERDIG) | 1
|
|
|
|
y = (1 << Bignum::BITSPERDIG*16) | 1
|
|
|
|
assert_equal(x.big_mul_normal(y), x.big_mul_balance(y))
|
|
|
|
end
|
|
|
|
|
|
|
|
def test_mul_balance_2x17
|
|
|
|
x = (1 << Bignum::BITSPERDIG) | 1
|
|
|
|
y = (1 << Bignum::BITSPERDIG*17) | 1
|
|
|
|
assert_equal(x.big_mul_normal(y), x.big_mul_balance(y))
|
|
|
|
end
|
|
|
|
|
2013-07-07 10:01:40 -04:00
|
|
|
def test_mul_karatsuba
|
|
|
|
x = (1 << BITSPERDIG) | 1
|
|
|
|
y = (1 << BITSPERDIG) | 1
|
|
|
|
z = (1 << (BITSPERDIG*2)) | (2 << BITSPERDIG) | 1
|
|
|
|
assert_equal(z, x.big_mul_karatsuba(y))
|
|
|
|
end
|
|
|
|
|
2013-07-07 11:21:26 -04:00
|
|
|
def test_mul_karatsuba_odd_y
|
|
|
|
x = (1 << BITSPERDIG) | 1
|
|
|
|
y = (1 << (2*BITSPERDIG)) | 1
|
|
|
|
assert_equal(x.big_mul_normal(y), x.big_mul_karatsuba(y))
|
|
|
|
end
|
|
|
|
|
|
|
|
def test_mul_karatsuba_odd_xy
|
|
|
|
x = (1 << (2*BITSPERDIG)) | 1
|
|
|
|
y = (1 << (2*BITSPERDIG)) | 1
|
|
|
|
assert_equal(x.big_mul_normal(y), x.big_mul_karatsuba(y))
|
|
|
|
end
|
|
|
|
|
|
|
|
def test_mul_karatsuba_x1_gt_x0
|
|
|
|
x = (2 << BITSPERDIG) | 1
|
|
|
|
y = (1 << BITSPERDIG) | 2
|
|
|
|
assert_equal(x.big_mul_normal(y), x.big_mul_karatsuba(y))
|
|
|
|
end
|
|
|
|
|
|
|
|
def test_mul_karatsuba_y1_gt_y0
|
|
|
|
x = (1 << BITSPERDIG) | 2
|
|
|
|
y = (2 << BITSPERDIG) | 1
|
|
|
|
assert_equal(x.big_mul_normal(y), x.big_mul_karatsuba(y))
|
|
|
|
end
|
|
|
|
|
|
|
|
def test_mul_karatsuba_x1_gt_x0_and_y1_gt_y0
|
|
|
|
x = (2 << BITSPERDIG) | 1
|
|
|
|
y = (2 << BITSPERDIG) | 1
|
|
|
|
assert_equal(x.big_mul_normal(y), x.big_mul_karatsuba(y))
|
|
|
|
end
|
|
|
|
|
|
|
|
def test_mul_karatsuba_carry2
|
|
|
|
x = (1 << BITSPERDIG) | BDIGMAX
|
|
|
|
y = (1 << BITSPERDIG) | BDIGMAX
|
|
|
|
assert_equal(x.big_mul_normal(y), x.big_mul_karatsuba(y))
|
|
|
|
end
|
|
|
|
|
|
|
|
def test_mul_karatsuba_borrow
|
|
|
|
x = (BDIGMAX << BITSPERDIG) | 1
|
|
|
|
y = (BDIGMAX << BITSPERDIG) | 1
|
|
|
|
assert_equal(x.big_mul_normal(y), x.big_mul_karatsuba(y))
|
|
|
|
end
|
|
|
|
|
2013-07-07 10:01:40 -04:00
|
|
|
end
|
|
|
|
end
|