54b6a45f30815358977b07ed8331745019fe478d
[openssl.git] / crypto / ec / ec_lcl.h
1 /* crypto/ec/ec_lcl.h */
2 /*
3  * Originally written by Bodo Moeller for the OpenSSL project.
4  */
5 /* ====================================================================
6  * Copyright (c) 1998-2003 The OpenSSL Project.  All rights reserved.
7  *
8  * Redistribution and use in source and binary forms, with or without
9  * modification, are permitted provided that the following conditions
10  * are met:
11  *
12  * 1. Redistributions of source code must retain the above copyright
13  *    notice, this list of conditions and the following disclaimer. 
14  *
15  * 2. Redistributions in binary form must reproduce the above copyright
16  *    notice, this list of conditions and the following disclaimer in
17  *    the documentation and/or other materials provided with the
18  *    distribution.
19  *
20  * 3. All advertising materials mentioning features or use of this
21  *    software must display the following acknowledgment:
22  *    "This product includes software developed by the OpenSSL Project
23  *    for use in the OpenSSL Toolkit. (http://www.openssl.org/)"
24  *
25  * 4. The names "OpenSSL Toolkit" and "OpenSSL Project" must not be used to
26  *    endorse or promote products derived from this software without
27  *    prior written permission. For written permission, please contact
28  *    openssl-core@openssl.org.
29  *
30  * 5. Products derived from this software may not be called "OpenSSL"
31  *    nor may "OpenSSL" appear in their names without prior written
32  *    permission of the OpenSSL Project.
33  *
34  * 6. Redistributions of any form whatsoever must retain the following
35  *    acknowledgment:
36  *    "This product includes software developed by the OpenSSL Project
37  *    for use in the OpenSSL Toolkit (http://www.openssl.org/)"
38  *
39  * THIS SOFTWARE IS PROVIDED BY THE OpenSSL PROJECT ``AS IS'' AND ANY
40  * EXPRESSED OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
41  * IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR
42  * PURPOSE ARE DISCLAIMED.  IN NO EVENT SHALL THE OpenSSL PROJECT OR
43  * ITS CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL,
44  * SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT
45  * NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES;
46  * LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION)
47  * HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT,
48  * STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE)
49  * ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED
50  * OF THE POSSIBILITY OF SUCH DAMAGE.
51  * ====================================================================
52  *
53  * This product includes cryptographic software written by Eric Young
54  * (eay@cryptsoft.com).  This product includes software written by Tim
55  * Hudson (tjh@cryptsoft.com).
56  *
57  */
58 /* ====================================================================
59  * Copyright 2002 Sun Microsystems, Inc. ALL RIGHTS RESERVED.
60  *
61  * Portions of the attached software ("Contribution") are developed by 
62  * SUN MICROSYSTEMS, INC., and are contributed to the OpenSSL project.
63  *
64  * The Contribution is licensed pursuant to the OpenSSL open source
65  * license provided above.
66  *
67  * The elliptic curve binary polynomial software is originally written by 
68  * Sheueling Chang Shantz and Douglas Stebila of Sun Microsystems Laboratories.
69  *
70  */
71
72
73 #include <stdlib.h>
74
75 #include <openssl/obj_mac.h>
76 #include <openssl/ec.h>
77
78
79 /* Structure details are not part of the exported interface,
80  * so all this may change in future versions. */
81
82 struct ec_method_st {
83         /* used by EC_METHOD_get_field_type: */
84         int field_type; /* a NID */
85
86         /* used by EC_GROUP_new, EC_GROUP_free, EC_GROUP_clear_free, EC_GROUP_copy: */
87         int (*group_init)(EC_GROUP *);
88         void (*group_finish)(EC_GROUP *);
89         void (*group_clear_finish)(EC_GROUP *);
90         int (*group_copy)(EC_GROUP *, const EC_GROUP *);
91
92         /* used by EC_GROUP_set_curve_GFp, EC_GROUP_get_curve_GFp, */
93         /* EC_GROUP_set_curve_GF2m, and EC_GROUP_get_curve_GF2m: */
94         int (*group_set_curve)(EC_GROUP *, const BIGNUM *p, const BIGNUM *a, const BIGNUM *b, BN_CTX *);
95         int (*group_get_curve)(const EC_GROUP *, BIGNUM *p, BIGNUM *a, BIGNUM *b, BN_CTX *);
96
97         /* used by EC_GROUP_get_degree: */
98         int (*group_get_degree)(const EC_GROUP *);
99
100         /* used by EC_GROUP_check: */
101         int (*group_check_discriminant)(const EC_GROUP *, BN_CTX *);
102
103         /* used by EC_POINT_new, EC_POINT_free, EC_POINT_clear_free, EC_POINT_copy: */
104         int (*point_init)(EC_POINT *);
105         void (*point_finish)(EC_POINT *);
106         void (*point_clear_finish)(EC_POINT *);
107         int (*point_copy)(EC_POINT *, const EC_POINT *);
108
109         /* used by EC_POINT_set_to_infinity,
110          * EC_POINT_set_Jprojective_coordinates_GFp,
111          * EC_POINT_get_Jprojective_coordinates_GFp,
112          * EC_POINT_set_affine_coordinates_GFp,     ..._GF2m,
113          * EC_POINT_get_affine_coordinates_GFp,     ..._GF2m,
114          * EC_POINT_set_compressed_coordinates_GFp, ..._GF2m:
115          */
116         int (*point_set_to_infinity)(const EC_GROUP *, EC_POINT *);
117         int (*point_set_Jprojective_coordinates_GFp)(const EC_GROUP *, EC_POINT *,
118                 const BIGNUM *x, const BIGNUM *y, const BIGNUM *z, BN_CTX *);
119         int (*point_get_Jprojective_coordinates_GFp)(const EC_GROUP *, const EC_POINT *,
120                 BIGNUM *x, BIGNUM *y, BIGNUM *z, BN_CTX *);
121         int (*point_set_affine_coordinates)(const EC_GROUP *, EC_POINT *,
122                 const BIGNUM *x, const BIGNUM *y, BN_CTX *);
123         int (*point_get_affine_coordinates)(const EC_GROUP *, const EC_POINT *,
124                 BIGNUM *x, BIGNUM *y, BN_CTX *);
125         int (*point_set_compressed_coordinates)(const EC_GROUP *, EC_POINT *,
126                 const BIGNUM *x, int y_bit, BN_CTX *);
127
128         /* used by EC_POINT_point2oct, EC_POINT_oct2point: */
129         size_t (*point2oct)(const EC_GROUP *, const EC_POINT *, point_conversion_form_t form,
130                 unsigned char *buf, size_t len, BN_CTX *);
131         int (*oct2point)(const EC_GROUP *, EC_POINT *,
132                 const unsigned char *buf, size_t len, BN_CTX *);
133
134         /* used by EC_POINT_add, EC_POINT_dbl, ECP_POINT_invert: */
135         int (*add)(const EC_GROUP *, EC_POINT *r, const EC_POINT *a, const EC_POINT *b, BN_CTX *);
136         int (*dbl)(const EC_GROUP *, EC_POINT *r, const EC_POINT *a, BN_CTX *);
137         int (*invert)(const EC_GROUP *, EC_POINT *, BN_CTX *);
138
139         /* used by EC_POINT_is_at_infinity, EC_POINT_is_on_curve, EC_POINT_cmp: */
140         int (*is_at_infinity)(const EC_GROUP *, const EC_POINT *);
141         int (*is_on_curve)(const EC_GROUP *, const EC_POINT *, BN_CTX *);
142         int (*point_cmp)(const EC_GROUP *, const EC_POINT *a, const EC_POINT *b, BN_CTX *);
143
144         /* used by EC_POINT_make_affine, EC_POINTs_make_affine: */
145         int (*make_affine)(const EC_GROUP *, EC_POINT *, BN_CTX *);
146         int (*points_make_affine)(const EC_GROUP *, size_t num, EC_POINT *[], BN_CTX *);
147
148         /* used by EC_POINTs_mul, EC_POINT_mul, EC_POINT_precompute_mult, EC_POINT_have_precompute_mult
149          * (default implementations are used if the 'mul' pointer is 0): */
150         int (*mul)(const EC_GROUP *group, EC_POINT *r, const BIGNUM *scalar,
151                 size_t num, const EC_POINT *points[], const BIGNUM *scalars[], BN_CTX *);
152         int (*precompute_mult)(EC_GROUP *group, BN_CTX *);
153         int (*have_precompute_mult)(const EC_GROUP *group);
154
155
156         /* internal functions */
157
158         /* 'field_mul', 'field_sqr', and 'field_div' can be used by 'add' and 'dbl' so that
159          * the same implementations of point operations can be used with different
160          * optimized implementations of expensive field operations: */
161         int (*field_mul)(const EC_GROUP *, BIGNUM *r, const BIGNUM *a, const BIGNUM *b, BN_CTX *);
162         int (*field_sqr)(const EC_GROUP *, BIGNUM *r, const BIGNUM *a, BN_CTX *);
163         int (*field_div)(const EC_GROUP *, BIGNUM *r, const BIGNUM *a, const BIGNUM *b, BN_CTX *);
164
165         int (*field_encode)(const EC_GROUP *, BIGNUM *r, const BIGNUM *a, BN_CTX *); /* e.g. to Montgomery */
166         int (*field_decode)(const EC_GROUP *, BIGNUM *r, const BIGNUM *a, BN_CTX *); /* e.g. from Montgomery */
167         int (*field_set_to_one)(const EC_GROUP *, BIGNUM *r, BN_CTX *);
168 } /* EC_METHOD */;
169
170
171 struct ec_group_st {
172         const EC_METHOD *meth;
173
174         EC_POINT *generator; /* optional */
175         BIGNUM order, cofactor;
176
177         int curve_name;/* optional NID for named curve */
178         int asn1_flag; /* flag to control the asn1 encoding */
179         point_conversion_form_t asn1_form;
180
181         unsigned char *seed; /* optional seed for parameters (appears in ASN1) */
182         size_t seed_len;
183
184         void *extra_data;
185         void *(*extra_data_dup_func)(void *);
186         void (*extra_data_free_func)(void *);
187         void (*extra_data_clear_free_func)(void *);
188
189         /* The following members are handled by the method functions,
190          * even if they appear generic */
191         
192         BIGNUM field; /* Field specification.
193                        * For curves over GF(p), this is the modulus;
194                        * for curves over GF(2^m), this is the 
195                        * irreducible polynomial defining the field.
196                        */
197
198         unsigned int poly[5]; /* Field specification for curves over GF(2^m).
199                                * The irreducible f(t) is then of the form:
200                                *     t^poly[0] + t^poly[1] + ... + t^poly[k]
201                                * where m = poly[0] > poly[1] > ... > poly[k] = 0.
202                                */
203
204         BIGNUM a, b; /* Curve coefficients.
205                       * (Here the assumption is that BIGNUMs can be used
206                       * or abused for all kinds of fields, not just GF(p).)
207                       * For characteristic  > 3,  the curve is defined
208                       * by a Weierstrass equation of the form
209                       *     y^2 = x^3 + a*x + b.
210                       * For characteristic  2,  the curve is defined by
211                       * an equation of the form
212                       *     y^2 + x*y = x^3 + a*x^2 + b.
213                       */
214
215         int a_is_minus3; /* enable optimized point arithmetics for special case */
216
217         void *field_data1; /* method-specific (e.g., Montgomery structure) */
218         void *field_data2; /* method-specific */
219         int (*field_mod_func)(BIGNUM *, const BIGNUM *, const BIGNUM *, BN_CTX *); /* method-specific */
220 } /* EC_GROUP */;
221
222
223 /* Basically a 'mixin' for extra data, but available for EC_GROUPs only
224  * (with visibility limited to 'package' level for now).
225  * We use the function pointers as index for retrieval; this obviates
226  * global ex_data-style index tables.
227  * (Currently, we have one slot only, but is is possible to extend this
228  * if necessary.) */
229 int EC_GROUP_set_extra_data(EC_GROUP *, void *extra_data, void *(*extra_data_dup_func)(void *),
230         void (*extra_data_free_func)(void *), void (*extra_data_clear_free_func)(void *));
231 void *EC_GROUP_get_extra_data(const EC_GROUP *, void *(*extra_data_dup_func)(void *),
232         void (*extra_data_free_func)(void *), void (*extra_data_clear_free_func)(void *));
233 void EC_GROUP_free_extra_data(EC_GROUP *);
234 void EC_GROUP_clear_free_extra_data(EC_GROUP *);
235
236
237
238 struct ec_point_st {
239         const EC_METHOD *meth;
240
241         /* All members except 'meth' are handled by the method functions,
242          * even if they appear generic */
243
244         BIGNUM X;
245         BIGNUM Y;
246         BIGNUM Z; /* Jacobian projective coordinates:
247                    * (X, Y, Z)  represents  (X/Z^2, Y/Z^3)  if  Z != 0 */
248         int Z_is_one; /* enable optimized point arithmetics for special case */
249 } /* EC_POINT */;
250
251
252
253 /* method functions in ec_mult.c
254  * (ec_lib.c uses these as defaults if group->method->mul is 0 */
255 int ec_wNAF_mul(const EC_GROUP *group, EC_POINT *r, const BIGNUM *scalar,
256         size_t num, const EC_POINT *points[], const BIGNUM *scalars[], BN_CTX *);
257 int ec_wNAF_precompute_mult(EC_GROUP *group, BN_CTX *);
258 int ec_wNAF_have_precompute_mult(const EC_GROUP *group);
259
260
261 /* method functions in ecp_smpl.c */
262 int ec_GFp_simple_group_init(EC_GROUP *);
263 void ec_GFp_simple_group_finish(EC_GROUP *);
264 void ec_GFp_simple_group_clear_finish(EC_GROUP *);
265 int ec_GFp_simple_group_copy(EC_GROUP *, const EC_GROUP *);
266 int ec_GFp_simple_group_set_curve(EC_GROUP *, const BIGNUM *p, const BIGNUM *a, const BIGNUM *b, BN_CTX *);
267 int ec_GFp_simple_group_get_curve(const EC_GROUP *, BIGNUM *p, BIGNUM *a, BIGNUM *b, BN_CTX *);
268 int ec_GFp_simple_group_get_degree(const EC_GROUP *);
269 int ec_GFp_simple_group_check_discriminant(const EC_GROUP *, BN_CTX *);
270 int ec_GFp_simple_point_init(EC_POINT *);
271 void ec_GFp_simple_point_finish(EC_POINT *);
272 void ec_GFp_simple_point_clear_finish(EC_POINT *);
273 int ec_GFp_simple_point_copy(EC_POINT *, const EC_POINT *);
274 int ec_GFp_simple_point_set_to_infinity(const EC_GROUP *, EC_POINT *);
275 int ec_GFp_simple_set_Jprojective_coordinates_GFp(const EC_GROUP *, EC_POINT *,
276         const BIGNUM *x, const BIGNUM *y, const BIGNUM *z, BN_CTX *);
277 int ec_GFp_simple_get_Jprojective_coordinates_GFp(const EC_GROUP *, const EC_POINT *,
278         BIGNUM *x, BIGNUM *y, BIGNUM *z, BN_CTX *);
279 int ec_GFp_simple_point_set_affine_coordinates(const EC_GROUP *, EC_POINT *,
280         const BIGNUM *x, const BIGNUM *y, BN_CTX *);
281 int ec_GFp_simple_point_get_affine_coordinates(const EC_GROUP *, const EC_POINT *,
282         BIGNUM *x, BIGNUM *y, BN_CTX *);
283 int ec_GFp_simple_set_compressed_coordinates(const EC_GROUP *, EC_POINT *,
284         const BIGNUM *x, int y_bit, BN_CTX *);
285 size_t ec_GFp_simple_point2oct(const EC_GROUP *, const EC_POINT *, point_conversion_form_t form,
286         unsigned char *buf, size_t len, BN_CTX *);
287 int ec_GFp_simple_oct2point(const EC_GROUP *, EC_POINT *,
288         const unsigned char *buf, size_t len, BN_CTX *);
289 int ec_GFp_simple_add(const EC_GROUP *, EC_POINT *r, const EC_POINT *a, const EC_POINT *b, BN_CTX *);
290 int ec_GFp_simple_dbl(const EC_GROUP *, EC_POINT *r, const EC_POINT *a, BN_CTX *);
291 int ec_GFp_simple_invert(const EC_GROUP *, EC_POINT *, BN_CTX *);
292 int ec_GFp_simple_is_at_infinity(const EC_GROUP *, const EC_POINT *);
293 int ec_GFp_simple_is_on_curve(const EC_GROUP *, const EC_POINT *, BN_CTX *);
294 int ec_GFp_simple_cmp(const EC_GROUP *, const EC_POINT *a, const EC_POINT *b, BN_CTX *);
295 int ec_GFp_simple_make_affine(const EC_GROUP *, EC_POINT *, BN_CTX *);
296 int ec_GFp_simple_points_make_affine(const EC_GROUP *, size_t num, EC_POINT *[], BN_CTX *);
297 int ec_GFp_simple_field_mul(const EC_GROUP *, BIGNUM *r, const BIGNUM *a, const BIGNUM *b, BN_CTX *);
298 int ec_GFp_simple_field_sqr(const EC_GROUP *, BIGNUM *r, const BIGNUM *a, BN_CTX *);
299
300
301 /* method functions in ecp_mont.c */
302 int ec_GFp_mont_group_init(EC_GROUP *);
303 int ec_GFp_mont_group_set_curve(EC_GROUP *, const BIGNUM *p, const BIGNUM *a, const BIGNUM *b, BN_CTX *);
304 void ec_GFp_mont_group_finish(EC_GROUP *);
305 void ec_GFp_mont_group_clear_finish(EC_GROUP *);
306 int ec_GFp_mont_group_copy(EC_GROUP *, const EC_GROUP *);
307 int ec_GFp_mont_field_mul(const EC_GROUP *, BIGNUM *r, const BIGNUM *a, const BIGNUM *b, BN_CTX *);
308 int ec_GFp_mont_field_sqr(const EC_GROUP *, BIGNUM *r, const BIGNUM *a, BN_CTX *);
309 int ec_GFp_mont_field_encode(const EC_GROUP *, BIGNUM *r, const BIGNUM *a, BN_CTX *);
310 int ec_GFp_mont_field_decode(const EC_GROUP *, BIGNUM *r, const BIGNUM *a, BN_CTX *);
311 int ec_GFp_mont_field_set_to_one(const EC_GROUP *, BIGNUM *r, BN_CTX *);
312
313
314 /* method functions in ecp_recp.c */
315 int ec_GFp_recp_group_init(EC_GROUP *);
316 int ec_GFp_recp_group_set_curve(EC_GROUP *, const BIGNUM *p, const BIGNUM *a, const BIGNUM *b, BN_CTX *);
317 void ec_GFp_recp_group_finish(EC_GROUP *);
318 void ec_GFp_recp_group_clear_finish(EC_GROUP *);
319 int ec_GFp_recp_group_copy(EC_GROUP *, const EC_GROUP *);
320 int ec_GFp_recp_field_mul(const EC_GROUP *, BIGNUM *r, const BIGNUM *a, const BIGNUM *b, BN_CTX *);
321 int ec_GFp_recp_field_sqr(const EC_GROUP *, BIGNUM *r, const BIGNUM *a, BN_CTX *);
322
323
324 /* method functions in ecp_nist.c */
325 int ec_GFp_nist_group_set_curve(EC_GROUP *, const BIGNUM *p, const BIGNUM *a, const BIGNUM *b, BN_CTX *);
326 int ec_GFp_nist_field_mul(const EC_GROUP *, BIGNUM *r, const BIGNUM *a, const BIGNUM *b, BN_CTX *);
327 int ec_GFp_nist_field_sqr(const EC_GROUP *, BIGNUM *r, const BIGNUM *a, BN_CTX *);
328
329
330 /* method functions in ec2_smpl.c */
331 int ec_GF2m_simple_group_init(EC_GROUP *);
332 void ec_GF2m_simple_group_finish(EC_GROUP *);
333 void ec_GF2m_simple_group_clear_finish(EC_GROUP *);
334 int ec_GF2m_simple_group_copy(EC_GROUP *, const EC_GROUP *);
335 int ec_GF2m_simple_group_set_curve(EC_GROUP *, const BIGNUM *p, const BIGNUM *a, const BIGNUM *b, BN_CTX *);
336 int ec_GF2m_simple_group_get_curve(const EC_GROUP *, BIGNUM *p, BIGNUM *a, BIGNUM *b, BN_CTX *);
337 int ec_GF2m_simple_group_get_degree(const EC_GROUP *);
338 int ec_GF2m_simple_group_check_discriminant(const EC_GROUP *, BN_CTX *);
339 int ec_GF2m_simple_point_init(EC_POINT *);
340 void ec_GF2m_simple_point_finish(EC_POINT *);
341 void ec_GF2m_simple_point_clear_finish(EC_POINT *);
342 int ec_GF2m_simple_point_copy(EC_POINT *, const EC_POINT *);
343 int ec_GF2m_simple_point_set_to_infinity(const EC_GROUP *, EC_POINT *);
344 int ec_GF2m_simple_point_set_affine_coordinates(const EC_GROUP *, EC_POINT *,
345         const BIGNUM *x, const BIGNUM *y, BN_CTX *);
346 int ec_GF2m_simple_point_get_affine_coordinates(const EC_GROUP *, const EC_POINT *,
347         BIGNUM *x, BIGNUM *y, BN_CTX *);
348 int ec_GF2m_simple_set_compressed_coordinates(const EC_GROUP *, EC_POINT *,
349         const BIGNUM *x, int y_bit, BN_CTX *);
350 size_t ec_GF2m_simple_point2oct(const EC_GROUP *, const EC_POINT *, point_conversion_form_t form,
351         unsigned char *buf, size_t len, BN_CTX *);
352 int ec_GF2m_simple_oct2point(const EC_GROUP *, EC_POINT *,
353         const unsigned char *buf, size_t len, BN_CTX *);
354 int ec_GF2m_simple_add(const EC_GROUP *, EC_POINT *r, const EC_POINT *a, const EC_POINT *b, BN_CTX *);
355 int ec_GF2m_simple_dbl(const EC_GROUP *, EC_POINT *r, const EC_POINT *a, BN_CTX *);
356 int ec_GF2m_simple_invert(const EC_GROUP *, EC_POINT *, BN_CTX *);
357 int ec_GF2m_simple_is_at_infinity(const EC_GROUP *, const EC_POINT *);
358 int ec_GF2m_simple_is_on_curve(const EC_GROUP *, const EC_POINT *, BN_CTX *);
359 int ec_GF2m_simple_cmp(const EC_GROUP *, const EC_POINT *a, const EC_POINT *b, BN_CTX *);
360 int ec_GF2m_simple_make_affine(const EC_GROUP *, EC_POINT *, BN_CTX *);
361 int ec_GF2m_simple_points_make_affine(const EC_GROUP *, size_t num, EC_POINT *[], BN_CTX *);
362 int ec_GF2m_simple_field_mul(const EC_GROUP *, BIGNUM *r, const BIGNUM *a, const BIGNUM *b, BN_CTX *);
363 int ec_GF2m_simple_field_sqr(const EC_GROUP *, BIGNUM *r, const BIGNUM *a, BN_CTX *);
364 int ec_GF2m_simple_field_div(const EC_GROUP *, BIGNUM *r, const BIGNUM *a, const BIGNUM *b, BN_CTX *);
365
366
367 /* method functions in ec2_mult.c */
368 int ec_GF2m_simple_mul(const EC_GROUP *group, EC_POINT *r, const BIGNUM *scalar,
369         size_t num, const EC_POINT *points[], const BIGNUM *scalars[], BN_CTX *);
370 int ec_GF2m_precompute_mult(EC_GROUP *group, BN_CTX *ctx);
371 int ec_GF2m_have_precompute_mult(const EC_GROUP *group);