Changeset d8e71a in git for libpolys/polys/pCoeff.cc

Timestamp:

Aug 25, 2022, 11:43:42 AM (21 months ago)

Author:

Hans Schoenemann <hannes@…>

Branches:

(u'fieker-DuVal', '117eb8c30fc9e991c4decca4832b1d19036c4c65')(u'spielwiese', 'd08f5f0bb3329b8ca19f23b74cb1473686415c3a')

Children:

aaa2126025a96fdc04e3b30a2afd09f23a3bd228

Parents:

2d1047928fc91699976942d5b18b0096edd8439e

Message:

add CoeffTerm(<vector>,<ideal>)

File:

• libpolys/polys/pCoeff.cc (modified) (1 diff)

libpolys/polys/pCoeff.cc

@@ -82 +82 @@
   return res;
 }
+
+/// find coeffs of a vector of a list of given monomials, n>=max_comp(v)
+poly p_CoeffTermId(poly v, ideal m, int n, const ring r)
+{
+  if ((n<=0)||(v==NULL)) return NULL;
+  poly q;
+  poly u=NULL;
+  int dummy;
+  const int ncols_m=IDELEMS(m);
+  v=p_Copy(v,r);
+  for(int i=1;i<=n;i++)
+  {
+    p_TakeOutComp(&v,i,&q,&dummy,r);
+    for(int j=0;j<ncols_m;j++)
+    {
+      number n=p_CoeffTerm(q,m->m[j],r);
+      poly uu=p_NSet(n,r);
+      if (uu!=NULL)
+      {
+        p_SetComp(uu,(i-1)*ncols_m+j+1,r);
+        u=p_Add_q(u,uu,r);
+      }
+    }
+    if (v==NULL) break;
+  }
+  return u;
+}
+