Singular
https://www.singular.uni-kl.de/forum/

primary decomposition in Q[i]
https://www.singular.uni-kl.de/forum/viewtopic.php?f=10&t=1390
Page 1 of 1

Author:  Alberto Damiano [ Thu Aug 11, 2005 5:32 pm ]
Post subject:  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

Author:  levandov [ Thu Aug 11, 2005 8:39 pm ]
Post subject:  Re: primary decomposition in Q[i]

Dear Alberto Damiano,
Quote:
> 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;

here you've forgotten a comma, it should be rather
Code:
ring r=(0,i),(x,y),dp;
minpoly=i^2+1;

Quote:
> 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?
>

There should be an error in you Singular script, because the correct code
Code:
ring r=(0,i),(x,y),dp;
minpoly=i^2+1;
ideal I=x2+y2;
LIB "primdec.lib";
list L=primdecGTZ(I);
L;

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):

Code:
[1]:
   [1]:
       _[1]=x+(-i)*y
   [2]:
       _[1]=x+(-i)*y
[2]:
   [1]:
       _[1]=x+(i)*y
   [2]:
       _[1]=x+(i)*y



Quote:
> Thank you


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.

Quote:
> Alberto Damiano


Best regards,

Page 1 of 1 All times are UTC + 1 hour [ DST ]
Powered by phpBB © 2000, 2002, 2005, 2007 phpBB Group
http://www.phpbb.com/