util/find-docs-nits: Recognise SPARSE_ARRAY_OF
[oweals/openssl.git] / test / exptest.c
1 /*
2  * Copyright 1995-2017 The OpenSSL Project Authors. All Rights Reserved.
3  *
4  * Licensed under the Apache License 2.0 (the "License").  You may not use
5  * this file except in compliance with the License.  You can obtain a copy
6  * in the file LICENSE in the source distribution or at
7  * https://www.openssl.org/source/license.html
8  */
9
10 #include <stdio.h>
11 #include <stdlib.h>
12 #include <string.h>
13
14 #include "internal/nelem.h"
15
16 #include <openssl/bio.h>
17 #include <openssl/bn.h>
18 #include <openssl/rand.h>
19 #include <openssl/err.h>
20
21 #include "testutil.h"
22
23 #define NUM_BITS        (BN_BITS2 * 4)
24
25 #define BN_print_var(v) test_output_bignum(#v, v)
26
27 /*
28  * Test that r == 0 in test_exp_mod_zero(). Returns one on success,
29  * returns zero and prints debug output otherwise.
30  */
31 static int a_is_zero_mod_one(const char *method, const BIGNUM *r,
32                              const BIGNUM *a)
33 {
34     if (!BN_is_zero(r)) {
35         TEST_error("%s failed: a ** 0 mod 1 = r (should be 0)", method);
36         BN_print_var(a);
37         BN_print_var(r);
38         return 0;
39     }
40     return 1;
41 }
42
43 /*
44  * test_mod_exp_zero tests that x**0 mod 1 == 0. It returns zero on success.
45  */
46 static int test_mod_exp_zero(void)
47 {
48     BIGNUM *a = NULL, *p = NULL, *m = NULL;
49     BIGNUM *r = NULL;
50     BN_ULONG one_word = 1;
51     BN_CTX *ctx = BN_CTX_new();
52     int ret = 1, failed = 0;
53
54     if (!TEST_ptr(m = BN_new())
55         || !TEST_ptr(a = BN_new())
56         || !TEST_ptr(p = BN_new())
57         || !TEST_ptr(r = BN_new()))
58         goto err;
59
60     BN_one(m);
61     BN_one(a);
62     BN_zero(p);
63
64     if (!TEST_true(BN_rand(a, 1024, BN_RAND_TOP_ONE, BN_RAND_BOTTOM_ANY)))
65         goto err;
66
67     if (!TEST_true(BN_mod_exp(r, a, p, m, ctx)))
68         goto err;
69
70     if (!TEST_true(a_is_zero_mod_one("BN_mod_exp", r, a)))
71         failed = 1;
72
73     if (!TEST_true(BN_mod_exp_recp(r, a, p, m, ctx)))
74         goto err;
75
76     if (!TEST_true(a_is_zero_mod_one("BN_mod_exp_recp", r, a)))
77         failed = 1;
78
79     if (!TEST_true(BN_mod_exp_simple(r, a, p, m, ctx)))
80         goto err;
81
82     if (!TEST_true(a_is_zero_mod_one("BN_mod_exp_simple", r, a)))
83         failed = 1;
84
85     if (!TEST_true(BN_mod_exp_mont(r, a, p, m, ctx, NULL)))
86         goto err;
87
88     if (!TEST_true(a_is_zero_mod_one("BN_mod_exp_mont", r, a)))
89         failed = 1;
90
91     if (!TEST_true(BN_mod_exp_mont_consttime(r, a, p, m, ctx, NULL)))
92         goto err;
93
94     if (!TEST_true(a_is_zero_mod_one("BN_mod_exp_mont_consttime", r, a)))
95         failed = 1;
96
97     /*
98      * A different codepath exists for single word multiplication
99      * in non-constant-time only.
100      */
101     if (!TEST_true(BN_mod_exp_mont_word(r, one_word, p, m, ctx, NULL)))
102         goto err;
103
104     if (!TEST_BN_eq_zero(r)) {
105         TEST_error("BN_mod_exp_mont_word failed: "
106                    "1 ** 0 mod 1 = r (should be 0)");
107         BN_print_var(r);
108         goto err;
109     }
110
111     ret = !failed;
112  err:
113     BN_free(r);
114     BN_free(a);
115     BN_free(p);
116     BN_free(m);
117     BN_CTX_free(ctx);
118
119     return ret;
120 }
121
122 static int test_mod_exp(int round)
123 {
124     BN_CTX *ctx;
125     unsigned char c;
126     int ret = 0;
127     BIGNUM *r_mont = NULL;
128     BIGNUM *r_mont_const = NULL;
129     BIGNUM *r_recp = NULL;
130     BIGNUM *r_simple = NULL;
131     BIGNUM *a = NULL;
132     BIGNUM *b = NULL;
133     BIGNUM *m = NULL;
134
135     if (!TEST_ptr(ctx = BN_CTX_new()))
136         goto err;
137
138     if (!TEST_ptr(r_mont = BN_new())
139         || !TEST_ptr(r_mont_const = BN_new())
140         || !TEST_ptr(r_recp = BN_new())
141         || !TEST_ptr(r_simple = BN_new())
142         || !TEST_ptr(a = BN_new())
143         || !TEST_ptr(b = BN_new())
144         || !TEST_ptr(m = BN_new()))
145         goto err;
146
147     RAND_bytes(&c, 1);
148     c = (c % BN_BITS) - BN_BITS2;
149     BN_rand(a, NUM_BITS + c, BN_RAND_TOP_ONE, BN_RAND_BOTTOM_ANY);
150
151     RAND_bytes(&c, 1);
152     c = (c % BN_BITS) - BN_BITS2;
153     BN_rand(b, NUM_BITS + c, BN_RAND_TOP_ONE, BN_RAND_BOTTOM_ANY);
154
155     RAND_bytes(&c, 1);
156     c = (c % BN_BITS) - BN_BITS2;
157     BN_rand(m, NUM_BITS + c, BN_RAND_TOP_ONE, BN_RAND_BOTTOM_ODD);
158
159     if (!TEST_true(BN_mod(a, a, m, ctx))
160         || !TEST_true(BN_mod(b, b, m, ctx))
161         || !TEST_true(BN_mod_exp_mont(r_mont, a, b, m, ctx, NULL))
162         || !TEST_true(BN_mod_exp_recp(r_recp, a, b, m, ctx))
163         || !TEST_true(BN_mod_exp_simple(r_simple, a, b, m, ctx))
164         || !TEST_true(BN_mod_exp_mont_consttime(r_mont_const, a, b, m, ctx, NULL)))
165         goto err;
166
167     if (!TEST_BN_eq(r_simple, r_mont)
168         || !TEST_BN_eq(r_simple, r_recp)
169         || !TEST_BN_eq(r_simple, r_mont_const)) {
170         if (BN_cmp(r_simple, r_mont) != 0)
171             TEST_info("simple and mont results differ");
172         if (BN_cmp(r_simple, r_mont_const) != 0)
173             TEST_info("simple and mont const time results differ");
174         if (BN_cmp(r_simple, r_recp) != 0)
175             TEST_info("simple and recp results differ");
176
177         BN_print_var(a);
178         BN_print_var(b);
179         BN_print_var(m);
180         BN_print_var(r_simple);
181         BN_print_var(r_recp);
182         BN_print_var(r_mont);
183         BN_print_var(r_mont_const);
184         goto err;
185     }
186
187     ret = 1;
188  err:
189     BN_free(r_mont);
190     BN_free(r_mont_const);
191     BN_free(r_recp);
192     BN_free(r_simple);
193     BN_free(a);
194     BN_free(b);
195     BN_free(m);
196     BN_CTX_free(ctx);
197
198     return ret;
199 }
200
201 int setup_tests(void)
202 {
203     ADD_TEST(test_mod_exp_zero);
204     ADD_ALL_TESTS(test_mod_exp, 200);
205     return 1;
206 }