primary decomposition in Q[i]

I am interested in computing the Primary decomposition of ideals in C[x,y]. For example i tried the following ideal:

I=(x^2+y^2)

in the polinomial ring Q[i][x,y] where i^2+1=0

the definition of the ring was as follows:

ring r=(0,i)(x,y),dp;minpoly=i^2+1;

PrimdecGTZ returns incorrectly the ideal itself, as if it was primary, while all the other primdec command from the primdec.lib and mprimdec.lib return rerror messages. My suispect is that we cannot compute the primary decomposistion with an algebraic extension of Q as coefficient ring. Anyone has any idea?

Thank you

Alberto Damiano

email: adamiano@gmu.edu
Posted in old Singular Forum on: 2004-05-05 17:17:36+02

Dear Alberto Damiano,

here you've forgotten a comma, it should be rather


There should be an error in you Singular script, because the correct code

gives always a correct answer (I've checked it extra with versions 2-0-4, 2-0-5 and even with the development version 2-1-2):






If the given code produces something different from the given answer, please let us know your Singular version and also the version of primdec.lib.



Best regards,