misc_code/ugly_computations/cm/test-cm-rank.sage


1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16

def test():
    attach("cm-rank+-all.sage")
    no = 0
    unk = 0
    for cf in data:
        E = EllipticCurve(cf)
        E_K = E.base_extend(NumberField(x^2-E.cm_discriminant(),'a'))
        try:
            if E_K.rank() != 2* E.rank():
                print "Non doubling:", cf, E.rank(), E_K.rank()
                no += 1
        except:
            print "Can't compute the rank for", cf
            unk += 1
    print "Total ok curves:", (len(data)-(no+unk)), "/", len(data), \
            "( unknown:", unk, ")"