Changeset 2f028d in git for kernel

Timestamp:

May 30, 2016, 5:03:40 PM (8 years ago)

Author:

Adi Popescu <adi_popescum@…>

Branches:

(u'spielwiese', '4a9821a93ffdc22a6696668bd4f6b8c9de3e6c5f')

Children:

e07fa918e46981e53aa69816b1ee0abceadaab88

Parents:

48cbb08675244a2938cd9318ad2f38b0c70e5172

Message:

add: posInSMonFirst

Location:

kernel/GBEngine

Files:

• kutil.cc (modified) (1 diff)
• kutil.h (modified) (1 diff)

kernel/GBEngine/kutil.cc

@@ -4648 +4648 @@
 }
 
+
+// sorts by degree and pLtCmp
+// but puts pure monomials at the beginning
+int posInSMonFirst (const kStrategy strat, const int length,const poly p,
+            const int ecart_p)
+{
+  if (length<0) return 0;
+  if(pNext(p) == NULL) return 0;
+  polyset set=strat->S;
+  int o = p_Deg(p,currRing);
+  int op = p_Deg(set[length],currRing);
+
+  if ((op < o)
+  || ((op == o) && (pLtCmp(set[length],p) == -1)))
+    return length+1;
+  int i;
+  int an = 0;
+  for(i=0;i<=length;i++)
+    if(set[i] != NULL && pNext(set[i]) == NULL)
+      an++;
+  int en= length;
+  loop
+  {
+    if (an >= en-1)
+    {
+      op = p_Deg(set[an],currRing);
+      if ((op < o)
+      || ((op == o) && (pLtCmp(set[an],p) == -1)))
+        return en;
+      return an;
+    }
+    i=(an+en) / 2;
+    op = p_Deg(set[i],currRing);
+    if ((op < o)
+    || ((op == o) && (pLtCmp(set[i],p) == -1)))
+      an=i;
+    else
+      en=i;
+  }
+}
 
 /*2

kernel/GBEngine/kutil.h

@@ -426 +426 @@
 int posInS (const kStrategy strat, const int length, const poly p,
             const int ecart_p);
+int posInSMonFirst (const kStrategy strat, const int length, const poly p,
+            const int ecart_p);
 int posInT0 (const TSet set,const int length,LObject &p);
 int posInT1 (const TSet set,const int length,LObject &p);