search

adipopescu / To-Do-List

My Singular Project
0 stars 0 forks source link

examples ok now: #27

Closed jakobkroeker closed 9 years ago

jakobkroeker commented 9 years ago

The patches in SignatureSB fixes several bugs, but now I hit a new one:

//options: redefine loadLib usage prompt

ring rng = (integer),(x(1),x(2),x(3),x(4),x(5)),(dp(5),C);
ideal I  = x(4)^2-x(3)*x(5),x(3)*x(4)-x(5),x(3)^2-x(4),x(1)^2-121*x(5)^2-11*x(5)-1,2*x(3)^3-3*x(3)*x(4)+x(5),-4*x(1)*x(4)^2-2*x(1)*x(3)*x(5),2*x(1)*x(3)*x(4)+2*x(1)*x(5),2*x(1)*x(3)^2-4*x(1)*x(4),484*x(4)^2*x(5)+242*x(3)*x(5)^2+22*x(4)^2+11*x(3)*x(5),-8*x(1)*x(3)*x(4)+2*x(1)*x(5),4*x(1)*x(3)^2,-2*x(1)*x(3),968*x(3)*x(4)*x(5)+44*x(3)*x(4)-242*x(5)^2-11*x(5),-4*x(1)*x(3)^2-2*x(1)*x(4),4*x(1)*x(3),-2*x(1),484*x(3)^2*x(5)+22*x(3)^2+242*x(4)*x(5)+11*x(4);
ideal J  = std(I); // error

update: seems ok now

ring rng = (integer),(x,y),(dp(2),C);
ideal I  = xy-2x-8y,y-2,x-8;
ideal J  = std(I);

finishes now

ring rng = (integer),(x,y),(dp(2),C);
poly  a = 2x;
ideal J = 3y,2x,xy;
ideal aJJ =  quotient( a*J, J );

finishes now

LIB("primdec.lib");
ring rng = (integer),(xs,xp,xm),(lp(2),lp(1),C);
ideal I = 4*xm^2-8,14*xm^2+1;
ideal J = -5*xs-14;
ideal A = intersect(std(I),std(J)); //ok
ideal B= intersect(I,J); // takes long time or hangs

seems ok now

problem with parser?


LIB("ring.lib");

proc intersect_via_eliminate(def i, def j) 
"USAGE:   intersect_via_eliminate(i,j);
@*        i,j = ideals
ASSUME:   i,j = ideals over the integers or a field
RETURN:   ideal, the intersection of i and j
EXAMPLE:  example intersect_via_eliminate;  shows an example" 
{

  i = std(i);
  j = std(j);

  int bIsQring =   isQuotientRing( basering );
  ideal q = ideal(basering);

  def BAS = basering;

  // add one new variable "t" to a copy of basering
  int numVariables = nvars(basering);  
  int pos_t = numVariables + 1; // position of new variable 't' is numVariables + 1

  // introduce new variable called "t" 
  // hardcoded, but does not hurt, since there are no parameters. and 'x' is not in conflict with 't'
  string tempstr = "ring pr2=" + charstr(BAS) + ",(x(1.."+  string(numVariables) + "), t), dp;"; 
  execute(tempstr);   

  if ( not bIsQring )
  {
     def r2 = pr2;
  }
  else 
  {
     ideal   q = imap(BAS,q);
     qring r2 = q;
  }

  /*
  list rl = ringlist(basering);

  rl[2] =  indexedVariables("x", numVariables) + list("t");;
  dbprint(printlevel-voice, "rl[2] = ", rl[2] );
  // change variable ordering to dp:
  rl =  changeRingListVarOrderTo_dp (rl);
  dbprint(printlevel-voice, "rl[2] = ", rl[2] );

  // difference to previous variant: if basering is a qring, then r2 is also a qring.
  def r2 = ring(rl);
  setring r2;
  // */
  //dbprint(printlevel-voice, "intersect_via_eliminate:: r2 = ", basering );
  //dbprint(printlevel-voice, "intersect_via_eliminate:: var(pos_t)= ", var(pos_t) );

  // for i=(f_1,...,f_s) and j=(g_1,...,g_t) compute the ideal
  // J=(tf_1,...,tf_s,(1-t)g_1,...,(1-t)g_t)
  // then we eliminate "t" from J and map the new J to our basering
  ideal J = 0;

  int k;
  ideal ii = fetch(BAS,i);
  ii = std(ii); //not necessary due to Anne
  for(k=1; k <= ncols(ii); k++) 
  {
     J = J,var(pos_t)*ii[k];
  }
  ideal jj = fetch(BAS,j);
  jj = std(jj); //not necessary due to Anne
  for(k=1; k <= ncols(jj); k++)
  {
     J = J,(1 - var(pos_t) )*jj[k];
  }

  ideal e = eliminate(J, var(pos_t) );
  setring BAS;

  ideal result1 = fetch(r2,e); 

  for(k=1; k <= ncols(result1); k++)
  {
     if(leadcoef(result1[k]) < 0)
     {
        //ERROR("intersect_via_eliminate::Warning: negative leadcoef occured");
        result1[k] = -result1[k];
     }
  }
  return(result1);
}

//system("random",1412595701);
ring rng = (integer),(xn,xl,xf),(lp(2),lp(1),C);
//int trial: =14;
// ** redefining I **
// ** redefining J **
ideal I = -9*xl-10*xf,-10*xl^2-xf;
ideal J = 15*xl*xf-10*xl,-4*xl*xf+9*xl;
 intersect(I,J);
 intersect(std(I),std(J) );
intersect_via_eliminate(std(I),std(J));

ASSUME(0, idealsEqual(      intersect(I,J), intersect( std(I),  std(J) )      ) );
ASSUME(0, idealsEqual( intersect(I,J), intersect_via_eliminate(I,J) ) );
ideal stdI = std(I);
ideal stdJ = std(J);
ASSUME(0, idealsEqual( intersect(stdI,stdJ), intersect_via_eliminate(stdI,stdJ) ) );
 intersect(std(I),std(J));
 intersect_via_eliminate(std(I),std(J));  // segfault

ok now

ring rng = (integer),(xm,xv,xa),(lp(2),lp(1),C);
//int trial: =2;
// ** redefining I **
// ** redefining J **
ideal I = 4*xm-4,-6*xm+11;
ideal J = -2*xv+11,-15*xm*xv;
 intersect(I,J);
 intersect(std(I),std(J) ); //segfault
 intersect(I,J);
 intersect(std(I),std(J) );

seems ok now

jakobkroeker commented 9 years ago

fifths example

ring rng = (integer),(x(1),x(2),x(3),x(4),y(1),y(2),y(3),y(4),y(5)),(dp(9),C);
ideal id = x(3)-y(2),x(2)-y(3),y(3)*y(5)-x(1)-1,y(3)^2-y(2),y(2)*y(3)-x(4),y(1)*y(3)-2,242*x(4)*y(3)+2*y(1)*y(5)-2*y(5)^2+22*y(3),x(1)*y(3)-y(2)*y(5)+y(3),y(2)^2-x(4)*y(3),y(1)*y(2)-2*y(3),121*x(4)*y(2)-x(1)*y(5)+11*y(2)+y(5),x(1)*y(2)-x(4)*y(5)+y(2),x(4)*y(1)-2*y(2),x(1)*y(1)+y(1)-2*y(5),x(1)^2-121*x(4)^2-11*x(4)-1,y(2)*y(5)^2-121*x(4)^2-2*x(1)-11*x(4)-2,2*y(1)*y(5)^2-2*y(5)^3+242*x(1)*x(4)+22*x(1)+242*x(4)+22,y(1)^2*y(5)-y(1)*y(5)^2+242*x(4)+22,121*x(4)^2*y(3)-x(4)*y(5)^2+11*x(4)*y(3)+2*y(2)*y(5),y(1)*y(5)^3-y(5)^4+121*x(1)*x(4)*y(5)+11*x(1)*y(5)+121*x(4)*y(5)+11*y(5),x(4)*y(5)^3-121*x(1)*x(4)^2-11*x(1)*x(4)-363*x(4)^2-4*x(1)-33*x(4)-4,x(1)*y(5)^3-14641*x(4)^3-y(5)^3-242*x(1)*x(4)-2662*x(4)^2-22*x(1)-363*x(4)-22;
ideal j = 2,y(3),y(2),x(4),x(3),x(2),x(1)+1,y(1)^2*y(5)-y(1)*y(5)^2,y(1)*y(5)^3-y(5)^4;
ideal i=quotient(id,j);

also ok now (did not check correctness)

ring rng =  (integer),(x,y,z),dp;
int trial= 3;
ideal I = 15x2-5,-15xy+5;
ideal J = -12y2-13yz,4x+7y;
intersect(I,J); //crash vs does not finish in debug mode
ring rng = (integer),(x,y,z),(dp(2),dp(1),C);
int trial= 11;
ideal I = 8y2z+10x+11,7xyz-4;
ideal J = -7x-10,-x2z+14xz-11x;
intersect(I,J); //crash vs does not finish in debug mode

also ok

here are more examples:

ring rng = (integer),(y(1..4),x(1..3)),dp;
ideal I  = 12*y(3)*y(4)+19*y(1)+6*y(4),2*x(3)*y(4)-2*y(1)+7,2*y(1)*y(3)+y(1)-7*y(3)+6,7*x(3)*y(3)-6*x(3)-6*y(3)-3,y(1)^2+126*y(3)*y(4)+196*y(1)+60*y(4),x(3)*y(1)-3,x(3)*y(3)*y(4)+30*y(3)*y(4)+48*y(1)+15*y(4)+3,48*y(3)*y(4)^2-76*y(1)*y(4)+24*y(4)^2+266*y(4),48*y(3)^2*y(4)+48*y(3)*y(4)+12*y(4),48*y(1)*y(3)*y(4)+24*y(1)*y(4)-168*y(3)*y(4)-84*y(4),48*x(3)*y(3)*y(4)-76*x(3)*y(1)-48*y(1)*y(3)+24*x(3)*y(4)+266*x(3)-24*y(1)+168*y(3)+84,-266*x(3)*y(4)+168*y(3)*y(4)+84*y(4),266*x(3)*y(3)+168*y(3)^2-228*x(3)-60*y(3)-72,336*x(3)*y(3)*y(4)-60*x(3)*y(4)-144*y(3)*y(4)-72*y(4),-266*x(3)^2-168*x(3)*y(3)+144*x(3)+144*y(3)+72,-168*y(3)^2*y(4)+266*y(1)*y(3)+60*y(3)*y(4)-228*y(1)-931*y(3)+72*y(4)+798,-168*y(3)^3-24*y(3)^2+102*y(3)+36,-168*y(1)*y(3)^2+60*y(1)*y(3)+588*y(3)^2+72*y(1)-210*y(3)-252,168*x(3)*y(3)^2+168*x(3)*y(3)-144*y(3)^2+42*x(3)-144*y(3)-36,28*x(3)*y(3)*y(4)-28*y(1)*y(3)+14*x(3)*y(4)-24*y(3)*y(4)+24*y(1)+98*y(3)-12*y(4)-84,-28*x(3)*y(3)^2+10*x(3)*y(3)+12*x(3),-28*x(3)*y(1)*y(3)+24*x(3)*y(1)+98*x(3)*y(3)-84*x(3),28*x(3)^2*y(3)+14*x(3)^2-24*x(3)*y(3)-12*x(3),48*y(1)*y(4)^2-84*y(4)^2,48*y(1)*y(3)*y(4)+24*y(1)*y(4)-84*y(3)*y(4)+72*y(4),72*y(4)^2,48*x(3)*y(1)*y(4)-84*x(3)*y(4)-72*y(4),-48*y(1)^2*y(3)-24*y(1)^2+252*y(1)*y(3)+72*y(3)*y(4)+12*y(1)-294*y(3)+36*y(4)+252,-8*y(1)^2*y(4)+504*y(3)*y(4)^2-756*y(1)*y(4)+252*y(4)^2+2744*y(4),504*y(3)^2*y(4)+492*y(3)*y(4)+120*y(4),504*y(1)*y(3)*y(4)+240*y(1)*y(4)-1764*y(3)*y(4)-840*y(4),-8*x(3)*y(1)^2+504*x(3)*y(3)*y(4)-756*x(3)*y(1)-504*y(1)*y(3)+252*x(3)*y(4)+2744*x(3)-240*y(1)+1764*y(3)+840,168*y(1)*y(3)*y(4)-144*y(1)*y(4)-294*y(3)*y(4)+252*y(4),168*y(1)*y(3)^2-60*y(1)*y(3)-294*y(3)^2-72*y(1)+504*y(3)-216,252*y(3)*y(4)-216*y(4),-168*x(3)*y(1)*y(3)-84*x(3)*y(1)+294*x(3)*y(3)+144*y(1)*y(3)-252*x(3)+72*y(1)-252*y(3)+216,-28*x(3)*y(1)*y(4)-2744*x(3)*y(4)+24*y(1)*y(4)+1764*y(3)*y(4)+840*y(4),28*x(3)*y(1)*y(3)-24*x(3)*y(1)+2744*x(3)*y(3)+1764*y(3)^2-2352*x(3)-672*y(3)-720,3528*x(3)*y(3)*y(4)-672*x(3)*y(4)-1512*y(3)*y(4)-720*y(4),-28*x(3)^2*y(1)-2744*x(3)^2+24*x(3)*y(1)-1764*x(3)*y(3)+1512*x(3)+1512*y(3)+720,28*y(1)^2*y(3)-1764*y(3)^2*y(4)-24*y(1)^2+2646*y(1)*y(3)+630*y(3)*y(4)-2268*y(1)-9604*y(3)+756*y(4)+8232,-1764*y(3)^3-210*y(3)^2+1056*y(3)+360,-1764*y(1)*y(3)^2+672*y(1)*y(3)+6174*y(3)^2+720*y(1)-2352*y(3)-2520,1764*x(3)*y(3)^2+1722*x(3)*y(3)-1512*y(3)^2+420*x(3)-1476*y(3)-360,24*x(3)*y(4)^2+24*y(1)*y(4),24*x(3)*y(3)*y(4)+38*x(3)*y(1)+24*y(1)*y(3)+12*x(3)*y(4)+12*y(1),24*x(3)*y(1)*y(4),24*x(3)^2*y(4),-24*y(1)*y(3)*y(4)+38*y(1)^2-12*y(1)*y(4)-133*y(1),-24*y(1)*y(3)^2-24*y(1)*y(3)-6*y(1),-24*y(1)^2*y(3)-12*y(1)^2+84*y(1)*y(3)+42*y(1),-24*x(3)*y(1)*y(3)-12*x(3)*y(1)+84*x(3)*y(3)+42*x(3),-4*x(3)*y(1)*y(4)-4*y(1)^2+14*x(3)*y(4)+14*y(1),-4*x(3)*y(1)*y(3)-2*x(3)*y(1),-4*x(3)*y(1)^2+14*x(3)*y(1),-4*x(3)^2*y(1)+14*x(3)^2,84*x(3)*y(3)*y(4)+133*x(3)*y(1)-72*x(3)*y(4)-114*y(1),84*x(3)*y(3)^2-30*x(3)*y(3)-36*x(3),-84*x(3)*y(1)*y(3)-42*x(3)*y(1)+72*y(1)*y(3)+36*y(1),-84*x(3)^2*y(3)-42*x(3)^2+72*x(3)*y(3)+36*x(3),-14*x(3)^2*y(4)-14*x(3)*y(1)+12*x(3)*y(4)+12*y(1),14*x(3)^2*y(3)-12*x(3)^2,-14*x(3)^2*y(1)+12*x(3)*y(1),-14*x(3)^3+12*x(3)^2,28*x(3)*y(1)*y(3)-5*x(3)*y(1)-49*x(3)*y(3)-12*y(1)*y(3)+42*x(3)-6*y(1),-24*y(1)^2*y(4)+42*y(1)*y(4),-24*y(1)^2*y(3)-12*y(1)^2+42*y(1)*y(3)-36*y(1),-36*y(1)*y(4),-36*x(3)*y(4),-252*x(3)*y(4)^2-252*y(1)*y(4),-4*x(3)*y(1)^2-252*x(3)*y(3)*y(4)-392*x(3)*y(1)-252*y(1)*y(3)-120*x(3)*y(4)-120*y(1),-252*x(3)*y(1)*y(4),-252*x(3)^2*y(4),-4*y(1)^3+252*y(1)*y(3)*y(4)-378*y(1)^2+126*y(1)*y(4)+1372*y(1),252*y(1)*y(3)^2+246*y(1)*y(3)+60*y(1),252*y(1)^2*y(3)+120*y(1)^2-882*y(1)*y(3)-420*y(1),252*x(3)*y(1)*y(3)+120*x(3)*y(1)-882*x(3)*y(3)-420*x(3),-14*x(3)*y(1)^2-882*x(3)*y(3)*y(4)-1372*x(3)*y(1)+12*y(1)^2+756*x(3)*y(4)+1176*y(1),-882*x(3)*y(3)^2+336*x(3)*y(3)+360*x(3),882*x(3)*y(1)*y(3)+420*x(3)*y(1)-756*y(1)*y(3)-360*y(1),882*x(3)^2*y(3)+420*x(3)^2-756*x(3)*y(3)-360*x(3),-38*x(3)*y(4)^2+24*y(3)*y(4)^2+12*y(4)^2,24*y(3)^2*y(4)+24*y(3)*y(4)+6*y(4),24*x(3)*y(3)*y(4)^2+12*x(3)*y(4)^2,-38*x(3)^2*y(4),-24*y(3)^2*y(4)^2+38*y(1)*y(3)*y(4)-12*y(3)*y(4)^2-133*y(3)*y(4),-24*y(3)^3*y(4)-24*y(3)^2*y(4)-6*y(3)*y(4),-24*y(1)*y(3)^2*y(4)-12*y(1)*y(3)*y(4)+84*y(3)^2*y(4)+42*y(3)*y(4),38*x(3)*y(1)*y(3)+12*x(3)*y(3)*y(4)-133*x(3)*y(3)-12*y(1)*y(3)+6*x(3)*y(4)-6*y(1)+42*y(3)+21,4*x(3)*y(3)*y(4)^2-4*y(1)*y(3)*y(4)+2*x(3)*y(4)^2+120*y(3)*y(4)^2-192*y(1)*y(4)+14*y(3)*y(4)+60*y(4)^2+672*y(4),120*y(3)^2*y(4)+120*y(3)*y(4)+30*y(4),120*y(1)*y(3)*y(4)+60*y(1)*y(4)-420*y(3)*y(4)-210*y(4),4*x(3)^2*y(3)*y(4)-4*x(3)*y(1)*y(3)+2*x(3)^2*y(4)+120*x(3)*y(3)*y(4)-192*x(3)*y(1)+14*x(3)*y(3)-120*y(1)*y(3)+60*x(3)*y(4)+672*x(3)-60*y(1)+420*y(3)+210,114*x(3)*y(4)-72*y(3)*y(4)-36*y(4),-133*x(3)*y(3)^2+114*x(3)*y(3)+42*y(3)^2-15*y(3)-18,-84*x(3)*y(3)^2*y(4)+72*y(3)^2*y(4)-36*x(3)*y(4)+36*y(3)*y(4),133*x(3)^2*y(3)-156*x(3)*y(3)-21*x(3)+36*y(3)+18,-684*x(3)*y(4)+432*y(3)*y(4)+216*y(4),14*x(3)*y(3)^2+660*x(3)*y(3)+420*y(3)^2-576*x(3)-150*y(3)-180,14*x(3)^2*y(3)*y(4)-12*x(3)^2*y(4)+840*x(3)*y(3)*y(4)-150*x(3)*y(4)-360*y(3)*y(4)-180*y(4),-14*x(3)^2*y(3)-672*x(3)^2-408*x(3)*y(3)+366*x(3)+360*y(3)+180,12*x(3)*y(3)*y(4)-432*y(3)^2*y(4)+672*y(1)*y(3)+6*x(3)*y(4)+144*y(3)*y(4)-576*y(1)-2352*y(3)+180*y(4)+2016,-14*x(3)*y(3)^3+5*x(3)*y(3)^2-420*y(3)^3+6*x(3)*y(3)-60*y(3)^2+255*y(3)+90,-14*x(3)*y(1)*y(3)^2+12*x(3)*y(1)*y(3)+49*x(3)*y(3)^2-420*y(1)*y(3)^2-42*x(3)*y(3)+150*y(1)*y(3)+1470*y(3)^2+180*y(1)-525*y(3)-630,14*x(3)^2*y(3)^2+7*x(3)^2*y(3)+408*x(3)*y(3)^2+414*x(3)*y(3)-360*y(3)^2+105*x(3)-360*y(3)-90,-24*y(1)*y(3)*y(4)^2+42*y(3)*y(4)^2,-24*y(1)*y(3)^2*y(4)-12*y(1)*y(3)*y(4)+42*y(3)^2*y(4)-36*y(3)*y(4),-36*y(3)*y(4)^2,12*x(3)*y(1)*y(4)+36*x(3)*y(4)-18*y(4),4*x(3)*y(1)*y(4)^2+392*x(3)*y(4)^2+120*y(1)*y(4)^2-252*y(3)*y(4)^2-336*y(4)^2,120*y(1)*y(3)*y(4)-252*y(3)^2*y(4)+60*y(1)*y(4)-456*y(3)*y(4)+120*y(4),-252*x(3)*y(3)*y(4)^2-120*x(3)*y(4)^2+180*y(4)^2,4*x(3)^2*y(1)*y(4)+392*x(3)^2*y(4)+120*x(3)*y(1)*y(4)-216*x(3)*y(4)-180*y(4),-4*y(1)^2*y(3)*y(4)+252*y(3)^2*y(4)^2-378*y(1)*y(3)*y(4)+126*y(3)*y(4)^2+1372*y(3)*y(4),252*y(3)^3*y(4)+246*y(3)^2*y(4)+60*y(3)*y(4),252*y(1)*y(3)^2*y(4)+120*y(1)*y(3)*y(4)-882*y(3)^2*y(4)-420*y(3)*y(4),-4*x(3)*y(1)^2*y(3)-378*x(3)*y(1)*y(3)-120*y(1)^2*y(3)-120*x(3)*y(3)*y(4)-60*y(1)^2+1372*x(3)*y(3)+756*y(1)*y(3)-60*x(3)*y(4)+180*y(3)*y(4)+90*y(1)-1176*y(3)+90*y(4)+420,-12*x(3)*y(1)*y(4)+432*y(1)*y(3)*y(4)-1176*x(3)*y(4)-360*y(1)*y(4)+1008*y(4),14*x(3)*y(1)*y(3)^2-12*x(3)*y(1)*y(3)+1372*x(3)*y(3)^2+420*y(1)*y(3)^2-1176*x(3)*y(3)-150*y(1)*y(3)-1176*y(3)^2-180*y(1)+1428*y(3)-360,882*x(3)*y(3)^2*y(4)-756*y(3)^2*y(4)+360*x(3)*y(4)+270*y(3)*y(4)-540*y(4),-14*x(3)^2*y(1)*y(3)-1372*x(3)^2*y(3)-408*x(3)*y(1)*y(3)-210*x(3)*y(1)+2352*x(3)*y(3)+360*y(1)*y(3)-420*x(3)+180*y(1)-1008*y(3)+360,-12*x(3)*y(3)*y(4)^2-19*x(3)*y(1)*y(4)+6*y(1)*y(4),-12*x(3)*y(3)^2*y(4)-19*x(3)*y(1)*y(3)-6*x(3)*y(3)*y(4)+6*y(1)*y(3)+3*y(1),6*x(3)*y(1)*y(4),6*x(3)^2*y(4),2*x(3)^2*y(4)^2+2*x(3)*y(1)*y(4)+60*x(3)*y(4)^2+60*y(1)*y(4),2*x(3)*y(1)*y(3)+60*x(3)*y(3)*y(4)+96*x(3)*y(1)+60*y(1)*y(3)+30*x(3)*y(4)+30*y(1),2*x(3)^2*y(1)*y(4)+60*x(3)*y(1)*y(4),2*x(3)^3*y(4)+60*x(3)^2*y(4),-4*x(3)*y(1)*y(3)*y(4)-x(3)*y(1)*y(4)+7*x(3)*y(3)*y(4)-60*y(1)*y(3)*y(4)+96*y(1)^2-30*y(1)*y(4)-336*y(1),-2*x(3)*y(1)*y(3)^2-x(3)*y(1)*y(3)-60*y(1)*y(3)^2-60*y(1)*y(3)-15*y(1),-2*x(3)*y(1)^2*y(3)+7*x(3)*y(1)*y(3)-60*y(1)^2*y(3)-30*y(1)^2+210*y(1)*y(3)+105*y(1),-2*x(3)^2*y(1)*y(3)+7*x(3)^2*y(3)-60*x(3)*y(1)*y(3)-30*x(3)*y(1)+210*x(3)*y(3)+105*x(3),-6*x(3)^2*y(4)+216*x(3)*y(3)*y(4)+336*x(3)*y(1)-180*x(3)*y(4)-288*y(1),7*x(3)^2*y(3)^2-6*x(3)^2*y(3)+210*x(3)*y(3)^2-75*x(3)*y(3)-90*x(3),-7*x(3)^2*y(1)*y(3)-204*x(3)*y(1)*y(3)-105*x(3)*y(1)+180*y(1)*y(3)+90*y(1),-7*x(3)^3*y(3)-204*x(3)^2*y(3)-105*x(3)^2+180*x(3)*y(3)+90*x(3),-2*x(3)*y(1)^2*y(4)-126*x(3)*y(3)*y(4)^2-196*x(3)*y(1)*y(4)-60*y(1)^2*y(4)+168*y(1)*y(4),-2*x(3)*y(1)^2*y(3)-126*x(3)*y(3)^2*y(4)-196*x(3)*y(1)*y(3)-60*y(1)^2*y(3)-60*x(3)*y(3)*y(4)-30*y(1)^2+168*y(1)*y(3)-60*y(1),60*x(3)*y(1)*y(4)-90*y(1)*y(4),60*x(3)^2*y(4)-90*x(3)*y(4);
ideal J  = std(I);

rng = ZZ[y1,y2,y3,y4,x1,x2,x3]
I = ideal(12*y3*y4+19*y1+6*y4,2*x3*y4-2*y1+7,2*y1*y3+y1-7*y3+6,7*x3*y3-6*x3-6*y3-3,y1^2+126*y3*y4+196*y1+60*y4,x3*y1-3,x3*y3*y4+30*y3*y4+48*y1+15*y4+3,48*y3*y4^2-76*y1*y4+24*y4^2+266*y4,48*y3^2*y4+48*y3*y4+12*y4,48*y1*y3*y4+24*y1*y4-168*y3*y4-84*y4,48*x3*y3*y4-76*x3*y1-48*y1*y3+24*x3*y4+266*x3-24*y1+168*y3+84,-266*x3*y4+168*y3*y4+84*y4,266*x3*y3+168*y3^2-228*x3-60*y3-72,336*x3*y3*y4-60*x3*y4-144*y3*y4-72*y4,-266*x3^2-168*x3*y3+144*x3+144*y3+72,-168*y3^2*y4+266*y1*y3+60*y3*y4-228*y1-931*y3+72*y4+798,-168*y3^3-24*y3^2+102*y3+36,-168*y1*y3^2+60*y1*y3+588*y3^2+72*y1-210*y3-252,168*x3*y3^2+168*x3*y3-144*y3^2+42*x3-144*y3-36,28*x3*y3*y4-28*y1*y3+14*x3*y4-24*y3*y4+24*y1+98*y3-12*y4-84,-28*x3*y3^2+10*x3*y3+12*x3,-28*x3*y1*y3+24*x3*y1+98*x3*y3-84*x3,28*x3^2*y3+14*x3^2-24*x3*y3-12*x3,48*y1*y4^2-84*y4^2,48*y1*y3*y4+24*y1*y4-84*y3*y4+72*y4,72*y4^2,48*x3*y1*y4-84*x3*y4-72*y4,-48*y1^2*y3-24*y1^2+252*y1*y3+72*y3*y4+12*y1-294*y3+36*y4+252,-8*y1^2*y4+504*y3*y4^2-756*y1*y4+252*y4^2+2744*y4,504*y3^2*y4+492*y3*y4+120*y4,504*y1*y3*y4+240*y1*y4-1764*y3*y4-840*y4,-8*x3*y1^2+504*x3*y3*y4-756*x3*y1-504*y1*y3+252*x3*y4+2744*x3-240*y1+1764*y3+840,168*y1*y3*y4-144*y1*y4-294*y3*y4+252*y4,168*y1*y3^2-60*y1*y3-294*y3^2-72*y1+504*y3-216,252*y3*y4-216*y4,-168*x3*y1*y3-84*x3*y1+294*x3*y3+144*y1*y3-252*x3+72*y1-252*y3+216,-28*x3*y1*y4-2744*x3*y4+24*y1*y4+1764*y3*y4+840*y4,28*x3*y1*y3-24*x3*y1+2744*x3*y3+1764*y3^2-2352*x3-672*y3-720,3528*x3*y3*y4-672*x3*y4-1512*y3*y4-720*y4,-28*x3^2*y1-2744*x3^2+24*x3*y1-1764*x3*y3+1512*x3+1512*y3+720,28*y1^2*y3-1764*y3^2*y4-24*y1^2+2646*y1*y3+630*y3*y4-2268*y1-9604*y3+756*y4+8232,-1764*y3^3-210*y3^2+1056*y3+360,-1764*y1*y3^2+672*y1*y3+6174*y3^2+720*y1-2352*y3-2520,1764*x3*y3^2+1722*x3*y3-1512*y3^2+420*x3-1476*y3-360,24*x3*y4^2+24*y1*y4,24*x3*y3*y4+38*x3*y1+24*y1*y3+12*x3*y4+12*y1,24*x3*y1*y4,24*x3^2*y4,-24*y1*y3*y4+38*y1^2-12*y1*y4-133*y1,-24*y1*y3^2-24*y1*y3-6*y1,-24*y1^2*y3-12*y1^2+84*y1*y3+42*y1,-24*x3*y1*y3-12*x3*y1+84*x3*y3+42*x3,-4*x3*y1*y4-4*y1^2+14*x3*y4+14*y1,-4*x3*y1*y3-2*x3*y1,-4*x3*y1^2+14*x3*y1,-4*x3^2*y1+14*x3^2,84*x3*y3*y4+133*x3*y1-72*x3*y4-114*y1,84*x3*y3^2-30*x3*y3-36*x3,-84*x3*y1*y3-42*x3*y1+72*y1*y3+36*y1,-84*x3^2*y3-42*x3^2+72*x3*y3+36*x3,-14*x3^2*y4-14*x3*y1+12*x3*y4+12*y1,14*x3^2*y3-12*x3^2,-14*x3^2*y1+12*x3*y1,-14*x3^3+12*x3^2,28*x3*y1*y3-5*x3*y1-49*x3*y3-12*y1*y3+42*x3-6*y1,-24*y1^2*y4+42*y1*y4,-24*y1^2*y3-12*y1^2+42*y1*y3-36*y1,-36*y1*y4,-36*x3*y4,-252*x3*y4^2-252*y1*y4,-4*x3*y1^2-252*x3*y3*y4-392*x3*y1-252*y1*y3-120*x3*y4-120*y1,-252*x3*y1*y4,-252*x3^2*y4,-4*y1^3+252*y1*y3*y4-378*y1^2+126*y1*y4+1372*y1,252*y1*y3^2+246*y1*y3+60*y1,252*y1^2*y3+120*y1^2-882*y1*y3-420*y1,252*x3*y1*y3+120*x3*y1-882*x3*y3-420*x3,-14*x3*y1^2-882*x3*y3*y4-1372*x3*y1+12*y1^2+756*x3*y4+1176*y1,-882*x3*y3^2+336*x3*y3+360*x3,882*x3*y1*y3+420*x3*y1-756*y1*y3-360*y1,882*x3^2*y3+420*x3^2-756*x3*y3-360*x3,-38*x3*y4^2+24*y3*y4^2+12*y4^2,24*y3^2*y4+24*y3*y4+6*y4,24*x3*y3*y4^2+12*x3*y4^2,-38*x3^2*y4,-24*y3^2*y4^2+38*y1*y3*y4-12*y3*y4^2-133*y3*y4,-24*y3^3*y4-24*y3^2*y4-6*y3*y4,-24*y1*y3^2*y4-12*y1*y3*y4+84*y3^2*y4+42*y3*y4,38*x3*y1*y3+12*x3*y3*y4-133*x3*y3-12*y1*y3+6*x3*y4-6*y1+42*y3+21,4*x3*y3*y4^2-4*y1*y3*y4+2*x3*y4^2+120*y3*y4^2-192*y1*y4+14*y3*y4+60*y4^2+672*y4,120*y3^2*y4+120*y3*y4+30*y4,120*y1*y3*y4+60*y1*y4-420*y3*y4-210*y4,4*x3^2*y3*y4-4*x3*y1*y3+2*x3^2*y4+120*x3*y3*y4-192*x3*y1+14*x3*y3-120*y1*y3+60*x3*y4+672*x3-60*y1+420*y3+210,114*x3*y4-72*y3*y4-36*y4,-133*x3*y3^2+114*x3*y3+42*y3^2-15*y3-18,-84*x3*y3^2*y4+72*y3^2*y4-36*x3*y4+36*y3*y4,133*x3^2*y3-156*x3*y3-21*x3+36*y3+18,-684*x3*y4+432*y3*y4+216*y4,14*x3*y3^2+660*x3*y3+420*y3^2-576*x3-150*y3-180,14*x3^2*y3*y4-12*x3^2*y4+840*x3*y3*y4-150*x3*y4-360*y3*y4-180*y4,-14*x3^2*y3-672*x3^2-408*x3*y3+366*x3+360*y3+180,12*x3*y3*y4-432*y3^2*y4+672*y1*y3+6*x3*y4+144*y3*y4-576*y1-2352*y3+180*y4+2016,-14*x3*y3^3+5*x3*y3^2-420*y3^3+6*x3*y3-60*y3^2+255*y3+90,-14*x3*y1*y3^2+12*x3*y1*y3+49*x3*y3^2-420*y1*y3^2-42*x3*y3+150*y1*y3+1470*y3^2+180*y1-525*y3-630,14*x3^2*y3^2+7*x3^2*y3+408*x3*y3^2+414*x3*y3-360*y3^2+105*x3-360*y3-90,-24*y1*y3*y4^2+42*y3*y4^2,-24*y1*y3^2*y4-12*y1*y3*y4+42*y3^2*y4-36*y3*y4,-36*y3*y4^2,12*x3*y1*y4+36*x3*y4-18*y4,4*x3*y1*y4^2+392*x3*y4^2+120*y1*y4^2-252*y3*y4^2-336*y4^2,120*y1*y3*y4-252*y3^2*y4+60*y1*y4-456*y3*y4+120*y4,-252*x3*y3*y4^2-120*x3*y4^2+180*y4^2,4*x3^2*y1*y4+392*x3^2*y4+120*x3*y1*y4-216*x3*y4-180*y4,-4*y1^2*y3*y4+252*y3^2*y4^2-378*y1*y3*y4+126*y3*y4^2+1372*y3*y4,252*y3^3*y4+246*y3^2*y4+60*y3*y4,252*y1*y3^2*y4+120*y1*y3*y4-882*y3^2*y4-420*y3*y4,-4*x3*y1^2*y3-378*x3*y1*y3-120*y1^2*y3-120*x3*y3*y4-60*y1^2+1372*x3*y3+756*y1*y3-60*x3*y4+180*y3*y4+90*y1-1176*y3+90*y4+420,-12*x3*y1*y4+432*y1*y3*y4-1176*x3*y4-360*y1*y4+1008*y4,14*x3*y1*y3^2-12*x3*y1*y3+1372*x3*y3^2+420*y1*y3^2-1176*x3*y3-150*y1*y3-1176*y3^2-180*y1+1428*y3-360,882*x3*y3^2*y4-756*y3^2*y4+360*x3*y4+270*y3*y4-540*y4,-14*x3^2*y1*y3-1372*x3^2*y3-408*x3*y1*y3-210*x3*y1+2352*x3*y3+360*y1*y3-420*x3+180*y1-1008*y3+360,-12*x3*y3*y4^2-19*x3*y1*y4+6*y1*y4,-12*x3*y3^2*y4-19*x3*y1*y3-6*x3*y3*y4+6*y1*y3+3*y1,6*x3*y1*y4,6*x3^2*y4,2*x3^2*y4^2+2*x3*y1*y4+60*x3*y4^2+60*y1*y4,2*x3*y1*y3+60*x3*y3*y4+96*x3*y1+60*y1*y3+30*x3*y4+30*y1,2*x3^2*y1*y4+60*x3*y1*y4,2*x3^3*y4+60*x3^2*y4,-4*x3*y1*y3*y4-x3*y1*y4+7*x3*y3*y4-60*y1*y3*y4+96*y1^2-30*y1*y4-336*y1,-2*x3*y1*y3^2-x3*y1*y3-60*y1*y3^2-60*y1*y3-15*y1,-2*x3*y1^2*y3+7*x3*y1*y3-60*y1^2*y3-30*y1^2+210*y1*y3+105*y1,-2*x3^2*y1*y3+7*x3^2*y3-60*x3*y1*y3-30*x3*y1+210*x3*y3+105*x3,-6*x3^2*y4+216*x3*y3*y4+336*x3*y1-180*x3*y4-288*y1,7*x3^2*y3^2-6*x3^2*y3+210*x3*y3^2-75*x3*y3-90*x3,-7*x3^2*y1*y3-204*x3*y1*y3-105*x3*y1+180*y1*y3+90*y1,-7*x3^3*y3-204*x3^2*y3-105*x3^2+180*x3*y3+90*x3,-2*x3*y1^2*y4-126*x3*y3*y4^2-196*x3*y1*y4-60*y1^2*y4+168*y1*y4,-2*x3*y1^2*y3-126*x3*y3^2*y4-196*x3*y1*y3-60*y1^2*y3-60*x3*y3*y4-30*y1^2+168*y1*y3-60*y1,60*x3*y1*y4-90*y1*y4,60*x3^2*y4-90*x3*y4);
)
J = groebnerBasis I
adipopescu commented 9 years ago

The most examples that don't work now are with the ws order (the ones from Yue)

jakobkroeker commented 9 years ago

Am 27.11.2014 11:22, schrieb Adi Popescu:

The most examples that don't work now are with the ws order (the ones from Yue)

— Reply to this email directly or view it on GitHub https://github.com/adipopescu/To-Do-List/issues/27#issuecomment-64771464.

I just moved working examples from

https://github.com/adipopescu/To-Do-List/issues/20

to

https://github.com/adipopescu/To-Do-List/issues/27

in order to shorten the issue list in

https://github.com/adipopescu/To-Do-List/issues/20