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