Check whether the element in Ideal

[kirtsar]

Qx, x = QQ["X"]
Qi, a = NumberField(x^2 + 1)
Zi = maximal_order(Qi)
I = (1 + i)^2 * Zi
in(2, I)                   # not working
in(ZZ(2), I)            # not working
in(Zi(2), I)             # not working
in(Qi(2), I)            # only this happen to work

2 // (1 + i)^2         # works
Zi(2) // (1 + i)^2   # not working

[thofma]

Yes, they should all work. Thanks for reporting this!

You probably know it, but // will always return something outside the ring and you can use divexact and divides to do division in the ring.

[kirtsar]

@thofma
No, I didn't know it, thanks!

[CarloSircana]

Fixed.