source: git/libpolys/coeffs/longrat.h @ 3216ec

jengelh-datetimespielwiese
Last change on this file since 3216ec was 3216ec, checked in by Hans Schoenemann <hannes@…>, 10 years ago
add: extgcd for bigint
  • Property mode set to 100644
File size: 3.9 KB
Line 
1#ifndef LONGRAT_H
2#define LONGRAT_H
3/****************************************
4*  Computer Algebra System SINGULAR     *
5****************************************/
6/*
7* ABSTRACT: computation with long rational numbers
8*/
9#include <misc/auxiliary.h>
10#include <omalloc/omalloc.h>
11
12#include <coeffs/si_gmp.h>
13#include <coeffs/coeffs.h>
14
15#define SR_HDL(A) ((long)(A))
16
17#define SR_INT    1L
18#define INT_TO_SR(INT)  ((number) (((long)INT << 2) + SR_INT))
19#define SR_TO_INT(SR)   (((long)SR) >> 2)
20
21#define MP_SMALL 1
22
23
24struct snumber;
25typedef struct snumber  *number;
26struct snumber
27{
28  mpz_t z;
29  mpz_t n;
30#if defined(LDEBUG)
31  int debug;
32#endif
33  BOOLEAN s;
34};
35
36// allow inlining only from p_Numbers.h and if ! LDEBUG
37
38#if defined(DO_LINLINE) && defined(P_NUMBERS_H) && !defined(LDEBUG)
39#define LINLINE static inline
40#else
41#define LINLINE
42#undef DO_LINLINE
43#endif // DO_LINLINE
44
45LINLINE BOOLEAN  nlEqual(number a, number b, const coeffs r);
46LINLINE number   nlInit(long i, const coeffs r);
47number nlRInit (long i);
48LINLINE BOOLEAN  nlIsOne(number a, const coeffs r);
49LINLINE BOOLEAN  nlIsZero(number za, const coeffs r);
50LINLINE number   nlCopy(number a, const coeffs r);
51LINLINE number   nl_Copy(number a, const coeffs r);
52LINLINE void     nlDelete(number *a, const coeffs r);
53LINLINE number   nlNeg(number za, const coeffs r);
54LINLINE number   nlAdd(number la, number li, const coeffs r);
55LINLINE number   nlSub(number la, number li, const coeffs r);
56LINLINE number   nlMult(number a, number b, const coeffs r);
57
58BOOLEAN nlInitChar(coeffs r, void*);
59
60number   nlInit2 (int i, int j, const coeffs r);
61number   nlInit2gmp (mpz_t i, mpz_t j);
62
63// number nlInitMPZ(mpz_t m, const coeffs r);
64// void nlMPZ(mpz_t m, number &n, const coeffs r);
65
66number   nlGcd(number a, number b, const coeffs r);
67number nlExtGcd(number a, number b, number *s, number *t, const coeffs);
68number   nlLcm(number a, number b, const coeffs r);   /*special routine !*/
69BOOLEAN  nlGreater(number a, number b, const coeffs r);
70BOOLEAN  nlIsMOne(number a, const coeffs r);
71int      nlInt(number &n, const coeffs r);
72number   nlBigInt(number &n);
73
74#ifdef HAVE_RINGS
75number nlMapGMP(number from, const coeffs src, const coeffs dst);
76void     nlGMP(number &i, number n, const coeffs r);
77#endif
78
79BOOLEAN  nlGreaterZero(number za, const coeffs r);
80number   nlInvers(number a, const coeffs r);
81void     nlNormalize(number &x, const coeffs r);
82number   nlDiv(number a, number b, const coeffs r);
83number   nlExactDiv(number a, number b, const coeffs r);
84number   nlIntDiv(number a, number b, const coeffs r);
85number   nlIntMod(number a, number b, const coeffs r);
86void     nlPower(number x, int exp, number *lu, const coeffs r);
87const char *   nlRead (const char *s, number *a, const coeffs r);
88void     nlWrite(number &a, const coeffs r);
89
90/// Map q \in QQ \to Zp
91number nlModP(number q, const coeffs Q, const coeffs Zp);
92
93int      nlSize(number n, const coeffs r);
94number   nlGetDenom(number &n, const coeffs r);
95number   nlGetNumerator(number &n, const coeffs r);
96void     nlCoeffWrite(const coeffs r, BOOLEAN details);
97number   nlChineseRemainder(number *x, number *q,int rl, const coeffs C);
98number   nlFarey(number nN, number nP, const coeffs CF);
99
100#ifdef LDEBUG
101BOOLEAN  nlDBTest(number a, const char *f, const int l);
102#endif
103extern number nlOne;
104
105nMapFunc nlSetMap(const coeffs src, const coeffs dst);
106
107extern omBin rnumber_bin;
108
109#define FREE_RNUMBER(x) omFreeBin((void *)x, rnumber_bin)
110#define ALLOC_RNUMBER() (number)omAllocBin(rnumber_bin)
111#define ALLOC0_RNUMBER() (number)omAlloc0Bin(rnumber_bin)
112
113// in-place operations
114void nlInpGcd(number &a, number b, const coeffs r);
115void nlInpIntDiv(number &a, number b, const coeffs r);
116
117LINLINE void nlInpAdd(number &a, number b, const coeffs r);
118LINLINE void nlInpMult(number &a, number b, const coeffs r);
119
120#ifdef LDEBUG
121#define nlTest(a, r) nlDBTest(a,__FILE__,__LINE__, r)
122BOOLEAN nlDBTest(number a, char *f,int l, const coeffs r);
123#else
124#define nlTest(a, r) ((void)0)
125#endif
126
127#endif
128
129
Note: See TracBrowser for help on using the repository browser.