[0e1846] | 1 | #ifndef MODULOP_H |
---|
| 2 | #define MODULOP_H |
---|
| 3 | /**************************************** |
---|
| 4 | * Computer Algebra System SINGULAR * |
---|
| 5 | ****************************************/ |
---|
[6ae4f5] | 6 | /* $Id: modulop.h,v 1.3 1997-04-09 12:20:01 Singular Exp $ */ |
---|
[0e1846] | 7 | /* |
---|
[6ae4f5] | 8 | * ABSTRACT: numbers modulo p (<=32003) |
---|
[0e1846] | 9 | */ |
---|
| 10 | #include "structs.h" |
---|
| 11 | |
---|
| 12 | extern int npPrimeM; |
---|
| 13 | extern int npGen; |
---|
| 14 | extern int npMapPrime; |
---|
| 15 | |
---|
| 16 | BOOLEAN npGreaterZero (number k); |
---|
| 17 | number npMult (number a, number b); |
---|
| 18 | number npInit (int i); |
---|
| 19 | int npInt (number &n); |
---|
| 20 | number npCopy (number k1); |
---|
| 21 | number npAdd (number a, number b); |
---|
| 22 | number npSub (number a, number b); |
---|
| 23 | void npPower (number a, int i, number * result); |
---|
| 24 | BOOLEAN npIsZero (number a); |
---|
| 25 | BOOLEAN npIsOne (number a); |
---|
| 26 | BOOLEAN npIsMOne (number a); |
---|
| 27 | number npDiv (number a, number b); |
---|
| 28 | number npNeg (number c); |
---|
| 29 | number npInvers (number c); |
---|
| 30 | BOOLEAN npGreater (number a, number b); |
---|
| 31 | BOOLEAN npEqual (number a, number b); |
---|
| 32 | void npWrite (number &a); |
---|
| 33 | char * npRead (char *s, number *a); |
---|
| 34 | number npIntMod (number a, number b); |
---|
| 35 | #ifdef LDEBUG |
---|
| 36 | BOOLEAN npDBTest (number a, char *f, int l); |
---|
| 37 | #endif |
---|
| 38 | void npSetChar(int c); |
---|
| 39 | //int npGetChar(); |
---|
| 40 | |
---|
| 41 | BOOLEAN npSetMap(int c, char ** par, int nop, number minpol); |
---|
| 42 | number npMapP(number from); |
---|
| 43 | /*-------specials for spolys, do NOT use otherwise--------------------------*/ |
---|
| 44 | /* for npMultM, npSubM, npNegM, npEqualM : */ |
---|
| 45 | extern int npPminus1M; |
---|
| 46 | extern CARDINAL *npExpTable; |
---|
| 47 | extern CARDINAL *npLogTable; |
---|
| 48 | |
---|
| 49 | inline number npMultM(number a, number b) |
---|
| 50 | { |
---|
| 51 | int x = npLogTable[(int)a]+npLogTable[(int)b]; |
---|
| 52 | return (number)npExpTable[x<npPminus1M ? x : x-npPminus1M]; |
---|
| 53 | } |
---|
| 54 | /* |
---|
| 55 | *inline number npMultM(number a, number b) |
---|
| 56 | *{ |
---|
| 57 | * return (number)(((int)a*(int)b) % npPrimeM); |
---|
| 58 | *} |
---|
| 59 | */ |
---|
| 60 | #define npSubM(a,b) (number)((int)a<(int)b ?\ |
---|
| 61 | npPrimeM-(int)b+(int)a : (int)a-(int)b) |
---|
| 62 | #define npNegM(A) (number)(npPrimeM-(int)(A)) |
---|
| 63 | #define npEqualM(A,B) ((int)A==(int)B) |
---|
| 64 | |
---|
| 65 | #endif |
---|
| 66 | |
---|