1 /*
2 * CDDL HEADER START
3 *
4 * The contents of this file are subject to the terms of the
5 * Common Development and Distribution License (the "License").
6 * You may not use this file except in compliance with the License.
7 *
8 * You can obtain a copy of the license at usr/src/OPENSOLARIS.LICENSE
9 * or http://www.opensolaris.org/os/licensing.
10 * See the License for the specific language governing permissions
11 * and limitations under the License.
12 *
13 * When distributing Covered Code, include this CDDL HEADER in each
14 * file and include the License file at usr/src/OPENSOLARIS.LICENSE.
15 * If applicable, add the following below this CDDL HEADER, with the
16 * fields enclosed by brackets "[]" replaced with your own identifying
17 * information: Portions Copyright [yyyy] [name of copyright owner]
18 *
19 * CDDL HEADER END
20 */
21 /*
22 * Copyright 2009 Sun Microsystems, Inc. All rights reserved.
23 * Use is subject to license terms.
24 */
25
26 #ifndef _BIGNUM_H
27 #define _BIGNUM_H
28
29 #ifdef __cplusplus
30 extern "C" {
31 #endif
32
33 #include <sys/types.h>
34
35 #if defined(__sparcv9) || defined(__amd64) /* 64-bit chunk size */
36 #ifndef UMUL64
37 #define UMUL64 /* 64-bit multiplication results are supported */
38 #endif
39 #else
40 #define BIGNUM_CHUNK_32
41 #endif
42
43
44 #define BITSINBYTE 8
45
46 /* Bignum "digits" (aka "chunks" or "words") are either 32- or 64-bits */
47 #ifdef BIGNUM_CHUNK_32
48 #define BIG_CHUNK_SIZE 32
49 #define BIG_CHUNK_TYPE uint32_t
50 #define BIG_CHUNK_TYPE_SIGNED int32_t
51 #define BIG_CHUNK_HIGHBIT 0x80000000
52 #define BIG_CHUNK_ALLBITS 0xffffffff
53 #define BIG_CHUNK_LOWHALFBITS 0xffff
54 #define BIG_CHUNK_HALF_HIGHBIT 0x8000
55
56 #else
57 #define BIG_CHUNK_SIZE 64
58 #define BIG_CHUNK_TYPE uint64_t
59 #define BIG_CHUNK_TYPE_SIGNED int64_t
60 #define BIG_CHUNK_HIGHBIT 0x8000000000000000ULL
61 #define BIG_CHUNK_ALLBITS 0xffffffffffffffffULL
62 #define BIG_CHUNK_LOWHALFBITS 0xffffffffULL
63 #define BIG_CHUNK_HALF_HIGHBIT 0x80000000ULL
64 #endif
65
66 #define BITLEN2BIGNUMLEN(x) (((x) + BIG_CHUNK_SIZE - 1) / BIG_CHUNK_SIZE)
67 #define CHARLEN2BIGNUMLEN(x) (((x) + sizeof (BIG_CHUNK_TYPE) - 1) / \
68 sizeof (BIG_CHUNK_TYPE))
69
70 #define BIGNUM_WORDSIZE (BIG_CHUNK_SIZE / BITSINBYTE) /* word size in bytes */
71 #define BIG_CHUNKS_FOR_160BITS ((160 + BIG_CHUNK_SIZE - 1) / BIG_CHUNK_SIZE)
72
73
74 /*
75 * leading 0's are permitted
76 * 0 should be represented by size>=1, size>=len>=1, sign=1,
77 * value[i]=0 for 0<i<len
78 */
79 typedef struct {
80 /* size and len in units of BIG_CHUNK_TYPE words */
81 uint32_t size; /* size of memory allocated for value */
82 uint32_t len; /* number of valid data words in value */
83 int sign; /* 1 for nonnegative, -1 for negative */
84 int malloced; /* 1 if value was malloced, 0 if not */
85 BIG_CHUNK_TYPE *value;
86 } BIGNUM;
87
88 #define BIGTMPSIZE 65
89
90 #define BIG_TRUE 1
91 #define BIG_FALSE 0
92
93 typedef int BIG_ERR_CODE;
94
95 /* error codes */
96 #define BIG_OK 0
97 #define BIG_NO_MEM -1
98 #define BIG_INVALID_ARGS -2
99 #define BIG_DIV_BY_0 -3
100 #define BIG_NO_RANDOM -4
101 #define BIG_GENERAL_ERR -5
102 #define BIG_TEST_FAILED -6
103 #define BIG_BUFFER_TOO_SMALL -7
104
105
106 #define arraysize(x) (sizeof (x) / sizeof (x[0]))
107
108 typedef BIG_ERR_CODE (*big_modexp_ncp_func_ptr)(BIGNUM *result,
109 BIGNUM *ma, BIGNUM *e, BIGNUM *n,
110 BIGNUM *tmp, BIG_CHUNK_TYPE n0, void *ncp, void *req);
111
112 typedef struct {
113 big_modexp_ncp_func_ptr func;
114 void *ncp;
115 void *reqp;
116 } big_modexp_ncp_info_t;
117
118
119 #ifdef USE_FLOATING_POINT
120 void conv_d16_to_i32(uint32_t *i32, double *d16, int64_t *tmp, int ilen);
121 void conv_i32_to_d32(double *d32, uint32_t *i32, int len);
122 void conv_i32_to_d16(double *d16, uint32_t *i32, int len);
123 void conv_i32_to_d32_and_d16(double *d32, double *d16,
124 uint32_t *i32, int len);
125 void mont_mulf_noconv(uint32_t *result, double *dm1, double *dm2, double *dt,
126 double *dn, uint32_t *nint, int nlen, double dn0);
127 #endif /* USE_FLOATING_POINT */
128
129 extern BIGNUM big_One;
130 extern BIGNUM big_Two;
131
132
133 void printbignum(char *aname, BIGNUM *a);
134
135 BIG_ERR_CODE big_init(BIGNUM *number, int size);
136 BIG_ERR_CODE big_extend(BIGNUM *number, int size);
137 void big_finish(BIGNUM *number);
138 void bytestring2bignum(BIGNUM *bn, uchar_t *kn, size_t len);
139 void bignum2bytestring(uchar_t *kn, BIGNUM *bn, size_t len);
140 BIG_ERR_CODE big_mont_rr(BIGNUM *result, BIGNUM *n);
141 BIG_ERR_CODE big_modexp(BIGNUM *result, BIGNUM *a, BIGNUM *e,
142 BIGNUM *n, BIGNUM *n_rr);
143 BIG_ERR_CODE big_modexp_ext(BIGNUM *result, BIGNUM *a, BIGNUM *e,
144 BIGNUM *n, BIGNUM *n_rr, big_modexp_ncp_info_t *info);
145 BIG_ERR_CODE big_modexp_crt(BIGNUM *result, BIGNUM *a, BIGNUM *dmodpminus1,
146 BIGNUM *dmodqminus1, BIGNUM *p, BIGNUM *q, BIGNUM *pinvmodq,
147 BIGNUM *p_rr, BIGNUM *q_rr);
148 BIG_ERR_CODE big_modexp_crt_ext(BIGNUM *result, BIGNUM *a, BIGNUM *dmodpminus1,
149 BIGNUM *dmodqminus1, BIGNUM *p, BIGNUM *q, BIGNUM *pinvmodq,
150 BIGNUM *p_rr, BIGNUM *q_rr, big_modexp_ncp_info_t *info);
151 int big_cmp_abs(BIGNUM *a, BIGNUM *b);
152 BIG_ERR_CODE randombignum(BIGNUM *r, int length);
153 BIG_ERR_CODE big_div_pos(BIGNUM *result, BIGNUM *remainder,
154 BIGNUM *aa, BIGNUM *bb);
155 BIG_ERR_CODE big_ext_gcd_pos(BIGNUM *gcd, BIGNUM *cm, BIGNUM *ce,
156 BIGNUM *m, BIGNUM *e);
157 BIG_ERR_CODE big_add(BIGNUM *result, BIGNUM *aa, BIGNUM *bb);
158 BIG_ERR_CODE big_add_abs(BIGNUM *result, BIGNUM *aa, BIGNUM *bb);
159 BIG_ERR_CODE big_mul(BIGNUM *result, BIGNUM *aa, BIGNUM *bb);
160 void big_shiftright(BIGNUM *result, BIGNUM *aa, int offs);
161 BIG_ERR_CODE big_nextprime_pos(BIGNUM *result, BIGNUM *n);
162 BIG_ERR_CODE big_nextprime_pos_ext(BIGNUM *result, BIGNUM *n,
163 big_modexp_ncp_info_t *info);
164 BIG_ERR_CODE big_sub_pos(BIGNUM *result, BIGNUM *aa, BIGNUM *bb);
165 BIG_ERR_CODE big_copy(BIGNUM *dest, BIGNUM *src);
166 BIG_ERR_CODE big_sub(BIGNUM *result, BIGNUM *aa, BIGNUM *bb);
167 int big_bitlength(BIGNUM *n);
168 BIG_ERR_CODE big_init1(BIGNUM *number, int size,
169 BIG_CHUNK_TYPE *buf, int bufsize);
170 BIG_ERR_CODE big_mont_mul(BIGNUM *ret,
171 BIGNUM *a, BIGNUM *b, BIGNUM *n, BIG_CHUNK_TYPE n0);
172 int big_is_zero(BIGNUM *n);
173 BIG_CHUNK_TYPE big_n0(BIG_CHUNK_TYPE n);
174
175
176 #if defined(HWCAP)
177
178 #if (BIG_CHUNK_SIZE != 32)
179 #error HWCAP works only with 32-bit bignum chunks
180 #endif
181
182 #define BIG_MUL_SET_VEC(r, a, len, digit) \
183 (*big_mul_set_vec_impl)(r, a, len, digit)
184 #define BIG_MUL_ADD_VEC(r, a, len, digit) \
185 (*big_mul_add_vec_impl)(r, a, len, digit)
186 #define BIG_MUL_VEC(r, a, alen, b, blen) \
187 (*big_mul_vec_impl)(r, a, alen, b, blen)
188 #define BIG_SQR_VEC(r, a, len) \
189 (*big_sqr_vec_impl)(r, a, len)
190
191 extern BIG_CHUNK_TYPE (*big_mul_set_vec_impl)
192 (BIG_CHUNK_TYPE *r, BIG_CHUNK_TYPE *a, int len, BIG_CHUNK_TYPE digit);
193 extern BIG_CHUNK_TYPE (*big_mul_add_vec_impl)
194 (BIG_CHUNK_TYPE *r, BIG_CHUNK_TYPE *a, int len, BIG_CHUNK_TYPE digit);
195 extern void (*big_mul_vec_impl)
196 (BIG_CHUNK_TYPE *r, BIG_CHUNK_TYPE *a, int alen, BIG_CHUNK_TYPE *b,
197 int blen);
198 extern void (*big_sqr_vec_impl)
199 (BIG_CHUNK_TYPE *r, BIG_CHUNK_TYPE *a, int len);
200
201 #else /* ! HWCAP */
202
203 #define BIG_MUL_SET_VEC(r, a, len, digit) big_mul_set_vec(r, a, len, digit)
204 #define BIG_MUL_ADD_VEC(r, a, len, digit) big_mul_add_vec(r, a, len, digit)
205 #define BIG_MUL_VEC(r, a, alen, b, blen) big_mul_vec(r, a, alen, b, blen)
206 #define BIG_SQR_VEC(r, a, len) big_sqr_vec(r, a, len)
207
208 extern BIG_CHUNK_TYPE big_mul_set_vec(BIG_CHUNK_TYPE *r, BIG_CHUNK_TYPE *a,
209 int len, BIG_CHUNK_TYPE d);
210 extern BIG_CHUNK_TYPE big_mul_add_vec(BIG_CHUNK_TYPE *r,
211 BIG_CHUNK_TYPE *a, int len, BIG_CHUNK_TYPE d);
212 extern void big_mul_vec(BIG_CHUNK_TYPE *r, BIG_CHUNK_TYPE *a, int alen,
213 BIG_CHUNK_TYPE *b, int blen);
214 extern void big_sqr_vec(BIG_CHUNK_TYPE *r, BIG_CHUNK_TYPE *a, int len);
215
216 #endif /* HWCAP */
217
218 #ifdef __cplusplus
219 }
220 #endif
221
222 #endif /* _BIGNUM_H */