Changeset 8c484e in git for coeffs/modulop.cc

Timestamp:

May 5, 2010, 3:05:12 PM (13 years ago)

Author:

Hans Schoenemann <hannes@…>

Branches:

(u'spielwiese', '8e0ad00ce244dfd0756200662572aef8402f13d5')

Children:

2336d0a674d92415a0e20e0b4abf6b1a96ab56e2

Parents:

63c3a8df898418cfb9484633a8d9eb78c299a95c

git-author:

Hans Schoenemann <hannes@mathematik.uni-kl.de>2010-05-05 15:05:12+02:00

git-committer:

Mohamed Barakat <mohamed.barakat@rwth-aachen.de>2011-11-09 11:51:39+01:00

Message:

init for Z/p

File:

• coeffs/modulop.cc (modified) (8 diffs)

coeffs/modulop.cc

@@ -8 +8 @@
 
 #include <string.h>
+#include "config.h"
+#include <omalloc.h>
 #include "coeffs.h"
+#include "output.h"
 #include "numbers.h"
 #include "longrat.h"
 #include "mpr_complex.h"
+#include "mylimits.h"
 #include "modulop.h"
 
 int npGen=0;
-long npMapPrime;
 
 #ifdef HAVE_DIV_MOD
@@ -243 +246 @@
       ii *= 10;
       ii += *s++ - '0';
-      if (ii >= (MAX_INT_VAL / 10)) ii = ii % r->npPrimeM;
+      if (ii >= (MAX_INT / 10)) ii = ii % r->npPrimeM;
     }
     while (((*s) >= '0') && ((*s) <= '9'));
@@ -304 +307 @@
 }
 
-void npInitChar(int c, coeffs r)
+void npKillChar(coeffs r)
+{
+  #ifdef HAVE_DIV_MOD
+  if (r->npInvTable!=NULL)
+  omFreeSize( (void *)r->npInvTable, r->npPrimeM*sizeof(unsigned short) );
+  r->npInvTable=NULL;
+  #else
+  if (r->npExpTable!=NULL)
+  {
+    omFreeSize( (void *)r->npExpTable, r->npPrimeM*sizeof(unsigned short) );
+    omFreeSize( (void *)r->npLogTable, r->npPrimeM*sizeof(unsigned short) );
+    r->npExpTable=NULL; r->npLogTable=NULL;
+  }
+  #endif
+}
+
+void npInitChar(coeffs r, int c)
 {
   int i, w;
@@ -353 +372 @@
 #endif
     }
+    r->cfKillChar=npKillChar;
+    r->cfSetChar=NULL;
+    r->cfInit = npInit;
+    r->n_Int  = npInt;
+    r->nAdd   = npAdd;
+    r->nSub   = npSub;
+    r->nMult  = npMult;
+    r->nDiv   = npDiv;
+    r->nExactDiv= npDiv;
+    r->nNeg   = npNeg;
+    r->nInvers= npInvers;
+    r->cfCopy  = ndCopy;
+    r->nGreater = npGreater;
+    r->nEqual = npEqual;
+    r->nIsZero = npIsZero;
+    r->nIsOne = npIsOne;
+    r->nIsMOne = npIsMOne;
+    r->nGreaterZero = npGreaterZero;
+    r->cfWrite = npWrite;
+    r->nRead = npRead;
+    r->nPower = npPower;
+    r->cfSetMap = npSetMap;
+    r->nName= ndName; 
+    r->nSize  = ndSize;
+#ifdef LDEBUG
+    r->nDBTest=npDBTest;
+#endif
+#ifdef NV_OPS
+    if (c>NV_MAX_PRIME)
+    {
+      r->nMult  = nvMult;
+      r->nDiv   = nvDiv;
+      r->nExactDiv= nvDiv;
+      r->nInvers= nvInvers;
+      r->nPower= nvPower;
+    }
+#endif
   }
   else
@@ -372 +428 @@
 #endif
 
-number npMap0(number from, const coeffs dst_r)
+number npMap0(number from, const coeffs src, const coeffs dst_r)
 {
   return npInit(nlModP(from,dst_r->npPrimeM),dst_r);
 }
 
-number npMapP(number from, const coeffs dst_r)
+number npMapP(number from, const coeffs src, const coeffs dst_r)
 {
   long i = (long)from;
-  if (i>npMapPrime/2)
-  {
-    i-=npMapPrime;
+  if (i>src->npPrimeM/2)
+  {
+    i-=src->npPrimeM;
     while (i < 0) i+=dst_r->npPrimeM;
   }
@@ -389 +445 @@
 }
 
-static number npMapLongR(number from, const coeffs dst_r)
+static number npMapLongR(number from, const coeffs src, const coeffs dst_r)
 {
   gmp_float *ff=(gmp_float*)from;
@@ -455 +511 @@
   mpz_clear(dest);
   if(res->s==0)
-    iz=(long)npDiv((number)iz,(number)in);
+    iz=(long)npDiv((number)iz,(number)in,dst_r);
   omFreeBin((void *)res, rnumber_bin);
   return (number)iz;
@@ -507 +563 @@
     if (n_GetChar(src) == n_GetChar(dst))
     {
-      return ndCopy;
+      return ndCopyMap;
     }
     else
     {
-      npMapPrime=n_GetChar(src);
       return npMapP;
     }