R.<x> = QQ[]
K.<theta_K> = NumberField(x^6 - 18*x^4 - 24*x^3 + 27*x^2 + 36*x - 6)
fp = K.fractional_ideal((2, -7/44*theta_K^5 + 19/44*theta_K^4 + 87/44*theta_K^3 - 87/44*theta_K^2 - 5/2*theta_K + 39/22))
print(fp.norm()) # yields 2
print(fp in K.primes_above(2)) # yields True
v = K.valuation(fp) # raises AssertionError
--> 303 assert(len(candidates_for_I) > 0) # This should not be possible, unless I contains a unit