Contact
CoCalc Logo Icon
StoreFeaturesDocsShareSupport News AboutSign UpSign In
| Download
Views: 2206
C = CremonaModularSymbols(100001, cuspidal=True, sign=-1) C.dimension()
9091
T = C.hecke_matrix(2) TZ = T.sage_matrix_over_ZZ()
from time import clock T2 = TZ.change_ring(GF(2)).dense_matrix() t = clock() T2.rank() clock() - t #runtime in seconds
9047 0.36098400000000197
T3 = TZ.change_ring(GF(3)).dense_matrix() t = clock() T3.rank() clock() - t
9085 49.191177
p = 2^30 - 35 Tp = TZ.change_ring(GF(p)).dense_matrix() #t = clock() #Tp.rank()#runs for days, don't attempt #print clock() - t
Tz = TZ.dense_matrix() t = clock() Tz.rank() print clock() - t
9091 81.011494