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 __call__(self, val): xself = long(self.false_self) yself = rpm.mpw(self.false_self) xm = long.__getattribute__(xself, self.method_name) ym = rpm.mpw.__getattribute__(yself, self.method_name) xa = xm(long(val)) ya = ym(rpm.mpw(val)) print " Comparing", self.method_name, type(xa), xa, type(ya), ya assert xa == ya return xa
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 zs = oct(za) print "__oct__:\t", oct(wa), "\t", zs assert oct(wa) == zs 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
import rpm #rpm.mpw().Debug(-1); for i in range(100): a = rpm.mpw(i) for j in range(100): c = a**j print i, '^', j, '=', c
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 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
import rpm rpm.mpw().Debug(1) Methods = ( '__neg__', '__pos__', '__abs__', '__invert__', '__add__', '__sub__', '__mul__', '__div__', '__mod__', '__divmod__', '__pow__', '__lshift__', '__rshift__', '__and__', '__xor__', '__or__') class Factory(object): def __init__(self, false_self, method_name): self.false_self = false_self self.method_name = method_name def __call__(self, val): xself = long(self.false_self) yself = rpm.mpw(self.false_self) xm = long.__getattribute__(xself, self.method_name) ym = rpm.mpw.__getattribute__(yself, self.method_name) xa = xm(long(val)) ya = ym(rpm.mpw(val)) print " Comparing", self.method_name, type(xa), xa, type(ya), ya assert xa == ya return xa class Long(long): def __getattribute__(self, name): print "__getattribute__ ~%s~" % name if name not in Methods: return long.getattr(self, name) return Factory(self, name) a1 = Long(0x987654321) a2 = Long(0x10)
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 zs = oct(za) print "__oct__:\t", oct(wa), "\t", zs assert oct(wa) == zs 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 setUp(self): mpw().Debug(0) pass
def tearDown(self): mpw().Debug(0) pass
import rpm rpm.mpw().Debug(1) Methods = ('__neg__', '__pos__', '__abs__', '__invert__', '__add__', '__sub__', '__mul__', '__div__', '__mod__', '__divmod__', '__pow__', '__lshift__', '__rshift__', '__and__', '__xor__', '__or__') class Factory(object): def __init__(self, false_self, method_name): self.false_self = false_self self.method_name = method_name def __call__(self, val): xself = long(self.false_self) yself = rpm.mpw(self.false_self) xm = long.__getattribute__(xself, self.method_name) ym = rpm.mpw.__getattribute__(yself, self.method_name) xa = xm(long(val)) ya = ym(rpm.mpw(val)) print " Comparing", self.method_name, type(xa), xa, type(ya), ya assert xa == ya return xa class Long(long): def __getattribute__(self, name): print "__getattribute__ ~%s~" % name if name not in Methods: return long.getattr(self, name)
import rpm rpm.rng().Debug(1) # rpm.mpw().Debug(-1); # Generate random intergers r = rpm.rng() for i in range(100): print r.next() # Generate random numbers with lots of bits m = rpm.mpw("800000000000000000000000000000000000000000000000000000000000000000000000") for i in range(100): r.next(m) # Generate 160 bit random primes bits = 160 trials = -1 a = r.prime(bits, trials) for i in range(100): b = r.prime(bits, trials) print a.gcd(a, b) a = b
import rpm rpm.rng().Debug(1) #rpm.mpw().Debug(-1); # Generate random intergers r = rpm.rng() for i in range(100): print r.next() # Generate random numbers with lots of bits m = rpm.mpw( "800000000000000000000000000000000000000000000000000000000000000000000000") for i in range(100): r.next(m) # Generate 160 bit random primes bits = 160 trials = -1 a = r.prime(bits, trials) for i in range(100): b = r.prime(bits, trials) print a.gcd(a, b) a = b