search

jakobkroeker / test_singular

manage a personal issue list for Singular CAS
0 stars 0 forks source link

bug in primdecSY #124

Open jakobkroeker opened 9 years ago

jakobkroeker commented 9 years ago

see also http://www.singular.uni-kl.de:8002/trac/ticket/678

first example

LIB("primdec.lib");
ring rng = (7,vv),(x,y,z),dp;
minpoly = (vv^2+vv+3);
ideal I =  (2*vv+3)*y^2*z+(2*vv+1)*y^2+(vv-2), x*z^2+(-2*vv+3)*x^2, (3*vv)*x^3+(-vv) ;
def L1 = primdecSY (I); // wrong?
def L2 = primdecSY ( std(I) );
size(L1) == size(L2); //
L1;
L2;
isSubModule( L2[1][2], L2[2][2] );
isSubModule( L2[2][2], L2[1][2] );
testPrimary(L2,I);

output


> L1;
[1]
   [1]
      _[1]=z^2+(-2*vv+3)*x
      _[2]=y^2*z+(-3*vv-1)*y^2-3
      _[3]=x^2*z+(3*vv+1)*x^2+(2*vv-1)*y^2+(2*vv)*x*z+(3*vv+3)*x+(3*vv+2)*z+3
      _[4]=x*y^2+(-2*vv)*y^2+(vv-1)*z+(2*vv-3)
      _[5]=x^3+2
      _[6]=y^4+(vv-2)*x^2+(-vv+2)*y^2+(-vv)*x*z+(3*vv+3)*x+(-3*vv-2)*z-2
   [2]
      _[1]=z^2+(-2*vv+3)*x
      _[2]=y^2*z+(-3*vv-1)*y^2-3
      _[3]=x^2*z+(3*vv+1)*x^2+(2*vv-1)*y^2+(2*vv)*x*z+(3*vv+3)*x+(3*vv+2)*z+3
      _[4]=x*y^2+(-2*vv)*y^2+(vv-1)*z+(2*vv-3)
      _[5]=x^3+2
      _[6]=y^4+(vv-2)*x^2+(-vv+2)*y^2+(-vv)*x*z+(3*vv+3)*x+(-3*vv-2)*z-2
> L2;
[1]
   [1]
      _[1]=z^2+(-2*vv+3)*x
      _[2]=y^2+(2*vv-3)*x*z+(-2*vv-1)*y*z+2*x+(3*vv+3)*y+(-3*vv+2)*z+(2*vv+3)
      _[3]=x*y+y^2+3*y*z+(vv+1)*x+(2*vv+1)*y+(3*vv-1)*z+(-3*vv+3)
      _[4]=x^2+(2*vv-1)*y^2+x*z+2*y*z+(-3*vv-2)*x+(-3*vv-3)*y+(3*vv)*z-2
   [2]:
      _[1]=z^2+(-2*vv+3)*x
      _[2]=y^2+(2*vv-3)*x*z+(-2*vv-1)*y*z+2*x+(3*vv+3)*y+(-3*vv+2)*z+(2*vv+3)
      _[3]=x*y+y^2+3*y*z+(vv+1)*x+(2*vv+1)*y+(3*vv-1)*z+(-3*vv+3)
      _[4]=x^2+(2*vv-1)*y^2+x*z+2*y*z+(-3*vv-2)*x+(-3*vv-3)*y+(3*vv)*z-2
[2]:
   [1]
      _[1]=z^2+(-2*vv+3)*x
      _[2]=y^2+(2*vv-3)*x*z+(2*vv+1)*y*z+2*x+(-3*vv-3)*y+(-3*vv+2)*z+(2*vv+3)
      _[3]=x*y-y^2+3*y*z+(-vv-1)*x+(2*vv+1)*y+(-3*vv+1)*z+(3*vv-3)
      _[4]=x^2+(2*vv-1)*y^2+x*z-2*y*z+(-3*vv-2)*x+(3*vv+3)*y+(3*vv)*z-2
   [2]:
      _[1]=z^2+(-2*vv+3)*x
      _[2]=y^2+(2*vv-3)*x*z+(2*vv+1)*y*z+2*x+(-3*vv-3)*y+(-3*vv+2)*z+(2*vv+3)
      _[3]=x*y-y^2+3*y*z+(-vv-1)*x+(2*vv+1)*y+(-3*vv+1)*z+(3*vv-3)
      _[4]=x^2+(2*vv-1)*y^2+x*z-2*y*z+(-3*vv-2)*x+(3*vv+3)*y+(3*vv)*z-2
> isSubModule( L2[1][2], L2[2][2] );
0
> isSubModule( L2[2][2], L2[1][2] );
0
> testPrimary(L2,I);
1

second example:

LIB("primdec.lib");
ring rng = (27271,vv),(xx,yy,zz,ww),dp;
minpoly = (vv^3+25);
ideal I = (12841*vv^2)*xx^2*ww+(12841*vv^2),(2883*vv)*xx^2+(vv)*xx*yy^2+(-8499*vv)*xx*yy,625*ww^2+6885*ww-7544;
list L1  = primdecGTZ (I);
list L2  = primdecSY (I);
isSubModule( L1[1][2], L1[2][2] );
isSubModule( L1[2][2], L1[1][2] );
testPrimary( L1, I );
size( L1 ); 
size( L2 );

output

> ring rng = (27271,vv),(xx,yy,zz,ww),dp;
// ** redefining rng **
> minpoly = (vv^3+25);
> ideal I = (12841*vv^2)*xx^2*ww+(12841*vv^2),(2883*vv)*xx^2+(vv)*xx*yy^2+(-8499*vv)*xx*yy,625*ww^2+6885*ww-7544;
> list L1  = primdecGTZ (I);
> list L2  = primdecSY (I);
> isSubModule( L1[1][2], L1[2][2] );
0
> isSubModule( L1[2][2], L1[1][2] );
0
> testPrimary( L1, I );
1
> size( L1 ); 
2
> size( L2 );
1
> L1;
[1]:
   [1]:
      _[1]=ww^2-2607*ww+6053
      _[2]=yy*ww-2040*xx-1083*yy+9076*ww-5748
      _[3]=xx*ww-482*yy*ww-3801*ww^2+13331*yy-11066*ww
      _[4]=yy^2+2883*xx-8499*yy
      _[5]=xx*yy+5174*yy^2-2846*xx*ww+13159*yy*ww-11066*yy
      _[6]=xx^2+5174*xx*yy-4552*yy^2-12023*xx*ww+1698*yy*ww-11066*xx+1613*yy
   [2]:
      _[1]=ww^2-2607*ww+6053
      _[2]=yy*ww-2040*xx-1083*yy+9076*ww-5748
      _[3]=xx*ww-482*yy*ww-3801*ww^2+13331*yy-11066*ww
      _[4]=yy^2+2883*xx-8499*yy
      _[5]=xx*yy+5174*yy^2-2846*xx*ww+13159*yy*ww-11066*yy
      _[6]=xx^2+5174*xx*yy-4552*yy^2-12023*xx*ww+1698*yy*ww-11066*xx+1613*yy
[2]:
   [1]:
      _[1]=ww^2-2607*ww+6053
      _[2]=yy*ww+2040*xx-1083*yy+9696*ww-7433
      _[3]=xx*ww+482*yy*ww+13159*ww^2-13331*yy+1908*ww
      _[4]=yy^2+2883*xx-8499*yy
      _[5]=xx*yy-5174*yy^2+2846*xx*ww-3801*yy*ww+1908*yy
      _[6]=xx^2-5174*xx*yy-4552*yy^2+4937*xx*ww-8278*yy*ww+1908*xx-11714*yy
   [2]:
      _[1]=ww^2-2607*ww+6053
      _[2]=yy*ww+2040*xx-1083*yy+9696*ww-7433
      _[3]=xx*ww+482*yy*ww+13159*ww^2-13331*yy+1908*ww
      _[4]=yy^2+2883*xx-8499*yy
      _[5]=xx*yy-5174*yy^2+2846*xx*ww-3801*yy*ww+1908*yy
      _[6]=xx^2-5174*xx*yy-4552*yy^2+4937*xx*ww-8278*yy*ww+1908*xx-11714*yy
> L2;
[1]:
   [1]:
      _[1]=ww^2-2607*ww+6053
      _[2]=yy^2+2883*xx-8499*yy
      _[3]=xx^2+6713*ww+7191
   [2]:
      _[1]=ww^2-2607*ww+6053
      _[2]=yy^2+2883*xx-8499*yy
      _[3]=xx^2+6713*ww+7191