[a6904c] | 1 | #ifndef LONGALG_H |
---|
| 2 | #define LONGALG_H |
---|
| 3 | /**************************************** |
---|
| 4 | * Computer Algebra System SINGULAR * |
---|
| 5 | ****************************************/ |
---|
| 6 | /* $Id: longalg.h 12749 2010-04-27 14:17:56Z hannes $ */ |
---|
| 7 | /* |
---|
| 8 | * ABSTRACT: algebraic numbers |
---|
| 9 | */ |
---|
[20b794] | 10 | #include <coeffs/coeffs.h> |
---|
| 11 | #include <coeffs/longrat.h> |
---|
[6fd69c] | 12 | #include <polys/monomials/ring.h> |
---|
[20b794] | 13 | #include <polys/ext_fields/longtrans.h> |
---|
[a6904c] | 14 | |
---|
[3dbee61] | 15 | #define ALLOC_LNUMBER() (lnumber)omAllocBin(rnumber_bin) |
---|
| 16 | #define ALLOC0_LNUMBER() (lnumber)omAlloc0Bin(rnumber_bin) |
---|
| 17 | #define FREE_LNUMBER(x) omFreeBin((ADDRESS)x, rnumber_bin) |
---|
| 18 | |
---|
| 19 | |
---|
| 20 | typedef polyrec * napoly; |
---|
[21972b] | 21 | |
---|
[3dbee61] | 22 | struct slnumber; |
---|
| 23 | typedef struct slnumber * lnumber; |
---|
| 24 | |
---|
| 25 | struct slnumber |
---|
[a6904c] | 26 | { |
---|
[6fd69c] | 27 | int anz; |
---|
| 28 | poly *liste; |
---|
[a6904c] | 29 | }; |
---|
[6fd69c] | 30 | typedef struct snaIdeal * naIdeal; |
---|
| 31 | extern omBin snaIdeal_bin; |
---|
| 32 | extern naIdeal naI; |
---|
| 33 | extern poly naMinimalPoly; |
---|
[9a07cd] | 34 | |
---|
[3dbee61] | 35 | void naSetChar(int p, ring r); |
---|
| 36 | void naDelete (number *p, const ring r); |
---|
| 37 | number naInit(int i, const ring r); /* z := i */ |
---|
| 38 | number naPar(int i); /* z := par(i) */ |
---|
| 39 | int naParDeg(number n); /* i := deg(n) */ |
---|
| 40 | int naSize(number n); /* size desc. */ |
---|
| 41 | int naInt(number &n, const ring r); |
---|
| 42 | BOOLEAN naIsZero(number za); /* za = 0 ? */ |
---|
| 43 | BOOLEAN naIsOne(number za); /* za = 1 ? */ |
---|
| 44 | BOOLEAN naIsMOne(number za); /* za = -1 ? */ |
---|
| 45 | BOOLEAN naEqual(number a, number b); /* a = b ? */ |
---|
| 46 | BOOLEAN naGreater(number a, number b); /* dummy */ |
---|
| 47 | number naNeg(number za); /* za := - za */ |
---|
| 48 | number naInvers(number a); |
---|
| 49 | void naPower(number x, int exp, number *lo); |
---|
| 50 | BOOLEAN naGreaterZero(number a); |
---|
| 51 | number naCopy(number p); /* erg:= p */ |
---|
| 52 | number na_Copy(number p, const ring r); /* erg:= p */ |
---|
| 53 | number naAdd(number la, number li); /* lu := la+li */ |
---|
| 54 | number naMult(number la, number li); /* lo := la*li */ |
---|
| 55 | number naDiv(number la, number li); /* lo := la/li */ |
---|
| 56 | number naIntDiv(number la, number li); /* lo := la/li */ |
---|
| 57 | //number naIntMod(number la, number li); /* lo := la/li */ |
---|
| 58 | number naSub(number la, number li); /* lu := la-li */ |
---|
| 59 | void naNormalize(number &p); |
---|
| 60 | number naGcd(number a, number b, const ring r); |
---|
| 61 | number naLcm(number a, number b, const ring r); |
---|
[a6904c] | 62 | const char * naRead(const char * s, number * p); |
---|
[3dbee61] | 63 | void naWrite(number &p, const ring r); |
---|
| 64 | char * naName(number n); |
---|
| 65 | nMapFunc naSetMap(const ring src, const ring dst); |
---|
| 66 | number naMap0P(number c); |
---|
| 67 | number naMap00(number c); |
---|
[a6904c] | 68 | #ifdef LDEBUG |
---|
[3dbee61] | 69 | BOOLEAN naDBTest(number a, const char *f,const int l); |
---|
[a6904c] | 70 | #endif |
---|
| 71 | |
---|
[3dbee61] | 72 | void naSetIdeal(ideal I); |
---|
| 73 | |
---|
| 74 | // external access to the interna |
---|
| 75 | poly naPermNumber(number z, int * par_perm, int P, ring r); |
---|
| 76 | #define napAddExp(p,i,e) (p_AddExp(p,i,e,currRing->algring)) |
---|
| 77 | #define napLength(p) pLength(p) |
---|
| 78 | #define napNeg(p) (p_Neg(p,currRing->algring)) |
---|
| 79 | #define napVariables naNumbOfPar |
---|
| 80 | #define napGetCoeff(p) pGetCoeff(p) |
---|
| 81 | #define napGetExpFrom(p,i,r) (p_GetExp(p,i,r->algring)) |
---|
| 82 | #define napSetExp(p,i,e) (p_SetExp(p,i,e,currRing->algring)) |
---|
| 83 | #define napNew() (p_Init(currRing->algring)) |
---|
| 84 | #define napAdd(p1,p2) (p_Add_q(p1,p2,currRing->algring)) |
---|
| 85 | #define napSetm(p) p_Setm(p,currRing->algring) |
---|
| 86 | #define nanumber lnumber |
---|
| 87 | napoly napRemainder(napoly f, const napoly g); |
---|
| 88 | number naGetDenom(number &n, const ring r); |
---|
| 89 | number naGetNumerator(number &n, const ring r); |
---|
| 90 | number napGetDenom(number &n, const ring r); |
---|
| 91 | number napGetNumerator(number &n, const ring r); |
---|
| 92 | poly napPermNumber(number z, int * par_perm, int P, ring r); |
---|
[a6904c] | 93 | #endif |
---|
| 94 | |
---|