There is an server with 48 cpus in my office. I want to use all of them to compute the following codes in sage. How can I do that?
%time
import itertools
def is_matrix_integral(M):
for i,j in itertools.product(range(2),range(2)):
if M[i][j] not in ZZ:
return False
return True
def Gamma(M):
if is_matrix_integral(M) and M[1][0]%243==0 and M.det()==1:
return True
return False
p=4
T=Matrix([[2,-1],[9,-4]])*Matrix([[4,0],[0,9]])
W=Matrix([[0,1],[-243,0]])
S=Matrix([[1,1/3],[0,1]])
A=Matrix([[1,0],[81,1]])
B=-1/(243)WA^(-1)WA
C=Matrix([[1,1/9],[-27,-2]])
e=Matrix([[1,0],[0,-1]])
def Q(a,c):
return T*Matrix([[a,a+c],[c,-a]])*T^(-1)*e
I=[-100000..-1]+[1..100000]
for i,j in itertools.product([1..100000],I):
if sqrt(i*i+i*j+j*j) in QQ:
for k in [0,1,2]:
L=Q(i,j)*W*Q(i,j)^(-1)*W^(-1)*B^k
if Gamma(L):
print x,Q(i,j),Q(i,j).det(),i,j
for l in [0,1,2]:
M=Q(i,j)*B*Q(i,j)^(-1)*B^l
if Gamma(M):
print x,x,M,M.det(),i,j,k,l
break
