Changeset d681e8 in git for kernel/polys.cc

Timestamp:

May 19, 2007, 3:22:23 PM (16 years ago)

Author:

Oliver Wienand <wienand@…>

Branches:

(u'jengelh-datetime', 'ceac47cbc86fe4a15902392bdbb9bd2ae0ea02c6')(u'spielwiese', 'a800fe4b3e9d37a38c5a10cc0ae9dfa0c15a4ee6')

Children:

aef8bb403e723cf3c64015379df7cbc448fbb3b3

Parents:

1591a873ae939d12cc28f2fbb1408a9e8fdad9b8

Message:

kInline.cc, kutil.h:
new: k_GetStringLeadTerms

kutil.cc:
lots of KDEBUG messages inside of HAVE_RINGS switches
added ntl xgcd import
new procs:
enterOneStrongPoly
initenterstrongPairs

rmodulon.cc:
nrnLcm small error
Question: Is x == modul Then return 0 faster than return x % modul?
nrnComp missing case corrected
intDiv(0,0) = 0, intDiv(0,1) = 1, intDiv(0, a) = modul DIV a

rmodulo2m.cc:
intDiv(0,0) = 0, intDiv(0,1) = 1, intDiv(0, a) = modul DIV a


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

File:

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

kernel/polys.cc

@@ -2 +2 @@
 *  Computer Algebra System SINGULAR     *
 ****************************************/
-/* $Id: polys.cc,v 1.17 2007-05-11 10:48:04 wienand Exp $ */
+/* $Id: polys.cc,v 1.18 2007-05-19 13:22:23 wienand Exp $ */
 
 /*
@@ -691 +691 @@
 #ifdef HAVE_RINGS  //HACK TODO Oliver
 number nGetUnit(number k) {
-  return (number) nIntDiv(k, nGcd(k, 0, currRing));
+  number unit = nIntDiv(k, nGcd(k, 0, currRing));
+  number gcd = nGcd(unit, 0, currRing);
+  if (!nIsOne(gcd))
+  {
+    number tmp = nMult(unit, unit);
+    number gcd_new = nGcd(tmp, 0, currRing);
+    while (gcd_new != gcd)
+    {
+      gcd = gcd_new;
+      tmp = nMult(tmp, unit);
+      gcd_new = nGcd(tmp, 0, currRing);
+    }
+    unit = nAdd(unit, nIntDiv(0, gcd_new));
+  }
+  return unit;
 }
 #endif