1:
2:
3:
4:
5:
6:
7:
8:
9:
10:
11:
12:
13:
14:
15:
16:
17:
18:
19:
20:
21:
22:
23:
24:
25:
26:
27:
28:
29:
30:
31:
32:
33:
34:
35:
36:
37:
38:
39:
40:
41:
42:
43:
44:
45:
46:
47:
48:
49:
50:
51:
52:
53:
54:
55:
56:
57:
58:
59: #include <stdio.h>
60: #include <stdlib.h>
61: #include <string.h>
62:
63: #include "../e_os.h"
64:
65: #include <openssl/bio.h>
66: #include <openssl/bn.h>
67: #include <openssl/rand.h>
68: #include <openssl/err.h>
69:
70: #define NUM_BITS (BN_BITS*2)
71:
72: static const char rnd_seed[] = "string to make the random number generator think it has entropy";
73:
74: int main(int argc, char *argv[])
75: {
76: BN_CTX *ctx;
77: BIO *out=NULL;
78: int i,ret;
79: unsigned char c;
80: BIGNUM *r_mont,*r_mont_const,*r_recp,*r_simple,*a,*b,*m;
81:
82: RAND_seed(rnd_seed, sizeof rnd_seed);
83:
84:
85:
86: ERR_load_BN_strings();
87:
88: ctx=BN_CTX_new();
89: if (ctx == NULL) EXIT(1);
90: r_mont=BN_new();
91: r_mont_const=BN_new();
92: r_recp=BN_new();
93: r_simple=BN_new();
94: a=BN_new();
95: b=BN_new();
96: m=BN_new();
97: if ( (r_mont == NULL) || (r_recp == NULL) ||
98: (a == NULL) || (b == NULL))
99: goto err;
100:
101: out=BIO_new(BIO_s_file());
102:
103: if (out == NULL) EXIT(1);
104: BIO_set_fp(out,stdout,BIO_NOCLOSE);
105:
106: for (i=0; i<200; i++)
107: {
108: RAND_bytes(&c,1);
109: c=(c%BN_BITS)-BN_BITS2;
110: BN_rand(a,NUM_BITS+c,0,0);
111:
112: RAND_bytes(&c,1);
113: c=(c%BN_BITS)-BN_BITS2;
114: BN_rand(b,NUM_BITS+c,0,0);
115:
116: RAND_bytes(&c,1);
117: c=(c%BN_BITS)-BN_BITS2;
118: BN_rand(m,NUM_BITS+c,0,1);
119:
120: BN_mod(a,a,m,ctx);
121: BN_mod(b,b,m,ctx);
122:
123: ret=BN_mod_exp_mont(r_mont,a,b,m,ctx,NULL);
124: if (ret <= 0)
125: {
126: printf("BN_mod_exp_mont() problems\n");
127: ERR_print_errors(out);
128: EXIT(1);
129: }
130:
131: ret=BN_mod_exp_recp(r_recp,a,b,m,ctx);
132: if (ret <= 0)
133: {
134: printf("BN_mod_exp_recp() problems\n");
135: ERR_print_errors(out);
136: EXIT(1);
137: }
138:
139: ret=BN_mod_exp_simple(r_simple,a,b,m,ctx);
140: if (ret <= 0)
141: {
142: printf("BN_mod_exp_simple() problems\n");
143: ERR_print_errors(out);
144: EXIT(1);
145: }
146:
147: ret=BN_mod_exp_mont_consttime(r_mont_const,a,b,m,ctx,NULL);
148: if (ret <= 0)
149: {
150: printf("BN_mod_exp_mont_consttime() problems\n");
151: ERR_print_errors(out);
152: EXIT(1);
153: }
154:
155: if (BN_cmp(r_simple, r_mont) == 0
156: && BN_cmp(r_simple,r_recp) == 0
157: && BN_cmp(r_simple,r_mont_const) == 0)
158: {
159: printf(".");
160: fflush(stdout);
161: }
162: else
163: {
164: if (BN_cmp(r_simple,r_mont) != 0)
165: printf("\nsimple and mont results differ\n");
166: if (BN_cmp(r_simple,r_mont) != 0)
167: printf("\nsimple and mont const time results differ\n");
168: if (BN_cmp(r_simple,r_recp) != 0)
169: printf("\nsimple and recp results differ\n");
170:
171: printf("a (%3d) = ",BN_num_bits(a)); BN_print(out,a);
172: printf("\nb (%3d) = ",BN_num_bits(b)); BN_print(out,b);
173: printf("\nm (%3d) = ",BN_num_bits(m)); BN_print(out,m);
174: printf("\nsimple ="); BN_print(out,r_simple);
175: printf("\nrecp ="); BN_print(out,r_recp);
176: printf("\nmont ="); BN_print(out,r_mont);
177: printf("\nmont_ct ="); BN_print(out,r_mont_const);
178: printf("\n");
179: EXIT(1);
180: }
181: }
182: BN_free(r_mont);
183: BN_free(r_mont_const);
184: BN_free(r_recp);
185: BN_free(r_simple);
186: BN_free(a);
187: BN_free(b);
188: BN_free(m);
189: BN_CTX_free(ctx);
190: ERR_remove_state(0);
191: CRYPTO_mem_leaks(out);
192: BIO_free(out);
193: printf(" done\n");
194: EXIT(0);
195: err:
196: ERR_load_crypto_strings();
197: ERR_print_errors(out);
198: #ifdef OPENSSL_SYS_NETWARE
199: printf("ERROR\n");
200: #endif
201: EXIT(1);
202: return(1);
203: }
204: