Changeset 02d009 in git

Timestamp:

Feb 15, 2021, 6:01:06 PM (2 years ago)

Author:

Hans Schoenemann <hannes@…>

Branches:

(u'jengelh-datetime', 'ceac47cbc86fe4a15902392bdbb9bd2ae0ea02c6')(u'spielwiese', '0d6b7fcd9813a1ca1ed4220cfa2b104b97a0a003')

Children:

7b5b095ed2434be9e1783ac23414e01c2a75fec4

Parents:

e6f77c64ad84a4393b83a375b4207c80fe495718

Message:

add: p_Content_n

Location:

libpolys/polys/monomials

Files:

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

libpolys/polys/monomials/p_polys.cc

@@ -2337 +2337 @@
 }
 
+void p_Content_n(poly ph, number &c,const ring r)
+{
+  if (ph==NULL)
+  {
+    c=n_Init(1,r->cf);
+    return;
+  }
+  const coeffs cf=r->cf;
+  if (pNext(ph)==NULL)
+  {
+    c=pGetCoeff(ph);
+    p_SetCoeff0(ph,n_Init(1,cf),r);
+  }
+  if ((cf->cfSubringGcd==ndGcd)
+  || (cf->cfGcd==ndGcd)) /* trivial gcd*/
+  {
+    c=n_Init(1,r->cf);
+    return;
+  }
+  number h;
+  if ((rField_is_Q(r))
+  || (rField_is_Q_a(r))
+  || (rField_is_Zp_a)(r)
+  || (rField_is_Z(r))
+  )
+  {
+    h=p_InitContent(ph,r); /* first guess of a gcd of all coeffs */
+  }
+  else
+  {
+    h=n_Copy(pGetCoeff(ph),cf);
+  }
+  poly p;
+  if(n_IsOne(h,cf))
+  {
+    goto content_finish;
+  }
+  p=ph;
+  // take the SubringGcd of all coeffs
+  while (p!=NULL)
+  {
+    n_Normalize(pGetCoeff(p),cf);
+    number d=n_SubringGcd(h,pGetCoeff(p),cf);
+    n_Delete(&h,cf);
+    h = d;
+    if(n_IsOne(h,cf))
+    {
+      goto content_finish;
+    }
+    pIter(p);
+  }
+  // if found<>1, divide by it
+  p = ph;
+  while (p!=NULL)
+  {
+    number d = n_ExactDiv(pGetCoeff(p),h,cf);
+    p_SetCoeff(p,d,r);
+    pIter(p);
+  }
+content_finish:
+  c=h;
+  // and last: check leading sign:
+  if(!n_GreaterZero(pGetCoeff(ph),r->cf))
+  {
+    c = n_InpNeg(c,r->cf);
+    ph = p_Neg(ph,r);
+  }
+}
+
 #define CLEARENUMERATORS 1
 

libpolys/polys/monomials/p_polys.h

@@ -212 +212 @@
 void      p_ContentForGB(poly p, const ring r);
 void      p_Content(poly p, const ring r);
+void      p_Content_n(poly p, number &c,const ring r);
 #if 1
 // currently only used by Singular/janet