Changeset 8316b0 in git for libpolys/polys

Timestamp:

Sep 6, 2022, 2:07:08 PM (20 months ago)

Author:

Hans Schoenemann <hannes@…>

Branches:

(u'spielwiese', '2a584933abf2a2d3082034c7586d38bb6de1a30a')

Children:

9487e9aedf3f8533bbf8ca861dae3f53b4dd9c4e

Parents:

61b23dbd8490ceddb7e3db6a82bba09eff76d347

Message:

add system("CoeffTerm",<vector>,<module>)

Location:

libpolys/polys

Files:

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

libpolys/polys/pCoeff.cc

@@ -110 +110 @@
 }
 
+/// find coeffs of a vector of a matrix(module) of given monomials, n>=max_comp(v)
+poly p_CoeffTermMo(poly v, ideal m, int n, const ring r)
+{
+  poly res=NULL;
+  for(int i=0;i<IDELEMS(m);i++)
+  {
+    ideal mm=id_Vec2Ideal(m->m[i],r);
+    poly tmp=p_CoeffTermId(v,mm,n,r);
+    if (i>0)
+    {
+      p_Shift(&tmp,m->rank,r); // m->rank >= p_MaxComp(m->m[i-1],r)
+      res=p_Add_q(res,tmp,r);
+    }
+    else
+      res=tmp;
+  }
+  return res;
+}

libpolys/polys/pCoeff.h

@@ -33 +33 @@
 /// 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);
+
+/// find coeffs of a vector of a matrix(as module) of given monomials, n>=max_comp(v)
+poly p_CoeffTermMo(poly v, ideal m, int n, const ring r);
 #endif