incomplete primdecSY() decomp example over field extension

[kroeker@…]

Here is another example from the incomplete primdecSY() decomposion hell:

LIB("primdec.lib");


proc ocassionlyFailingPrimdecSY()
{
 system("random",491308186);

 ring rng = (7,vv),(xc,xt,xp,xq),dp;
 option(noredSB);
 minpoly = (vv^2+vv+3);
 ideal I = (-2*vv+2)*xc^2*xt+(vv+1)*xc*xt*xq+(-2*vv+1)*xc*xt+(vv-3)*xc*xp,(-3*vv+3)*xc^2*xp+(3*vv)*xc^2;
 list L = primdecSY(I);
 ideal c2 = L[2][1];
 c2;
 list pdc2 = primdecGTZ(c2);
 pdc2;
 size(pdc2);
 ASSUME(0,  size(pdc2)==1);
}

int i;
for (i=1;i<10;i++)
{ 
    dbprint(1,"i = "+string(i) );
    ocassionlyFailingPrimdecSY(); 
}

output:

i = 1
c2[1]=xc^2
c2[2]=xc*xt*xq+(-vv-2)*xc*xt+(-vv+1)*xc*xp
[1]:
   [1]:
      _[1]=xc
   [2]:
      _[1]=xc
[2]:
   [1]:
      _[1]=xt*xq+(-vv-2)*xt+(-vv+1)*xp
      _[2]=xc^2
   [2]:
      _[1]=xc
      _[2]=xt*xq+(-vv-2)*xt+(-vv+1)*xp
2
   ? ASSUME failed: ASSUME(0,  size(pdc2)==1);

[hannes]

not reproducible with current version (dda219ee370f1a2dc241951de15186d1c02c0141, Tue Oct 28 12:12:29 2014)