xref: /openbsd-src/regress/lib/libcrypto/bn/bn_primes.c (revision 472a54ea59ddef54c0ec79bcc032e3253f2b4d50)
1*472a54eaStb /*	$OpenBSD: bn_primes.c,v 1.3 2023/04/25 15:30:03 tb Exp $ */
2b8d22d11Stb /*
3b8d22d11Stb  * Copyright (c) 2022 Theo Buehler <tb@openbsd.org>
4b8d22d11Stb  *
5b8d22d11Stb  * Permission to use, copy, modify, and distribute this software for any
6b8d22d11Stb  * purpose with or without fee is hereby granted, provided that the above
7b8d22d11Stb  * copyright notice and this permission notice appear in all copies.
8b8d22d11Stb  *
9b8d22d11Stb  * THE SOFTWARE IS PROVIDED "AS IS" AND THE AUTHOR DISCLAIMS ALL WARRANTIES
10b8d22d11Stb  * WITH REGARD TO THIS SOFTWARE INCLUDING ALL IMPLIED WARRANTIES OF
11b8d22d11Stb  * MERCHANTABILITY AND FITNESS. IN NO EVENT SHALL THE AUTHOR BE LIABLE FOR
12b8d22d11Stb  * ANY SPECIAL, DIRECT, INDIRECT, OR CONSEQUENTIAL DAMAGES OR ANY DAMAGES
13b8d22d11Stb  * WHATSOEVER RESULTING FROM LOSS OF USE, DATA OR PROFITS, WHETHER IN AN
14b8d22d11Stb  * ACTION OF CONTRACT, NEGLIGENCE OR OTHER TORTIOUS ACTION, ARISING OUT OF
15b8d22d11Stb  * OR IN CONNECTION WITH THE USE OR PERFORMANCE OF THIS SOFTWARE.
16b8d22d11Stb  */
17b8d22d11Stb 
1820c8b67dStb #include <stdint.h>
1920c8b67dStb #include <stdio.h>
2020c8b67dStb #include <stdlib.h>
2120c8b67dStb 
22b8d22d11Stb #include <openssl/bn.h>
23b8d22d11Stb 
24b8d22d11Stb #include "bn_prime.h"
25b8d22d11Stb 
26b8d22d11Stb static int
test_bn_is_prime_fasttest(int do_trial_division)27b8d22d11Stb test_bn_is_prime_fasttest(int do_trial_division)
28b8d22d11Stb {
29b8d22d11Stb 	BIGNUM *n = NULL;
30b8d22d11Stb 	char *descr = NULL;
31b8d22d11Stb 	uint16_t i, j, max;
32b8d22d11Stb 	int is_prime, ret;
33b8d22d11Stb 	int failed = 1;
34b8d22d11Stb 
35b8d22d11Stb 	if (asprintf(&descr, "with%s trial divisions",
36b8d22d11Stb 	    do_trial_division ? "" : "out") == -1) {
37b8d22d11Stb 		descr = NULL;
38b8d22d11Stb 		fprintf(stderr, "asprintf failed\n");
39b8d22d11Stb 		goto err;
40b8d22d11Stb 	}
41b8d22d11Stb 
42b8d22d11Stb 	if ((n = BN_new()) == NULL) {
43b8d22d11Stb 		fprintf(stderr, "BN_new failed\n");
44b8d22d11Stb 		goto err;
45b8d22d11Stb 	}
46b8d22d11Stb 
47b8d22d11Stb 	max = primes[NUMPRIMES - 1] + 1;
48b8d22d11Stb 
49b8d22d11Stb 	failed = 0;
50b8d22d11Stb 	for (i = 1, j = 0; i < max && j < NUMPRIMES; i++) {
51b8d22d11Stb 		if (!BN_set_word(n, i)) {
52b8d22d11Stb 			fprintf(stderr, "BN_set_word(%d) failed", i);
53b8d22d11Stb 			failed = 1;
54b8d22d11Stb 			goto err;
55b8d22d11Stb 		}
56b8d22d11Stb 
57b8d22d11Stb 		is_prime = i == primes[j];
58b8d22d11Stb 		if (is_prime)
59b8d22d11Stb 			j++;
60b8d22d11Stb 
61b8d22d11Stb 		ret = BN_is_prime_fasttest_ex(n, BN_prime_checks, NULL,
62b8d22d11Stb 		    do_trial_division, NULL);
63b8d22d11Stb 		if (ret != is_prime) {
64b8d22d11Stb 			fprintf(stderr,
65b8d22d11Stb 			    "BN_is_prime_fasttest_ex(%d) %s: want %d, got %d\n",
66b8d22d11Stb 			    i, descr, is_prime, ret);
67b8d22d11Stb 			failed = 1;
68b8d22d11Stb 		}
69b8d22d11Stb 	}
70b8d22d11Stb 
71b8d22d11Stb 	if (i < max || j < NUMPRIMES) {
72b8d22d11Stb 		fprintf(stderr, "%s: %d < %d or %d < %d\n", descr, i, max, j,
73b8d22d11Stb 		    NUMPRIMES);
74b8d22d11Stb 		failed = 1;
75b8d22d11Stb 	}
76b8d22d11Stb 
77b8d22d11Stb  err:
78b8d22d11Stb 	BN_free(n);
79b8d22d11Stb 	free(descr);
80b8d22d11Stb 	return failed;
81b8d22d11Stb }
82b8d22d11Stb 
83b8d22d11Stb #define BN_PRIME_FN_INIT(a) { .fn = a, .name = #a }
84b8d22d11Stb 
85b8d22d11Stb static const struct test_dynamic_api {
86b8d22d11Stb 	BIGNUM *(*fn)(BIGNUM *);
87b8d22d11Stb 	const char *name;
88b8d22d11Stb } dynamic_api_data[] = {
89b8d22d11Stb 	BN_PRIME_FN_INIT(BN_get_rfc2409_prime_1024),
90b8d22d11Stb 	BN_PRIME_FN_INIT(BN_get_rfc2409_prime_768),
91b8d22d11Stb 	BN_PRIME_FN_INIT(BN_get_rfc3526_prime_1536),
92b8d22d11Stb 	BN_PRIME_FN_INIT(BN_get_rfc3526_prime_2048),
93b8d22d11Stb 	BN_PRIME_FN_INIT(BN_get_rfc3526_prime_3072),
94b8d22d11Stb 	BN_PRIME_FN_INIT(BN_get_rfc3526_prime_4096),
95b8d22d11Stb 	BN_PRIME_FN_INIT(BN_get_rfc3526_prime_6144),
96b8d22d11Stb 	BN_PRIME_FN_INIT(BN_get_rfc3526_prime_8192),
97b8d22d11Stb };
98b8d22d11Stb 
99b8d22d11Stb #define N_DYNAMIC_TESTS (sizeof(dynamic_api_data) / sizeof(dynamic_api_data[0]))
100b8d22d11Stb 
101b8d22d11Stb static int
test_prime_dynamic_api(const struct test_dynamic_api * tc)102b8d22d11Stb test_prime_dynamic_api(const struct test_dynamic_api *tc)
103b8d22d11Stb {
104b8d22d11Stb 	BIGNUM *prime;
105b8d22d11Stb 	int ret;
106b8d22d11Stb 	int failed = 1;
107b8d22d11Stb 
108b8d22d11Stb 	if ((prime = tc->fn(NULL)) == NULL) {
109b8d22d11Stb 		fprintf(stderr, "%s failed\n", tc->name);
110b8d22d11Stb 		goto err;
111b8d22d11Stb 	}
112b8d22d11Stb 
113b8d22d11Stb 	if ((ret = BN_is_prime_fasttest_ex(prime, 1, NULL, 1, NULL)) != 1) {
114b8d22d11Stb 		fprintf(stderr, "%s: %s want 1, got %d\n", tc->name,
115b8d22d11Stb 		    "BN_is_prime_fasttest_ex", ret);
116b8d22d11Stb 		goto err;
117b8d22d11Stb 	}
118b8d22d11Stb 
119b8d22d11Stb 	failed = 0;
120b8d22d11Stb 
121b8d22d11Stb  err:
122b8d22d11Stb 	BN_free(prime);
123b8d22d11Stb 	return failed;
124b8d22d11Stb }
125b8d22d11Stb 
126b8d22d11Stb static int
test_prime_constants(void)127b8d22d11Stb test_prime_constants(void)
128b8d22d11Stb {
129b8d22d11Stb 	size_t i;
130b8d22d11Stb 	int failed = 0;
131b8d22d11Stb 
132b8d22d11Stb 	for (i = 0; i < N_DYNAMIC_TESTS; i++)
133b8d22d11Stb 		failed |= test_prime_dynamic_api(&dynamic_api_data[i]);
134b8d22d11Stb 
135b8d22d11Stb 	return failed;
136b8d22d11Stb }
137b8d22d11Stb 
138b8d22d11Stb int
main(void)139b8d22d11Stb main(void)
140b8d22d11Stb {
141b8d22d11Stb 	int failed = 0;
142b8d22d11Stb 
143b8d22d11Stb 	failed |= test_bn_is_prime_fasttest(0);
144b8d22d11Stb 	failed |= test_bn_is_prime_fasttest(1);
145b8d22d11Stb 	failed |= test_prime_constants();
146b8d22d11Stb 
147b8d22d11Stb 	return failed;
148b8d22d11Stb }
149