gap> Q:=QuadraticNumberField(-5);;
gap> OQ:=RingOfIntegers(Q);;
gap> I:=QuadraticIdeal(OQ,7);;
gap> IsPrime(I);
false

gap> Factors(Discriminant(OQ));
[ -2, 2, 5 ]
