+# Simple tests of modmul.
+for ai in range(20, 200, 60):
+ a = sqrt(3<<(2*ai-1))
+ for bi in range(20, 200, 60):
+ b = sqrt(5<<(2*bi-1))
+ for m in range(20, 600, 32):
+ m = sqrt(2**(m+1))
+ print "modmul", hexstr(a), hexstr(b), hexstr(m), hexstr((a*b) % m)
+