crypto/ec/ec_lib.c

   1 /* crypto/ec/ec_lib.c */
   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  * Binary polynomial ECC support in OpenSSL originally developed by
  61  * SUN MICROSYSTEMS, INC., and contributed to the OpenSSL project.
  62  */
  63
  64 #include <string.h>
  65
  66 #include <openssl/err.h>
  67 #include <openssl/opensslv.h>
  68
  69 #include "ec_lcl.h"
  70
  71 const char EC_version[] = "EC" OPENSSL_VERSION_PTEXT;
  72
  73 /* functions for EC_GROUP objects */
  74
  75 EC_GROUP *EC_GROUP_new(const EC_METHOD *meth)
  76 {
  77     EC_GROUP *ret;
  78
  79     if (meth == NULL) {
  80         ECerr(EC_F_EC_GROUP_NEW, EC_R_SLOT_FULL);
  81         return NULL;
  82     }
  83     if (meth->group_init == 0) {
  84         ECerr(EC_F_EC_GROUP_NEW, ERR_R_SHOULD_NOT_HAVE_BEEN_CALLED);
  85         return NULL;
  86     }
  87
  88     ret = OPENSSL_malloc(sizeof *ret);
  89     if (ret == NULL) {
  90         ECerr(EC_F_EC_GROUP_NEW, ERR_R_MALLOC_FAILURE);
  91         return NULL;
  92     }
  93
  94     ret->meth = meth;
  95
  96     ret->extra_data = NULL;
  97     ret->mont_data = NULL;
  98
  99     ret->generator = NULL;
 100     ret->order = BN_new();
 101     ret->cofactor = NULL;
 102     if (!ret->order)
 103         goto err;
 104     ret->cofactor = BN_new();
 105     if (!ret->cofactor)
 106         goto err;
 107
 108     ret->curve_name = 0;
 109     ret->asn1_flag = OPENSSL_EC_NAMED_CURVE;
 110     ret->asn1_form = POINT_CONVERSION_UNCOMPRESSED;
 111
 112     ret->seed = NULL;
 113     ret->seed_len = 0;
 114
 115     if (!meth->group_init(ret))
 116         goto err;
 117
 118     return ret;
 119  err:
 120     if (ret->order)
 121         BN_free(ret->order);
 122     if (ret->cofactor)
 123         BN_free(ret->cofactor);
 124     OPENSSL_free(ret);
 125     return NULL;
 126 }
 127
 128 void EC_GROUP_free(EC_GROUP *group)
 129 {
 130     if (!group)
 131         return;
 132
 133     if (group->meth->group_finish != 0)
 134         group->meth->group_finish(group);
 135
 136     EC_EX_DATA_free_all_data(&group->extra_data);
 137
 138     if (group->mont_data)
 139         BN_MONT_CTX_free(group->mont_data);
 140
 141     if (group->generator != NULL)
 142         EC_POINT_free(group->generator);
 143     BN_free(group->order);
 144     BN_free(group->cofactor);
 145
 146     if (group->seed)
 147         OPENSSL_free(group->seed);
 148
 149     OPENSSL_free(group);
 150 }
 151
 152 void EC_GROUP_clear_free(EC_GROUP *group)
 153 {
 154     if (!group)
 155         return;
 156
 157     if (group->meth->group_clear_finish != 0)
 158         group->meth->group_clear_finish(group);
 159     else if (group->meth->group_finish != 0)
 160         group->meth->group_finish(group);
 161
 162     EC_EX_DATA_clear_free_all_data(&group->extra_data);
 163
 164     if (group->mont_data)
 165         BN_MONT_CTX_free(group->mont_data);
 166
 167     if (group->generator != NULL)
 168         EC_POINT_clear_free(group->generator);
 169     BN_clear_free(group->order);
 170     BN_clear_free(group->cofactor);
 171
 172     if (group->seed) {
 173         OPENSSL_cleanse(group->seed, group->seed_len);
 174         OPENSSL_free(group->seed);
 175     }
 176
 177     OPENSSL_cleanse(group, sizeof *group);
 178     OPENSSL_free(group);
 179 }
 180
 181 int EC_GROUP_copy(EC_GROUP *dest, const EC_GROUP *src)
 182 {
 183     EC_EXTRA_DATA *d;
 184
 185     if (dest->meth->group_copy == 0) {
 186         ECerr(EC_F_EC_GROUP_COPY, ERR_R_SHOULD_NOT_HAVE_BEEN_CALLED);
 187         return 0;
 188     }
 189     if (dest->meth != src->meth) {
 190         ECerr(EC_F_EC_GROUP_COPY, EC_R_INCOMPATIBLE_OBJECTS);
 191         return 0;
 192     }
 193     if (dest == src)
 194         return 1;
 195
 196     EC_EX_DATA_free_all_data(&dest->extra_data);
 197
 198     for (d = src->extra_data; d != NULL; d = d->next) {
 199         void *t = d->dup_func(d->data);
 200
 201         if (t == NULL)
 202             return 0;
 203         if (!EC_EX_DATA_set_data
 204             (&dest->extra_data, t, d->dup_func, d->free_func,
 205              d->clear_free_func))
 206             return 0;
 207     }
 208
 209     if (src->mont_data != NULL) {
 210         if (dest->mont_data == NULL) {
 211             dest->mont_data = BN_MONT_CTX_new();
 212             if (dest->mont_data == NULL)
 213                 return 0;
 214         }
 215         if (!BN_MONT_CTX_copy(dest->mont_data, src->mont_data))
 216             return 0;
 217     } else {
 218         /* src->generator == NULL */
 219         if (dest->mont_data != NULL) {
 220             BN_MONT_CTX_free(dest->mont_data);
 221             dest->mont_data = NULL;
 222         }
 223     }
 224
 225     if (src->generator != NULL) {
 226         if (dest->generator == NULL) {
 227             dest->generator = EC_POINT_new(dest);
 228             if (dest->generator == NULL)
 229                 return 0;
 230         }
 231         if (!EC_POINT_copy(dest->generator, src->generator))
 232             return 0;
 233     } else {
 234         /* src->generator == NULL */
 235         if (dest->generator != NULL) {
 236             EC_POINT_clear_free(dest->generator);
 237             dest->generator = NULL;
 238         }
 239     }
 240
 241     if (!BN_copy(dest->order, src->order))
 242         return 0;
 243     if (!BN_copy(dest->cofactor, src->cofactor))
 244         return 0;
 245
 246     dest->curve_name = src->curve_name;
 247     dest->asn1_flag = src->asn1_flag;
 248     dest->asn1_form = src->asn1_form;
 249
 250     if (src->seed) {
 251         if (dest->seed)
 252             OPENSSL_free(dest->seed);
 253         dest->seed = OPENSSL_malloc(src->seed_len);
 254         if (dest->seed == NULL)
 255             return 0;
 256         if (!memcpy(dest->seed, src->seed, src->seed_len))
 257             return 0;
 258         dest->seed_len = src->seed_len;
 259     } else {
 260         if (dest->seed)
 261             OPENSSL_free(dest->seed);
 262         dest->seed = NULL;
 263         dest->seed_len = 0;
 264     }
 265
 266     return dest->meth->group_copy(dest, src);
 267 }
 268
 269 EC_GROUP *EC_GROUP_dup(const EC_GROUP *a)
 270 {
 271     EC_GROUP *t = NULL;
 272     int ok = 0;
 273
 274     if (a == NULL)
 275         return NULL;
 276
 277     if ((t = EC_GROUP_new(a->meth)) == NULL)
 278         return (NULL);
 279     if (!EC_GROUP_copy(t, a))
 280         goto err;
 281
 282     ok = 1;
 283
 284  err:
 285     if (!ok) {
 286         if (t)
 287             EC_GROUP_free(t);
 288         return NULL;
 289     } else
 290         return t;
 291 }
 292
 293 const EC_METHOD *EC_GROUP_method_of(const EC_GROUP *group)
 294 {
 295     return group->meth;
 296 }
 297
 298 int EC_METHOD_get_field_type(const EC_METHOD *meth)
 299 {
 300     return meth->field_type;
 301 }
 302
 303 int EC_GROUP_set_generator(EC_GROUP *group, const EC_POINT *generator,
 304                            const BIGNUM *order, const BIGNUM *cofactor)
 305 {
 306     if (generator == NULL) {
 307         ECerr(EC_F_EC_GROUP_SET_GENERATOR, ERR_R_PASSED_NULL_PARAMETER);
 308         return 0;
 309     }
 310
 311     if (group->generator == NULL) {
 312         group->generator = EC_POINT_new(group);
 313         if (group->generator == NULL)
 314             return 0;
 315     }
 316     if (!EC_POINT_copy(group->generator, generator))
 317         return 0;
 318
 319     if (order != NULL) {
 320         if (!BN_copy(group->order, order))
 321             return 0;
 322     } else
 323         BN_zero(group->order);
 324
 325     if (cofactor != NULL) {
 326         if (!BN_copy(group->cofactor, cofactor))
 327             return 0;
 328     } else
 329         BN_zero(group->cofactor);
 330
 331     /*
 332      * We ignore the return value because some groups have an order with
 333      * factors of two, which makes the Montgomery setup fail.
 334      * |group->mont_data| will be NULL in this case.
 335      */
 336     ec_precompute_mont_data(group);
 337
 338     return 1;
 339 }
 340
 341 const EC_POINT *EC_GROUP_get0_generator(const EC_GROUP *group)
 342 {
 343     return group->generator;
 344 }
 345
 346 BN_MONT_CTX *EC_GROUP_get_mont_data(const EC_GROUP *group)
 347 {
 348     return group->mont_data;
 349 }
 350
 351 int EC_GROUP_get_order(const EC_GROUP *group, BIGNUM *order, BN_CTX *ctx)
 352 {
 353     if (!BN_copy(order, group->order))
 354         return 0;
 355
 356     return !BN_is_zero(order);
 357 }
 358
 359 int EC_GROUP_get_cofactor(const EC_GROUP *group, BIGNUM *cofactor,
 360                           BN_CTX *ctx)
 361 {
 362     if (!BN_copy(cofactor, group->cofactor))
 363         return 0;
 364
 365     return !BN_is_zero(group->cofactor);
 366 }
 367
 368 void EC_GROUP_set_curve_name(EC_GROUP *group, int nid)
 369 {
 370     group->curve_name = nid;
 371 }
 372
 373 int EC_GROUP_get_curve_name(const EC_GROUP *group)
 374 {
 375     return group->curve_name;
 376 }
 377
 378 void EC_GROUP_set_asn1_flag(EC_GROUP *group, int flag)
 379 {
 380     group->asn1_flag = flag;
 381 }
 382
 383 int EC_GROUP_get_asn1_flag(const EC_GROUP *group)
 384 {
 385     return group->asn1_flag;
 386 }
 387
 388 void EC_GROUP_set_point_conversion_form(EC_GROUP *group,
 389                                         point_conversion_form_t form)
 390 {
 391     group->asn1_form = form;
 392 }
 393
 394 point_conversion_form_t EC_GROUP_get_point_conversion_form(const EC_GROUP
 395                                                            *group)
 396 {
 397     return group->asn1_form;
 398 }
 399
 400 size_t EC_GROUP_set_seed(EC_GROUP *group, const unsigned char *p, size_t len)
 401 {
 402     if (group->seed) {
 403         OPENSSL_free(group->seed);
 404         group->seed = NULL;
 405         group->seed_len = 0;
 406     }
 407
 408     if (!len || !p)
 409         return 1;
 410
 411     if ((group->seed = OPENSSL_malloc(len)) == NULL)
 412         return 0;
 413     memcpy(group->seed, p, len);
 414     group->seed_len = len;
 415
 416     return len;
 417 }
 418
 419 unsigned char *EC_GROUP_get0_seed(const EC_GROUP *group)
 420 {
 421     return group->seed;
 422 }
 423
 424 size_t EC_GROUP_get_seed_len(const EC_GROUP *group)
 425 {
 426     return group->seed_len;
 427 }
 428
 429 int EC_GROUP_set_curve_GFp(EC_GROUP *group, const BIGNUM *p, const BIGNUM *a,
 430                            const BIGNUM *b, BN_CTX *ctx)
 431 {
 432     if (group->meth->group_set_curve == 0) {
 433         ECerr(EC_F_EC_GROUP_SET_CURVE_GFP, ERR_R_SHOULD_NOT_HAVE_BEEN_CALLED);
 434         return 0;
 435     }
 436     return group->meth->group_set_curve(group, p, a, b, ctx);
 437 }
 438
 439 int EC_GROUP_get_curve_GFp(const EC_GROUP *group, BIGNUM *p, BIGNUM *a,
 440                            BIGNUM *b, BN_CTX *ctx)
 441 {
 442     if (group->meth->group_get_curve == 0) {
 443         ECerr(EC_F_EC_GROUP_GET_CURVE_GFP, ERR_R_SHOULD_NOT_HAVE_BEEN_CALLED);
 444         return 0;
 445     }
 446     return group->meth->group_get_curve(group, p, a, b, ctx);
 447 }
 448
 449 #ifndef OPENSSL_NO_EC2M
 450 int EC_GROUP_set_curve_GF2m(EC_GROUP *group, const BIGNUM *p, const BIGNUM *a,
 451                             const BIGNUM *b, BN_CTX *ctx)
 452 {
 453     if (group->meth->group_set_curve == 0) {
 454         ECerr(EC_F_EC_GROUP_SET_CURVE_GF2M,
 455               ERR_R_SHOULD_NOT_HAVE_BEEN_CALLED);
 456         return 0;
 457     }
 458     return group->meth->group_set_curve(group, p, a, b, ctx);
 459 }
 460
 461 int EC_GROUP_get_curve_GF2m(const EC_GROUP *group, BIGNUM *p, BIGNUM *a,
 462                             BIGNUM *b, BN_CTX *ctx)
 463 {
 464     if (group->meth->group_get_curve == 0) {
 465         ECerr(EC_F_EC_GROUP_GET_CURVE_GF2M,
 466               ERR_R_SHOULD_NOT_HAVE_BEEN_CALLED);
 467         return 0;
 468     }
 469     return group->meth->group_get_curve(group, p, a, b, ctx);
 470 }
 471 #endif
 472
 473 int EC_GROUP_get_degree(const EC_GROUP *group)
 474 {
 475     if (group->meth->group_get_degree == 0) {
 476         ECerr(EC_F_EC_GROUP_GET_DEGREE, ERR_R_SHOULD_NOT_HAVE_BEEN_CALLED);
 477         return 0;
 478     }
 479     return group->meth->group_get_degree(group);
 480 }
 481
 482 int EC_GROUP_check_discriminant(const EC_GROUP *group, BN_CTX *ctx)
 483 {
 484     if (group->meth->group_check_discriminant == 0) {
 485         ECerr(EC_F_EC_GROUP_CHECK_DISCRIMINANT,
 486               ERR_R_SHOULD_NOT_HAVE_BEEN_CALLED);
 487         return 0;
 488     }
 489     return group->meth->group_check_discriminant(group, ctx);
 490 }
 491
 492 int EC_GROUP_cmp(const EC_GROUP *a, const EC_GROUP *b, BN_CTX *ctx)
 493 {
 494     int r = 0;
 495     BIGNUM *a1, *a2, *a3, *b1, *b2, *b3;
 496     BN_CTX *ctx_new = NULL;
 497
 498     /* compare the field types */
 499     if (EC_METHOD_get_field_type(EC_GROUP_method_of(a)) !=
 500         EC_METHOD_get_field_type(EC_GROUP_method_of(b)))
 501         return 1;
 502     /* compare the curve name (if present in both) */
 503     if (EC_GROUP_get_curve_name(a) && EC_GROUP_get_curve_name(b) &&
 504         EC_GROUP_get_curve_name(a) != EC_GROUP_get_curve_name(b))
 505         return 1;
 506
 507     if (!ctx)
 508         ctx_new = ctx = BN_CTX_new();
 509     if (!ctx)
 510         return -1;
 511
 512     BN_CTX_start(ctx);
 513     a1 = BN_CTX_get(ctx);
 514     a2 = BN_CTX_get(ctx);
 515     a3 = BN_CTX_get(ctx);
 516     b1 = BN_CTX_get(ctx);
 517     b2 = BN_CTX_get(ctx);
 518     b3 = BN_CTX_get(ctx);
 519     if (!b3) {
 520         BN_CTX_end(ctx);
 521         if (ctx_new)
 522             BN_CTX_free(ctx);
 523         return -1;
 524     }
 525
 526     /*
 527      * XXX This approach assumes that the external representation of curves
 528      * over the same field type is the same.
 529      */
 530     if (!a->meth->group_get_curve(a, a1, a2, a3, ctx) ||
 531         !b->meth->group_get_curve(b, b1, b2, b3, ctx))
 532         r = 1;
 533
 534     if (r || BN_cmp(a1, b1) || BN_cmp(a2, b2) || BN_cmp(a3, b3))
 535         r = 1;
 536
 537     /* XXX EC_POINT_cmp() assumes that the methods are equal */
 538     if (r || EC_POINT_cmp(a, EC_GROUP_get0_generator(a),
 539                           EC_GROUP_get0_generator(b), ctx))
 540         r = 1;
 541
 542     if (!r) {
 543         /* compare the order and cofactor */
 544         if (!EC_GROUP_get_order(a, a1, ctx) ||
 545             !EC_GROUP_get_order(b, b1, ctx) ||
 546             !EC_GROUP_get_cofactor(a, a2, ctx) ||
 547             !EC_GROUP_get_cofactor(b, b2, ctx)) {
 548             BN_CTX_end(ctx);
 549             if (ctx_new)
 550                 BN_CTX_free(ctx);
 551             return -1;
 552         }
 553         if (BN_cmp(a1, b1) || BN_cmp(a2, b2))
 554             r = 1;
 555     }
 556
 557     BN_CTX_end(ctx);
 558     if (ctx_new)
 559         BN_CTX_free(ctx);
 560
 561     return r;
 562 }
 563
 564 /* this has 'package' visibility */
 565 int EC_EX_DATA_set_data(EC_EXTRA_DATA **ex_data, void *data,
 566                         void *(*dup_func) (void *),
 567                         void (*free_func) (void *),
 568                         void (*clear_free_func) (void *))
 569 {
 570     EC_EXTRA_DATA *d;
 571
 572     if (ex_data == NULL)
 573         return 0;
 574
 575     for (d = *ex_data; d != NULL; d = d->next) {
 576         if (d->dup_func == dup_func && d->free_func == free_func
 577             && d->clear_free_func == clear_free_func) {
 578             ECerr(EC_F_EC_EX_DATA_SET_DATA, EC_R_SLOT_FULL);
 579             return 0;
 580         }
 581     }
 582
 583     if (data == NULL)
 584         /* no explicit entry needed */
 585         return 1;
 586
 587     d = OPENSSL_malloc(sizeof *d);
 588     if (d == NULL)
 589         return 0;
 590
 591     d->data = data;
 592     d->dup_func = dup_func;
 593     d->free_func = free_func;
 594     d->clear_free_func = clear_free_func;
 595
 596     d->next = *ex_data;
 597     *ex_data = d;
 598
 599     return 1;
 600 }
 601
 602 /* this has 'package' visibility */
 603 void *EC_EX_DATA_get_data(const EC_EXTRA_DATA *ex_data,
 604                           void *(*dup_func) (void *),
 605                           void (*free_func) (void *),
 606                           void (*clear_free_func) (void *))
 607 {
 608     const EC_EXTRA_DATA *d;
 609
 610     for (d = ex_data; d != NULL; d = d->next) {
 611         if (d->dup_func == dup_func && d->free_func == free_func
 612             && d->clear_free_func == clear_free_func)
 613             return d->data;
 614     }
 615
 616     return NULL;
 617 }
 618
 619 /* this has 'package' visibility */
 620 void EC_EX_DATA_free_data(EC_EXTRA_DATA **ex_data,
 621                           void *(*dup_func) (void *),
 622                           void (*free_func) (void *),
 623                           void (*clear_free_func) (void *))
 624 {
 625     EC_EXTRA_DATA **p;
 626
 627     if (ex_data == NULL)
 628         return;
 629
 630     for (p = ex_data; *p != NULL; p = &((*p)->next)) {
 631         if ((*p)->dup_func == dup_func && (*p)->free_func == free_func
 632             && (*p)->clear_free_func == clear_free_func) {
 633             EC_EXTRA_DATA *next = (*p)->next;
 634
 635             (*p)->free_func((*p)->data);
 636             OPENSSL_free(*p);
 637
 638             *p = next;
 639             return;
 640         }
 641     }
 642 }
 643
 644 /* this has 'package' visibility */
 645 void EC_EX_DATA_clear_free_data(EC_EXTRA_DATA **ex_data,
 646                                 void *(*dup_func) (void *),
 647                                 void (*free_func) (void *),
 648                                 void (*clear_free_func) (void *))
 649 {
 650     EC_EXTRA_DATA **p;
 651
 652     if (ex_data == NULL)
 653         return;
 654
 655     for (p = ex_data; *p != NULL; p = &((*p)->next)) {
 656         if ((*p)->dup_func == dup_func && (*p)->free_func == free_func
 657             && (*p)->clear_free_func == clear_free_func) {
 658             EC_EXTRA_DATA *next = (*p)->next;
 659
 660             (*p)->clear_free_func((*p)->data);
 661             OPENSSL_free(*p);
 662
 663             *p = next;
 664             return;
 665         }
 666     }
 667 }
 668
 669 /* this has 'package' visibility */
 670 void EC_EX_DATA_free_all_data(EC_EXTRA_DATA **ex_data)
 671 {
 672     EC_EXTRA_DATA *d;
 673
 674     if (ex_data == NULL)
 675         return;
 676
 677     d = *ex_data;
 678     while (d) {
 679         EC_EXTRA_DATA *next = d->next;
 680
 681         d->free_func(d->data);
 682         OPENSSL_free(d);
 683
 684         d = next;
 685     }
 686     *ex_data = NULL;
 687 }
 688
 689 /* this has 'package' visibility */
 690 void EC_EX_DATA_clear_free_all_data(EC_EXTRA_DATA **ex_data)
 691 {
 692     EC_EXTRA_DATA *d;
 693
 694     if (ex_data == NULL)
 695         return;
 696
 697     d = *ex_data;
 698     while (d) {
 699         EC_EXTRA_DATA *next = d->next;
 700
 701         d->clear_free_func(d->data);
 702         OPENSSL_free(d);
 703
 704         d = next;
 705     }
 706     *ex_data = NULL;
 707 }
 708
 709 /* functions for EC_POINT objects */
 710
 711 EC_POINT *EC_POINT_new(const EC_GROUP *group)
 712 {
 713     EC_POINT *ret;
 714
 715     if (group == NULL) {
 716         ECerr(EC_F_EC_POINT_NEW, ERR_R_PASSED_NULL_PARAMETER);
 717         return NULL;
 718     }
 719     if (group->meth->point_init == 0) {
 720         ECerr(EC_F_EC_POINT_NEW, ERR_R_SHOULD_NOT_HAVE_BEEN_CALLED);
 721         return NULL;
 722     }
 723
 724     ret = OPENSSL_malloc(sizeof *ret);
 725     if (ret == NULL) {
 726         ECerr(EC_F_EC_POINT_NEW, ERR_R_MALLOC_FAILURE);
 727         return NULL;
 728     }
 729
 730     ret->meth = group->meth;
 731
 732     if (!ret->meth->point_init(ret)) {
 733         OPENSSL_free(ret);
 734         return NULL;
 735     }
 736
 737     return ret;
 738 }
 739
 740 void EC_POINT_free(EC_POINT *point)
 741 {
 742     if (!point)
 743         return;
 744
 745     if (point->meth->point_finish != 0)
 746         point->meth->point_finish(point);
 747     OPENSSL_free(point);
 748 }
 749
 750 void EC_POINT_clear_free(EC_POINT *point)
 751 {
 752     if (!point)
 753         return;
 754
 755     if (point->meth->point_clear_finish != 0)
 756         point->meth->point_clear_finish(point);
 757     else if (point->meth->point_finish != 0)
 758         point->meth->point_finish(point);
 759     OPENSSL_cleanse(point, sizeof *point);
 760     OPENSSL_free(point);
 761 }
 762
 763 int EC_POINT_copy(EC_POINT *dest, const EC_POINT *src)
 764 {
 765     if (dest->meth->point_copy == 0) {
 766         ECerr(EC_F_EC_POINT_COPY, ERR_R_SHOULD_NOT_HAVE_BEEN_CALLED);
 767         return 0;
 768     }
 769     if (dest->meth != src->meth) {
 770         ECerr(EC_F_EC_POINT_COPY, EC_R_INCOMPATIBLE_OBJECTS);
 771         return 0;
 772     }
 773     if (dest == src)
 774         return 1;
 775     return dest->meth->point_copy(dest, src);
 776 }
 777
 778 EC_POINT *EC_POINT_dup(const EC_POINT *a, const EC_GROUP *group)
 779 {
 780     EC_POINT *t;
 781     int r;
 782
 783     if (a == NULL)
 784         return NULL;
 785
 786     t = EC_POINT_new(group);
 787     if (t == NULL)
 788         return (NULL);
 789     r = EC_POINT_copy(t, a);
 790     if (!r) {
 791         EC_POINT_free(t);
 792         return NULL;
 793     } else
 794         return t;
 795 }
 796
 797 const EC_METHOD *EC_POINT_method_of(const EC_POINT *point)
 798 {
 799     return point->meth;
 800 }
 801
 802 int EC_POINT_set_to_infinity(const EC_GROUP *group, EC_POINT *point)
 803 {
 804     if (group->meth->point_set_to_infinity == 0) {
 805         ECerr(EC_F_EC_POINT_SET_TO_INFINITY,
 806               ERR_R_SHOULD_NOT_HAVE_BEEN_CALLED);
 807         return 0;
 808     }
 809     if (group->meth != point->meth) {
 810         ECerr(EC_F_EC_POINT_SET_TO_INFINITY, EC_R_INCOMPATIBLE_OBJECTS);
 811         return 0;
 812     }
 813     return group->meth->point_set_to_infinity(group, point);
 814 }
 815
 816 int EC_POINT_set_Jprojective_coordinates_GFp(const EC_GROUP *group,
 817                                              EC_POINT *point, const BIGNUM *x,
 818                                              const BIGNUM *y, const BIGNUM *z,
 819                                              BN_CTX *ctx)
 820 {
 821     if (group->meth->point_set_Jprojective_coordinates_GFp == 0) {
 822         ECerr(EC_F_EC_POINT_SET_JPROJECTIVE_COORDINATES_GFP,
 823               ERR_R_SHOULD_NOT_HAVE_BEEN_CALLED);
 824         return 0;
 825     }
 826     if (group->meth != point->meth) {
 827         ECerr(EC_F_EC_POINT_SET_JPROJECTIVE_COORDINATES_GFP,
 828               EC_R_INCOMPATIBLE_OBJECTS);
 829         return 0;
 830     }
 831     return group->meth->point_set_Jprojective_coordinates_GFp(group, point, x,
 832                                                               y, z, ctx);
 833 }
 834
 835 int EC_POINT_get_Jprojective_coordinates_GFp(const EC_GROUP *group,
 836                                              const EC_POINT *point, BIGNUM *x,
 837                                              BIGNUM *y, BIGNUM *z,
 838                                              BN_CTX *ctx)
 839 {
 840     if (group->meth->point_get_Jprojective_coordinates_GFp == 0) {
 841         ECerr(EC_F_EC_POINT_GET_JPROJECTIVE_COORDINATES_GFP,
 842               ERR_R_SHOULD_NOT_HAVE_BEEN_CALLED);
 843         return 0;
 844     }
 845     if (group->meth != point->meth) {
 846         ECerr(EC_F_EC_POINT_GET_JPROJECTIVE_COORDINATES_GFP,
 847               EC_R_INCOMPATIBLE_OBJECTS);
 848         return 0;
 849     }
 850     return group->meth->point_get_Jprojective_coordinates_GFp(group, point, x,
 851                                                               y, z, ctx);
 852 }
 853
 854 int EC_POINT_set_affine_coordinates_GFp(const EC_GROUP *group,
 855                                         EC_POINT *point, const BIGNUM *x,
 856                                         const BIGNUM *y, BN_CTX *ctx)
 857 {
 858     if (group->meth->point_set_affine_coordinates == 0) {
 859         ECerr(EC_F_EC_POINT_SET_AFFINE_COORDINATES_GFP,
 860               ERR_R_SHOULD_NOT_HAVE_BEEN_CALLED);
 861         return 0;
 862     }
 863     if (group->meth != point->meth) {
 864         ECerr(EC_F_EC_POINT_SET_AFFINE_COORDINATES_GFP,
 865               EC_R_INCOMPATIBLE_OBJECTS);
 866         return 0;
 867     }
 868     return group->meth->point_set_affine_coordinates(group, point, x, y, ctx);
 869 }
 870
 871 #ifndef OPENSSL_NO_EC2M
 872 int EC_POINT_set_affine_coordinates_GF2m(const EC_GROUP *group,
 873                                          EC_POINT *point, const BIGNUM *x,
 874                                          const BIGNUM *y, BN_CTX *ctx)
 875 {
 876     if (group->meth->point_set_affine_coordinates == 0) {
 877         ECerr(EC_F_EC_POINT_SET_AFFINE_COORDINATES_GF2M,
 878               ERR_R_SHOULD_NOT_HAVE_BEEN_CALLED);
 879         return 0;
 880     }
 881     if (group->meth != point->meth) {
 882         ECerr(EC_F_EC_POINT_SET_AFFINE_COORDINATES_GF2M,
 883               EC_R_INCOMPATIBLE_OBJECTS);
 884         return 0;
 885     }
 886     return group->meth->point_set_affine_coordinates(group, point, x, y, ctx);
 887 }
 888 #endif
 889
 890 int EC_POINT_get_affine_coordinates_GFp(const EC_GROUP *group,
 891                                         const EC_POINT *point, BIGNUM *x,
 892                                         BIGNUM *y, BN_CTX *ctx)
 893 {
 894     if (group->meth->point_get_affine_coordinates == 0) {
 895         ECerr(EC_F_EC_POINT_GET_AFFINE_COORDINATES_GFP,
 896               ERR_R_SHOULD_NOT_HAVE_BEEN_CALLED);
 897         return 0;
 898     }
 899     if (group->meth != point->meth) {
 900         ECerr(EC_F_EC_POINT_GET_AFFINE_COORDINATES_GFP,
 901               EC_R_INCOMPATIBLE_OBJECTS);
 902         return 0;
 903     }
 904     return group->meth->point_get_affine_coordinates(group, point, x, y, ctx);
 905 }
 906
 907 #ifndef OPENSSL_NO_EC2M
 908 int EC_POINT_get_affine_coordinates_GF2m(const EC_GROUP *group,
 909                                          const EC_POINT *point, BIGNUM *x,
 910                                          BIGNUM *y, BN_CTX *ctx)
 911 {
 912     if (group->meth->point_get_affine_coordinates == 0) {
 913         ECerr(EC_F_EC_POINT_GET_AFFINE_COORDINATES_GF2M,
 914               ERR_R_SHOULD_NOT_HAVE_BEEN_CALLED);
 915         return 0;
 916     }
 917     if (group->meth != point->meth) {
 918         ECerr(EC_F_EC_POINT_GET_AFFINE_COORDINATES_GF2M,
 919               EC_R_INCOMPATIBLE_OBJECTS);
 920         return 0;
 921     }
 922     return group->meth->point_get_affine_coordinates(group, point, x, y, ctx);
 923 }
 924 #endif
 925
 926 int EC_POINT_add(const EC_GROUP *group, EC_POINT *r, const EC_POINT *a,
 927                  const EC_POINT *b, BN_CTX *ctx)
 928 {
 929     if (group->meth->add == 0) {
 930         ECerr(EC_F_EC_POINT_ADD, ERR_R_SHOULD_NOT_HAVE_BEEN_CALLED);
 931         return 0;
 932     }
 933     if ((group->meth != r->meth) || (r->meth != a->meth)
 934         || (a->meth != b->meth)) {
 935         ECerr(EC_F_EC_POINT_ADD, EC_R_INCOMPATIBLE_OBJECTS);
 936         return 0;
 937     }
 938     return group->meth->add(group, r, a, b, ctx);
 939 }
 940
 941 int EC_POINT_dbl(const EC_GROUP *group, EC_POINT *r, const EC_POINT *a,
 942                  BN_CTX *ctx)
 943 {
 944     if (group->meth->dbl == 0) {
 945         ECerr(EC_F_EC_POINT_DBL, ERR_R_SHOULD_NOT_HAVE_BEEN_CALLED);
 946         return 0;
 947     }
 948     if ((group->meth != r->meth) || (r->meth != a->meth)) {
 949         ECerr(EC_F_EC_POINT_DBL, EC_R_INCOMPATIBLE_OBJECTS);
 950         return 0;
 951     }
 952     return group->meth->dbl(group, r, a, ctx);
 953 }
 954
 955 int EC_POINT_invert(const EC_GROUP *group, EC_POINT *a, BN_CTX *ctx)
 956 {
 957     if (group->meth->invert == 0) {
 958         ECerr(EC_F_EC_POINT_INVERT, ERR_R_SHOULD_NOT_HAVE_BEEN_CALLED);
 959         return 0;
 960     }
 961     if (group->meth != a->meth) {
 962         ECerr(EC_F_EC_POINT_INVERT, EC_R_INCOMPATIBLE_OBJECTS);
 963         return 0;
 964     }
 965     return group->meth->invert(group, a, ctx);
 966 }
 967
 968 int EC_POINT_is_at_infinity(const EC_GROUP *group, const EC_POINT *point)
 969 {
 970     if (group->meth->is_at_infinity == 0) {
 971         ECerr(EC_F_EC_POINT_IS_AT_INFINITY,
 972               ERR_R_SHOULD_NOT_HAVE_BEEN_CALLED);
 973         return 0;
 974     }
 975     if (group->meth != point->meth) {
 976         ECerr(EC_F_EC_POINT_IS_AT_INFINITY, EC_R_INCOMPATIBLE_OBJECTS);
 977         return 0;
 978     }
 979     return group->meth->is_at_infinity(group, point);
 980 }
 981
 982 int EC_POINT_is_on_curve(const EC_GROUP *group, const EC_POINT *point,
 983                          BN_CTX *ctx)
 984 {
 985     if (group->meth->is_on_curve == 0) {
 986         ECerr(EC_F_EC_POINT_IS_ON_CURVE, ERR_R_SHOULD_NOT_HAVE_BEEN_CALLED);
 987         return 0;
 988     }
 989     if (group->meth != point->meth) {
 990         ECerr(EC_F_EC_POINT_IS_ON_CURVE, EC_R_INCOMPATIBLE_OBJECTS);
 991         return 0;
 992     }
 993     return group->meth->is_on_curve(group, point, ctx);
 994 }
 995
 996 int EC_POINT_cmp(const EC_GROUP *group, const EC_POINT *a, const EC_POINT *b,
 997                  BN_CTX *ctx)
 998 {
 999     if (group->meth->point_cmp == 0) {
1000         ECerr(EC_F_EC_POINT_CMP, ERR_R_SHOULD_NOT_HAVE_BEEN_CALLED);
1001         return -1;
1002     }
1003     if ((group->meth != a->meth) || (a->meth != b->meth)) {
1004         ECerr(EC_F_EC_POINT_CMP, EC_R_INCOMPATIBLE_OBJECTS);
1005         return -1;
1006     }
1007     return group->meth->point_cmp(group, a, b, ctx);
1008 }
1009
1010 int EC_POINT_make_affine(const EC_GROUP *group, EC_POINT *point, BN_CTX *ctx)
1011 {
1012     if (group->meth->make_affine == 0) {
1013         ECerr(EC_F_EC_POINT_MAKE_AFFINE, ERR_R_SHOULD_NOT_HAVE_BEEN_CALLED);
1014         return 0;
1015     }
1016     if (group->meth != point->meth) {
1017         ECerr(EC_F_EC_POINT_MAKE_AFFINE, EC_R_INCOMPATIBLE_OBJECTS);
1018         return 0;
1019     }
1020     return group->meth->make_affine(group, point, ctx);
1021 }
1022
1023 int EC_POINTs_make_affine(const EC_GROUP *group, size_t num,
1024                           EC_POINT *points[], BN_CTX *ctx)
1025 {
1026     size_t i;
1027
1028     if (group->meth->points_make_affine == 0) {
1029         ECerr(EC_F_EC_POINTS_MAKE_AFFINE, ERR_R_SHOULD_NOT_HAVE_BEEN_CALLED);
1030         return 0;
1031     }
1032     for (i = 0; i < num; i++) {
1033         if (group->meth != points[i]->meth) {
1034             ECerr(EC_F_EC_POINTS_MAKE_AFFINE, EC_R_INCOMPATIBLE_OBJECTS);
1035             return 0;
1036         }
1037     }
1038     return group->meth->points_make_affine(group, num, points, ctx);
1039 }
1040
1041 /*
1042  * Functions for point multiplication. If group->meth->mul is 0, we use the
1043  * wNAF-based implementations in ec_mult.c; otherwise we dispatch through
1044  * methods.
1045  */
1046
1047 int EC_POINTs_mul(const EC_GROUP *group, EC_POINT *r, const BIGNUM *scalar,
1048                   size_t num, const EC_POINT *points[],
1049                   const BIGNUM *scalars[], BN_CTX *ctx)
1050 {
1051     if (group->meth->mul == 0)
1052         /* use default */
1053         return ec_wNAF_mul(group, r, scalar, num, points, scalars, ctx);
1054
1055     return group->meth->mul(group, r, scalar, num, points, scalars, ctx);
1056 }
1057
1058 int EC_POINT_mul(const EC_GROUP *group, EC_POINT *r, const BIGNUM *g_scalar,
1059                  const EC_POINT *point, const BIGNUM *p_scalar, BN_CTX *ctx)
1060 {
1061     /* just a convenient interface to EC_POINTs_mul() */
1062
1063     const EC_POINT *points[1];
1064     const BIGNUM *scalars[1];
1065
1066     points[0] = point;
1067     scalars[0] = p_scalar;
1068
1069     return EC_POINTs_mul(group, r, g_scalar,
1070                          (point != NULL
1071                           && p_scalar != NULL), points, scalars, ctx);
1072 }
1073
1074 int EC_GROUP_precompute_mult(EC_GROUP *group, BN_CTX *ctx)
1075 {
1076     if (group->meth->mul == 0)
1077         /* use default */
1078         return ec_wNAF_precompute_mult(group, ctx);
1079
1080     if (group->meth->precompute_mult != 0)
1081         return group->meth->precompute_mult(group, ctx);
1082     else
1083         return 1;               /* nothing to do, so report success */
1084 }
1085
1086 int EC_GROUP_have_precompute_mult(const EC_GROUP *group)
1087 {
1088     if (group->meth->mul == 0)
1089         /* use default */
1090         return ec_wNAF_have_precompute_mult(group);
1091
1092     if (group->meth->have_precompute_mult != 0)
1093         return group->meth->have_precompute_mult(group);
1094     else
1095         return 0;               /* cannot tell whether precomputation has
1096                                  * been performed */
1097 }
1098
1099 /*
1100  * ec_precompute_mont_data sets |group->mont_data| from |group->order| and
1101  * returns one on success. On error it returns zero.
1102  */
1103 int ec_precompute_mont_data(EC_GROUP *group)
1104 {
1105     BN_CTX *ctx = BN_CTX_new();
1106     int ret = 0;
1107
1108     if (group->mont_data) {
1109         BN_MONT_CTX_free(group->mont_data);
1110         group->mont_data = NULL;
1111     }
1112
1113     if (ctx == NULL)
1114         goto err;
1115
1116     group->mont_data = BN_MONT_CTX_new();
1117     if (!group->mont_data)
1118         goto err;
1119
1120     if (!BN_MONT_CTX_set(group->mont_data, group->order, ctx)) {
1121         BN_MONT_CTX_free(group->mont_data);
1122         group->mont_data = NULL;
1123         goto err;
1124     }
1125
1126     ret = 1;
1127
1128  err:
1129
1130     if (ctx)
1131         BN_CTX_free(ctx);
1132     return ret;
1133 }