xref: /openbsd-src/regress/lib/libcrypto/ecdh/ecdhtest.c (revision 50b7afb2c2c0993b0894d4e34bf857cb13ed9c80)
1 /* crypto/ecdh/ecdhtest.c */
2 /* ====================================================================
3  * Copyright 2002 Sun Microsystems, Inc. ALL RIGHTS RESERVED.
4  *
5  * The Elliptic Curve Public-Key Crypto Library (ECC Code) included
6  * herein is developed by SUN MICROSYSTEMS, INC., and is contributed
7  * to the OpenSSL project.
8  *
9  * The ECC Code is licensed pursuant to the OpenSSL open source
10  * license provided below.
11  *
12  * The ECDH software is originally written by Douglas Stebila of
13  * Sun Microsystems Laboratories.
14  *
15  */
16 /* ====================================================================
17  * Copyright (c) 1998-2003 The OpenSSL Project.  All rights reserved.
18  *
19  * Redistribution and use in source and binary forms, with or without
20  * modification, are permitted provided that the following conditions
21  * are met:
22  *
23  * 1. Redistributions of source code must retain the above copyright
24  *    notice, this list of conditions and the following disclaimer.
25  *
26  * 2. Redistributions in binary form must reproduce the above copyright
27  *    notice, this list of conditions and the following disclaimer in
28  *    the documentation and/or other materials provided with the
29  *    distribution.
30  *
31  * 3. All advertising materials mentioning features or use of this
32  *    software must display the following acknowledgment:
33  *    "This product includes software developed by the OpenSSL Project
34  *    for use in the OpenSSL Toolkit. (http://www.openssl.org/)"
35  *
36  * 4. The names "OpenSSL Toolkit" and "OpenSSL Project" must not be used to
37  *    endorse or promote products derived from this software without
38  *    prior written permission. For written permission, please contact
39  *    openssl-core@openssl.org.
40  *
41  * 5. Products derived from this software may not be called "OpenSSL"
42  *    nor may "OpenSSL" appear in their names without prior written
43  *    permission of the OpenSSL Project.
44  *
45  * 6. Redistributions of any form whatsoever must retain the following
46  *    acknowledgment:
47  *    "This product includes software developed by the OpenSSL Project
48  *    for use in the OpenSSL Toolkit (http://www.openssl.org/)"
49  *
50  * THIS SOFTWARE IS PROVIDED BY THE OpenSSL PROJECT ``AS IS'' AND ANY
51  * EXPRESSED OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
52  * IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR
53  * PURPOSE ARE DISCLAIMED.  IN NO EVENT SHALL THE OpenSSL PROJECT OR
54  * ITS CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL,
55  * SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT
56  * NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES;
57  * LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION)
58  * HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT,
59  * STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE)
60  * ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED
61  * OF THE POSSIBILITY OF SUCH DAMAGE.
62  * ====================================================================
63  *
64  * This product includes cryptographic software written by Eric Young
65  * (eay@cryptsoft.com).  This product includes software written by Tim
66  * Hudson (tjh@cryptsoft.com).
67  *
68  */
69 
70 #include <stdio.h>
71 #include <stdlib.h>
72 #include <string.h>
73 
74 #include <openssl/crypto.h>
75 #include <openssl/bio.h>
76 #include <openssl/bn.h>
77 #include <openssl/objects.h>
78 #include <openssl/rand.h>
79 #include <openssl/sha.h>
80 #include <openssl/err.h>
81 
82 #include <openssl/ec.h>
83 #include <openssl/ecdh.h>
84 
85 static const int KDF1_SHA1_len = 20;
86 static void *KDF1_SHA1(const void *in, size_t inlen, void *out, size_t *outlen)
87 	{
88 #ifndef OPENSSL_NO_SHA
89 	if (*outlen < SHA_DIGEST_LENGTH)
90 		return NULL;
91 	else
92 		*outlen = SHA_DIGEST_LENGTH;
93 	return SHA1(in, inlen, out);
94 #else
95 	return NULL;
96 #endif
97 	}
98 
99 
100 static int test_ecdh_curve(int nid, const char *text, BN_CTX *ctx, BIO *out)
101 	{
102 	EC_KEY *a=NULL;
103 	EC_KEY *b=NULL;
104 	BIGNUM *x_a=NULL, *y_a=NULL,
105 	       *x_b=NULL, *y_b=NULL;
106 	char buf[12];
107 	unsigned char *abuf=NULL,*bbuf=NULL;
108 	int i,alen,blen,aout,bout,ret=0;
109 	const EC_GROUP *group;
110 
111 	a = EC_KEY_new_by_curve_name(nid);
112 	b = EC_KEY_new_by_curve_name(nid);
113 	if (a == NULL || b == NULL)
114 		goto err;
115 
116 	group = EC_KEY_get0_group(a);
117 
118 	if ((x_a=BN_new()) == NULL) goto err;
119 	if ((y_a=BN_new()) == NULL) goto err;
120 	if ((x_b=BN_new()) == NULL) goto err;
121 	if ((y_b=BN_new()) == NULL) goto err;
122 
123 	BIO_puts(out,"Testing key generation with ");
124 	BIO_puts(out,text);
125 	(void)BIO_flush(out);
126 
127 	if (!EC_KEY_generate_key(a)) goto err;
128 
129 	if (EC_METHOD_get_field_type(EC_GROUP_method_of(group)) == NID_X9_62_prime_field)
130 		{
131 		if (!EC_POINT_get_affine_coordinates_GFp(group,
132 			EC_KEY_get0_public_key(a), x_a, y_a, ctx)) goto err;
133 		}
134 #ifndef OPENSSL_NO_EC2M
135 	else
136 		{
137 		if (!EC_POINT_get_affine_coordinates_GF2m(group,
138 			EC_KEY_get0_public_key(a), x_a, y_a, ctx)) goto err;
139 		}
140 #endif
141 	BIO_printf(out," .");
142 	(void)BIO_flush(out);
143 
144 	if (!EC_KEY_generate_key(b)) goto err;
145 
146 	if (EC_METHOD_get_field_type(EC_GROUP_method_of(group)) == NID_X9_62_prime_field)
147 		{
148 		if (!EC_POINT_get_affine_coordinates_GFp(group,
149 			EC_KEY_get0_public_key(b), x_b, y_b, ctx)) goto err;
150 		}
151 #ifndef OPENSSL_NO_EC2M
152 	else
153 		{
154 		if (!EC_POINT_get_affine_coordinates_GF2m(group,
155 			EC_KEY_get0_public_key(b), x_b, y_b, ctx)) goto err;
156 		}
157 #endif
158 
159 	BIO_printf(out,".");
160 	(void)BIO_flush(out);
161 
162 	alen=KDF1_SHA1_len;
163 	abuf=malloc(alen);
164 	aout=ECDH_compute_key(abuf,alen,EC_KEY_get0_public_key(b),a,KDF1_SHA1);
165 
166 	BIO_printf(out,".");
167 	(void)BIO_flush(out);
168 
169 	blen=KDF1_SHA1_len;
170 	bbuf=malloc(blen);
171 	bout=ECDH_compute_key(bbuf,blen,EC_KEY_get0_public_key(a),b,KDF1_SHA1);
172 
173 	BIO_printf(out,".");
174 	(void)BIO_flush(out);
175 
176 	if ((aout < 4) || (bout != aout) || (memcmp(abuf,bbuf,aout) != 0))
177 		{
178 		BIO_printf(out, " failed\n\n");
179 		BIO_printf(out, "key a:\n");
180 		BIO_printf(out, "private key: ");
181 		BN_print(out, EC_KEY_get0_private_key(a));
182 		BIO_printf(out, "\n");
183 		BIO_printf(out, "public key (x,y): ");
184 		BN_print(out, x_a);
185 		BIO_printf(out, ",");
186 		BN_print(out, y_a);
187 		BIO_printf(out, "\nkey b:\n");
188 		BIO_printf(out, "private key: ");
189 		BN_print(out, EC_KEY_get0_private_key(b));
190 		BIO_printf(out, "\n");
191 		BIO_printf(out, "public key (x,y): ");
192 		BN_print(out, x_b);
193 		BIO_printf(out, ",");
194 		BN_print(out, y_b);
195 		BIO_printf(out, "\n");
196 		BIO_printf(out, "generated key a: ");
197 		for (i=0; i<bout; i++)
198 			{
199 			snprintf(buf, sizeof buf, "%02X", bbuf[i]);
200 			BIO_puts(out, buf);
201 			}
202 		BIO_printf(out, "\n");
203 		BIO_printf(out, "generated key b: ");
204 		for (i=0; i<aout; i++)
205 			{
206 			snprintf(buf, sizeof buf, "%02X", abuf[i]);
207 			BIO_puts(out,buf);
208 			}
209 		BIO_printf(out, "\n");
210 		fprintf(stderr,"Error in ECDH routines\n");
211 		ret=0;
212 		}
213 	else
214 		{
215 		BIO_printf(out, " ok\n");
216 		ret=1;
217 		}
218 err:
219 	ERR_print_errors_fp(stderr);
220 
221 	free(abuf);
222 	free(bbuf);
223 	if (x_a) BN_free(x_a);
224 	if (y_a) BN_free(y_a);
225 	if (x_b) BN_free(x_b);
226 	if (y_b) BN_free(y_b);
227 	if (b) EC_KEY_free(b);
228 	if (a) EC_KEY_free(a);
229 	return(ret);
230 	}
231 
232 /* Keys and shared secrets from RFC 7027 */
233 
234 static const unsigned char bp256_da[] = {
235 	0x81,0xDB,0x1E,0xE1,0x00,0x15,0x0F,0xF2,0xEA,0x33,0x8D,0x70,
236 	0x82,0x71,0xBE,0x38,0x30,0x0C,0xB5,0x42,0x41,0xD7,0x99,0x50,
237 	0xF7,0x7B,0x06,0x30,0x39,0x80,0x4F,0x1D
238 };
239 
240 static const unsigned char bp256_db[] = {
241 	0x55,0xE4,0x0B,0xC4,0x1E,0x37,0xE3,0xE2,0xAD,0x25,0xC3,0xC6,
242 	0x65,0x45,0x11,0xFF,0xA8,0x47,0x4A,0x91,0xA0,0x03,0x20,0x87,
243 	0x59,0x38,0x52,0xD3,0xE7,0xD7,0x6B,0xD3
244 };
245 
246 static const unsigned char bp256_Z[] = {
247 	0x89,0xAF,0xC3,0x9D,0x41,0xD3,0xB3,0x27,0x81,0x4B,0x80,0x94,
248 	0x0B,0x04,0x25,0x90,0xF9,0x65,0x56,0xEC,0x91,0xE6,0xAE,0x79,
249 	0x39,0xBC,0xE3,0x1F,0x3A,0x18,0xBF,0x2B
250 };
251 
252 static const unsigned char bp384_da[] = {
253 	0x1E,0x20,0xF5,0xE0,0x48,0xA5,0x88,0x6F,0x1F,0x15,0x7C,0x74,
254 	0xE9,0x1B,0xDE,0x2B,0x98,0xC8,0xB5,0x2D,0x58,0xE5,0x00,0x3D,
255 	0x57,0x05,0x3F,0xC4,0xB0,0xBD,0x65,0xD6,0xF1,0x5E,0xB5,0xD1,
256 	0xEE,0x16,0x10,0xDF,0x87,0x07,0x95,0x14,0x36,0x27,0xD0,0x42
257 };
258 
259 static const unsigned char bp384_db[] = {
260 	0x03,0x26,0x40,0xBC,0x60,0x03,0xC5,0x92,0x60,0xF7,0x25,0x0C,
261 	0x3D,0xB5,0x8C,0xE6,0x47,0xF9,0x8E,0x12,0x60,0xAC,0xCE,0x4A,
262 	0xCD,0xA3,0xDD,0x86,0x9F,0x74,0xE0,0x1F,0x8B,0xA5,0xE0,0x32,
263 	0x43,0x09,0xDB,0x6A,0x98,0x31,0x49,0x7A,0xBA,0xC9,0x66,0x70
264 };
265 
266 static const unsigned char bp384_Z[] = {
267 	0x0B,0xD9,0xD3,0xA7,0xEA,0x0B,0x3D,0x51,0x9D,0x09,0xD8,0xE4,
268 	0x8D,0x07,0x85,0xFB,0x74,0x4A,0x6B,0x35,0x5E,0x63,0x04,0xBC,
269 	0x51,0xC2,0x29,0xFB,0xBC,0xE2,0x39,0xBB,0xAD,0xF6,0x40,0x37,
270 	0x15,0xC3,0x5D,0x4F,0xB2,0xA5,0x44,0x4F,0x57,0x5D,0x4F,0x42
271 };
272 
273 static const unsigned char bp512_da[] = {
274 	0x16,0x30,0x2F,0xF0,0xDB,0xBB,0x5A,0x8D,0x73,0x3D,0xAB,0x71,
275 	0x41,0xC1,0xB4,0x5A,0xCB,0xC8,0x71,0x59,0x39,0x67,0x7F,0x6A,
276 	0x56,0x85,0x0A,0x38,0xBD,0x87,0xBD,0x59,0xB0,0x9E,0x80,0x27,
277 	0x96,0x09,0xFF,0x33,0x3E,0xB9,0xD4,0xC0,0x61,0x23,0x1F,0xB2,
278 	0x6F,0x92,0xEE,0xB0,0x49,0x82,0xA5,0xF1,0xD1,0x76,0x4C,0xAD,
279 	0x57,0x66,0x54,0x22
280 };
281 
282 static const unsigned char bp512_db[] = {
283 	0x23,0x0E,0x18,0xE1,0xBC,0xC8,0x8A,0x36,0x2F,0xA5,0x4E,0x4E,
284 	0xA3,0x90,0x20,0x09,0x29,0x2F,0x7F,0x80,0x33,0x62,0x4F,0xD4,
285 	0x71,0xB5,0xD8,0xAC,0xE4,0x9D,0x12,0xCF,0xAB,0xBC,0x19,0x96,
286 	0x3D,0xAB,0x8E,0x2F,0x1E,0xBA,0x00,0xBF,0xFB,0x29,0xE4,0xD7,
287 	0x2D,0x13,0xF2,0x22,0x45,0x62,0xF4,0x05,0xCB,0x80,0x50,0x36,
288 	0x66,0xB2,0x54,0x29
289 };
290 
291 
292 static const unsigned char bp512_Z[] = {
293 	0xA7,0x92,0x70,0x98,0x65,0x5F,0x1F,0x99,0x76,0xFA,0x50,0xA9,
294 	0xD5,0x66,0x86,0x5D,0xC5,0x30,0x33,0x18,0x46,0x38,0x1C,0x87,
295 	0x25,0x6B,0xAF,0x32,0x26,0x24,0x4B,0x76,0xD3,0x64,0x03,0xC0,
296 	0x24,0xD7,0xBB,0xF0,0xAA,0x08,0x03,0xEA,0xFF,0x40,0x5D,0x3D,
297 	0x24,0xF1,0x1A,0x9B,0x5C,0x0B,0xEF,0x67,0x9F,0xE1,0x45,0x4B,
298 	0x21,0xC4,0xCD,0x1F
299 };
300 
301 /* Given private value and NID, create EC_KEY structure */
302 
303 static EC_KEY *mk_eckey(int nid, const unsigned char *p, size_t plen)
304 	{
305 	int ok = 0;
306 	EC_KEY *k = NULL;
307 	BIGNUM *priv = NULL;
308 	EC_POINT *pub = NULL;
309 	const EC_GROUP *grp;
310 	k = EC_KEY_new_by_curve_name(nid);
311 	if (!k)
312 		goto err;
313 	priv = BN_bin2bn(p, plen, NULL);
314 	if (!priv)
315 		goto err;
316 	if (!EC_KEY_set_private_key(k, priv))
317 		goto err;
318 	grp = EC_KEY_get0_group(k);
319 	pub = EC_POINT_new(grp);
320 	if (!pub)
321 		goto err;
322 	if (!EC_POINT_mul(grp, pub, priv, NULL, NULL, NULL))
323 		goto err;
324 	if (!EC_KEY_set_public_key(k, pub))
325 		goto err;
326 	ok = 1;
327 	err:
328 	if (priv)
329 		BN_clear_free(priv);
330 	if (pub)
331 		EC_POINT_free(pub);
332 	if (ok)
333 		return k;
334 	else if (k)
335 		EC_KEY_free(k);
336 	return NULL;
337 	}
338 
339 /* Known answer test: compute shared secret and check it matches
340  * expected value.
341  */
342 
343 static int ecdh_kat(BIO *out, const char *cname, int nid,
344 		const unsigned char *k1, size_t k1_len,
345 		const unsigned char *k2, size_t k2_len,
346 		const unsigned char *Z, size_t Zlen)
347 	{
348 	int rv = 0;
349 	EC_KEY *key1 = NULL, *key2 = NULL;
350 	unsigned char *Ztmp = NULL;
351 	size_t Ztmplen;
352 	BIO_puts(out, "Testing ECDH shared secret with ");
353 	BIO_puts(out, cname);
354 	key1 = mk_eckey(nid, k1, k1_len);
355 	key2 = mk_eckey(nid, k2, k2_len);
356 	if (!key1 || !key2)
357 		goto err;
358 	Ztmplen = (EC_GROUP_get_degree(EC_KEY_get0_group(key1)) + 7)/8;
359 	if (Ztmplen != Zlen)
360 		goto err;
361 	Ztmp = malloc(Ztmplen);
362 	if (!ECDH_compute_key(Ztmp, Ztmplen,
363 				EC_KEY_get0_public_key(key2), key1, 0))
364 		goto err;
365 	if (memcmp(Ztmp, Z, Zlen))
366 		goto err;
367 	memset(Ztmp, 0, Zlen);
368 	if (!ECDH_compute_key(Ztmp, Ztmplen,
369 				EC_KEY_get0_public_key(key1), key2, 0))
370 		goto err;
371 	if (memcmp(Ztmp, Z, Zlen))
372 		goto err;
373 	rv = 1;
374 	err:
375 	if (key1)
376 		EC_KEY_free(key1);
377 	if (key2)
378 		EC_KEY_free(key2);
379 	free(Ztmp);
380 	if (rv)
381 		BIO_puts(out, " ok\n");
382 	else
383 		{
384 		fprintf(stderr, "Error in ECDH routines\n");
385 		ERR_print_errors_fp(stderr);
386 		}
387 	return rv;
388 	}
389 
390 #define test_ecdh_kat(bio, curve, bits) \
391 	ecdh_kat(bio, curve, NID_brainpoolP##bits##r1, \
392 		bp##bits##_da, sizeof(bp##bits##_da), \
393 		bp##bits##_db, sizeof(bp##bits##_db), \
394 		bp##bits##_Z, sizeof(bp##bits##_Z))
395 
396 int main(int argc, char *argv[])
397 	{
398 	BN_CTX *ctx=NULL;
399 	int ret=1;
400 	BIO *out;
401 
402 	out=BIO_new(BIO_s_file());
403 	if (out == NULL) exit(1);
404 	BIO_set_fp(out,stdout,BIO_NOCLOSE);
405 
406 	if ((ctx=BN_CTX_new()) == NULL) goto err;
407 
408 	/* NIST PRIME CURVES TESTS */
409 	if (!test_ecdh_curve(NID_X9_62_prime192v1, "NIST Prime-Curve P-192", ctx, out)) goto err;
410 	if (!test_ecdh_curve(NID_secp224r1, "NIST Prime-Curve P-224", ctx, out)) goto err;
411 	if (!test_ecdh_curve(NID_X9_62_prime256v1, "NIST Prime-Curve P-256", ctx, out)) goto err;
412 	if (!test_ecdh_curve(NID_secp384r1, "NIST Prime-Curve P-384", ctx, out)) goto err;
413 	if (!test_ecdh_curve(NID_secp521r1, "NIST Prime-Curve P-521", ctx, out)) goto err;
414 #ifndef OPENSSL_NO_EC2M
415 	/* NIST BINARY CURVES TESTS */
416 	if (!test_ecdh_curve(NID_sect163k1, "NIST Binary-Curve K-163", ctx, out)) goto err;
417 	if (!test_ecdh_curve(NID_sect163r2, "NIST Binary-Curve B-163", ctx, out)) goto err;
418 	if (!test_ecdh_curve(NID_sect233k1, "NIST Binary-Curve K-233", ctx, out)) goto err;
419 	if (!test_ecdh_curve(NID_sect233r1, "NIST Binary-Curve B-233", ctx, out)) goto err;
420 	if (!test_ecdh_curve(NID_sect283k1, "NIST Binary-Curve K-283", ctx, out)) goto err;
421 	if (!test_ecdh_curve(NID_sect283r1, "NIST Binary-Curve B-283", ctx, out)) goto err;
422 	if (!test_ecdh_curve(NID_sect409k1, "NIST Binary-Curve K-409", ctx, out)) goto err;
423 	if (!test_ecdh_curve(NID_sect409r1, "NIST Binary-Curve B-409", ctx, out)) goto err;
424 	if (!test_ecdh_curve(NID_sect571k1, "NIST Binary-Curve K-571", ctx, out)) goto err;
425 	if (!test_ecdh_curve(NID_sect571r1, "NIST Binary-Curve B-571", ctx, out)) goto err;
426 #endif
427 	if (!test_ecdh_kat(out, "Brainpool Prime-Curve brainpoolP256r1", 256))
428 		goto err;
429 	if (!test_ecdh_kat(out, "Brainpool Prime-Curve brainpoolP384r1", 384))
430 		goto err;
431 	if (!test_ecdh_kat(out, "Brainpool Prime-Curve brainpoolP512r1", 512))
432 		goto err;
433 
434 	ret = 0;
435 
436 err:
437 	ERR_print_errors_fp(stderr);
438 	if (ctx) BN_CTX_free(ctx);
439 	BIO_free(out);
440 	CRYPTO_cleanup_all_ex_data();
441 	ERR_remove_thread_state(NULL);
442 	CRYPTO_mem_leaks_fp(stderr);
443 	exit(ret);
444 	}
445