Changeset 5dad8b3 in git

Timestamp:

Nov 10, 2010, 3:16:13 PM (13 years ago)

Author:

Hans Schoenemann <hannes@…>

Branches:

(u'jengelh-datetime', 'ceac47cbc86fe4a15902392bdbb9bd2ae0ea02c6')(u'spielwiese', '00e2e9c41af3fde1273eb3633f4c0c7c3db2579d')

Children:

1446e7f6b858dd03816a2e1b9a0c6022428f8b67

Parents:

40157303341732339fd6bc4f53baf18c09715ae0

Message:

more static routines

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

Location:

kernel

Files:

• polys.cc (modified) (3 diffs)
• polys.h (modified) (1 diff)

kernel/polys.cc

@@ -275 +275 @@
     }
 #endif
-    
     p_Setm(rc,r);
   }
-finish:  
+finish:
   return s;
 }
@@ -354 +353 @@
 }
 
+/*4
+*Returns the exponent of the maximal power of the leading monomial of
+*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;
+  }
+}
 /*2
 *replaces the maximal powers of the leading monomial of p2 in p1 by
@@ -389 +467 @@
 }
 
-/*4
-*Returns the exponent of the maximal power of the leading monomial of
-*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 as i-th entry of P the coefficient of the (i-1) power of
-*the leading monomial of p2 in p1
-*/
-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;
-  }
-}
-
-/*2
-*returns the leading monomial of p1 divided by the maximal power of that
-*of p2
-*/
-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;
-  }
-}
+
 /*----------utilities for syzygies--------------*/
 poly pTakeOutComp(poly * p, int k)

kernel/polys.h

@@ -365 +365 @@
 poly      pDehomogen (poly p1,poly p2,number n);
 BOOLEAN   pIsHomogeneous (poly p);
-
-// returns the leading monomial of p1 divided by the maximal power of
-// that of p2
-poly      pDivByMonom (poly p1,poly p2);
-
-// Returns as i-th entry of P the coefficient of the (i-1) power of
-// the leading monomial of p2 in p1
-void      pCancelPolyByMonom (poly p1,poly p2,polyset * P,int * SizeOfSet);
 
 poly      pPermPoly (poly p, int * perm,const ring OldRing, nMapFunc nMap,