Changeset 8c484e in git

Timestamp:

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

Author:

Hans Schoenemann <hannes@…>

Branches:

(u'fieker-DuVal', '117eb8c30fc9e991c4decca4832b1d19036c4c65')(u'spielwiese', 'b4f17ed1d25f93d46dbe29e4b499baecc2fd51bb')

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

Location:

coeffs

Files:

• coeffs.h (modified) (1 diff)
• modulop.cc (modified) (8 diffs)
• modulop.h (modified) (1 diff)
• numbers.cc (modified) (5 diffs)
• numbers.h (modified) (1 diff)

coeffs/coeffs.h

@@ -128 +128 @@
    // ?
    // initialisation:
-   void (*cfInitChar)(coeffs r); // do one-time initialisations
+   void (*cfInitChar)(coeffs r, int parameter); // do one-time initialisations
    void (*cfKillChar)(coeffs r); //  undo all initialisations
                                 // or NULL

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;
     }

coeffs/modulop.h

@@ -24 +24 @@
 extern int npGen;
 extern long npMapPrime;
+
+void npInitChar(coeffs r, int ch);
 
 BOOLEAN npGreaterZero (number k, const coeffs r);

coeffs/numbers.cc

@@ -396 +396 @@
   /*----------------------char. p----------------*/
   {
+    /* never again:
     npInitChar(c,r);
     n->cfInit = npInit;
@@ -417 +418 @@
     n->nPower = npPower;
     n->cfSetMap = npSetMap;
-    /* nName= ndName; */
-    /*nSize  = ndSize;*/
+    // nName= ndName;
+    // nSize  = ndSize;
 #ifdef LDEBUG
     n->nDBTest=npDBTest;
@@ -433 +434 @@
     }
 #endif
+    */
+    r->cfInitChar=npInitChar;
+    npInitChar(r,c);
   }
   /* -------------- GF(p^m) -----------------------*/
@@ -594 +598 @@
         if (cf_root==r) cf_root=n->next;
         r->cfDelete(&(r->nNULL),r);
+        if (r->cfKillChar!=NULL) r->cfKillChar(r);
+        /* was:
         switch(r->type)
         {
-          case n_Zp:
-               #ifdef HAVE_DIV_MOD
-               if (r->npInvTable!=NULL)
-               omFreeSize( (void *)r->npInvTable,
-                           r->npPrimeM*sizeof(unsigned short) );
-               #else
-               if (r->npExpTable!=NULL)
-               {
-                 omFreeSize( (void *)r->npExpTable,
-                             r->npPrimeM*sizeof(unsigned short) );
-                 omFreeSize( (void *)r->npLogTable,
-                             r->npPrimeM*sizeof(unsigned short) );
-               }
-               #endif
-               break;
           case n_Zp_a:
           case n_Q_a:
@@ -626 +617 @@
                  break;
           }
+          */
           omFreeSize((void *)r, sizeof(n_Procs_s));
           r=NULL;

coeffs/numbers.h

@@ -66 +66 @@
 number ndCopy(number a, const coeffs r);
 number ndCopyMap(number a, const coeffs r, const coeffs aRing);
+int ndSize(number a, const coeffs r);
+char * ndName(number n, const coeffs r);
 
 void   ndInpMult(number &a, number b, const coeffs r);