Changeset 8251cc in git for kernel/polys.cc

Timestamp:

Apr 26, 2007, 11:22:34 AM (17 years ago)

Author:

Oliver Wienand <wienand@…>

Branches:

(u'spielwiese', '17f1d200f27c5bd38f5dfc6e8a0879242279d1d8')

Children:

a2589fbeca8f1226cf84a107ae899a974fbe264b

Parents:

c6eecb9d3c2515ce9b40ef65205d28266ba90999

Message:

division of monomials in the case of coefficient rings


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

File:

• kernel/polys.cc (modified) (4 diffs)

kernel/polys.cc

@@ -2 +2 @@
 *  Computer Algebra System SINGULAR     *
 ****************************************/
-/* $Id: polys.cc,v 1.13 2006-12-15 17:16:07 Singular Exp $ */
+/* $Id: polys.cc,v 1.14 2007-04-26 09:22:34 wienand Exp $ */
 
 /*
@@ -110 +110 @@
 }
 
+#ifdef HAVE_RING2TOM
+#define pDiv_nn(p, n)              p_Div_nn(p, n, currRing)
+
+poly p_Div_nn(poly p, const number n, const ring r)
+{
+  pAssume(!n_IsZero(n,r));
+  p_Test(p, r);
+
+  poly q = p;
+  while (p != NULL)
+  {
+    number nc = pGetCoeff(p);
+    pSetCoeff0(p, n_Div(nc, n, r));
+    n_Delete(&nc, r);
+    pIter(p);
+  }
+  p_Test(q, r);
+  return q;
+}
+#endif
+
 /*2
 * divides a by the monomial b, ignores monomials which are not divisible
@@ -120 +141 @@
   poly prev=NULL;
   int i;
+#ifdef HAVE_RING2TOM
+  bool unit = true;
+  number inv = pGetCoeff(b);
+  if ((currRing->cring == 1) && ((long) pGetCoeff(b) % 2 == 0))
+  {
+    unit = false;
+  }
+  if (unit) inv=nInvers(inv);
+#else
   number inv=nInvers(pGetCoeff(b));
+#endif
 
   while (a!=NULL)
@@ -147 +178 @@
     }
   }
+#ifdef HAVE_RING2TOM
+  if (!unit)
+  {
+    pDiv_nn(result,inv);
+  }
+  else
+#endif
   pMult_nn(result,inv);
   nDelete(&inv);