E = EllipticCurve([1,0,1,-1,-2]) F.<cbrt5> = NumberField(x^3-5) EF = E.change_ring(F) p = F.fractional_ideal(cbrt5) EF.reduction(p).cardinality()