source: git/libpolys/polys/ext_fields/longtrans.h @ 6ccdd3a

Last change on this file since 6ccdd3a was 6ccdd3a, checked in by Frank Seelisch <seelisch@…>, 12 years ago
renamed algring; comments; transext.* added
  • Property mode set to 100644
File size: 7.2 KB
1#ifndef LONGTRANS_H
2#define LONGTRANS_H
4*  Computer Algebra System SINGULAR     *
6/* $Id: longtrans.h 12469 2011-02-25 13:38:49Z seelisch $ */
8* ABSTRACT:   numbers in transcendental field extensions,
9              i.e., in rational function fields
11#include <coeffs/longrat.h>
12#include <polys/monomials/polys-impl.h>
16   Instantiation of an algebraic field extension in SINGULAR
17   works by first creating a transcendental field extension
18   and then providing a minimal polynomial / minimal ideal.
19   Consequently, first the code for transcendental field
20   extensions will be activated; see longtrans.*.
21   When providing a minimal polynomial / minimal ideal, all
22   function pointers will be re-defined, using the below method
23   redefineFunctionPointers(). After that, the code for algebraic
24   field extensions is active; see longalg.*.
26   NOTE:
27   Elements of algebraic and transcendental field extensions
28   are polynomials or quotients of two polynomials, respectively.
29   All these polynomials, as well as the minimal polynomial (in
30   the case of an algebraic extension) live in the globally
31   accessible ring 'nacRing', defined in longtrans.*.
34   (not true for types / structs / global variables):
35   nap*   macros and methods operating on polynomials living in
36          nacRing (defined in longtrans.*),
37   na*    (but not nap*) methods and macros for algebraic field
38          extensions (defined in longalg.*),
39   nt*    methods and macros for transcendental field extensions,
40          (defined in longtrans.*)
41   nac*   function pointers for computing with the coefficients of
42          polynomials living in nacRing (defined in longtrans.*)
45struct slnumber;
46typedef struct slnumber * lnumber;
47struct slnumber
49  poly z;
50  poly n;
51  BOOLEAN s;
54extern int  ntNumbOfPar;
55#define     ntParNames (currRing->parameter)
56extern int  ntIsChar0;                                  /* == 1 iff char = 0,
57                                                           == 0 otherwise */ 
58extern ring ntMapRing;
59extern int  ntParsToCopy;
61/* specific methods / macros for transcendental field extensions */
62void        ntSetChar(int p, ring r);
63void        ntDelete (number *p, const ring r);
64number      ntInit(int i, const ring r);                /* z := i      */
65number      ntPar(int i);                               /* z := par(i) */
66int         ntParDeg(number n);                         /* i := deg(n) */
67int         ntSize(number n);                           /* size desc.  */
68int         ntInt(number &n, const ring r);
69BOOLEAN     ntIsZero(number za);                        /* za = 0 ?    */
70BOOLEAN     ntIsOne(number  za);                        /* za = 1 ?    */
71BOOLEAN     ntIsMOne(number  za);                       /* za = -1 ?   */
72BOOLEAN     ntEqual(number a, number b);                /* a = b ?     */
73BOOLEAN     ntGreater(number a, number b);              /* dummy       */
74number      ntNeg(number za);                           /* za := - za  */
75number      ntInvers(number a);
76void        ntPower(number x, int exp, number *lo);
77BOOLEAN     ntGreaterZero(number a);
78number      ntCopy(number p);                           /* erg:= p     */
79number      nt_Copy(number p, const ring r);            /* erg:= p     */
80number      ntAdd(number la, number li);                /* lu := la+li */
81number      ntMult(number la, number li);               /* lo := la*li */
82number      ntDiv(number la, number li);                /* lo := la/li */
83number      ntIntDiv(number la, number li);             /* lo := la/li */
84//number      ntIntMod(number la, number li);           /* lo := la/li */
85number      ntSub(number la, number li);                /* lu := la-li */
86void        ntNormalize(number &p);
87number      ntGcd(number a, number b, const ring r);
88number      ntLcm(number a, number b, const ring r);
89const char* ntRead(const char * s, number * p);
90void        ntWrite(number &p, const ring r);
91char*       ntName(number n);
92nMapFunc    ntSetMap(const ring src, const ring dst);
93number      ntMap0P(number c);
94number      ntMap00(number c);
95#ifdef LDEBUG
96BOOLEAN     ntDBTest(number a, const char *f,const int l);
98poly      ntRemainder(poly f, const poly  g);
99void        ntSetIdeal(ideal I);
100void        ntCoefNormalize(number pp);
101extern number (*ntMap)(number from);
103/* procedure variables for computing with the coefficients of
104   polynomials living in nacRing */
105extern ring        nacRing;
106extern numberfunc  nacMult, nacSub, nacAdd, nacDiv, nacIntDiv;
107extern number      (*nacGcd)(number a, number b, const ring r);
108extern number      (*nacLcm)(number a, number b, const ring r);
109extern number      (*nacInit)(int i, const ring r);
110extern int         (*nacInt)(number &n, const ring r);
111extern void        (*nacNormalize)(number &a);
112extern number      (*nacNeg)(number a);
113extern number      (*nacCopy)(number a);
114extern number      (*nacInvers)(number a);
115extern BOOLEAN     (*nacIsZero)(number a);
116extern BOOLEAN     (*nacIsOne)(number a);
117extern BOOLEAN     (*nacGreaterZero)(number a);
118extern BOOLEAN     (*nacGreater)(number a, number b);
119extern number      (*nacMap)(number);
121/* for computing with polynomials living in nacRing */
122poly    napPermNumber(number z, int * par_perm, int P, ring r);
123#define napAddExp(p,i,e)       (p_AddExp(p,i,e,currRing->extRing))
124#define napLength(p)           pLength(p)
125#define napNeg(p)              (p_Neg(p,currRing->extRing))
126#define napVariables           naNumbOfPar
127#define napGetCoeff(p)         pGetCoeff(p)
128#define napGetExpFrom(p,i,r)   (p_GetExp(p,i,r->extRing))
129#define napSetExp(p,i,e)       (p_SetExp(p,i,e,currRing->extRing))
130#define napNew()               (p_Init(currRing->extRing))
131#define napAdd(p1,p2)          (p_Add_q(p1,p2,currRing->extRing))
132#define napSetm(p)             p_Setm(p,currRing->extRing)
133#define napCopy(p)             p_Copy(p,nacRing)
134#define napSetCoeff(p,n)       {n_Delete(&pGetCoeff(p),nacRing);pGetCoeff(p)=n;}
135#define napComp(p,q)           p_LmCmp((poly)p,(poly)q, nacRing)
136#define napMultT(A,E)          A=(poly)p_Mult_mm((poly)A,(poly)E,nacRing)
137#define napDeg(p)              (int)p_Totaldegree(p, nacRing)
138number  napGetDenom(number &n, const ring r);
139number  napGetNumerator(number &n, const ring r);
140void    napTest(poly p);
141poly  napInitz(number z);
142poly  napCopyNeg(const poly p);
143void    napMultN(poly p, const number z);
144void    napDivMod(poly f, const poly g, poly *q, poly *r);
145poly  napRemainder(poly f, const poly g);
146poly  napInvers(poly x, const poly c);
147int     napMaxDeg(poly p);
148int     napMaxDegLen(poly p, int &l);
149void    napWrite(poly p,const BOOLEAN has_denom, const ring r);
150int     napExp(poly a, poly b);
151int     napExpi(int i, poly a, poly b);
152void    napContent(poly ph);
153void    napCleardenom(poly ph);
154poly  napGcd0(poly a, poly b);
155poly  napGcd(poly a, poly b);
156number  napLcm(poly a);
157BOOLEAN napDivPoly (poly p, poly q);
158poly  napRedp (poly q);
159poly  napTailred (poly q);
160poly  napMap(poly p);
161poly  napPerm(poly p, const int *par_perm, const ring src_ring,
162                const nMapFunc nMap);
163const char* napHandleMons(const char *s, int i, poly ex);
164const char* napHandlePars(const char *s, int i, poly ex);
165const char* napRead(const char *s, poly *b);
Note: See TracBrowser for help on using the repository browser.