source: git/coeffs/longrat.h @ 95d8df

spielwiese
Last change on this file since 95d8df was 95d8df, checked in by Mohamed Barakat <mohamed.barakat@…>, 14 years ago
make omalloc, misc, reporter, resources, coeffs + make test in coeffs
  • Property mode set to 100644
File size: 3.7 KB
Line 
1#ifndef LONGRAT_H
2#define LONGRAT_H
3/****************************************
4*  Computer Algebra System SINGULAR     *
5****************************************/
6/* $Id$ */
7/*
8* ABSTRACT: computation with long rational numbers
9*/
10#include "coeffs.h"
11// #include "config.h"
12#include "si_gmp.h"
13#include <auxiliary.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#define mpz_size1(A) (ABS((A)->_mp_size))
24//#define mpz_size1(A) mpz_size(A)
25
26struct snumber;
27typedef struct snumber  *number;
28struct snumber
29{
30  mpz_t z;
31  mpz_t n;
32#if defined(LDEBUG)
33  int debug;
34#endif
35  BOOLEAN s;
36};
37
38// allow inlining only from p_Numbers.h and if ! LDEBUG
39
40#if defined(DO_LINLINE) && defined(P_NUMBERS_H) && !defined(LDEBUG)
41#define LINLINE static inline
42#else
43#define LINLINE
44#undef DO_LINLINE
45#endif // DO_LINLINE
46
47LINLINE BOOLEAN  nlEqual(number a, number b, const coeffs r);
48LINLINE number   nlInit(int i, const coeffs r);
49number nlRInit (int i);
50LINLINE BOOLEAN  nlIsOne(number a, const coeffs r);
51LINLINE BOOLEAN  nlIsZero(number za, const coeffs r);
52LINLINE number   nlCopy(number a, const coeffs r);
53LINLINE number   nl_Copy(number a, const coeffs r);
54LINLINE void     nlDelete(number *a, const coeffs r);
55LINLINE number   nlNeg(number za, const coeffs r);
56LINLINE number   nlAdd(number la, number li, const coeffs r);
57LINLINE number   nlSub(number la, number li, const coeffs r);
58LINLINE number   nlMult(number a, number b, const coeffs r);
59
60void     nlInitChar(coeffs r, void*);
61static BOOLEAN nlCoeffsEqual(const coeffs r, n_coeffType n, void * parameter);
62
63number   nlInit2 (int i, int j, const coeffs r);
64number   nlInit2gmp (mpz_t i, mpz_t j);
65number   nlGcd(number a, number b, const coeffs r);
66number   nlLcm(number a, number b, const coeffs r);   /*special routine !*/
67BOOLEAN  nlGreater(number a, number b, const coeffs r);
68BOOLEAN  nlIsMOne(number a, const coeffs r);
69int      nlInt(number &n, const coeffs r);
70number   nlBigInt(number &n);
71
72#ifdef HAVE_RINGS
73void     nlGMP(number &i, number n, const coeffs r);
74#endif
75
76BOOLEAN  nlGreaterZero(number za, const coeffs r);
77number   nlInvers(number a, const coeffs r);
78void     nlNormalize(number &x, const coeffs r);
79number   nlDiv(number a, number b, const coeffs r);
80number   nlExactDiv(number a, number b, const coeffs r);
81number   nlIntDiv(number a, number b, const coeffs r);
82number   nlIntMod(number a, number b, const coeffs r);
83void     nlPower(number x, int exp, number *lu, const coeffs r);
84const char *   nlRead (const char *s, number *a, const coeffs r);
85void     nlWrite(number &a, const coeffs r);
86int      nlModP(number n, int p);
87int      nlSize(number n, const coeffs r);
88number   nlGetDenom(number &n, const coeffs r);
89number   nlGetNumerator(number &n, const coeffs r);
90number   nlChineseRemainder(number *x, number *q,int rl);
91#ifdef LDEBUG
92BOOLEAN  nlDBTest(number a, const char *f, const int l);
93#endif
94extern number nlOne;
95
96nMapFunc nlSetMap(const coeffs src, const coeffs dst);
97
98#ifndef OM_ALLOC_H
99struct omBin_s;
100#endif
101extern omBin_s* rnumber_bin;
102
103
104#define FREE_RNUMBER(x) omFreeBin((void *)x, rnumber_bin)
105#define ALLOC_RNUMBER() (number)omAllocBin(rnumber_bin)
106#define ALLOC0_RNUMBER() (number)omAlloc0Bin(rnumber_bin)
107
108// in-place operations
109void nlInpGcd(number &a, number b, const coeffs r);
110void nlInpIntDiv(number &a, number b, const coeffs r);
111LINLINE number nlInpAdd(number la, number li, const coeffs r);
112LINLINE void nlInpMult(number &a, number b, const coeffs r);
113
114number nlFarey(number nP, number nN);
115
116#ifdef LDEBUG
117#define nlTest(a, r) nlDBTest(a,__FILE__,__LINE__, r)
118BOOLEAN nlDBTest(number a, char *f,int l, const coeffs r);
119#else
120#define nlTest(a, r) ((void)0)
121#endif
122
123#endif
124
125
Note: See TracBrowser for help on using the repository browser.