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