Changeset cf5fc11 in git

Timestamp:

Jan 27, 2005, 5:44:13 PM (18 years ago)

Author:

Hans Schönemann <hannes@…>

Branches:

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

Children:

e21afbaa54be68a2741df57da6a4156aba4b81d7

Parents:

c95632f6751087d2f62be544ec7e07d86913a78a

Message:

*hannes: pCleardenom_n


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

Location:

kernel

Files:

• polys.h (modified) (2 diffs)
• polys1.cc (modified) (2 diffs)

kernel/polys.h

@@ -4 +4 @@
 *  Computer Algebra System SINGULAR     *
 ****************************************/
-/* $Id: polys.h,v 1.2 2004-04-23 13:56:40 Singular Exp $ */
+/* $Id: polys.h,v 1.3 2005-01-27 16:42:15 Singular Exp $ */
 /*
 * ABSTRACT - all basic methods to manipulate polynomials of the
@@ -360 +360 @@
 void      pSimpleContent(poly p, int s);
 void      pCleardenom(poly p);
+void      pCleardenom_n(poly p,number &c);
 void      p_Normalize(poly p, ring r);
 #define   pNormalize(p) p_Normalize(p,currRing)

kernel/polys1.cc

@@ -2 +2 @@
 *  Computer Algebra System SINGULAR     *
 ****************************************/
-/* $Id: polys1.cc,v 1.6 2004-08-27 12:22:34 Singular Exp $ */
+/* $Id: polys1.cc,v 1.7 2005-01-27 16:44:13 Singular Exp $ */
 
 /*
@@ -921 +921 @@
 }
 
+void pCleardenom_n(poly ph,number &c)
+{
+  number d, h;
+  poly p;
+
+  p = ph;
+  if(pNext(p)==NULL)
+  {
+    c=nInvers(pGetCoeff(p));
+    pSetCoeff(p,nInit(1));
+  }
+  else
+  {
+    h = nInit(1);
+    while (p!=NULL)
+    {
+      nNormalize(pGetCoeff(p));
+      d=nLcm(h,pGetCoeff(p),currRing);
+      nDelete(&h);
+      h=d;
+      pIter(p);
+    }
+    c=h;
+    /* contains the 1/lcm of all denominators */
+    if(!nIsOne(h))
+    {
+      p = ph;
+      while (p!=NULL)
+      {
+        /* should be:
+        * number hh;
+        * nGetDenom(p->coef,&hh);
+        * nMult(&h,&hh,&d);
+        * nNormalize(d);
+        * nDelete(&hh);
+        * nMult(d,p->coef,&hh);
+        * nDelete(&d);
+        * nDelete(&(p->coef));
+        * p->coef =hh;
+        */
+        d=nMult(h,pGetCoeff(p));
+        nNormalize(d);
+        pSetCoeff(p,d);
+        pIter(p);
+      }
+      if (nGetChar()==1)
+      {
+        loop
+        {
+          h = nInit(1);
+          p=ph;
+          while (p!=NULL)
+          {
+            d=nLcm(h,pGetCoeff(p),currRing);
+            nDelete(&h);
+            h=d;
+            pIter(p);
+          }
+          /* contains the 1/lcm of all denominators */
+          if(!nIsOne(h))
+          {
+            p = ph;
+            while (p!=NULL)
+            {
+              /* should be:
+              * number hh;
+              * nGetDenom(p->coef,&hh);
+              * nMult(&h,&hh,&d);
+              * nNormalize(d);
+              * nDelete(&hh);
+              * nMult(d,p->coef,&hh);
+              * nDelete(&d);
+              * nDelete(&(p->coef));
+              * p->coef =hh;
+              */
+              d=nMult(h,pGetCoeff(p));
+              nNormalize(d);
+              pSetCoeff(p,d);
+              pIter(p);
+            }
+            number t=nMult(c,h);
+            nDelete(&c);
+            c=t;
+          }
+          else
+          {
+            break;
+          }
+          nDelete(&h);
+        }
+      }
+    }
+  }
+}
+
 /*2
 *tests if p is homogeneous with respect to the actual weigths