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#ifndef _h_PSTMATH
45#define _h_PSTMATH
46#ifndef DISABLE_PSTM
47
48
49#ifndef CHAR_BIT
50#define CHAR_BIT 8
51#endif
52
53
54
55
56
57
58#ifndef HAVE_NATIVE_INT64
59 #define PSTM_16BIT
60#endif
61
62
63
64
65
66
67
68
69
70
71#ifdef PSTM_8BIT
72
73 typedef unsigned char pstm_digit;
74 typedef unsigned short pstm_word;
75 #define DIGIT_BIT 8
76
77#elif defined(PSTM_16BIT)
78
79 typedef unsigned short pstm_digit;
80 typedef unsigned long pstm_word;
81 #define DIGIT_BIT 16
82
83#elif defined(PSTM_64BIT)
84
85 #ifndef __GNUC__
86 #error "64bit digits requires GCC"
87 #endif
88 typedef unsigned long pstm_digit;
89 typedef unsigned long pstm_word __attribute__ ((mode(TI)));
90 #define DIGIT_BIT 64
91
92#else
93
94 typedef uint32 pstm_digit;
95 typedef uint64 pstm_word;
96 #define DIGIT_BIT 32
97 #define PSTM_32BIT
98#endif
99
100#define PSTM_MASK (pstm_digit)(-1)
101#define PSTM_DIGIT_MAX PSTM_MASK
102
103
104
105
106
107#define PSTM_LT -1
108#define PSTM_EQ 0
109#define PSTM_GT 1
110
111#define PSTM_ZPOS 0
112#define PSTM_NEG 1
113
114#define PSTM_OKAY PS_SUCCESS
115#define PSTM_MEM PS_MEM_FAIL
116
117
118
119
120
121#define PSTM_DEFAULT_INIT 64
122#define PSTM_MAX_SIZE 4096
123
124typedef struct {
125 int used, alloc, sign;
126 pstm_digit *dp;
127
128} pstm_int;
129
130
131
132
133
134#define pstm_iszero(a) (((a)->used == 0) ? PS_TRUE : PS_FALSE)
135#define pstm_iseven(a) (((a)->used > 0 && (((a)->dp[0] & 1) == 0)) ? PS_TRUE : PS_FALSE)
136#define pstm_isodd(a) (((a)->used > 0 && (((a)->dp[0] & 1) == 1)) ? PS_TRUE : PS_FALSE)
137#define pstm_abs(a, b) { pstm_copy(a, b); (b)->sign = 0; }
138
139extern void pstm_set(pstm_int *a, pstm_digit b);
140
141extern void pstm_zero(pstm_int * a);
142
143
144#define pstm_init(pool, a) \
145 pstm_init( a)
146extern int32 pstm_init(psPool_t *pool, pstm_int * a);
147
148
149#define pstm_init_size(pool, a, size) \
150 pstm_init_size( a, size)
151extern int32 pstm_init_size(psPool_t *pool, pstm_int * a, uint32 size);
152
153
154#define pstm_init_copy(pool, a, b, toSqr) \
155 pstm_init_copy( a, b, toSqr)
156extern int32 pstm_init_copy(psPool_t *pool, pstm_int * a, pstm_int * b,
157 int toSqr);
158
159extern int pstm_count_bits (pstm_int * a);
160
161
162#define pstm_init_for_read_unsigned_bin(pool, a, len) \
163 pstm_init_for_read_unsigned_bin( a, len)
164extern int32 pstm_init_for_read_unsigned_bin(psPool_t *pool, pstm_int *a,
165 uint32 len);
166
167extern int32 pstm_read_unsigned_bin(pstm_int *a, unsigned char *b, int32 c);
168
169extern int32 pstm_unsigned_bin_size(pstm_int *a);
170
171extern int32 pstm_copy(pstm_int * a, pstm_int * b);
172
173extern void pstm_exch(pstm_int * a, pstm_int * b);
174
175extern void pstm_clear(pstm_int * a);
176
177extern void pstm_clear_multi(pstm_int *mp0, pstm_int *mp1, pstm_int *mp2,
178 pstm_int *mp3, pstm_int *mp4, pstm_int *mp5, pstm_int *mp6,
179 pstm_int *mp7);
180
181extern int32 pstm_grow(pstm_int * a, int size);
182
183extern void pstm_clamp(pstm_int * a);
184
185extern int32 pstm_cmp(pstm_int * a, pstm_int * b);
186
187extern int32 pstm_cmp_mag(pstm_int * a, pstm_int * b);
188
189extern void pstm_rshd(pstm_int *a, int x);
190
191extern int32 pstm_lshd(pstm_int * a, int b);
192
193
194#define pstm_div(pool, a, b, c, d) \
195 pstm_div( a, b, c, d)
196extern int32 pstm_div(psPool_t *pool, pstm_int *a, pstm_int *b, pstm_int *c,
197 pstm_int *d);
198
199
200#define pstm_div_2d(pool, a, b, c, d) \
201 pstm_div_2d( a, b, c, d)
202extern int32 pstm_div_2d(psPool_t *pool, pstm_int *a, int b, pstm_int *c,
203 pstm_int *d);
204
205extern int32 pstm_div_2(pstm_int * a, pstm_int * b);
206
207extern int32 s_pstm_sub(pstm_int *a, pstm_int *b, pstm_int *c);
208
209extern int32 pstm_sub(pstm_int *a, pstm_int *b, pstm_int *c);
210
211
212#define pstm_sub_d(pool, a, b, c) \
213 pstm_sub_d( a, b, c)
214extern int32 pstm_sub_d(psPool_t *pool, pstm_int *a, pstm_digit b, pstm_int *c);
215
216extern int32 pstm_mul_2(pstm_int * a, pstm_int * b);
217
218
219#define pstm_mod(pool, a, b, c) \
220 pstm_mod( a, b, c)
221extern int32 pstm_mod(psPool_t *pool, pstm_int *a, pstm_int *b, pstm_int *c);
222
223
224#define pstm_mulmod(pool, a, b, c, d) \
225 pstm_mulmod( a, b, c, d)
226extern int32 pstm_mulmod(psPool_t *pool, pstm_int *a, pstm_int *b, pstm_int *c,
227 pstm_int *d);
228
229
230#define pstm_exptmod(pool, G, X, P, Y) \
231 pstm_exptmod( G, X, P, Y)
232extern int32 pstm_exptmod(psPool_t *pool, pstm_int *G, pstm_int *X, pstm_int *P,
233 pstm_int *Y);
234
235extern int32 pstm_2expt(pstm_int *a, int b);
236
237extern int32 pstm_add(pstm_int *a, pstm_int *b, pstm_int *c);
238
239
240#define pstm_to_unsigned_bin(pool, a, b) \
241 pstm_to_unsigned_bin( a, b)
242extern int32 pstm_to_unsigned_bin(psPool_t *pool, pstm_int *a,
243 unsigned char *b);
244
245
246#define pstm_to_unsigned_bin_nr(pool, a, b) \
247 pstm_to_unsigned_bin_nr( a, b)
248extern int32 pstm_to_unsigned_bin_nr(psPool_t *pool, pstm_int *a,
249 unsigned char *b);
250
251extern int32 pstm_montgomery_setup(pstm_int *a, pstm_digit *rho);
252
253
254#define pstm_montgomery_reduce(pool, a, m, mp, paD, paDlen) \
255 pstm_montgomery_reduce( a, m, mp, paD, paDlen)
256extern int32 pstm_montgomery_reduce(psPool_t *pool, pstm_int *a, pstm_int *m,
257 pstm_digit mp, pstm_digit *paD, uint32 paDlen);
258
259#define pstm_mul_comba(pool, A, B, C, paD, paDlen) \
260 pstm_mul_comba( A, B, C, paD, paDlen)
261extern int32 pstm_mul_comba(psPool_t *pool, pstm_int *A, pstm_int *B,
262 pstm_int *C, pstm_digit *paD, uint32 paDlen);
263
264
265#define pstm_sqr_comba(pool, A, B, paD, paDlen) \
266 pstm_sqr_comba( A, B, paD, paDlen)
267extern int32 pstm_sqr_comba(psPool_t *pool, pstm_int *A, pstm_int *B,
268 pstm_digit *paD, uint32 paDlen);
269
270extern int32 pstm_cmp_d(pstm_int *a, pstm_digit b);
271
272extern int32 pstm_montgomery_calc_normalization(pstm_int *a, pstm_int *b);
273
274extern int32 pstm_mul_d(pstm_int *a, pstm_digit b, pstm_int *c);
275
276
277#define pstm_invmod(pool, a, b, c) \
278 pstm_invmod( a, b, c)
279extern int32 pstm_invmod(psPool_t *pool, pstm_int * a, pstm_int * b,
280 pstm_int * c);
281
282#else
283 typedef int32 pstm_int;
284#endif
285#endif
286
287