forked from flintlib/flint
-
Notifications
You must be signed in to change notification settings - Fork 0
/
flint.h
169 lines (136 loc) · 3.86 KB
/
flint.h
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
/*=============================================================================
This file is part of FLINT.
FLINT is free software; you can redistribute it and/or modify
it under the terms of the GNU General Public License as published by
the Free Software Foundation; either version 2 of the License, or
(at your option) any later version.
FLINT is distributed in the hope that it will be useful,
but WITHOUT ANY WARRANTY; without even the implied warranty of
MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
GNU General Public License for more details.
You should have received a copy of the GNU General Public License
along with FLINT; if not, write to the Free Software
Foundation, Inc., 51 Franklin St, Fifth Floor, Boston, MA 02110-1301 USA
=============================================================================*/
/******************************************************************************
Copyright (C) 2009 William Hart
******************************************************************************/
#ifndef FLINT_H
#define FLINT_H
#include <mpir.h>
#include <mpfr.h>
#include "longlong.h"
/*
We define alternative key words for "asm" and "inline", allowing
the code to be compiled with the "-ansi" flag under GCC
*/
#ifndef __GNUC__
#define __asm__ asm
#define __inline__ inline
#endif
extern char version[];
#define ulong unsigned long
#if __GMP_BITS_PER_MP_LIMB == 64
#define FLINT_BITS 64
#define FLINT_D_BITS 53
#define FLINT64 1
#else
#define FLINT_BITS 32
#define FLINT_D_BITS 31
#endif
#define mp_bitcnt_t unsigned long
typedef struct
{
gmp_randstate_t gmp_state;
int gmp_init;
mp_limb_t __randval;
mp_limb_t __randval2;
} flint_rand_s;
typedef flint_rand_s flint_rand_t[1];
static __inline__
void flint_randinit(flint_rand_t state)
{
state->gmp_init = 0;
#if FLINT64
state->__randval = 4035456057UL;
state->__randval2 = 6748392731UL;
#else
state->__randval = 3119766748UL;
state->__randval2 = 4225528843UL;
#endif
}
static __inline__
void _flint_rand_init_gmp(flint_rand_t state)
{
if (!state->gmp_init)
{
gmp_randinit_default(state->gmp_state);
state->gmp_init = 1;
}
}
static __inline__
void flint_randclear(flint_rand_t state)
{
if (state->gmp_init)
gmp_randclear(state->gmp_state);
}
/*
We define this here as there is no mpfr.h
*/
typedef __mpfr_struct mpfr;
#define FLINT_ASSERT(param)
#define FLINT_MAX(x, y) ((x) > (y) ? (x) : (y))
#define FLINT_MIN(x, y) ((x) > (y) ? (y) : (x))
#define FLINT_ABS(x) ((long)(x) < 0 ? (-x) : (x))
#define MP_PTR_SWAP(x, y) \
do { \
mp_limb_t * __txxx; \
__txxx = x; \
x = y; \
y = __txxx; \
} while (0)
#define r_shift(in, shift) \
((shift == FLINT_BITS) ? 0L : ((in) >> (shift)))
#define l_shift(in, shift) \
((shift == FLINT_BITS) ? 0L : ((in) << (shift)))
#ifdef NEED_CLZ_TAB
extern unsigned char __flint_clz_tab[128];
#endif
static __inline__
unsigned int FLINT_BIT_COUNT(mp_limb_t x)
{
unsigned int zeros = FLINT_BITS;
if (x) count_leading_zeros(zeros, x);
return FLINT_BITS - zeros;
}
#undef mpn_zero
#define mpn_zero(xxx, nnn) \
do \
{ \
long ixxx; \
for (ixxx = 0; ixxx < nnn; ixxx++) \
(xxx)[ixxx] = 0UL; \
} while (0)
#undef mpn_copyi
#define mpn_copyi(xxx, yyy, nnn) \
do { \
long ixxx; \
for (ixxx = 0; ixxx < nnn; ixxx++) \
(xxx)[ixxx] = (yyy)[ixxx]; \
} while (0)
#undef mpn_copyd
#define mpn_copyd(xxx, yyy, nnn) \
do { \
long ixxx; \
for (ixxx = nnn - 1; ixxx >= 0; ixxx--) \
(xxx)[ixxx] = (yyy)[ixxx]; \
} while (0)
#undef mpn_store
#define mpn_store(xxx, nnn, yyy) \
do \
{ \
long ixxx; \
for (ixxx = 0; ixxx < nnn; ixxx++) \
(xxx)[ixxx] = yyy; \
} while (0)
#endif