Changeset b90479 in git

Timestamp:

Aug 7, 2019, 5:03:38 PM (4 years ago)

Author:

Hans Schoenemann <hannes@…>

Branches:

(u'spielwiese', '828514cf6e480e4bafc26df99217bf2a1ed1ef45')

Children:

fcc6e71bf5a9978aca272c257309c1804ba136c0

Parents:

4afedc31e655bd54e5f31232d2e6044d42600468

Message:

bug: create_ring in tropicalLifting

File:

• Singular/LIB/tropical.lib (modified) (20 diffs)

Singular/LIB/tropical.lib

@@ -340 +340 @@
   if (nvars(basering)==2)
   {
-    ring BASERING = create_ring(ringlist(ALTERRING)[1], "(t,x(1))", "("+ordstr(ALTERRING)+")", "no_minpoly");
+    execute("ring BASERING=("+charstr(ALTERRING)+"),(t,x(1)),("+ordstr(ALTERRING)+");");
   }
   else
@@ -470 +470 @@
     variablen=variablen+var(j);
   }
-  ring GRUNDRING = create_ring(ringlist(basering)[1], "("+string(variablen)+",t)", "(dp("+string(nvars(basering)-1)+"),lp(1))", "no_minpoly");
+  execute("ring GRUNDRING=("+charstr(basering)+"),("+string(variablen)+",t),(dp("+string(nvars(basering)-1)+"),lp(1));");
   ideal variablen;
   for (j=1;j<=nvars(basering)-1;j++)
@@ -573 +573 @@
         // from LIFTRING are present,
         // and where also the variables of CUTDOWNRING live
-        ring REPLACEMENTRING = create_ring(ringlist(LIFTRING)[1], "("+varstr(CUTDOWNRING)+")", "dp", "no_minpoly");
+        execute("ring REPLACEMENTRING=("+charstr(LIFTRING)+"),("+varstr(CUTDOWNRING)+"),dp;");
         list repl=imap(CUTDOWNRING,repl); // get the replacement rules
                                           // from CUTDOWNRING
@@ -1069 +1069 @@
       poly mp=minpoly;
       def OLDRING=basering;
-      ring NEWRING = create_ring(0, "("+varstr(basering)+","+parstr(basering)+")", "ds");
+      execute("ring NEWRING=0,("+varstr(basering)+","+parstr(basering)+"),ds;");
       ideal I=imap(OLDRING,mp),imap(OLDRING,f);
     }
@@ -1717 +1717 @@
     {
       string polynomstring=string(f);
-      ring drawring = create_ring("(0,"+parstr(basering)+")", "("+varstr(basering)+")", "dp");
+      execute("ring drawring=(0,"+parstr(basering)+"),("+varstr(basering)+"),dp;");
       execute("poly f="+polynomstring+";");
     }
@@ -2652 +2652 @@
   }
   // we want to homogenise the ideal i ....
-  ring HOMOGRING = create_ring(ringlist(basering)[1], "(@s,"+string(variablen)+")", "dp", "no_minpoly");
+  execute("ring HOMOGRING=("+charstr(basering)+"),(@s,"+string(variablen)+"),dp;");
   ideal i=homog(std(imap(BASERING,i)),@s);
   // ... and compute a standard basis with
@@ -2667 +2667 @@
   }
   intmat O=weightVectorToOrderMatrix(whomog);
-  ring WEIGHTRING = create_ring(ringlist(basering)[1], "("+varstr(basering)+")", "(M("+string(O)+"))", "no_minpoly");
+  execute("ring WEIGHTRING=("+charstr(basering)+"),("+varstr(basering)+"),(M("+string(O)+"));");
   // map i to the new ring and compute a GB of i, then dehomogenise i,
   // so that we can be sure, that the
@@ -2674 +2674 @@
   // compute the w-initial ideal with the help of the procedure tInitialForm;
   setring BASERING;
-  ring COMPINIRING = create_ring(ringlist(basering)[1], "("+string(variablen)+")", "dp", "no_minpoly");
+  execute("ring COMPINIRING=("+charstr(basering)+"),("+string(variablen)+"),dp;");
   ideal i=imap(WEIGHTRING,i);
   ideal ini;
@@ -2711 +2711 @@
   def BASERING=basering;
   intmat O=weightVectorToOrderMatrix(w);
-  ring INITIALRING = create_ring(ringlist(BASERING)[1], "("+varstr(basering)+")", "M("+string(O)+")", "no_minpoly");
+  execute("ring INITIALRING=("+charstr(BASERING)+"),("+varstr(basering)+"),M("+string(O)+");");
   poly f=imap(BASERING,f);
   int GRAD=deg(f);
@@ -2784 +2784 @@
   def BASERING=basering;
   intmat O=weightVectorToOrderMatrix(w);
-  ring INITIALRING = create_ring(ringlist(BASERING)[1], "("+varstr(basering)+")", "M("+string(O)+")", "no_minpoly");
+  execute("ring INITIALRING=("+charstr(BASERING)+"),("+varstr(basering)+"),M("+string(O)+");");
   ideal i=imap(BASERING,i);
   i=std(i);
@@ -2834 +2834 @@
   }
   def BASERING=basering;
-  ring PARAMETERRING = create_ring("("+string(char(basering))+")", "("+parstr(basering)+")", "ds");
+  execute("ring PARAMETERRING=("+string(char(basering))+"),("+parstr(basering)+"),ds;");
   poly den=imap(BASERING,den);
   poly num=imap(BASERING,num);
@@ -3493 +3493 @@
 {
   def BASERING=basering;
-  ring RADRING = create_ring(ringlist(basering)[1], "(@T,"+varstr(basering)+")", "(dp(1),"+ordstr(basering)+")", "no_minpoly");
+  execute("ring RADRING=("+charstr(basering)+"),(@T,"+varstr(basering)+"),(dp(1),"+ordstr(basering)+");");
   ideal I=ideal(imap(BASERING,i))+ideal(1-@T*imap(BASERING,f));
   if (reduce(1,std(I))==0)
@@ -4218 +4218 @@
         // return pideal, the initial and the list ergl which tells us
         // which variables we replaced by which form
-        ring BASERINGLESS1 = create_ring(ringlist(BASERING)[1], "("+string(variablen)+",t)", "(dp("+string(ncols(variablen))+"),lp(1))", "no_minpoly");
+        execute("ring BASERINGLESS1=("+charstr(BASERING)+"),("+string(variablen)+",t),(dp("+string(ncols(variablen))+"),lp(1));");
         ideal i=imap(BASERING,pideal);
         ideal ini=imap(BASERING,pini); //ideal ini2=tInitialIdeal(i,wvecp,1);
@@ -4295 +4295 @@
   }
   setring BASERING;
-  ring BASERINGLESS2 = create_ring(ringlist(BASERING)[1], "("+string(variablen)+",t)", "(dp("+string(ncols(variablen))+"),lp(1))", "no_minpoly");
+  execute("ring BASERINGLESS2=("+charstr(BASERING)+"),("+string(variablen)+",t),(dp("+string(ncols(variablen))+"),lp(1));");
   // using the 0/1-vector which tells us which variables belong
   // to the set of smallest entries of wvec
@@ -4643 +4643 @@
   if ((numberdeletedvariables>0) and (anzahlvariablen>1)) // if only t remains,
   { // all true variables are gone
-    ring NEURING = create_ring(ringlist(basering)[1], "("+string(variablen)+")", "(dp("+string(size(variablen)-1)+"),lp(1))", "no_minpoly");
+    execute("ring NEURING=("+charstr(basering)+"),("+string(variablen)+"),(dp("+string(size(variablen)-1)+"),lp(1));");
     ideal i=imap(BASERING,i);
     ideal gesamt_m=imap(BASERING,gesamt_m);
@@ -4659 +4659 @@
     {
       // pass to a ring which has variables which are suitable for gfan
-      ring GFANRING = create_ring(ringlist(basering)[1], "(a,b,c,d,e,f,g,h,i,j,k,l,m,n,o,p,q,r,s,t,u,v,w,x,y,z,A,B,C,D,E,F,G,H,I,J,K,L,M,N,O,P,Q,R,S,T,U,V,W,X,Y,Z)", "dp", "no_minpoly");
+      execute("ring GFANRING=("+charstr(basering)+"),(a,b,c,d,e,f,g,h,i,j,k,l,m,n,o,p,q,r,s,t,u,v,w,x,y,z,A,B,C,D,E,F,G,H,I,J,K,L,M,N,O,P,Q,R,S,T,U,V,W,X,Y,Z),dp;");
       ideal phiideal=b,c,d,e,f,g,h,i,j,k,l,m,n,o,p,q,r,s,t,u,v,w,x,y,z,A,B,C,D,E,F,G,H,I,J,K,L,M,N,O,P,Q,R,S,T,U,V,W,X,Y,Z;
       phiideal[nvars(PREGFANRING)]=a; // map t to a
@@ -4829 +4829 @@
     else
     {
-      ring PARARing = create_ring(ringlist(basering)[1], "t", "ls", "no_minpoly");
+      execute("ring PARARing=("+charstr(basering)+"),t,ls;");
     }
     ideal PARA; // will contain the parametrisation
@@ -4917 +4917 @@
     variablen=variablen+var(j);
   }
-  ring INITIALRING = create_ring(ringlist(basering)[1], "("+string(variablen)+")", "dp", "no_minpoly");
+  execute("ring INITIALRING=("+charstr(basering)+"),("+string(variablen)+"),dp;");
   ideal ini=imap(BASERING,ini);
   // compute the associated primes of the initialideal
@@ -5164 +5164 @@
   else
   {
-    ring TINRING = create_ring(ringlist(basering)[1], "("+string(variablen)+")", "dp", "no_minpoly");
+    execute("ring TINRING=("+charstr(basering)+"),("+string(variablen)+"),dp;");
   }
   ideal tin=imap(BASERING,tin);
@@ -5543 +5543 @@
     if (size(#)>0) // noAbs was used
     {
-      ring NOQRing = create_ring("("+string(char(LIFTRing))+")", "("+varstr(basering)+","+parstr(LIFTRing)+")", "dp");
+      execute("ring NOQRing=("+string(char(LIFTRing))+"),("+varstr(basering)+","+parstr(LIFTRing)+"),dp;");
       execute("qring TESTRing=std("+#[1]+");");
       ideal i=imap(BASERING,i);