Changeset d8e71a in git

Timestamp:

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

Author:

Hans Schoenemann <hannes@…>

Branches:

(u'spielwiese', 'fe61d9c35bf7c61f2b6cbf1b56e25e2f08d536cc')

Children:

aaa2126025a96fdc04e3b30a2afd09f23a3bd228

Parents:

2d1047928fc91699976942d5b18b0096edd8439e

Message:

add CoeffTerm(<vector>,<ideal>)

Files:

• Singular/extra.cc (modified) (2 diffs)
• libpolys/polys/pCoeff.cc (modified) (1 diff)
• libpolys/polys/pCoeff.h (modified) (1 diff)

Singular/extra.cc

@@ -3971 +3971 @@
       const short t5[]={2,VECTOR_CMD,POLY_CMD};
       const short t6[]={2,MODUL_CMD,POLY_CMD};
+      const short t7[]={2,VECTOR_CMD,IDEAL_CMD};
       if (iiCheckTypes(h,t1,0)
       || iiCheckTypes(h,t2,0))
@@ -4003 +4004 @@
         res->data=id_CoeffTermV(p,q,currRing);
         res->rtyp=MODUL_CMD;
+        return FALSE;
+      }
+      else if (iiCheckTypes(h,t7,0)) /* vector,ideal*/
+      {
+        poly p=(poly)h->Data();
+        ideal q=(ideal)h->next->Data();
+        res->data=p_CoeffTermId(p,q,p_MaxComp(p,currRing),currRing);
+        res->rtyp=VECTOR_CMD;
         return FALSE;
       }

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

libpolys/polys/pCoeff.h

@@ -30 +30 @@
 /// find coeffs of (polynomial) m in all vectors from I
 ideal id_CoeffTermV(ideal M, poly m, const ring r);
+
+/// 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);
 #endif