crypto/ec/ec_lib.c

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