Changeset d181ee in git

Timestamp:

Jan 22, 2015, 4:56:36 PM (9 years ago)

Author:

Hans Schoenemann <hannes@…>

Branches:

(u'spielwiese', '8e0ad00ce244dfd0756200662572aef8402f13d5')

Children:

76178cc2bf8ca154192fe714f17fee049c88365c

Parents:

29b5e9d336d980cc76f0fb6c9d67818a144f3ed34b278ca27388e13838313dcd3b914338f9e5dbf3

Message:

Merge pull request #688 from tigistdereje/algemodstd

chg: use check_leadmonom_and_size() in Algemodstd::LiftPolyCRT()

File:

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

Singular/LIB/algemodstd.lib

@@ -293 +293 @@
 ////////////////////////////////////////////////////////////////////////////////
 
+static proc check_leadmonom_and_size(list L)
+{
+    /*
+     * compare the size of ideals in the list and 
+     * check the corresponding leading monomials
+     * size(L)>=2
+     */
+    ideal J=L[1];
+    int i=size(L);
+    int sc=ncols(J);
+    int j,k;
+    poly g=leadmonom(J[1]);
+    for(j=1;j<=i;j++)
+    {
+        if(ncols(L[j])!=sc)
+        {
+            return(0);
+        }
+    }
+    for(k=2;k<=i;k++)
+    {
+        for(j=1;j<=sc;j++)
+        {
+            if(leadmonom(J[j])!=leadmonom(L[k][j]))
+            {
+                return(0);
+            }
+        }
+    }
+    return(1);
+}
+
+////////////////////////////////////////////////////////////////////////////////
+
 static proc LiftPolyCRT(ideal I)
 {
@@ -324 +358 @@
         }
     }
-    // apply CRT for polynomials
-    II =chinrempoly(Lk,LL),f;
+    if(check_leadmonom_and_size(Lk))
+    {
+        // apply CRT for polynomials
+        II =chinrempoly(Lk,LL),f;
+    }
+    else
+    {
+        II=0;
+    }
     return(II);
 }