Changeset d116de5 in git

Timestamp:

Sep 17, 2022, 5:05:30 PM (19 months ago)

Author:

Hans Schoenemann <hannes@…>

Branches:

(u'spielwiese', '17f1d200f27c5bd38f5dfc6e8a0879242279d1d8')

Children:

64af4704e7419dc5bf9280adb417f0d34f0b6ca2

Parents:

f7c2b38c7fee4700691783a2d80bee341805ebc5

Message:

add id_DelDiv_Sorted0

Location:

libpolys/polys

Files:

• simpleideals.cc (modified) (2 diffs)
• simpleideals.h (modified) (1 diff)

libpolys/polys/simpleideals.cc

@@ -385 +385 @@
 }
 
+/// delete id[j], if LT(j) == coeff*mon*LT(i) (j>i)
+/// assume LC(j)==NULL or n_Delete==ndDelete
+static void id_DelDiv_SEV0(ideal id, int k,const ring r)
+{
+  int kk = k+1;
+  long *sev=(long*)omAlloc0(kk*sizeof(long));
+  while(id->m[k]==NULL) k--;
+  for (int i=k; i>=0; i--)
+  {
+    if(id->m[i]!=NULL)
+      sev[i]=p_GetShortExpVector(id->m[i],r);
+  }
+  for (int i=0; i<k; i++)
+  {
+    if (id->m[i] != NULL)
+    {
+      poly m_i=id->m[i];
+      for (int j=i+1; j<=k; j++)
+      {
+        if (id->m[j]!=NULL)
+        {
+          if (p_LmShortDivisibleBy(m_i, sev[i],r, id->m[j],~sev[j],r))
+          {
+            p_LmFree(&id->m[j],r);
+          }
+          else if (p_LmShortDivisibleBy(id->m[j],sev[j],r, m_i,~sev[i],r))
+          {
+            p_LmFree(&id->m[i],r);
+            break;
+          }
+        }
+      }
+    }
+  }
+  omFreeSize(sev,kk*sizeof(long));
+}
+
+
 /// delete id[j], if LT(j) == coeff*mon*LT(i) and vice versa, i.e.,
 /// delete id[i], if LT(i) == coeff*mon*LT(j)
@@ -456 +494 @@
   int k = IDELEMS(id)-1;
   id_DelDiv_SEV(id,k,r);
+}
+void id_DelDiv_Sorted0(ideal id, const ring r)
+{
+  int k = IDELEMS(id)-1;
+  id_DelDiv_SEV0(id,k,r);
 }
 

libpolys/polys/simpleideals.h

@@ -111 +111 @@
 void id_DelDiv(ideal id, const ring r);
 void id_DelDiv_Sorted(ideal id, const ring r);
+void id_DelDiv_Sorted0(ideal id, const ring r);
 BOOLEAN id_IsConstant(ideal id, const ring r);