test-cm-rank.sage (523B)
1 def test(): 2 attach("cm-rank+-all.sage") 3 no = 0 4 unk = 0 5 for cf in data: 6 E = EllipticCurve(cf) 7 E_K = E.base_extend(NumberField(x^2-E.cm_discriminant(),'a')) 8 try: 9 if E_K.rank() != 2* E.rank(): 10 print "Non doubling:", cf, E.rank(), E_K.rank() 11 no += 1 12 except: 13 print "Can't compute the rank for", cf 14 unk += 1 15 print "Total ok curves:", (len(data)-(no+unk)), "/", len(data), \ 16 "( unknown:", unk, ")"