crypto/ecdh/ecdhtest.c - third_party/openssl - Git at Google

 /* crypto/ecdh/ecdhtest.c */
 /* ====================================================================
  * Copyright 2002 Sun Microsystems, Inc. ALL RIGHTS RESERVED.
  *
  * The Elliptic Curve Public-Key Crypto Library (ECC Code) included
  * herein is developed by SUN MICROSYSTEMS, INC., and is contributed
  * to the OpenSSL project.
  *
  * The ECC Code is licensed pursuant to the OpenSSL open source
  * license provided below.
  *
  * The ECDH software is originally written by Douglas Stebila of
  * Sun Microsystems Laboratories.
  *
  */
 /* ====================================================================
  * Copyright (c) 1998-2003 The OpenSSL Project.  All rights reserved.
  *
  * Redistribution and use in source and binary forms, with or without
  * modification, are permitted provided that the following conditions
  * are met:
  *
  * 1. Redistributions of source code must retain the above copyright
  *    notice, this list of conditions and the following disclaimer.
  *
  * 2. Redistributions in binary form must reproduce the above copyright
  *    notice, this list of conditions and the following disclaimer in
  *    the documentation and/or other materials provided with the
  *    distribution.
  *
  * 3. All advertising materials mentioning features or use of this
  *    software must display the following acknowledgment:
  *    "This product includes software developed by the OpenSSL Project
  *    for use in the OpenSSL Toolkit. (http://www.openssl.org/)"
  *
  * 4. The names "OpenSSL Toolkit" and "OpenSSL Project" must not be used to
  *    endorse or promote products derived from this software without
  *    prior written permission. For written permission, please contact
  *    openssl-core@openssl.org.
  *
  * 5. Products derived from this software may not be called "OpenSSL"
  *    nor may "OpenSSL" appear in their names without prior written
  *    permission of the OpenSSL Project.
  *
  * 6. Redistributions of any form whatsoever must retain the following
  *    acknowledgment:
  *    "This product includes software developed by the OpenSSL Project
  *    for use in the OpenSSL Toolkit (http://www.openssl.org/)"
  *
  * THIS SOFTWARE IS PROVIDED BY THE OpenSSL PROJECT ``AS IS'' AND ANY
  * EXPRESSED OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
  * IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR
  * PURPOSE ARE DISCLAIMED.  IN NO EVENT SHALL THE OpenSSL PROJECT OR
  * ITS CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL,
  * SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT
  * NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES;
  * LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION)
  * HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT,
  * STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE)
  * ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED
  * OF THE POSSIBILITY OF SUCH DAMAGE.
  * ====================================================================
  *
  * This product includes cryptographic software written by Eric Young
  * (eay@cryptsoft.com).  This product includes software written by Tim
  * Hudson (tjh@cryptsoft.com).
  *
  */

 #include <stdio.h>
 #include <stdlib.h>
 #include <string.h>

 #include "../e_os.h"

 #include <openssl/opensslconf.h> /* for OPENSSL_NO_ECDH */
 #include <openssl/crypto.h>
 #include <openssl/bio.h>
 #include <openssl/bn.h>
 #include <openssl/objects.h>
 #include <openssl/rand.h>
 #include <openssl/sha.h>
 #include <openssl/err.h>

 #ifdef OPENSSL_NO_ECDH
 int main(int argc, char *argv[])
 {
     printf("No ECDH support\n");
     return (0);
 }
 #else
 # include <openssl/ec.h>
 # include <openssl/ecdh.h>

 static const char rnd_seed[] =
     "string to make the random number generator think it has entropy";

 static const int KDF1_SHA1_len = 20;
 static void *KDF1_SHA1(const void *in, size_t inlen, void *out,
                        size_t *outlen)
 {
     if (*outlen < SHA_DIGEST_LENGTH)
         return NULL;
     *outlen = SHA_DIGEST_LENGTH;
     return SHA1(in, inlen, out);
 }

 static int test_ecdh_curve(int nid, const char *text, BN_CTX *ctx, BIO *out)
 {
     EC_KEY *a = NULL;
     EC_KEY *b = NULL;
     BIGNUM *x_a = NULL, *y_a = NULL, *x_b = NULL, *y_b = NULL;
     char buf[12];
     unsigned char *abuf = NULL, *bbuf = NULL;
     int i, alen, blen, aout, bout, ret = 0;
     const EC_GROUP *group;

     a = EC_KEY_new_by_curve_name(nid);
     b = EC_KEY_new_by_curve_name(nid);
     if (a == NULL || b == NULL)
         goto err;

     group = EC_KEY_get0_group(a);

     if ((x_a = BN_new()) == NULL)
         goto err;
     if ((y_a = BN_new()) == NULL)
         goto err;
     if ((x_b = BN_new()) == NULL)
         goto err;
     if ((y_b = BN_new()) == NULL)
         goto err;

     BIO_puts(out, "Testing key generation with ");
     BIO_puts(out, text);
 # ifdef NOISY
     BIO_puts(out, "\n");
 # else
     (void)BIO_flush(out);
 # endif

     if (!EC_KEY_generate_key(a))
         goto err;

     if (EC_METHOD_get_field_type(EC_GROUP_method_of(group)) ==
         NID_X9_62_prime_field) {
         if (!EC_POINT_get_affine_coordinates_GFp
             (group, EC_KEY_get0_public_key(a), x_a, y_a, ctx))
             goto err;
     }
 # ifndef OPENSSL_NO_EC2M
     else {
         if (!EC_POINT_get_affine_coordinates_GF2m(group,
                                                   EC_KEY_get0_public_key(a),
                                                   x_a, y_a, ctx))
             goto err;
     }
 # endif
 # ifdef NOISY
     BIO_puts(out, "  pri 1=");
     BN_print(out, a->priv_key);
     BIO_puts(out, "\n  pub 1=");
     BN_print(out, x_a);
     BIO_puts(out, ",");
     BN_print(out, y_a);
     BIO_puts(out, "\n");
 # else
     BIO_printf(out, " .");
     (void)BIO_flush(out);
 # endif

     if (!EC_KEY_generate_key(b))
         goto err;

     if (EC_METHOD_get_field_type(EC_GROUP_method_of(group)) ==
         NID_X9_62_prime_field) {
         if (!EC_POINT_get_affine_coordinates_GFp
             (group, EC_KEY_get0_public_key(b), x_b, y_b, ctx))
             goto err;
     }
 # ifndef OPENSSL_NO_EC2M
     else {
         if (!EC_POINT_get_affine_coordinates_GF2m(group,
                                                   EC_KEY_get0_public_key(b),
                                                   x_b, y_b, ctx))
             goto err;
     }
 # endif

 # ifdef NOISY
     BIO_puts(out, "  pri 2=");
     BN_print(out, b->priv_key);
     BIO_puts(out, "\n  pub 2=");
     BN_print(out, x_b);
     BIO_puts(out, ",");
     BN_print(out, y_b);
     BIO_puts(out, "\n");
 # else
     BIO_printf(out, ".");
     (void)BIO_flush(out);
 # endif

     alen = KDF1_SHA1_len;
     abuf = (unsigned char *)OPENSSL_malloc(alen);
     aout =
         ECDH_compute_key(abuf, alen, EC_KEY_get0_public_key(b), a, KDF1_SHA1);

 # ifdef NOISY
     BIO_puts(out, "  key1 =");
     for (i = 0; i < aout; i++) {
         sprintf(buf, "%02X", abuf[i]);
         BIO_puts(out, buf);
     }
     BIO_puts(out, "\n");
 # else
     BIO_printf(out, ".");
     (void)BIO_flush(out);
 # endif

     blen = KDF1_SHA1_len;
     bbuf = (unsigned char *)OPENSSL_malloc(blen);
     bout =
         ECDH_compute_key(bbuf, blen, EC_KEY_get0_public_key(a), b, KDF1_SHA1);

 # ifdef NOISY
     BIO_puts(out, "  key2 =");
     for (i = 0; i < bout; i++) {
         sprintf(buf, "%02X", bbuf[i]);
         BIO_puts(out, buf);
     }
     BIO_puts(out, "\n");
 # else
     BIO_printf(out, ".");
     (void)BIO_flush(out);
 # endif

     if ((aout < 4) || (bout != aout) || (memcmp(abuf, bbuf, aout) != 0)) {
 # ifndef NOISY
         BIO_printf(out, " failed\n\n");
         BIO_printf(out, "key a:\n");
         BIO_printf(out, "private key: ");
         BN_print(out, EC_KEY_get0_private_key(a));
         BIO_printf(out, "\n");
         BIO_printf(out, "public key (x,y): ");
         BN_print(out, x_a);
         BIO_printf(out, ",");
         BN_print(out, y_a);
         BIO_printf(out, "\nkey b:\n");
         BIO_printf(out, "private key: ");
         BN_print(out, EC_KEY_get0_private_key(b));
         BIO_printf(out, "\n");
         BIO_printf(out, "public key (x,y): ");
         BN_print(out, x_b);
         BIO_printf(out, ",");
         BN_print(out, y_b);
         BIO_printf(out, "\n");
         BIO_printf(out, "generated key a: ");
         for (i = 0; i < bout; i++) {
             sprintf(buf, "%02X", bbuf[i]);
             BIO_puts(out, buf);
         }
         BIO_printf(out, "\n");
         BIO_printf(out, "generated key b: ");
         for (i = 0; i < aout; i++) {
             sprintf(buf, "%02X", abuf[i]);
             BIO_puts(out, buf);
         }
         BIO_printf(out, "\n");
 # endif
         fprintf(stderr, "Error in ECDH routines\n");
         ret = 0;
     } else {
 # ifndef NOISY
         BIO_printf(out, " ok\n");
 # endif
         ret = 1;
     }
  err:
     ERR_print_errors_fp(stderr);

     if (abuf != NULL)
         OPENSSL_free(abuf);
     if (bbuf != NULL)
         OPENSSL_free(bbuf);
     if (x_a)
         BN_free(x_a);
     if (y_a)
         BN_free(y_a);
     if (x_b)
         BN_free(x_b);
     if (y_b)
         BN_free(y_b);
     if (b)
         EC_KEY_free(b);
     if (a)
         EC_KEY_free(a);
     return (ret);
 }

 /* Keys and shared secrets from RFC 7027 */

 static const unsigned char bp256_da[] = {
     0x81, 0xDB, 0x1E, 0xE1, 0x00, 0x15, 0x0F, 0xF2, 0xEA, 0x33, 0x8D, 0x70,
     0x82, 0x71, 0xBE, 0x38, 0x30, 0x0C, 0xB5, 0x42, 0x41, 0xD7, 0x99, 0x50,
     0xF7, 0x7B, 0x06, 0x30, 0x39, 0x80, 0x4F, 0x1D
 };

 static const unsigned char bp256_db[] = {
     0x55, 0xE4, 0x0B, 0xC4, 0x1E, 0x37, 0xE3, 0xE2, 0xAD, 0x25, 0xC3, 0xC6,
     0x65, 0x45, 0x11, 0xFF, 0xA8, 0x47, 0x4A, 0x91, 0xA0, 0x03, 0x20, 0x87,
     0x59, 0x38, 0x52, 0xD3, 0xE7, 0xD7, 0x6B, 0xD3
 };

 static const unsigned char bp256_Z[] = {
     0x89, 0xAF, 0xC3, 0x9D, 0x41, 0xD3, 0xB3, 0x27, 0x81, 0x4B, 0x80, 0x94,
     0x0B, 0x04, 0x25, 0x90, 0xF9, 0x65, 0x56, 0xEC, 0x91, 0xE6, 0xAE, 0x79,
     0x39, 0xBC, 0xE3, 0x1F, 0x3A, 0x18, 0xBF, 0x2B
 };

 static const unsigned char bp384_da[] = {
     0x1E, 0x20, 0xF5, 0xE0, 0x48, 0xA5, 0x88, 0x6F, 0x1F, 0x15, 0x7C, 0x74,
     0xE9, 0x1B, 0xDE, 0x2B, 0x98, 0xC8, 0xB5, 0x2D, 0x58, 0xE5, 0x00, 0x3D,
     0x57, 0x05, 0x3F, 0xC4, 0xB0, 0xBD, 0x65, 0xD6, 0xF1, 0x5E, 0xB5, 0xD1,
     0xEE, 0x16, 0x10, 0xDF, 0x87, 0x07, 0x95, 0x14, 0x36, 0x27, 0xD0, 0x42
 };

 static const unsigned char bp384_db[] = {
     0x03, 0x26, 0x40, 0xBC, 0x60, 0x03, 0xC5, 0x92, 0x60, 0xF7, 0x25, 0x0C,
     0x3D, 0xB5, 0x8C, 0xE6, 0x47, 0xF9, 0x8E, 0x12, 0x60, 0xAC, 0xCE, 0x4A,
     0xCD, 0xA3, 0xDD, 0x86, 0x9F, 0x74, 0xE0, 0x1F, 0x8B, 0xA5, 0xE0, 0x32,
     0x43, 0x09, 0xDB, 0x6A, 0x98, 0x31, 0x49, 0x7A, 0xBA, 0xC9, 0x66, 0x70
 };

 static const unsigned char bp384_Z[] = {
     0x0B, 0xD9, 0xD3, 0xA7, 0xEA, 0x0B, 0x3D, 0x51, 0x9D, 0x09, 0xD8, 0xE4,
     0x8D, 0x07, 0x85, 0xFB, 0x74, 0x4A, 0x6B, 0x35, 0x5E, 0x63, 0x04, 0xBC,
     0x51, 0xC2, 0x29, 0xFB, 0xBC, 0xE2, 0x39, 0xBB, 0xAD, 0xF6, 0x40, 0x37,
     0x15, 0xC3, 0x5D, 0x4F, 0xB2, 0xA5, 0x44, 0x4F, 0x57, 0x5D, 0x4F, 0x42
 };

 static const unsigned char bp512_da[] = {
     0x16, 0x30, 0x2F, 0xF0, 0xDB, 0xBB, 0x5A, 0x8D, 0x73, 0x3D, 0xAB, 0x71,
     0x41, 0xC1, 0xB4, 0x5A, 0xCB, 0xC8, 0x71, 0x59, 0x39, 0x67, 0x7F, 0x6A,
     0x56, 0x85, 0x0A, 0x38, 0xBD, 0x87, 0xBD, 0x59, 0xB0, 0x9E, 0x80, 0x27,
     0x96, 0x09, 0xFF, 0x33, 0x3E, 0xB9, 0xD4, 0xC0, 0x61, 0x23, 0x1F, 0xB2,
     0x6F, 0x92, 0xEE, 0xB0, 0x49, 0x82, 0xA5, 0xF1, 0xD1, 0x76, 0x4C, 0xAD,
     0x57, 0x66, 0x54, 0x22
 };

 static const unsigned char bp512_db[] = {
     0x23, 0x0E, 0x18, 0xE1, 0xBC, 0xC8, 0x8A, 0x36, 0x2F, 0xA5, 0x4E, 0x4E,
     0xA3, 0x90, 0x20, 0x09, 0x29, 0x2F, 0x7F, 0x80, 0x33, 0x62, 0x4F, 0xD4,
     0x71, 0xB5, 0xD8, 0xAC, 0xE4, 0x9D, 0x12, 0xCF, 0xAB, 0xBC, 0x19, 0x96,
     0x3D, 0xAB, 0x8E, 0x2F, 0x1E, 0xBA, 0x00, 0xBF, 0xFB, 0x29, 0xE4, 0xD7,
     0x2D, 0x13, 0xF2, 0x22, 0x45, 0x62, 0xF4, 0x05, 0xCB, 0x80, 0x50, 0x36,
     0x66, 0xB2, 0x54, 0x29
 };

 static const unsigned char bp512_Z[] = {
     0xA7, 0x92, 0x70, 0x98, 0x65, 0x5F, 0x1F, 0x99, 0x76, 0xFA, 0x50, 0xA9,
     0xD5, 0x66, 0x86, 0x5D, 0xC5, 0x30, 0x33, 0x18, 0x46, 0x38, 0x1C, 0x87,
     0x25, 0x6B, 0xAF, 0x32, 0x26, 0x24, 0x4B, 0x76, 0xD3, 0x64, 0x03, 0xC0,
     0x24, 0xD7, 0xBB, 0xF0, 0xAA, 0x08, 0x03, 0xEA, 0xFF, 0x40, 0x5D, 0x3D,
     0x24, 0xF1, 0x1A, 0x9B, 0x5C, 0x0B, 0xEF, 0x67, 0x9F, 0xE1, 0x45, 0x4B,
     0x21, 0xC4, 0xCD, 0x1F
 };

 /* Given private value and NID, create EC_KEY structure */

 static EC_KEY *mk_eckey(int nid, const unsigned char *p, size_t plen)
 {
     int ok = 0;
     EC_KEY *k = NULL;
     BIGNUM *priv = NULL;
     EC_POINT *pub = NULL;
     const EC_GROUP *grp;
     k = EC_KEY_new_by_curve_name(nid);
     if (!k)
         goto err;
     priv = BN_bin2bn(p, plen, NULL);
     if (!priv)
         goto err;
     if (!EC_KEY_set_private_key(k, priv))
         goto err;
     grp = EC_KEY_get0_group(k);
     pub = EC_POINT_new(grp);
     if (!pub)
         goto err;
     if (!EC_POINT_mul(grp, pub, priv, NULL, NULL, NULL))
         goto err;
     if (!EC_KEY_set_public_key(k, pub))
         goto err;
     ok = 1;
  err:
     if (priv)
         BN_clear_free(priv);
     if (pub)
         EC_POINT_free(pub);
     if (ok)
         return k;
     else if (k)
         EC_KEY_free(k);
     return NULL;
 }

 /*
  * Known answer test: compute shared secret and check it matches expected
  * value.
  */

 static int ecdh_kat(BIO *out, const char *cname, int nid,
                     const unsigned char *k1, size_t k1_len,
                     const unsigned char *k2, size_t k2_len,
                     const unsigned char *Z, size_t Zlen)
 {
     int rv = 0;
     EC_KEY *key1 = NULL, *key2 = NULL;
     unsigned char *Ztmp = NULL;
     size_t Ztmplen;
     BIO_puts(out, "Testing ECDH shared secret with ");
     BIO_puts(out, cname);
     key1 = mk_eckey(nid, k1, k1_len);
     key2 = mk_eckey(nid, k2, k2_len);
     if (!key1 || !key2)
         goto err;
     Ztmplen = (EC_GROUP_get_degree(EC_KEY_get0_group(key1)) + 7) / 8;
     if (Ztmplen != Zlen)
         goto err;
     Ztmp = OPENSSL_malloc(Ztmplen);
     if (!ECDH_compute_key(Ztmp, Ztmplen,
                           EC_KEY_get0_public_key(key2), key1, 0))
         goto err;
     if (memcmp(Ztmp, Z, Zlen))
         goto err;
     memset(Ztmp, 0, Zlen);
     if (!ECDH_compute_key(Ztmp, Ztmplen,
                           EC_KEY_get0_public_key(key1), key2, 0))
         goto err;
     if (memcmp(Ztmp, Z, Zlen))
         goto err;
     rv = 1;
  err:
     if (key1)
         EC_KEY_free(key1);
     if (key2)
         EC_KEY_free(key2);
     if (Ztmp)
         OPENSSL_free(Ztmp);
     if (rv)
         BIO_puts(out, " ok\n");
     else {
         fprintf(stderr, "Error in ECDH routines\n");
         ERR_print_errors_fp(stderr);
     }
     return rv;
 }

 # define test_ecdh_kat(bio, curve, bits) \
         ecdh_kat(bio, curve, NID_brainpoolP##bits##r1, \
                 bp##bits##_da, sizeof(bp##bits##_da), \
                 bp##bits##_db, sizeof(bp##bits##_db), \
                 bp##bits##_Z, sizeof(bp##bits##_Z))

 int main(int argc, char *argv[])
 {
     BN_CTX *ctx = NULL;
     int ret = 1;
     BIO *out;

     CRYPTO_malloc_debug_init();
     CRYPTO_dbg_set_options(V_CRYPTO_MDEBUG_ALL);
     CRYPTO_mem_ctrl(CRYPTO_MEM_CHECK_ON);

 # ifdef OPENSSL_SYS_WIN32
     CRYPTO_malloc_init();
 # endif

     RAND_seed(rnd_seed, sizeof rnd_seed);

     out = BIO_new(BIO_s_file());
     if (out == NULL)
         EXIT(1);
     BIO_set_fp(out, stdout, BIO_NOCLOSE);

     if ((ctx = BN_CTX_new()) == NULL)
         goto err;

     /* NIST PRIME CURVES TESTS */
     if (!test_ecdh_curve
         (NID_X9_62_prime192v1, "NIST Prime-Curve P-192", ctx, out))
         goto err;
     if (!test_ecdh_curve(NID_secp224r1, "NIST Prime-Curve P-224", ctx, out))
         goto err;
     if (!test_ecdh_curve
         (NID_X9_62_prime256v1, "NIST Prime-Curve P-256", ctx, out))
         goto err;
     if (!test_ecdh_curve(NID_secp384r1, "NIST Prime-Curve P-384", ctx, out))
         goto err;
     if (!test_ecdh_curve(NID_secp521r1, "NIST Prime-Curve P-521", ctx, out))
         goto err;
 # ifndef OPENSSL_NO_EC2M
     /* NIST BINARY CURVES TESTS */
     if (!test_ecdh_curve(NID_sect163k1, "NIST Binary-Curve K-163", ctx, out))
         goto err;
     if (!test_ecdh_curve(NID_sect163r2, "NIST Binary-Curve B-163", ctx, out))
         goto err;
     if (!test_ecdh_curve(NID_sect233k1, "NIST Binary-Curve K-233", ctx, out))
         goto err;
     if (!test_ecdh_curve(NID_sect233r1, "NIST Binary-Curve B-233", ctx, out))
         goto err;
     if (!test_ecdh_curve(NID_sect283k1, "NIST Binary-Curve K-283", ctx, out))
         goto err;
     if (!test_ecdh_curve(NID_sect283r1, "NIST Binary-Curve B-283", ctx, out))
         goto err;
     if (!test_ecdh_curve(NID_sect409k1, "NIST Binary-Curve K-409", ctx, out))
         goto err;
     if (!test_ecdh_curve(NID_sect409r1, "NIST Binary-Curve B-409", ctx, out))
         goto err;
     if (!test_ecdh_curve(NID_sect571k1, "NIST Binary-Curve K-571", ctx, out))
         goto err;
     if (!test_ecdh_curve(NID_sect571r1, "NIST Binary-Curve B-571", ctx, out))
         goto err;
 # endif
     if (!test_ecdh_kat(out, "Brainpool Prime-Curve brainpoolP256r1", 256))
         goto err;
     if (!test_ecdh_kat(out, "Brainpool Prime-Curve brainpoolP384r1", 384))
         goto err;
     if (!test_ecdh_kat(out, "Brainpool Prime-Curve brainpoolP512r1", 512))
         goto err;

     ret = 0;

  err:
     ERR_print_errors_fp(stderr);
     if (ctx)
         BN_CTX_free(ctx);
     BIO_free(out);
     CRYPTO_cleanup_all_ex_data();
     ERR_remove_thread_state(NULL);
     CRYPTO_mem_leaks_fp(stderr);
     EXIT(ret);
 }
 #endif
	/* crypto/ecdh/ecdhtest.c */
	/* ====================================================================
	* Copyright 2002 Sun Microsystems, Inc. ALL RIGHTS RESERVED.
	*
	* The Elliptic Curve Public-Key Crypto Library (ECC Code) included
	* herein is developed by SUN MICROSYSTEMS, INC., and is contributed
	* to the OpenSSL project.
	*
	* The ECC Code is licensed pursuant to the OpenSSL open source
	* license provided below.
	*
	* The ECDH software is originally written by Douglas Stebila of
	* Sun Microsystems Laboratories.
	*
	*/
	/* ====================================================================
	* Copyright (c) 1998-2003 The OpenSSL Project. All rights reserved.
	*
	* Redistribution and use in source and binary forms, with or without
	* modification, are permitted provided that the following conditions
	* are met:
	*
	* 1. Redistributions of source code must retain the above copyright
	* notice, this list of conditions and the following disclaimer.
	*
	* 2. Redistributions in binary form must reproduce the above copyright
	* notice, this list of conditions and the following disclaimer in
	* the documentation and/or other materials provided with the
	* distribution.
	*
	* 3. All advertising materials mentioning features or use of this
	* software must display the following acknowledgment:
	* "This product includes software developed by the OpenSSL Project
	* for use in the OpenSSL Toolkit. (http://www.openssl.org/)"
	*
	* 4. The names "OpenSSL Toolkit" and "OpenSSL Project" must not be used to
	* endorse or promote products derived from this software without
	* prior written permission. For written permission, please contact
	* openssl-core@openssl.org.
	*
	* 5. Products derived from this software may not be called "OpenSSL"
	* nor may "OpenSSL" appear in their names without prior written
	* permission of the OpenSSL Project.
	*
	* 6. Redistributions of any form whatsoever must retain the following
	* acknowledgment:
	* "This product includes software developed by the OpenSSL Project
	* for use in the OpenSSL Toolkit (http://www.openssl.org/)"
	*
	* THIS SOFTWARE IS PROVIDED BY THE OpenSSL PROJECT ``AS IS'' AND ANY
	* EXPRESSED OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
	* IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR
	* PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE OpenSSL PROJECT OR
	* ITS CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL,
	* SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT
	* NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES;
	* LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION)
	* HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT,
	* STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE)
	* ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED
	* OF THE POSSIBILITY OF SUCH DAMAGE.
	* ====================================================================
	*
	* This product includes cryptographic software written by Eric Young
	* (eay@cryptsoft.com). This product includes software written by Tim
	* Hudson (tjh@cryptsoft.com).
	*
	*/

	#include <stdio.h>
	#include <stdlib.h>
	#include <string.h>

	#include "../e_os.h"

	#include <openssl/opensslconf.h> /* for OPENSSL_NO_ECDH */
	#include <openssl/crypto.h>
	#include <openssl/bio.h>
	#include <openssl/bn.h>
	#include <openssl/objects.h>
	#include <openssl/rand.h>
	#include <openssl/sha.h>
	#include <openssl/err.h>

	#ifdef OPENSSL_NO_ECDH
	int main(int argc, char *argv[])
	{
	printf("No ECDH support\n");
	return (0);
	}
	#else
	# include <openssl/ec.h>
	# include <openssl/ecdh.h>

	static const char rnd_seed[] =
	"string to make the random number generator think it has entropy";

	static const int KDF1_SHA1_len = 20;
	static void KDF1_SHA1(const void in, size_t inlen, void *out,
	size_t *outlen)
	{
	if (*outlen < SHA_DIGEST_LENGTH)
	return NULL;
	*outlen = SHA_DIGEST_LENGTH;
	return SHA1(in, inlen, out);
	}

	static int test_ecdh_curve(int nid, const char text, BN_CTX ctx, BIO *out)
	{
	EC_KEY *a = NULL;
	EC_KEY *b = NULL;
	BIGNUM x_a = NULL, y_a = NULL, x_b = NULL, y_b = NULL;
	char buf[12];
	unsigned char abuf = NULL, bbuf = NULL;
	int i, alen, blen, aout, bout, ret = 0;
	const EC_GROUP *group;

	a = EC_KEY_new_by_curve_name(nid);
	b = EC_KEY_new_by_curve_name(nid);
	if (a == NULL \|\| b == NULL)
	goto err;

	group = EC_KEY_get0_group(a);

	if ((x_a = BN_new()) == NULL)
	goto err;
	if ((y_a = BN_new()) == NULL)
	goto err;
	if ((x_b = BN_new()) == NULL)
	goto err;
	if ((y_b = BN_new()) == NULL)
	goto err;

	BIO_puts(out, "Testing key generation with ");
	BIO_puts(out, text);
	# ifdef NOISY
	BIO_puts(out, "\n");
	# else
	(void)BIO_flush(out);
	# endif

	if (!EC_KEY_generate_key(a))
	goto err;

	if (EC_METHOD_get_field_type(EC_GROUP_method_of(group)) ==
	NID_X9_62_prime_field) {
	if (!EC_POINT_get_affine_coordinates_GFp
	(group, EC_KEY_get0_public_key(a), x_a, y_a, ctx))
	goto err;
	}
	# ifndef OPENSSL_NO_EC2M
	else {
	if (!EC_POINT_get_affine_coordinates_GF2m(group,
	EC_KEY_get0_public_key(a),
	x_a, y_a, ctx))
	goto err;
	}
	# endif
	# ifdef NOISY
	BIO_puts(out, " pri 1=");
	BN_print(out, a->priv_key);
	BIO_puts(out, "\n pub 1=");
	BN_print(out, x_a);
	BIO_puts(out, ",");
	BN_print(out, y_a);
	BIO_puts(out, "\n");
	# else
	BIO_printf(out, " .");
	(void)BIO_flush(out);
	# endif

	if (!EC_KEY_generate_key(b))
	goto err;

	if (EC_METHOD_get_field_type(EC_GROUP_method_of(group)) ==
	NID_X9_62_prime_field) {
	if (!EC_POINT_get_affine_coordinates_GFp
	(group, EC_KEY_get0_public_key(b), x_b, y_b, ctx))
	goto err;
	}
	# ifndef OPENSSL_NO_EC2M
	else {
	if (!EC_POINT_get_affine_coordinates_GF2m(group,
	EC_KEY_get0_public_key(b),
	x_b, y_b, ctx))
	goto err;
	}
	# endif

	# ifdef NOISY
	BIO_puts(out, " pri 2=");
	BN_print(out, b->priv_key);
	BIO_puts(out, "\n pub 2=");
	BN_print(out, x_b);
	BIO_puts(out, ",");
	BN_print(out, y_b);
	BIO_puts(out, "\n");
	# else
	BIO_printf(out, ".");
	(void)BIO_flush(out);
	# endif

	alen = KDF1_SHA1_len;
	abuf = (unsigned char *)OPENSSL_malloc(alen);
	aout =
	ECDH_compute_key(abuf, alen, EC_KEY_get0_public_key(b), a, KDF1_SHA1);

	# ifdef NOISY
	BIO_puts(out, " key1 =");
	for (i = 0; i < aout; i++) {
	sprintf(buf, "%02X", abuf[i]);
	BIO_puts(out, buf);
	}
	BIO_puts(out, "\n");
	# else
	BIO_printf(out, ".");
	(void)BIO_flush(out);
	# endif

	blen = KDF1_SHA1_len;
	bbuf = (unsigned char *)OPENSSL_malloc(blen);
	bout =
	ECDH_compute_key(bbuf, blen, EC_KEY_get0_public_key(a), b, KDF1_SHA1);

	# ifdef NOISY
	BIO_puts(out, " key2 =");
	for (i = 0; i < bout; i++) {
	sprintf(buf, "%02X", bbuf[i]);
	BIO_puts(out, buf);
	}
	BIO_puts(out, "\n");
	# else
	BIO_printf(out, ".");
	(void)BIO_flush(out);
	# endif

	if ((aout < 4) \|\| (bout != aout) \|\| (memcmp(abuf, bbuf, aout) != 0)) {
	# ifndef NOISY
	BIO_printf(out, " failed\n\n");
	BIO_printf(out, "key a:\n");
	BIO_printf(out, "private key: ");
	BN_print(out, EC_KEY_get0_private_key(a));
	BIO_printf(out, "\n");
	BIO_printf(out, "public key (x,y): ");
	BN_print(out, x_a);
	BIO_printf(out, ",");
	BN_print(out, y_a);
	BIO_printf(out, "\nkey b:\n");
	BIO_printf(out, "private key: ");
	BN_print(out, EC_KEY_get0_private_key(b));
	BIO_printf(out, "\n");
	BIO_printf(out, "public key (x,y): ");
	BN_print(out, x_b);
	BIO_printf(out, ",");
	BN_print(out, y_b);
	BIO_printf(out, "\n");
	BIO_printf(out, "generated key a: ");
	for (i = 0; i < bout; i++) {
	sprintf(buf, "%02X", bbuf[i]);
	BIO_puts(out, buf);
	}
	BIO_printf(out, "\n");
	BIO_printf(out, "generated key b: ");
	for (i = 0; i < aout; i++) {
	sprintf(buf, "%02X", abuf[i]);
	BIO_puts(out, buf);
	}
	BIO_printf(out, "\n");
	# endif
	fprintf(stderr, "Error in ECDH routines\n");
	ret = 0;
	} else {
	# ifndef NOISY
	BIO_printf(out, " ok\n");
	# endif
	ret = 1;
	}
	err:
	ERR_print_errors_fp(stderr);

	if (abuf != NULL)
	OPENSSL_free(abuf);
	if (bbuf != NULL)
	OPENSSL_free(bbuf);
	if (x_a)
	BN_free(x_a);
	if (y_a)
	BN_free(y_a);
	if (x_b)
	BN_free(x_b);
	if (y_b)
	BN_free(y_b);
	if (b)
	EC_KEY_free(b);
	if (a)
	EC_KEY_free(a);
	return (ret);
	}

	/* Keys and shared secrets from RFC 7027 */

	static const unsigned char bp256_da[] = {
	0x81, 0xDB, 0x1E, 0xE1, 0x00, 0x15, 0x0F, 0xF2, 0xEA, 0x33, 0x8D, 0x70,
	0x82, 0x71, 0xBE, 0x38, 0x30, 0x0C, 0xB5, 0x42, 0x41, 0xD7, 0x99, 0x50,
	0xF7, 0x7B, 0x06, 0x30, 0x39, 0x80, 0x4F, 0x1D
	};

	static const unsigned char bp256_db[] = {
	0x55, 0xE4, 0x0B, 0xC4, 0x1E, 0x37, 0xE3, 0xE2, 0xAD, 0x25, 0xC3, 0xC6,
	0x65, 0x45, 0x11, 0xFF, 0xA8, 0x47, 0x4A, 0x91, 0xA0, 0x03, 0x20, 0x87,
	0x59, 0x38, 0x52, 0xD3, 0xE7, 0xD7, 0x6B, 0xD3
	};

	static const unsigned char bp256_Z[] = {
	0x89, 0xAF, 0xC3, 0x9D, 0x41, 0xD3, 0xB3, 0x27, 0x81, 0x4B, 0x80, 0x94,
	0x0B, 0x04, 0x25, 0x90, 0xF9, 0x65, 0x56, 0xEC, 0x91, 0xE6, 0xAE, 0x79,
	0x39, 0xBC, 0xE3, 0x1F, 0x3A, 0x18, 0xBF, 0x2B
	};

	static const unsigned char bp384_da[] = {
	0x1E, 0x20, 0xF5, 0xE0, 0x48, 0xA5, 0x88, 0x6F, 0x1F, 0x15, 0x7C, 0x74,
	0xE9, 0x1B, 0xDE, 0x2B, 0x98, 0xC8, 0xB5, 0x2D, 0x58, 0xE5, 0x00, 0x3D,
	0x57, 0x05, 0x3F, 0xC4, 0xB0, 0xBD, 0x65, 0xD6, 0xF1, 0x5E, 0xB5, 0xD1,
	0xEE, 0x16, 0x10, 0xDF, 0x87, 0x07, 0x95, 0x14, 0x36, 0x27, 0xD0, 0x42
	};

	static const unsigned char bp384_db[] = {
	0x03, 0x26, 0x40, 0xBC, 0x60, 0x03, 0xC5, 0x92, 0x60, 0xF7, 0x25, 0x0C,
	0x3D, 0xB5, 0x8C, 0xE6, 0x47, 0xF9, 0x8E, 0x12, 0x60, 0xAC, 0xCE, 0x4A,
	0xCD, 0xA3, 0xDD, 0x86, 0x9F, 0x74, 0xE0, 0x1F, 0x8B, 0xA5, 0xE0, 0x32,
	0x43, 0x09, 0xDB, 0x6A, 0x98, 0x31, 0x49, 0x7A, 0xBA, 0xC9, 0x66, 0x70
	};

	static const unsigned char bp384_Z[] = {
	0x0B, 0xD9, 0xD3, 0xA7, 0xEA, 0x0B, 0x3D, 0x51, 0x9D, 0x09, 0xD8, 0xE4,
	0x8D, 0x07, 0x85, 0xFB, 0x74, 0x4A, 0x6B, 0x35, 0x5E, 0x63, 0x04, 0xBC,
	0x51, 0xC2, 0x29, 0xFB, 0xBC, 0xE2, 0x39, 0xBB, 0xAD, 0xF6, 0x40, 0x37,
	0x15, 0xC3, 0x5D, 0x4F, 0xB2, 0xA5, 0x44, 0x4F, 0x57, 0x5D, 0x4F, 0x42
	};

	static const unsigned char bp512_da[] = {
	0x16, 0x30, 0x2F, 0xF0, 0xDB, 0xBB, 0x5A, 0x8D, 0x73, 0x3D, 0xAB, 0x71,
	0x41, 0xC1, 0xB4, 0x5A, 0xCB, 0xC8, 0x71, 0x59, 0x39, 0x67, 0x7F, 0x6A,
	0x56, 0x85, 0x0A, 0x38, 0xBD, 0x87, 0xBD, 0x59, 0xB0, 0x9E, 0x80, 0x27,
	0x96, 0x09, 0xFF, 0x33, 0x3E, 0xB9, 0xD4, 0xC0, 0x61, 0x23, 0x1F, 0xB2,
	0x6F, 0x92, 0xEE, 0xB0, 0x49, 0x82, 0xA5, 0xF1, 0xD1, 0x76, 0x4C, 0xAD,
	0x57, 0x66, 0x54, 0x22
	};

	static const unsigned char bp512_db[] = {
	0x23, 0x0E, 0x18, 0xE1, 0xBC, 0xC8, 0x8A, 0x36, 0x2F, 0xA5, 0x4E, 0x4E,
	0xA3, 0x90, 0x20, 0x09, 0x29, 0x2F, 0x7F, 0x80, 0x33, 0x62, 0x4F, 0xD4,
	0x71, 0xB5, 0xD8, 0xAC, 0xE4, 0x9D, 0x12, 0xCF, 0xAB, 0xBC, 0x19, 0x96,
	0x3D, 0xAB, 0x8E, 0x2F, 0x1E, 0xBA, 0x00, 0xBF, 0xFB, 0x29, 0xE4, 0xD7,
	0x2D, 0x13, 0xF2, 0x22, 0x45, 0x62, 0xF4, 0x05, 0xCB, 0x80, 0x50, 0x36,
	0x66, 0xB2, 0x54, 0x29
	};

	static const unsigned char bp512_Z[] = {
	0xA7, 0x92, 0x70, 0x98, 0x65, 0x5F, 0x1F, 0x99, 0x76, 0xFA, 0x50, 0xA9,
	0xD5, 0x66, 0x86, 0x5D, 0xC5, 0x30, 0x33, 0x18, 0x46, 0x38, 0x1C, 0x87,
	0x25, 0x6B, 0xAF, 0x32, 0x26, 0x24, 0x4B, 0x76, 0xD3, 0x64, 0x03, 0xC0,
	0x24, 0xD7, 0xBB, 0xF0, 0xAA, 0x08, 0x03, 0xEA, 0xFF, 0x40, 0x5D, 0x3D,
	0x24, 0xF1, 0x1A, 0x9B, 0x5C, 0x0B, 0xEF, 0x67, 0x9F, 0xE1, 0x45, 0x4B,
	0x21, 0xC4, 0xCD, 0x1F
	};

	/* Given private value and NID, create EC_KEY structure */

	static EC_KEY mk_eckey(int nid, const unsigned char p, size_t plen)
	{
	int ok = 0;
	EC_KEY *k = NULL;
	BIGNUM *priv = NULL;
	EC_POINT *pub = NULL;
	const EC_GROUP *grp;
	k = EC_KEY_new_by_curve_name(nid);
	if (!k)
	goto err;
	priv = BN_bin2bn(p, plen, NULL);
	if (!priv)
	goto err;
	if (!EC_KEY_set_private_key(k, priv))
	goto err;
	grp = EC_KEY_get0_group(k);
	pub = EC_POINT_new(grp);
	if (!pub)
	goto err;
	if (!EC_POINT_mul(grp, pub, priv, NULL, NULL, NULL))
	goto err;
	if (!EC_KEY_set_public_key(k, pub))
	goto err;
	ok = 1;
	err:
	if (priv)
	BN_clear_free(priv);
	if (pub)
	EC_POINT_free(pub);
	if (ok)
	return k;
	else if (k)
	EC_KEY_free(k);
	return NULL;
	}

	/*
	* Known answer test: compute shared secret and check it matches expected
	* value.
	*/

	static int ecdh_kat(BIO out, const char cname, int nid,
	const unsigned char *k1, size_t k1_len,
	const unsigned char *k2, size_t k2_len,
	const unsigned char *Z, size_t Zlen)
	{
	int rv = 0;
	EC_KEY key1 = NULL, key2 = NULL;
	unsigned char *Ztmp = NULL;
	size_t Ztmplen;
	BIO_puts(out, "Testing ECDH shared secret with ");
	BIO_puts(out, cname);
	key1 = mk_eckey(nid, k1, k1_len);
	key2 = mk_eckey(nid, k2, k2_len);
	if (!key1 \|\| !key2)
	goto err;
	Ztmplen = (EC_GROUP_get_degree(EC_KEY_get0_group(key1)) + 7) / 8;
	if (Ztmplen != Zlen)
	goto err;
	Ztmp = OPENSSL_malloc(Ztmplen);
	if (!ECDH_compute_key(Ztmp, Ztmplen,
	EC_KEY_get0_public_key(key2), key1, 0))
	goto err;
	if (memcmp(Ztmp, Z, Zlen))
	goto err;
	memset(Ztmp, 0, Zlen);
	if (!ECDH_compute_key(Ztmp, Ztmplen,
	EC_KEY_get0_public_key(key1), key2, 0))
	goto err;
	if (memcmp(Ztmp, Z, Zlen))
	goto err;
	rv = 1;
	err:
	if (key1)
	EC_KEY_free(key1);
	if (key2)
	EC_KEY_free(key2);
	if (Ztmp)
	OPENSSL_free(Ztmp);
	if (rv)
	BIO_puts(out, " ok\n");
	else {
	fprintf(stderr, "Error in ECDH routines\n");
	ERR_print_errors_fp(stderr);
	}
	return rv;
	}

	# define test_ecdh_kat(bio, curve, bits) \
	ecdh_kat(bio, curve, NID_brainpoolP##bits##r1, \
	bp##bits##_da, sizeof(bp##bits##_da), \
	bp##bits##_db, sizeof(bp##bits##_db), \
	bp##bits##_Z, sizeof(bp##bits##_Z))

	int main(int argc, char *argv[])
	{
	BN_CTX *ctx = NULL;
	int ret = 1;
	BIO *out;

	CRYPTO_malloc_debug_init();
	CRYPTO_dbg_set_options(V_CRYPTO_MDEBUG_ALL);
	CRYPTO_mem_ctrl(CRYPTO_MEM_CHECK_ON);

	# ifdef OPENSSL_SYS_WIN32
	CRYPTO_malloc_init();
	# endif

	RAND_seed(rnd_seed, sizeof rnd_seed);

	out = BIO_new(BIO_s_file());
	if (out == NULL)
	EXIT(1);
	BIO_set_fp(out, stdout, BIO_NOCLOSE);

	if ((ctx = BN_CTX_new()) == NULL)
	goto err;

	/* NIST PRIME CURVES TESTS */
	if (!test_ecdh_curve
	(NID_X9_62_prime192v1, "NIST Prime-Curve P-192", ctx, out))
	goto err;
	if (!test_ecdh_curve(NID_secp224r1, "NIST Prime-Curve P-224", ctx, out))
	goto err;
	if (!test_ecdh_curve
	(NID_X9_62_prime256v1, "NIST Prime-Curve P-256", ctx, out))
	goto err;
	if (!test_ecdh_curve(NID_secp384r1, "NIST Prime-Curve P-384", ctx, out))
	goto err;
	if (!test_ecdh_curve(NID_secp521r1, "NIST Prime-Curve P-521", ctx, out))
	goto err;
	# ifndef OPENSSL_NO_EC2M
	/* NIST BINARY CURVES TESTS */
	if (!test_ecdh_curve(NID_sect163k1, "NIST Binary-Curve K-163", ctx, out))
	goto err;
	if (!test_ecdh_curve(NID_sect163r2, "NIST Binary-Curve B-163", ctx, out))
	goto err;
	if (!test_ecdh_curve(NID_sect233k1, "NIST Binary-Curve K-233", ctx, out))
	goto err;
	if (!test_ecdh_curve(NID_sect233r1, "NIST Binary-Curve B-233", ctx, out))
	goto err;
	if (!test_ecdh_curve(NID_sect283k1, "NIST Binary-Curve K-283", ctx, out))
	goto err;
	if (!test_ecdh_curve(NID_sect283r1, "NIST Binary-Curve B-283", ctx, out))
	goto err;
	if (!test_ecdh_curve(NID_sect409k1, "NIST Binary-Curve K-409", ctx, out))
	goto err;
	if (!test_ecdh_curve(NID_sect409r1, "NIST Binary-Curve B-409", ctx, out))
	goto err;
	if (!test_ecdh_curve(NID_sect571k1, "NIST Binary-Curve K-571", ctx, out))
	goto err;
	if (!test_ecdh_curve(NID_sect571r1, "NIST Binary-Curve B-571", ctx, out))
	goto err;
	# endif
	if (!test_ecdh_kat(out, "Brainpool Prime-Curve brainpoolP256r1", 256))
	goto err;
	if (!test_ecdh_kat(out, "Brainpool Prime-Curve brainpoolP384r1", 384))
	goto err;
	if (!test_ecdh_kat(out, "Brainpool Prime-Curve brainpoolP512r1", 512))
	goto err;

	ret = 0;

	err:
	ERR_print_errors_fp(stderr);
	if (ctx)
	BN_CTX_free(ctx);
	BIO_free(out);
	CRYPTO_cleanup_all_ex_data();
	ERR_remove_thread_state(NULL);
	CRYPTO_mem_leaks_fp(stderr);
	EXIT(ret);
	}
	#endif