Changeset 4b278c in git

Timestamp:

Jan 22, 2015, 1:19:41 PM (9 years ago)

Author:

dere <dkifle16@…>

Branches:

(u'spielwiese', '8e0ad00ce244dfd0756200662572aef8402f13d5')

Children:

2a5eb8ef50edb654ddd5bbc8e061d551ea44cc74d181eea62bf528f5777c224499fbf26cccac5ff7

Parents:

29b5e9d336d980cc76f0fb6c9d67818a144f3ed3

Message:

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);
 }