Changeset 90df7c in git

Timestamp:

Aug 20, 2013, 1:30:50 PM (10 years ago)

Author:

Hans Schoenemann <hannes@…>

Branches:

(u'jengelh-datetime', 'ceac47cbc86fe4a15902392bdbb9bd2ae0ea02c6')(u'spielwiese', 'a800fe4b3e9d37a38c5a10cc0ae9dfa0c15a4ee6')

Children:

066132a004c5c341957275f8c4b32c2090e4b532

Parents:

ab4cb22fc65acd82da63f7582fb678a52d0f755f

Message:

fix: by gpfister: extractZ

File:

• Singular/LIB/primdecint.lib (modified) (2 diffs)

Singular/LIB/primdecint.lib

@@ -1064 +1064 @@
 }
 
-////////////////////////////////////////////////////////////////////////////////
-
+//////////////////////////////////////////////////////////////////////////////
 static proc extractZ(ideal J, int j, list L, list B)
 {
@@ -1098 +1097 @@
          I[1]=0;
          I=simplify(I,2);
-         //=== this is our way to obtain the coefficients in Z[u] of the
-         //=== leading terms of the Groebner basis above
-         string quotring=Primdec::prepareQuotientring(nvars(basering)-L[1][3]);
-         execute(quotring);
-         ideal I=imap(Shelp,I);
-         list C;
-         int i;
-         for(i=1;i<=size(I);i++)
-         {
-            C[i]=leadcoef(I[i]);
-         }
-         setring Shelp;
-         list C=imap(quring,C);
-
+         if(L[1][3]==nvars(basering))
+         {
+            list C;
+            int i;
+            for(i=1;i<=size(I);i++)
+            {
+               C[i]=I[i];
+            }
+         }
+         else
+         {
+            //=== this is our way to obtain the coefficients in Z[u] of the
+            //=== leading terms of the Groebner basis above
+            string quotring=Primdec::prepareQuotientring(nvars(basering)-L[1][3]);
+            execute(quotring);
+            ideal I=imap(Shelp,I);
+            list C;
+            int i;
+            for(i=1;i<=size(I);i++)
+            {
+               C[i]=leadcoef(I[i]);
+            }
+            setring Shelp;
+            list C=imap(quring,C);
+         }
          setring R;
          list C=imap(Shelp,C);