def test02_CarryBorrow(self):
if verbose:
print '\n', '-=' * 30
print "Running %s.test02_CarryBorrow..." % \
self.__class__.__name__
a = 0x7fffffff
wa = -mpw(a)
wa = wa + wa
za = -long(a)
za = za + za
wb = -mpw(1)
zb = -long(1)
wc = mpw(1)
zc = long(1)
wd = mpw(a)
wd = wd + wd
zd = long(a)
zd = zd + zd
print "add --:\t", (wa + wa), "\t", (za + za)
print "add -+:\t", (wb + wd), "\t", (zb + zd)
print "add +-:\t", (wc + wa), "\t", (zc + za)
print "add ++:\t", (wd + wd), "\t", (zd + zd)
print "sub --:\t", (wb - wa), "\t", (zb - za)
# print "sub -+:\t", (wb-wd), "\t", (zb-zd)
# print "sub +-:\t", (wc-wa), "\t", (zc-za)
print "sub ++:\t", (wc - wd), "\t", (zc - zd)
pass
def test04_KnuthPoly(self):
self.t = 8
tfmt = "%o"
if verbose:
print '\n', '-=' * 30
print "Running %s.test04_KnuthPoly..." % \
self.__class__.__name__
print "\t(%d**m - 1) * (%d**n - 1), m,n in [%d,%d)" % (
self.t, self.t, self.lo, self.hi)
tm1 = tfmt % (self.t - 1)
tm2 = tfmt % (self.t - 2)
for m in range(self.lo, self.hi):
for n in range(m + 1, self.hi + 1):
wt = mpw(self.t)
wa = (wt**m - 1) * (wt**n - 1)
ws = tfmt % long(wa)
zs = tm1 * (m - 1) + tm2 + tm1 * (n - m) + "0" * (m - 1) + "1"
if ws != zs:
print "(%d**%d - 1) * (%d**%d - 1)\t%s" % (self.t, m,
self.t, n, ws)
assert ws == zs
self.t = 10
tfmt = "%d"
if verbose:
print "\t(%d**m - 1) * (%d**n - 1), m,n in [%d,%d)" % (
self.t, self.t, self.lo, self.hi)
tm1 = tfmt % (self.t - 1)
tm2 = tfmt % (self.t - 2)
for m in range(self.lo, self.hi):
for n in range(m + 1, self.hi + 1):
wt = mpw(self.t)
wa = (wt**m - 1) * (wt**n - 1)
ws = tfmt % long(wa)
zs = tm1 * (m - 1) + tm2 + tm1 * (n - m) + "0" * (m - 1) + "1"
if ws != zs:
print "(%d**%d - 1) * (%d**%d - 1)\t%s" % (self.t, m,
self.t, n, ws)
assert ws == zs
self.t = 16
tfmt = "%x"
if verbose:
print "\t(%d**m - 1) * (%d**n - 1), m,n in [%d,%d)" % (
self.t, self.t, self.lo, self.hi)
tm1 = tfmt % (self.t - 1)
tm2 = tfmt % (self.t - 2)
for m in range(self.lo, self.hi):
for n in range(m + 1, self.hi + 1):
wt = mpw(self.t)
wa = (wt**m - 1) * (wt**n - 1)
ws = tfmt % long(wa)
zs = tm1 * (m - 1) + tm2 + tm1 * (n - m) + "0" * (m - 1) + "1"
if ws != zs:
print "(%d**%d - 1) * (%d**%d - 1)\t%s" % (self.t, m,
self.t, n, ws)
assert ws == zs
pass
def test04_KnuthPoly(self):
self.t = 8
tfmt = "%o"
if verbose:
print '\n', '-=' * 30
print "Running %s.test04_KnuthPoly..." % \
self.__class__.__name__
print "\t(%d**m - 1) * (%d**n - 1), m,n in [%d,%d)" % (self.t,self.t,self.lo,self.hi)
tm1 = tfmt % (self.t - 1)
tm2 = tfmt % (self.t - 2)
for m in range(self.lo,self.hi):
for n in range(m+1,self.hi+1):
wt = mpw(self.t)
wa = (wt**m - 1) * (wt**n - 1)
ws = tfmt % long(wa)
zs = tm1 * (m - 1) + tm2 + tm1 * (n - m) + "0" * (m - 1) + "1"
if ws != zs:
print "(%d**%d - 1) * (%d**%d - 1)\t%s" % (self.t,m,self.t,n,ws)
assert ws == zs
self.t = 10
tfmt = "%d"
if verbose:
print "\t(%d**m - 1) * (%d**n - 1), m,n in [%d,%d)" % (self.t,self.t,self.lo,self.hi)
tm1 = tfmt % (self.t - 1)
tm2 = tfmt % (self.t - 2)
for m in range(self.lo,self.hi):
for n in range(m+1,self.hi+1):
wt = mpw(self.t)
wa = (wt**m - 1) * (wt**n - 1)
ws = tfmt % long(wa)
zs = tm1 * (m - 1) + tm2 + tm1 * (n - m) + "0" * (m - 1) + "1"
if ws != zs:
print "(%d**%d - 1) * (%d**%d - 1)\t%s" % (self.t,m,self.t,n,ws)
assert ws == zs
self.t = 16
tfmt = "%x"
if verbose:
print "\t(%d**m - 1) * (%d**n - 1), m,n in [%d,%d)" % (self.t,self.t,self.lo,self.hi)
tm1 = tfmt % (self.t - 1)
tm2 = tfmt % (self.t - 2)
for m in range(self.lo,self.hi):
for n in range(m+1,self.hi+1):
wt = mpw(self.t)
wa = (wt**m - 1) * (wt**n - 1)
ws = tfmt % long(wa)
zs = tm1 * (m - 1) + tm2 + tm1 * (n - m) + "0" * (m - 1) + "1"
if ws != zs:
print "(%d**%d - 1) * (%d**%d - 1)\t%s" % (self.t,m,self.t,n,ws)
assert ws == zs
pass
def test02_CarryBorrow(self):
if verbose:
print '\n', '-=' * 30
print "Running %s.test02_CarryBorrow..." % \
self.__class__.__name__
a = 0x7fffffff
wa = -mpw(a); wa = wa+wa
za = -long(a); za = za+za
wb = -mpw(1)
zb = -long(1)
wc = mpw(1)
zc = long(1)
wd = mpw(a); wd = wd+wd
zd = long(a); zd = zd+zd
print "add --:\t", (wa+wa), "\t", (za+za)
print "add -+:\t", (wb+wd), "\t", (zb+zd)
print "add +-:\t", (wc+wa), "\t", (zc+za)
print "add ++:\t", (wd+wd), "\t", (zd+zd)
print "sub --:\t", (wb-wa), "\t", (zb-za)
# print "sub -+:\t", (wb-wd), "\t", (zb-zd)
# print "sub +-:\t", (wc-wa), "\t", (zc-za)
print "sub ++:\t", (wc-wd), "\t", (zc-zd)
pass
def test05_IterativePowers(self):
if verbose:
print '\n', '-=' * 30
print "Running %s.test05_IterativePowers..." % \
self.__class__.__name__
print "\t(m**n)/(m**(n-1)) == m for m,n in [%d,%d)" % (self.lo,self.hi)
for m in range(self.lo,self.hi):
wa = mpw(m)
wd = wa
for n in range(self.lo,self.hi):
wc = wa**n
we = wc/wd
if we != m:
print m, '^', n, '=', we
assert we == m
if wc != 0:
wd = wc
pass
def test05_IterativePowers(self):
if verbose:
print '\n', '-=' * 30
print "Running %s.test05_IterativePowers..." % \
self.__class__.__name__
print "\t(m**n)/(m**(n-1)) == m for m,n in [%d,%d)" % (self.lo,
self.hi)
for m in range(self.lo, self.hi):
wa = mpw(m)
wd = wa
for n in range(self.lo, self.hi):
wc = wa**n
we = wc / wd
if we != m:
print m, '^', n, '=', we
assert we == m
if wc != 0:
wd = wc
pass
def test01_SimpleMethods(self):
if verbose:
print '\n', '-=' * 30
print "Running %s.test01_SimpleMethods..." % \
self.__class__.__name__
print "\ta:\t%s\t%s\t0x%x" % (type(self.a), self.a, self.a)
print "\tb:\t%s\t%s\t0x%x" % (type(self.b), self.b, self.b)
print "\tc:\t%s\t%s\t0x%x" % (type(self.c), self.c, self.c)
wa = mpw(self.a)
wb = mpw(self.b)
wc = mpw(self.c)
# xa - Long(self.a)
# xb = Long(self.b)
# xc = Long(self.c)
za = long(self.a)
zb = long(self.b)
zc = long(self.c)
print "__int__:\t", int(wb), "\t", int(zb)
assert int(wb) == int(zb)
print "__long__:\t", long(wa), "\t", long(za)
assert long(wb) == long(zb)
print "__float__:\t", float(wa), "\t", float(za)
assert float(wb) == float(zb)
zs = hex(za)
print "__hex__:\t", hex(wa), "\t", zs
# assert hex(wa) == zs # fails because of missing trailing L
zs = oct(za)
print "__oct__:\t", oct(wa), "\t", zs
# assert oct(wa) == zs # fails because of missing trailing L
print "__neg__:\t", (-wa), "\t", long(-za)
print "__pos__:\t", (+wa), "\t", long(+za)
print "__abs__:\t", abs(wa), "\t", long(abs(za))
print "__invert__:\t", (~wa), "\t", long(~za)
print "__add__:\t", (wa + wb), "\t", long(za + zb)
print "__sub__:\t", (wa - wb), "\t", long(za - zb)
print "__mul__:\t", (wa * wb), "\t", long(za * zb)
print "__div__:\t", (wa / wb), "\t", long(za / zb)
print "__mod__:\t", (wa % wb), "\t", long(za % zb)
wq, wr = divmod(wa, wb)
zq, zr = divmod(za, zb)
print "__divmod__ q:\t", wq, "\t", long(zq)
print "__divmod__ r:\t", wr, "\t", long(zr)
print "__pow__:\t", (wb ** wb), "\t", long(zb ** zb)
print "__lshift__:\t", (wa << wb), "\t", long(za << zb)
print "__rshift__:\t", (wa >> wb), "\t", long(za >> zb)
print "__and__:\t", (wa & wc), "\t", long(za & zc)
print "__xor__:\t", (wa ^ wa), "\t", long(za ^ za)
print "__or__:\t", (wa | wc), "\t", long(za | zc)
# print mpw.__complex__(b)
# print mpw.__coerce__(b, i)
del wa
del wb
del wc
del za
del zb
del zc
pass
def tearDown(self): mpw().Debug(0) pass
def setUp(self): mpw().Debug(0) pass
def test03_Signs(self):
if verbose:
print '\n', '-=' * 30
print "Running %s.test03_Signs..." % \
self.__class__.__name__
wpa = mpw(13)
wma = -wpa
wpb = wpa - 3
wmb = -wpb
zpa = long(13)
zma = -zpa
zpb = zpa - 3
zmb = -zpb
print "add --:\t", (wma+wmb), "\t", (zma+zmb)
print "add -+:\t", (wma+wpb), "\t", (zma+zpb)
print "add +-:\t", (wpa+wmb), "\t", (zpa+zmb)
print "add ++:\t", (wpa+wpb), "\t", (zpa+zpb)
print "sub --:\t", (wma-wmb), "\t", (zma-zmb)
print "sub -+:\t", (wma-wpb), "\t", (zma-zpb)
print "sub +-:\t", (wpa-wmb), "\t", (zpa-zmb)
print "sub ++:\t", (wpa-wpb), "\t", (zpa-zpb)
print "sub --:\t", (wmb-wma), "\t", (zmb-zma)
print "sub -+:\t", (wmb-wpa), "\t", (zmb-zpa)
print "sub +-:\t", (wpb-wma), "\t", (zpb-zma)
print "sub ++:\t", (wpb-wpa), "\t", (zpb-zpa)
print "mul --:\t", (wma*wmb), "\t", (zma*zmb)
print "mul -+:\t", (wma*wpb), "\t", (zma*zpb)
print "mul +-:\t", (wpa*wmb), "\t", (zpa*zmb)
print "mul ++:\t", (wpa*wpb), "\t", (zpa*zpb)
print "div --:\t", (wma/wmb), "\t", (zma/zmb)
print "div -+:\t", (wma/wpb), "\t", (zma/zpb)
print "div +-:\t", (wpa/wmb), "\t", (zpa/zmb)
print "div ++:\t", (wpa/wpb), "\t", (zpa/zpb)
print "div --:\t", (wmb/wma), "\t", (zmb/zma)
print "div -+:\t", (wmb/wpa), "\t", (zmb/zpa)
print "div +-:\t", (wpb/wma), "\t", (zpb/zma)
print "div ++:\t", (wpb/wpa), "\t", (zpb/zpa)
print "pow --:\t", (wma**wmb), "\t", (zma**zmb)
print "pow -+:\t", (wma**wpb), "\t", (zma**zpb)
print "pow +-:\t", (wpa**wmb), "\t", (zpa**zmb)
print "pow ++:\t", (wpa**wpb), "\t", (zpa**zpb)
print "pow --:\t", (wmb**wma), "\t", (zmb**zma)
print "pow -+:\t", (wmb**wpa), "\t", (zmb**zpa)
print "pow +-:\t", (wpb**wma), "\t", (zpb**zma)
print "pow ++:\t", (wpb**wpa), "\t", (zpb**zpa)
# wpa = mpw(13)
# wma = -wpa
# wpb = wpa - 3
# wmb = -wpb
# zpa = long(13)
# zma = -zpa
# zpb = zpa - 3
# zmb = -zpb
print "mod --:\t", (wma%wmb), "\t", (zma%zmb)
print "mod -+:\t", (wma%wpb), "\t", (zma%zpb)
print "mod +-:\t", (wpa%wmb), "\t", (zpa%zmb)
print "mod ++:\t", (wpa%wpb), "\t", (zpa%zpb)
print "mod --:\t", (wmb%wma), "\t", (zmb%zma)
print "mod -+:\t", (wmb%wpa), "\t", (zmb%zpa)
print "mod +-:\t", (wpb%wma), "\t", (zpb%zma)
print "mod ++:\t", (wpb%wpa), "\t", (zpb%zpa)
print "rem --:\t", divmod(wma, wmb), "\t", divmod(zma, zmb)
print "rem -+:\t", divmod(wma, wpb), "\t", divmod(zma, zpb)
print "rem +-:\t", divmod(wpa, wmb), "\t", divmod(zpa, zmb)
print "rem ++:\t", divmod(wpa, wpb), "\t", divmod(zpa, zpb)
print "rem --:\t", divmod(wmb, wma), "\t", divmod(zmb, zma)
print "rem -+:\t", divmod(wmb, wpa), "\t", divmod(zmb, zpa)
print "rem +-:\t", divmod(wpb, wma), "\t", divmod(zpb, zma)
print "rem ++:\t", divmod(wpb, wpa), "\t", divmod(zpb, zpa)
pass
def test01_SimpleMethods(self):
if verbose:
print '\n', '-=' * 30
print "Running %s.test01_SimpleMethods..." % \
self.__class__.__name__
print "\ta:\t%s\t%s\t0x%x" % (type(self.a), self.a, self.a)
print "\tb:\t%s\t%s\t0x%x" % (type(self.b), self.b, self.b)
print "\tc:\t%s\t%s\t0x%x" % (type(self.c), self.c, self.c)
wa = mpw(self.a)
wb = mpw(self.b)
wc = mpw(self.c)
# xa - Long(self.a)
# xb = Long(self.b)
# xc = Long(self.c)
za = long(self.a)
zb = long(self.b)
zc = long(self.c)
print "__int__:\t", int(wb), "\t", int(zb)
assert int(wb) == int(zb)
print "__long__:\t", long(wa), "\t", long(za)
assert long(wb) == long(zb)
print "__float__:\t", float(wa), "\t", float(za)
assert float(wb) == float(zb)
zs = hex(za)
print "__hex__:\t", hex(wa), "\t", zs
# assert hex(wa) == zs # fails because of missing trailing L
zs = oct(za)
print "__oct__:\t", oct(wa), "\t", zs
# assert oct(wa) == zs # fails because of missing trailing L
print "__neg__:\t", (-wa), "\t", long(-za)
print "__pos__:\t", (+wa), "\t", long(+za)
print "__abs__:\t", abs(wa), "\t", long(abs(za))
print "__invert__:\t", (~wa), "\t", long(~za)
print "__add__:\t", (wa + wb), "\t", long(za + zb)
print "__sub__:\t", (wa - wb), "\t", long(za - zb)
print "__mul__:\t", (wa * wb), "\t", long(za * zb)
print "__div__:\t", (wa / wb), "\t", long(za / zb)
print "__mod__:\t", (wa % wb), "\t", long(za % zb)
wq, wr = divmod(wa, wb)
zq, zr = divmod(za, zb)
print "__divmod__ q:\t", wq, "\t", long(zq)
print "__divmod__ r:\t", wr, "\t", long(zr)
print "__pow__:\t", (wb**wb), "\t", long(zb**zb)
print "__lshift__:\t", (wa << wb), "\t", long(za << zb)
print "__rshift__:\t", (wa >> wb), "\t", long(za >> zb)
print "__and__:\t", (wa & wc), "\t", long(za & zc)
print "__xor__:\t", (wa ^ wa), "\t", long(za ^ za)
print "__or__:\t", (wa | wc), "\t", long(za | zc)
# print mpw.__complex__(b)
# print mpw.__coerce__(b, i)
del wa
del wb
del wc
del za
del zb
del zc
pass
def tearDown(self):
mpw().Debug(0)
pass
def setUp(self):
mpw().Debug(0)
pass
def test03_Signs(self):
if verbose:
print '\n', '-=' * 30
print "Running %s.test03_Signs..." % \
self.__class__.__name__
wpa = mpw(13)
wma = -wpa
wpb = wpa - 3
wmb = -wpb
zpa = long(13)
zma = -zpa
zpb = zpa - 3
zmb = -zpb
print "add --:\t", (wma + wmb), "\t", (zma + zmb)
print "add -+:\t", (wma + wpb), "\t", (zma + zpb)
print "add +-:\t", (wpa + wmb), "\t", (zpa + zmb)
print "add ++:\t", (wpa + wpb), "\t", (zpa + zpb)
print "sub --:\t", (wma - wmb), "\t", (zma - zmb)
print "sub -+:\t", (wma - wpb), "\t", (zma - zpb)
print "sub +-:\t", (wpa - wmb), "\t", (zpa - zmb)
print "sub ++:\t", (wpa - wpb), "\t", (zpa - zpb)
print "sub --:\t", (wmb - wma), "\t", (zmb - zma)
print "sub -+:\t", (wmb - wpa), "\t", (zmb - zpa)
print "sub +-:\t", (wpb - wma), "\t", (zpb - zma)
print "sub ++:\t", (wpb - wpa), "\t", (zpb - zpa)
print "mul --:\t", (wma * wmb), "\t", (zma * zmb)
print "mul -+:\t", (wma * wpb), "\t", (zma * zpb)
print "mul +-:\t", (wpa * wmb), "\t", (zpa * zmb)
print "mul ++:\t", (wpa * wpb), "\t", (zpa * zpb)
print "div --:\t", (wma / wmb), "\t", (zma / zmb)
print "div -+:\t", (wma / wpb), "\t", (zma / zpb)
print "div +-:\t", (wpa / wmb), "\t", (zpa / zmb)
print "div ++:\t", (wpa / wpb), "\t", (zpa / zpb)
print "div --:\t", (wmb / wma), "\t", (zmb / zma)
print "div -+:\t", (wmb / wpa), "\t", (zmb / zpa)
print "div +-:\t", (wpb / wma), "\t", (zpb / zma)
print "div ++:\t", (wpb / wpa), "\t", (zpb / zpa)
print "pow --:\t", (wma**wmb), "\t", (zma**zmb)
print "pow -+:\t", (wma**wpb), "\t", (zma**zpb)
print "pow +-:\t", (wpa**wmb), "\t", (zpa**zmb)
print "pow ++:\t", (wpa**wpb), "\t", (zpa**zpb)
print "pow --:\t", (wmb**wma), "\t", (zmb**zma)
print "pow -+:\t", (wmb**wpa), "\t", (zmb**zpa)
print "pow +-:\t", (wpb**wma), "\t", (zpb**zma)
print "pow ++:\t", (wpb**wpa), "\t", (zpb**zpa)
# wpa = mpw(13)
# wma = -wpa
# wpb = wpa - 3
# wmb = -wpb
# zpa = long(13)
# zma = -zpa
# zpb = zpa - 3
# zmb = -zpb
print "mod --:\t", (wma % wmb), "\t", (zma % zmb)
print "mod -+:\t", (wma % wpb), "\t", (zma % zpb)
print "mod +-:\t", (wpa % wmb), "\t", (zpa % zmb)
print "mod ++:\t", (wpa % wpb), "\t", (zpa % zpb)
print "mod --:\t", (wmb % wma), "\t", (zmb % zma)
print "mod -+:\t", (wmb % wpa), "\t", (zmb % zpa)
print "mod +-:\t", (wpb % wma), "\t", (zpb % zma)
print "mod ++:\t", (wpb % wpa), "\t", (zpb % zpa)
print "rem --:\t", divmod(wma, wmb), "\t", divmod(zma, zmb)
print "rem -+:\t", divmod(wma, wpb), "\t", divmod(zma, zpb)
print "rem +-:\t", divmod(wpa, wmb), "\t", divmod(zpa, zmb)
print "rem ++:\t", divmod(wpa, wpb), "\t", divmod(zpa, zpb)
print "rem --:\t", divmod(wmb, wma), "\t", divmod(zmb, zma)
print "rem -+:\t", divmod(wmb, wpa), "\t", divmod(zmb, zpa)
print "rem +-:\t", divmod(wpb, wma), "\t", divmod(zpb, zma)
print "rem ++:\t", divmod(wpb, wpa), "\t", divmod(zpb, zpa)
pass