Changeset 66bff8 in git

Timestamp:

Dec 2, 2010, 4:53:31 PM (12 years ago)

Author:

Hans Schoenemann <hannes@…>

Branches:

(u'jengelh-datetime', 'ceac47cbc86fe4a15902392bdbb9bd2ae0ea02c6')(u'spielwiese', '0604212ebb110535022efecad887940825b97c3f')

Children:

722db68fa52e40537a49b654c41d22424285d092

Parents:

c3005ea46d9ce4dea15139c365c248317f216023

Message:

code cleanup: pGetMaxPower, pDivByMonom, pCancelPolyByMonom

git-svn-id: file:///usr/local/Singular/svn/trunk@13701 2c84dea3-7e68-4137-9b89-c4e89433aadc

File:

• kernel/polys.cc (modified) (1 diff)

kernel/polys.cc

@@ -357 +357 @@
 *p2 in that of p1
 */
-static int pGetMaxPower (poly p1,poly p2)
-{
-  int     i,k,res = 32000; /*a very large integer*/
-
-  if (p1 == NULL) return 0;
-  for (i=1; i<=pVariables; i++)
-  {
-    if ( pGetExp(p2,i) != 0)
-    {
-      k =  pGetExp(p1,i) /  pGetExp(p2,i);
-      if (k < res) res = k;
-    }
-  }
-  return res;
-}
-
-/*2
-*returns the leading monomial of p1 divided by the maximal power of that
-*of p2
-*/
-static poly pDivByMonom (poly p1,poly p2)
-{
-  int     k, i;
-
-  if (p1 == NULL) return NULL;
-  k = pGetMaxPower(p1,p2);
-  if (k == 0)
-    return pHead(p1);
-  else
-  {
-    number n;
-    poly p = pInit();
-
-    p->next = NULL;
-    for (i=1; i<=pVariables; i++)
-    {
-       pSetExp(p,i, pGetExp(p1,i)-k* pGetExp(p2,i));
-    }
-    nPower(p2->coef,k,&n);
-    pSetCoeff0(p,nDiv(p1->coef,n));
-    nDelete(&n);
-    pSetm(p);
-    return p;
-  }
-}
-// Returns as i-th entry of P the coefficient of the (i-1) power of
-// the leading monomial of p2 in p1
-//
-static void pCancelPolyByMonom (poly p1,poly p2,polyset * P,int * SizeOfSet)
-{
-  int   maxPow;
-  poly  p,qp,Coeff;
-
-  if (*P == NULL)
-  {
-    *P = (polyset) omAlloc(5*sizeof(poly));
-    *SizeOfSet = 5;
-  }
-  p = pCopy(p1);
-  while (p != NULL)
-  {
-    qp = p->next;
-    p->next = NULL;
-    maxPow = pGetMaxPower(p,p2);
-    Coeff = pDivByMonom(p,p2);
-    if (maxPow > *SizeOfSet)
-    {
-      pEnlargeSet(P,*SizeOfSet,maxPow+1-*SizeOfSet);
-      *SizeOfSet = maxPow+1;
-    }
-    (*P)[maxPow] = pAdd((*P)[maxPow],Coeff);
-    pDelete(&p);
-    p = qp;
-  }
-}
-
 /*----------utilities for syzygies--------------*/
 poly pTakeOutComp(poly * p, int k)