Changeset 52a01d in git for Singular/LIB/normal.lib

Timestamp:

Jan 28, 2015, 11:25:49 AM (9 years ago)

Author:

Hans Schoenemann <hannes@…>

Branches:

(u'spielwiese', 'a719bcf0b8dbc648b128303a49777a094b57592c')

Children:

6bd53ba88f906b9abf68e6983aac4cf495d7631c

Parents:

d8f9c942ee4bfbe1e3f486950f095b373d881a281386125efea2b96cc1b1eb9c36caa5004c0f47f6

Message:

Merge pull request #689 from surface-smoothers/fix.libcddgmp.check

use g++ instead of gcc for libcddgmp check (for ubuntu)
was somehow lost...

File:

• Singular/LIB/normal.lib (modified) (6 diffs)

Singular/LIB/normal.lib

@@ -2286 +2286 @@
    if(size(N) < nvars(R0))
    {
-     string newvar=string(N);
-     execute("ring R=("+charstr(R0)+"),("+newvar+"),dp;");
+     list newvar=N[1..ncols(N)];
+     list @rl0=ringlist(R0);
+     @rl0[2]=list(N[1..ncols(N)]);
+     @rl0[3]=list(list("dp",1:ncols(N)),list("C",0));
+     def R=ring(@rl0); kill @rl0;
+     setring R;
      ideal I =imap(R0,I);
-     if((ringlist(R0)[3][1][1]=="dp")&&(size(ringlist(R0)[3])==2))
+     list @rl0=ringlist(R0);
+     if((@rl0[3][1][1]=="dp")&&(size(@rl0[3])==2))
      {
        attrib(I,"isSB",1);
@@ -2297 +2302 @@
        I=std(I);
      }
+     kill @rl0;
    }
    else
@@ -2307 +2313 @@
    if(dim(I)==2)
    {
-      def newR=basering;
+     def newR=basering;
    }
    else
    {
-      short=0;
-      string smp = string(minpoly);
-      if(dim(I)==0)
-      {
-         execute("ring Rhelp=("+charstr(R0)+"),(@s,@t),dp;");
-      }
-      else
-      {
-         execute("ring Rhelp=("+charstr(R0)+"),(@s),dp;");
-      }
-      if (smp!="0")
-      { execute("minpoly = "+smp+";");}
-      def newR=R+Rhelp;
-      setring newR;
-      ideal I=imap(R,I);
-      I=homog(I,@s);
-      attrib(I,"isSB",1);
+     list @rl0=ringlist(R0);
+     if(dim(I)==0)
+     {
+       @rl0[2]=list("@s","@t");
+       @rl0[3]=list(list("dp",1:2),list("C",0));
+     }
+     else
+     {
+       @rl0[2]=list("@s");
+       @rl0[3]=list(list("dp",1),list("C",0));
+     }
+     def Rhelp=ring(@rl0);
+     kill @rl0;
+     def newR=R+Rhelp;
+     setring newR;
+     ideal I=imap(R,I);
+     I=homog(I,@s);
+     attrib(I,"isSB",1);
    }
 
    if((nvars(basering)<=3)&&(size(I)>1))
    {
-       ERROR("This is not equidimensional");
+     ERROR("This is not equidimensional");
    }
 
@@ -2652 +2659 @@
    option(redSB);
    def R=basering;
-   execute("ring S=("+charstr(R)+"),(x,y),lp;");
+   list @rl=ringlist(R);
+   @rl[2]=list("x","y");
+   @rl[3]=list(list("lp",1:2),list("C",0));
+   def S=ring(@rl); setring S;
    map phi=R,x,y;
    ideal singL=phi(singL);
@@ -2665 +2675 @@
       map alpha=S,var(1)-singL[2][2],var(2)-singL[1][2];
       f=alpha(f);
-      execute("ring C=("+charstr(S)+"),("+varstr(S)+"),ds;");
+      def C=changeordTo(S,"ds"); setring C;
       poly f=imap(S,f);
       ideal singL=imap(S,singL);
@@ -2712 +2722 @@
       poly p=beta(p);
 
-      execute("ring C=("+charstr(S)+",a),("+varstr(S)+"),ds;");
+      execute("ring C=("+charstr(S)+",a),(x,y),ds;");
       number p=number(imap(B,p));