Optimize sha/asm/keccak1600-avx2.pl.
[openssl.git] / crypto / ec / ec_lib.c
1 /*
2  * Copyright 2001-2016 The OpenSSL Project Authors. All Rights Reserved.
3  * Copyright (c) 2002, Oracle and/or its affiliates. All rights reserved
4  *
5  * Licensed under the OpenSSL license (the "License").  You may not use
6  * this file except in compliance with the License.  You can obtain a copy
7  * in the file LICENSE in the source distribution or at
8  * https://www.openssl.org/source/license.html
9  */
10
11 #include <string.h>
12
13 #include <openssl/err.h>
14 #include <openssl/opensslv.h>
15
16 #include "ec_lcl.h"
17
18 /* functions for EC_GROUP objects */
19
20 EC_GROUP *EC_GROUP_new(const EC_METHOD *meth)
21 {
22     EC_GROUP *ret;
23
24     if (meth == NULL) {
25         ECerr(EC_F_EC_GROUP_NEW, EC_R_SLOT_FULL);
26         return NULL;
27     }
28     if (meth->group_init == 0) {
29         ECerr(EC_F_EC_GROUP_NEW, ERR_R_SHOULD_NOT_HAVE_BEEN_CALLED);
30         return NULL;
31     }
32
33     ret = OPENSSL_zalloc(sizeof(*ret));
34     if (ret == NULL) {
35         ECerr(EC_F_EC_GROUP_NEW, ERR_R_MALLOC_FAILURE);
36         return NULL;
37     }
38
39     ret->meth = meth;
40     if ((ret->meth->flags & EC_FLAGS_CUSTOM_CURVE) == 0) {
41         ret->order = BN_new();
42         if (ret->order == NULL)
43             goto err;
44         ret->cofactor = BN_new();
45         if (ret->cofactor == NULL)
46             goto err;
47     }
48     ret->asn1_flag = OPENSSL_EC_NAMED_CURVE;
49     ret->asn1_form = POINT_CONVERSION_UNCOMPRESSED;
50     if (!meth->group_init(ret))
51         goto err;
52     return ret;
53
54  err:
55     BN_free(ret->order);
56     BN_free(ret->cofactor);
57     OPENSSL_free(ret);
58     return NULL;
59 }
60
61 void EC_pre_comp_free(EC_GROUP *group)
62 {
63     switch (group->pre_comp_type) {
64     case PCT_none:
65         break;
66     case PCT_nistz256:
67 #ifdef ECP_NISTZ256_ASM
68         EC_nistz256_pre_comp_free(group->pre_comp.nistz256);
69 #endif
70         break;
71 #ifndef OPENSSL_NO_EC_NISTP_64_GCC_128
72     case PCT_nistp224:
73         EC_nistp224_pre_comp_free(group->pre_comp.nistp224);
74         break;
75     case PCT_nistp256:
76         EC_nistp256_pre_comp_free(group->pre_comp.nistp256);
77         break;
78     case PCT_nistp521:
79         EC_nistp521_pre_comp_free(group->pre_comp.nistp521);
80         break;
81 #else
82     case PCT_nistp224:
83     case PCT_nistp256:
84     case PCT_nistp521:
85         break;
86 #endif
87     case PCT_ec:
88         EC_ec_pre_comp_free(group->pre_comp.ec);
89         break;
90     }
91     group->pre_comp.ec = NULL;
92 }
93
94 void EC_GROUP_free(EC_GROUP *group)
95 {
96     if (!group)
97         return;
98
99     if (group->meth->group_finish != 0)
100         group->meth->group_finish(group);
101
102     EC_pre_comp_free(group);
103     BN_MONT_CTX_free(group->mont_data);
104     EC_POINT_free(group->generator);
105     BN_free(group->order);
106     BN_free(group->cofactor);
107     OPENSSL_free(group->seed);
108     OPENSSL_free(group);
109 }
110
111 void EC_GROUP_clear_free(EC_GROUP *group)
112 {
113     if (!group)
114         return;
115
116     if (group->meth->group_clear_finish != 0)
117         group->meth->group_clear_finish(group);
118     else if (group->meth->group_finish != 0)
119         group->meth->group_finish(group);
120
121     EC_pre_comp_free(group);
122     BN_MONT_CTX_free(group->mont_data);
123     EC_POINT_clear_free(group->generator);
124     BN_clear_free(group->order);
125     BN_clear_free(group->cofactor);
126     OPENSSL_clear_free(group->seed, group->seed_len);
127     OPENSSL_clear_free(group, sizeof(*group));
128 }
129
130 int EC_GROUP_copy(EC_GROUP *dest, const EC_GROUP *src)
131 {
132     if (dest->meth->group_copy == 0) {
133         ECerr(EC_F_EC_GROUP_COPY, ERR_R_SHOULD_NOT_HAVE_BEEN_CALLED);
134         return 0;
135     }
136     if (dest->meth != src->meth) {
137         ECerr(EC_F_EC_GROUP_COPY, EC_R_INCOMPATIBLE_OBJECTS);
138         return 0;
139     }
140     if (dest == src)
141         return 1;
142
143     /* Copy precomputed */
144     dest->pre_comp_type = src->pre_comp_type;
145     switch (src->pre_comp_type) {
146     case PCT_none:
147         dest->pre_comp.ec = NULL;
148         break;
149     case PCT_nistz256:
150 #ifdef ECP_NISTZ256_ASM
151         dest->pre_comp.nistz256 = EC_nistz256_pre_comp_dup(src->pre_comp.nistz256);
152 #endif
153         break;
154 #ifndef OPENSSL_NO_EC_NISTP_64_GCC_128
155     case PCT_nistp224:
156         dest->pre_comp.nistp224 = EC_nistp224_pre_comp_dup(src->pre_comp.nistp224);
157         break;
158     case PCT_nistp256:
159         dest->pre_comp.nistp256 = EC_nistp256_pre_comp_dup(src->pre_comp.nistp256);
160         break;
161     case PCT_nistp521:
162         dest->pre_comp.nistp521 = EC_nistp521_pre_comp_dup(src->pre_comp.nistp521);
163         break;
164 #else
165     case PCT_nistp224:
166     case PCT_nistp256:
167     case PCT_nistp521:
168         break;
169 #endif
170     case PCT_ec:
171         dest->pre_comp.ec = EC_ec_pre_comp_dup(src->pre_comp.ec);
172         break;
173     }
174
175     if (src->mont_data != NULL) {
176         if (dest->mont_data == NULL) {
177             dest->mont_data = BN_MONT_CTX_new();
178             if (dest->mont_data == NULL)
179                 return 0;
180         }
181         if (!BN_MONT_CTX_copy(dest->mont_data, src->mont_data))
182             return 0;
183     } else {
184         /* src->generator == NULL */
185         BN_MONT_CTX_free(dest->mont_data);
186         dest->mont_data = NULL;
187     }
188
189     if (src->generator != NULL) {
190         if (dest->generator == NULL) {
191             dest->generator = EC_POINT_new(dest);
192             if (dest->generator == NULL)
193                 return 0;
194         }
195         if (!EC_POINT_copy(dest->generator, src->generator))
196             return 0;
197     } else {
198         /* src->generator == NULL */
199         EC_POINT_clear_free(dest->generator);
200         dest->generator = NULL;
201     }
202
203     if ((src->meth->flags & EC_FLAGS_CUSTOM_CURVE) == 0) {
204         if (!BN_copy(dest->order, src->order))
205             return 0;
206         if (!BN_copy(dest->cofactor, src->cofactor))
207             return 0;
208     }
209
210     dest->curve_name = src->curve_name;
211     dest->asn1_flag = src->asn1_flag;
212     dest->asn1_form = src->asn1_form;
213
214     if (src->seed) {
215         OPENSSL_free(dest->seed);
216         dest->seed = OPENSSL_malloc(src->seed_len);
217         if (dest->seed == NULL)
218             return 0;
219         if (!memcpy(dest->seed, src->seed, src->seed_len))
220             return 0;
221         dest->seed_len = src->seed_len;
222     } else {
223         OPENSSL_free(dest->seed);
224         dest->seed = NULL;
225         dest->seed_len = 0;
226     }
227
228     return dest->meth->group_copy(dest, src);
229 }
230
231 EC_GROUP *EC_GROUP_dup(const EC_GROUP *a)
232 {
233     EC_GROUP *t = NULL;
234     int ok = 0;
235
236     if (a == NULL)
237         return NULL;
238
239     if ((t = EC_GROUP_new(a->meth)) == NULL)
240         return (NULL);
241     if (!EC_GROUP_copy(t, a))
242         goto err;
243
244     ok = 1;
245
246  err:
247     if (!ok) {
248         EC_GROUP_free(t);
249         return NULL;
250     }
251         return t;
252 }
253
254 const EC_METHOD *EC_GROUP_method_of(const EC_GROUP *group)
255 {
256     return group->meth;
257 }
258
259 int EC_METHOD_get_field_type(const EC_METHOD *meth)
260 {
261     return meth->field_type;
262 }
263
264 int EC_GROUP_set_generator(EC_GROUP *group, const EC_POINT *generator,
265                            const BIGNUM *order, const BIGNUM *cofactor)
266 {
267     if (generator == NULL) {
268         ECerr(EC_F_EC_GROUP_SET_GENERATOR, ERR_R_PASSED_NULL_PARAMETER);
269         return 0;
270     }
271
272     if (group->generator == NULL) {
273         group->generator = EC_POINT_new(group);
274         if (group->generator == NULL)
275             return 0;
276     }
277     if (!EC_POINT_copy(group->generator, generator))
278         return 0;
279
280     if (order != NULL) {
281         if (!BN_copy(group->order, order))
282             return 0;
283     } else
284         BN_zero(group->order);
285
286     if (cofactor != NULL) {
287         if (!BN_copy(group->cofactor, cofactor))
288             return 0;
289     } else
290         BN_zero(group->cofactor);
291
292     /*
293      * Some groups have an order with
294      * factors of two, which makes the Montgomery setup fail.
295      * |group->mont_data| will be NULL in this case.
296      */
297     if (BN_is_odd(group->order)) {
298         return ec_precompute_mont_data(group);
299     }
300
301     BN_MONT_CTX_free(group->mont_data);
302     group->mont_data = NULL;
303     return 1;
304 }
305
306 const EC_POINT *EC_GROUP_get0_generator(const EC_GROUP *group)
307 {
308     return group->generator;
309 }
310
311 BN_MONT_CTX *EC_GROUP_get_mont_data(const EC_GROUP *group)
312 {
313     return group->mont_data;
314 }
315
316 int EC_GROUP_get_order(const EC_GROUP *group, BIGNUM *order, BN_CTX *ctx)
317 {
318     if (group->order == NULL)
319         return 0;
320     if (!BN_copy(order, group->order))
321         return 0;
322
323     return !BN_is_zero(order);
324 }
325
326 const BIGNUM *EC_GROUP_get0_order(const EC_GROUP *group)
327 {
328     return group->order;
329 }
330
331 int EC_GROUP_order_bits(const EC_GROUP *group)
332 {
333     OPENSSL_assert(group->meth->group_order_bits != NULL);
334     return group->meth->group_order_bits(group);
335 }
336
337 int EC_GROUP_get_cofactor(const EC_GROUP *group, BIGNUM *cofactor,
338                           BN_CTX *ctx)
339 {
340
341     if (group->cofactor == NULL)
342         return 0;
343     if (!BN_copy(cofactor, group->cofactor))
344         return 0;
345
346     return !BN_is_zero(group->cofactor);
347 }
348
349 const BIGNUM *EC_GROUP_get0_cofactor(const EC_GROUP *group)
350 {
351     return group->cofactor;
352 }
353
354 void EC_GROUP_set_curve_name(EC_GROUP *group, int nid)
355 {
356     group->curve_name = nid;
357 }
358
359 int EC_GROUP_get_curve_name(const EC_GROUP *group)
360 {
361     return group->curve_name;
362 }
363
364 void EC_GROUP_set_asn1_flag(EC_GROUP *group, int flag)
365 {
366     group->asn1_flag = flag;
367 }
368
369 int EC_GROUP_get_asn1_flag(const EC_GROUP *group)
370 {
371     return group->asn1_flag;
372 }
373
374 void EC_GROUP_set_point_conversion_form(EC_GROUP *group,
375                                         point_conversion_form_t form)
376 {
377     group->asn1_form = form;
378 }
379
380 point_conversion_form_t EC_GROUP_get_point_conversion_form(const EC_GROUP
381                                                            *group)
382 {
383     return group->asn1_form;
384 }
385
386 size_t EC_GROUP_set_seed(EC_GROUP *group, const unsigned char *p, size_t len)
387 {
388     OPENSSL_free(group->seed);
389     group->seed = NULL;
390     group->seed_len = 0;
391
392     if (!len || !p)
393         return 1;
394
395     if ((group->seed = OPENSSL_malloc(len)) == NULL)
396         return 0;
397     memcpy(group->seed, p, len);
398     group->seed_len = len;
399
400     return len;
401 }
402
403 unsigned char *EC_GROUP_get0_seed(const EC_GROUP *group)
404 {
405     return group->seed;
406 }
407
408 size_t EC_GROUP_get_seed_len(const EC_GROUP *group)
409 {
410     return group->seed_len;
411 }
412
413 int EC_GROUP_set_curve_GFp(EC_GROUP *group, const BIGNUM *p, const BIGNUM *a,
414                            const BIGNUM *b, BN_CTX *ctx)
415 {
416     if (group->meth->group_set_curve == 0) {
417         ECerr(EC_F_EC_GROUP_SET_CURVE_GFP, ERR_R_SHOULD_NOT_HAVE_BEEN_CALLED);
418         return 0;
419     }
420     return group->meth->group_set_curve(group, p, a, b, ctx);
421 }
422
423 int EC_GROUP_get_curve_GFp(const EC_GROUP *group, BIGNUM *p, BIGNUM *a,
424                            BIGNUM *b, BN_CTX *ctx)
425 {
426     if (group->meth->group_get_curve == 0) {
427         ECerr(EC_F_EC_GROUP_GET_CURVE_GFP, ERR_R_SHOULD_NOT_HAVE_BEEN_CALLED);
428         return 0;
429     }
430     return group->meth->group_get_curve(group, p, a, b, ctx);
431 }
432
433 #ifndef OPENSSL_NO_EC2M
434 int EC_GROUP_set_curve_GF2m(EC_GROUP *group, const BIGNUM *p, const BIGNUM *a,
435                             const BIGNUM *b, BN_CTX *ctx)
436 {
437     if (group->meth->group_set_curve == 0) {
438         ECerr(EC_F_EC_GROUP_SET_CURVE_GF2M,
439               ERR_R_SHOULD_NOT_HAVE_BEEN_CALLED);
440         return 0;
441     }
442     return group->meth->group_set_curve(group, p, a, b, ctx);
443 }
444
445 int EC_GROUP_get_curve_GF2m(const EC_GROUP *group, BIGNUM *p, BIGNUM *a,
446                             BIGNUM *b, BN_CTX *ctx)
447 {
448     if (group->meth->group_get_curve == 0) {
449         ECerr(EC_F_EC_GROUP_GET_CURVE_GF2M,
450               ERR_R_SHOULD_NOT_HAVE_BEEN_CALLED);
451         return 0;
452     }
453     return group->meth->group_get_curve(group, p, a, b, ctx);
454 }
455 #endif
456
457 int EC_GROUP_get_degree(const EC_GROUP *group)
458 {
459     if (group->meth->group_get_degree == 0) {
460         ECerr(EC_F_EC_GROUP_GET_DEGREE, ERR_R_SHOULD_NOT_HAVE_BEEN_CALLED);
461         return 0;
462     }
463     return group->meth->group_get_degree(group);
464 }
465
466 int EC_GROUP_check_discriminant(const EC_GROUP *group, BN_CTX *ctx)
467 {
468     if (group->meth->group_check_discriminant == 0) {
469         ECerr(EC_F_EC_GROUP_CHECK_DISCRIMINANT,
470               ERR_R_SHOULD_NOT_HAVE_BEEN_CALLED);
471         return 0;
472     }
473     return group->meth->group_check_discriminant(group, ctx);
474 }
475
476 int EC_GROUP_cmp(const EC_GROUP *a, const EC_GROUP *b, BN_CTX *ctx)
477 {
478     int r = 0;
479     BIGNUM *a1, *a2, *a3, *b1, *b2, *b3;
480     BN_CTX *ctx_new = NULL;
481
482     /* compare the field types */
483     if (EC_METHOD_get_field_type(EC_GROUP_method_of(a)) !=
484         EC_METHOD_get_field_type(EC_GROUP_method_of(b)))
485         return 1;
486     /* compare the curve name (if present in both) */
487     if (EC_GROUP_get_curve_name(a) && EC_GROUP_get_curve_name(b) &&
488         EC_GROUP_get_curve_name(a) != EC_GROUP_get_curve_name(b))
489         return 1;
490     if (a->meth->flags & EC_FLAGS_CUSTOM_CURVE)
491         return 0;
492
493     if (ctx == NULL)
494         ctx_new = ctx = BN_CTX_new();
495     if (ctx == NULL)
496         return -1;
497
498     BN_CTX_start(ctx);
499     a1 = BN_CTX_get(ctx);
500     a2 = BN_CTX_get(ctx);
501     a3 = BN_CTX_get(ctx);
502     b1 = BN_CTX_get(ctx);
503     b2 = BN_CTX_get(ctx);
504     b3 = BN_CTX_get(ctx);
505     if (b3 == NULL) {
506         BN_CTX_end(ctx);
507         BN_CTX_free(ctx_new);
508         return -1;
509     }
510
511     /*
512      * XXX This approach assumes that the external representation of curves
513      * over the same field type is the same.
514      */
515     if (!a->meth->group_get_curve(a, a1, a2, a3, ctx) ||
516         !b->meth->group_get_curve(b, b1, b2, b3, ctx))
517         r = 1;
518
519     if (r || BN_cmp(a1, b1) || BN_cmp(a2, b2) || BN_cmp(a3, b3))
520         r = 1;
521
522     /* XXX EC_POINT_cmp() assumes that the methods are equal */
523     if (r || EC_POINT_cmp(a, EC_GROUP_get0_generator(a),
524                           EC_GROUP_get0_generator(b), ctx))
525         r = 1;
526
527     if (!r) {
528         const BIGNUM *ao, *bo, *ac, *bc;
529         /* compare the order and cofactor */
530         ao = EC_GROUP_get0_order(a);
531         bo = EC_GROUP_get0_order(b);
532         ac = EC_GROUP_get0_cofactor(a);
533         bc = EC_GROUP_get0_cofactor(b);
534         if (ao == NULL || bo == NULL) {
535             BN_CTX_end(ctx);
536             BN_CTX_free(ctx_new);
537             return -1;
538         }
539         if (BN_cmp(ao, bo) || BN_cmp(ac, bc))
540             r = 1;
541     }
542
543     BN_CTX_end(ctx);
544     BN_CTX_free(ctx_new);
545
546     return r;
547 }
548
549 /* functions for EC_POINT objects */
550
551 EC_POINT *EC_POINT_new(const EC_GROUP *group)
552 {
553     EC_POINT *ret;
554
555     if (group == NULL) {
556         ECerr(EC_F_EC_POINT_NEW, ERR_R_PASSED_NULL_PARAMETER);
557         return NULL;
558     }
559     if (group->meth->point_init == 0) {
560         ECerr(EC_F_EC_POINT_NEW, ERR_R_SHOULD_NOT_HAVE_BEEN_CALLED);
561         return NULL;
562     }
563
564     ret = OPENSSL_zalloc(sizeof(*ret));
565     if (ret == NULL) {
566         ECerr(EC_F_EC_POINT_NEW, ERR_R_MALLOC_FAILURE);
567         return NULL;
568     }
569
570     ret->meth = group->meth;
571
572     if (!ret->meth->point_init(ret)) {
573         OPENSSL_free(ret);
574         return NULL;
575     }
576
577     return ret;
578 }
579
580 void EC_POINT_free(EC_POINT *point)
581 {
582     if (!point)
583         return;
584
585     if (point->meth->point_finish != 0)
586         point->meth->point_finish(point);
587     OPENSSL_free(point);
588 }
589
590 void EC_POINT_clear_free(EC_POINT *point)
591 {
592     if (!point)
593         return;
594
595     if (point->meth->point_clear_finish != 0)
596         point->meth->point_clear_finish(point);
597     else if (point->meth->point_finish != 0)
598         point->meth->point_finish(point);
599     OPENSSL_clear_free(point, sizeof(*point));
600 }
601
602 int EC_POINT_copy(EC_POINT *dest, const EC_POINT *src)
603 {
604     if (dest->meth->point_copy == 0) {
605         ECerr(EC_F_EC_POINT_COPY, ERR_R_SHOULD_NOT_HAVE_BEEN_CALLED);
606         return 0;
607     }
608     if (dest->meth != src->meth) {
609         ECerr(EC_F_EC_POINT_COPY, EC_R_INCOMPATIBLE_OBJECTS);
610         return 0;
611     }
612     if (dest == src)
613         return 1;
614     return dest->meth->point_copy(dest, src);
615 }
616
617 EC_POINT *EC_POINT_dup(const EC_POINT *a, const EC_GROUP *group)
618 {
619     EC_POINT *t;
620     int r;
621
622     if (a == NULL)
623         return NULL;
624
625     t = EC_POINT_new(group);
626     if (t == NULL)
627         return (NULL);
628     r = EC_POINT_copy(t, a);
629     if (!r) {
630         EC_POINT_free(t);
631         return NULL;
632     }
633     return t;
634 }
635
636 const EC_METHOD *EC_POINT_method_of(const EC_POINT *point)
637 {
638     return point->meth;
639 }
640
641 int EC_POINT_set_to_infinity(const EC_GROUP *group, EC_POINT *point)
642 {
643     if (group->meth->point_set_to_infinity == 0) {
644         ECerr(EC_F_EC_POINT_SET_TO_INFINITY,
645               ERR_R_SHOULD_NOT_HAVE_BEEN_CALLED);
646         return 0;
647     }
648     if (group->meth != point->meth) {
649         ECerr(EC_F_EC_POINT_SET_TO_INFINITY, EC_R_INCOMPATIBLE_OBJECTS);
650         return 0;
651     }
652     return group->meth->point_set_to_infinity(group, point);
653 }
654
655 int EC_POINT_set_Jprojective_coordinates_GFp(const EC_GROUP *group,
656                                              EC_POINT *point, const BIGNUM *x,
657                                              const BIGNUM *y, const BIGNUM *z,
658                                              BN_CTX *ctx)
659 {
660     if (group->meth->point_set_Jprojective_coordinates_GFp == 0) {
661         ECerr(EC_F_EC_POINT_SET_JPROJECTIVE_COORDINATES_GFP,
662               ERR_R_SHOULD_NOT_HAVE_BEEN_CALLED);
663         return 0;
664     }
665     if (group->meth != point->meth) {
666         ECerr(EC_F_EC_POINT_SET_JPROJECTIVE_COORDINATES_GFP,
667               EC_R_INCOMPATIBLE_OBJECTS);
668         return 0;
669     }
670     return group->meth->point_set_Jprojective_coordinates_GFp(group, point, x,
671                                                               y, z, ctx);
672 }
673
674 int EC_POINT_get_Jprojective_coordinates_GFp(const EC_GROUP *group,
675                                              const EC_POINT *point, BIGNUM *x,
676                                              BIGNUM *y, BIGNUM *z,
677                                              BN_CTX *ctx)
678 {
679     if (group->meth->point_get_Jprojective_coordinates_GFp == 0) {
680         ECerr(EC_F_EC_POINT_GET_JPROJECTIVE_COORDINATES_GFP,
681               ERR_R_SHOULD_NOT_HAVE_BEEN_CALLED);
682         return 0;
683     }
684     if (group->meth != point->meth) {
685         ECerr(EC_F_EC_POINT_GET_JPROJECTIVE_COORDINATES_GFP,
686               EC_R_INCOMPATIBLE_OBJECTS);
687         return 0;
688     }
689     return group->meth->point_get_Jprojective_coordinates_GFp(group, point, x,
690                                                               y, z, ctx);
691 }
692
693 int EC_POINT_set_affine_coordinates_GFp(const EC_GROUP *group,
694                                         EC_POINT *point, const BIGNUM *x,
695                                         const BIGNUM *y, BN_CTX *ctx)
696 {
697     if (group->meth->point_set_affine_coordinates == 0) {
698         ECerr(EC_F_EC_POINT_SET_AFFINE_COORDINATES_GFP,
699               ERR_R_SHOULD_NOT_HAVE_BEEN_CALLED);
700         return 0;
701     }
702     if (group->meth != point->meth) {
703         ECerr(EC_F_EC_POINT_SET_AFFINE_COORDINATES_GFP,
704               EC_R_INCOMPATIBLE_OBJECTS);
705         return 0;
706     }
707     if (!group->meth->point_set_affine_coordinates(group, point, x, y, ctx))
708         return 0;
709
710     if (EC_POINT_is_on_curve(group, point, ctx) <= 0) {
711         ECerr(EC_F_EC_POINT_SET_AFFINE_COORDINATES_GFP,
712               EC_R_POINT_IS_NOT_ON_CURVE);
713         return 0;
714     }
715     return 1;
716 }
717
718 #ifndef OPENSSL_NO_EC2M
719 int EC_POINT_set_affine_coordinates_GF2m(const EC_GROUP *group,
720                                          EC_POINT *point, const BIGNUM *x,
721                                          const BIGNUM *y, BN_CTX *ctx)
722 {
723     if (group->meth->point_set_affine_coordinates == 0) {
724         ECerr(EC_F_EC_POINT_SET_AFFINE_COORDINATES_GF2M,
725               ERR_R_SHOULD_NOT_HAVE_BEEN_CALLED);
726         return 0;
727     }
728     if (group->meth != point->meth) {
729         ECerr(EC_F_EC_POINT_SET_AFFINE_COORDINATES_GF2M,
730               EC_R_INCOMPATIBLE_OBJECTS);
731         return 0;
732     }
733     if (!group->meth->point_set_affine_coordinates(group, point, x, y, ctx))
734         return 0;
735
736     if (EC_POINT_is_on_curve(group, point, ctx) <= 0) {
737         ECerr(EC_F_EC_POINT_SET_AFFINE_COORDINATES_GF2M,
738               EC_R_POINT_IS_NOT_ON_CURVE);
739         return 0;
740     }
741     return 1;
742 }
743 #endif
744
745 int EC_POINT_get_affine_coordinates_GFp(const EC_GROUP *group,
746                                         const EC_POINT *point, BIGNUM *x,
747                                         BIGNUM *y, BN_CTX *ctx)
748 {
749     if (group->meth->point_get_affine_coordinates == 0) {
750         ECerr(EC_F_EC_POINT_GET_AFFINE_COORDINATES_GFP,
751               ERR_R_SHOULD_NOT_HAVE_BEEN_CALLED);
752         return 0;
753     }
754     if (group->meth != point->meth) {
755         ECerr(EC_F_EC_POINT_GET_AFFINE_COORDINATES_GFP,
756               EC_R_INCOMPATIBLE_OBJECTS);
757         return 0;
758     }
759     return group->meth->point_get_affine_coordinates(group, point, x, y, ctx);
760 }
761
762 #ifndef OPENSSL_NO_EC2M
763 int EC_POINT_get_affine_coordinates_GF2m(const EC_GROUP *group,
764                                          const EC_POINT *point, BIGNUM *x,
765                                          BIGNUM *y, BN_CTX *ctx)
766 {
767     if (group->meth->point_get_affine_coordinates == 0) {
768         ECerr(EC_F_EC_POINT_GET_AFFINE_COORDINATES_GF2M,
769               ERR_R_SHOULD_NOT_HAVE_BEEN_CALLED);
770         return 0;
771     }
772     if (group->meth != point->meth) {
773         ECerr(EC_F_EC_POINT_GET_AFFINE_COORDINATES_GF2M,
774               EC_R_INCOMPATIBLE_OBJECTS);
775         return 0;
776     }
777     return group->meth->point_get_affine_coordinates(group, point, x, y, ctx);
778 }
779 #endif
780
781 int EC_POINT_add(const EC_GROUP *group, EC_POINT *r, const EC_POINT *a,
782                  const EC_POINT *b, BN_CTX *ctx)
783 {
784     if (group->meth->add == 0) {
785         ECerr(EC_F_EC_POINT_ADD, ERR_R_SHOULD_NOT_HAVE_BEEN_CALLED);
786         return 0;
787     }
788     if ((group->meth != r->meth) || (r->meth != a->meth)
789         || (a->meth != b->meth)) {
790         ECerr(EC_F_EC_POINT_ADD, EC_R_INCOMPATIBLE_OBJECTS);
791         return 0;
792     }
793     return group->meth->add(group, r, a, b, ctx);
794 }
795
796 int EC_POINT_dbl(const EC_GROUP *group, EC_POINT *r, const EC_POINT *a,
797                  BN_CTX *ctx)
798 {
799     if (group->meth->dbl == 0) {
800         ECerr(EC_F_EC_POINT_DBL, ERR_R_SHOULD_NOT_HAVE_BEEN_CALLED);
801         return 0;
802     }
803     if ((group->meth != r->meth) || (r->meth != a->meth)) {
804         ECerr(EC_F_EC_POINT_DBL, EC_R_INCOMPATIBLE_OBJECTS);
805         return 0;
806     }
807     return group->meth->dbl(group, r, a, ctx);
808 }
809
810 int EC_POINT_invert(const EC_GROUP *group, EC_POINT *a, BN_CTX *ctx)
811 {
812     if (group->meth->invert == 0) {
813         ECerr(EC_F_EC_POINT_INVERT, ERR_R_SHOULD_NOT_HAVE_BEEN_CALLED);
814         return 0;
815     }
816     if (group->meth != a->meth) {
817         ECerr(EC_F_EC_POINT_INVERT, EC_R_INCOMPATIBLE_OBJECTS);
818         return 0;
819     }
820     return group->meth->invert(group, a, ctx);
821 }
822
823 int EC_POINT_is_at_infinity(const EC_GROUP *group, const EC_POINT *point)
824 {
825     if (group->meth->is_at_infinity == 0) {
826         ECerr(EC_F_EC_POINT_IS_AT_INFINITY,
827               ERR_R_SHOULD_NOT_HAVE_BEEN_CALLED);
828         return 0;
829     }
830     if (group->meth != point->meth) {
831         ECerr(EC_F_EC_POINT_IS_AT_INFINITY, EC_R_INCOMPATIBLE_OBJECTS);
832         return 0;
833     }
834     return group->meth->is_at_infinity(group, point);
835 }
836
837 /*
838  * Check whether an EC_POINT is on the curve or not. Note that the return
839  * value for this function should NOT be treated as a boolean. Return values:
840  *  1: The point is on the curve
841  *  0: The point is not on the curve
842  * -1: An error occurred
843  */
844 int EC_POINT_is_on_curve(const EC_GROUP *group, const EC_POINT *point,
845                          BN_CTX *ctx)
846 {
847     if (group->meth->is_on_curve == 0) {
848         ECerr(EC_F_EC_POINT_IS_ON_CURVE, ERR_R_SHOULD_NOT_HAVE_BEEN_CALLED);
849         return 0;
850     }
851     if (group->meth != point->meth) {
852         ECerr(EC_F_EC_POINT_IS_ON_CURVE, EC_R_INCOMPATIBLE_OBJECTS);
853         return 0;
854     }
855     return group->meth->is_on_curve(group, point, ctx);
856 }
857
858 int EC_POINT_cmp(const EC_GROUP *group, const EC_POINT *a, const EC_POINT *b,
859                  BN_CTX *ctx)
860 {
861     if (group->meth->point_cmp == 0) {
862         ECerr(EC_F_EC_POINT_CMP, ERR_R_SHOULD_NOT_HAVE_BEEN_CALLED);
863         return -1;
864     }
865     if ((group->meth != a->meth) || (a->meth != b->meth)) {
866         ECerr(EC_F_EC_POINT_CMP, EC_R_INCOMPATIBLE_OBJECTS);
867         return -1;
868     }
869     return group->meth->point_cmp(group, a, b, ctx);
870 }
871
872 int EC_POINT_make_affine(const EC_GROUP *group, EC_POINT *point, BN_CTX *ctx)
873 {
874     if (group->meth->make_affine == 0) {
875         ECerr(EC_F_EC_POINT_MAKE_AFFINE, ERR_R_SHOULD_NOT_HAVE_BEEN_CALLED);
876         return 0;
877     }
878     if (group->meth != point->meth) {
879         ECerr(EC_F_EC_POINT_MAKE_AFFINE, EC_R_INCOMPATIBLE_OBJECTS);
880         return 0;
881     }
882     return group->meth->make_affine(group, point, ctx);
883 }
884
885 int EC_POINTs_make_affine(const EC_GROUP *group, size_t num,
886                           EC_POINT *points[], BN_CTX *ctx)
887 {
888     size_t i;
889
890     if (group->meth->points_make_affine == 0) {
891         ECerr(EC_F_EC_POINTS_MAKE_AFFINE, ERR_R_SHOULD_NOT_HAVE_BEEN_CALLED);
892         return 0;
893     }
894     for (i = 0; i < num; i++) {
895         if (group->meth != points[i]->meth) {
896             ECerr(EC_F_EC_POINTS_MAKE_AFFINE, EC_R_INCOMPATIBLE_OBJECTS);
897             return 0;
898         }
899     }
900     return group->meth->points_make_affine(group, num, points, ctx);
901 }
902
903 /*
904  * Functions for point multiplication. If group->meth->mul is 0, we use the
905  * wNAF-based implementations in ec_mult.c; otherwise we dispatch through
906  * methods.
907  */
908
909 int EC_POINTs_mul(const EC_GROUP *group, EC_POINT *r, const BIGNUM *scalar,
910                   size_t num, const EC_POINT *points[],
911                   const BIGNUM *scalars[], BN_CTX *ctx)
912 {
913     if (group->meth->mul == 0)
914         /* use default */
915         return ec_wNAF_mul(group, r, scalar, num, points, scalars, ctx);
916
917     return group->meth->mul(group, r, scalar, num, points, scalars, ctx);
918 }
919
920 int EC_POINT_mul(const EC_GROUP *group, EC_POINT *r, const BIGNUM *g_scalar,
921                  const EC_POINT *point, const BIGNUM *p_scalar, BN_CTX *ctx)
922 {
923     /* just a convenient interface to EC_POINTs_mul() */
924
925     const EC_POINT *points[1];
926     const BIGNUM *scalars[1];
927
928     points[0] = point;
929     scalars[0] = p_scalar;
930
931     return EC_POINTs_mul(group, r, g_scalar,
932                          (point != NULL
933                           && p_scalar != NULL), points, scalars, ctx);
934 }
935
936 int EC_GROUP_precompute_mult(EC_GROUP *group, BN_CTX *ctx)
937 {
938     if (group->meth->mul == 0)
939         /* use default */
940         return ec_wNAF_precompute_mult(group, ctx);
941
942     if (group->meth->precompute_mult != 0)
943         return group->meth->precompute_mult(group, ctx);
944     else
945         return 1;               /* nothing to do, so report success */
946 }
947
948 int EC_GROUP_have_precompute_mult(const EC_GROUP *group)
949 {
950     if (group->meth->mul == 0)
951         /* use default */
952         return ec_wNAF_have_precompute_mult(group);
953
954     if (group->meth->have_precompute_mult != 0)
955         return group->meth->have_precompute_mult(group);
956     else
957         return 0;               /* cannot tell whether precomputation has
958                                  * been performed */
959 }
960
961 /*
962  * ec_precompute_mont_data sets |group->mont_data| from |group->order| and
963  * returns one on success. On error it returns zero.
964  */
965 int ec_precompute_mont_data(EC_GROUP *group)
966 {
967     BN_CTX *ctx = BN_CTX_new();
968     int ret = 0;
969
970     BN_MONT_CTX_free(group->mont_data);
971     group->mont_data = NULL;
972
973     if (ctx == NULL)
974         goto err;
975
976     group->mont_data = BN_MONT_CTX_new();
977     if (group->mont_data == NULL)
978         goto err;
979
980     if (!BN_MONT_CTX_set(group->mont_data, group->order, ctx)) {
981         BN_MONT_CTX_free(group->mont_data);
982         group->mont_data = NULL;
983         goto err;
984     }
985
986     ret = 1;
987
988  err:
989
990     BN_CTX_free(ctx);
991     return ret;
992 }
993
994 int EC_KEY_set_ex_data(EC_KEY *key, int idx, void *arg)
995 {
996     return CRYPTO_set_ex_data(&key->ex_data, idx, arg);
997 }
998
999 void *EC_KEY_get_ex_data(const EC_KEY *key, int idx)
1000 {
1001     return CRYPTO_get_ex_data(&key->ex_data, idx);
1002 }
1003
1004 int ec_group_simple_order_bits(const EC_GROUP *group)
1005 {
1006     if (group->order == NULL)
1007         return 0;
1008     return BN_num_bits(group->order);
1009 }