Thursday, 24 November 2011

ac.commutative algebra - Any implemented algorithm to compute the closure of an affine variety in a product of projective spaces?

This can be done in a few steps in probably any computer algebra package. You take the generators of your original ideal I, and bi-homogenize them, as described in the question. Then saturate with respect to the two hyperplanes at infinity, which are defined by the equation x0y0.



For example, the diagonal in mathbbA3timesmathbbA3 is defined by x1y1, x2y2, and x3y3. If I wanted to use this to compute the ideal of the diagonal in mathbbP3timesmathbbP3, I would use the following commands in Macaulay2:



 r = QQ[x0,x1,x2,x3,y0,y1,y2,y3]
i = ideal(x1*y0-y1*x0, x2*y0-y2*x0, x3*y0-y3*x0)
saturate(i, x0*y0)


The code in Singular would be:



 ring r = 0, (x0,x1,x2,x3,y0,y1,y2,y3), dp;
ideal i = x1*y0-y1*x0, x2*y0-y2*x0, x3*y0-y3*x0;
LIB "elim.lib";
sat(i, x0*y0);

No comments:

Post a Comment