OLD | NEW |
| (Empty) |
1 /* unused */ | |
2 | |
3 /* crypto/bn/expspeed.c */ | |
4 /* Copyright (C) 1995-1998 Eric Young (eay@cryptsoft.com) | |
5 * All rights reserved. | |
6 * | |
7 * This package is an SSL implementation written | |
8 * by Eric Young (eay@cryptsoft.com). | |
9 * The implementation was written so as to conform with Netscapes SSL. | |
10 * | |
11 * This library is free for commercial and non-commercial use as long as | |
12 * the following conditions are aheared to. The following conditions | |
13 * apply to all code found in this distribution, be it the RC4, RSA, | |
14 * lhash, DES, etc., code; not just the SSL code. The SSL documentation | |
15 * included with this distribution is covered by the same copyright terms | |
16 * except that the holder is Tim Hudson (tjh@cryptsoft.com). | |
17 * | |
18 * Copyright remains Eric Young's, and as such any Copyright notices in | |
19 * the code are not to be removed. | |
20 * If this package is used in a product, Eric Young should be given attribution | |
21 * as the author of the parts of the library used. | |
22 * This can be in the form of a textual message at program startup or | |
23 * in documentation (online or textual) provided with the package. | |
24 * | |
25 * Redistribution and use in source and binary forms, with or without | |
26 * modification, are permitted provided that the following conditions | |
27 * are met: | |
28 * 1. Redistributions of source code must retain the copyright | |
29 * notice, this list of conditions and the following disclaimer. | |
30 * 2. Redistributions in binary form must reproduce the above copyright | |
31 * notice, this list of conditions and the following disclaimer in the | |
32 * documentation and/or other materials provided with the distribution. | |
33 * 3. All advertising materials mentioning features or use of this software | |
34 * must display the following acknowledgement: | |
35 * "This product includes cryptographic software written by | |
36 * Eric Young (eay@cryptsoft.com)" | |
37 * The word 'cryptographic' can be left out if the rouines from the library | |
38 * being used are not cryptographic related :-). | |
39 * 4. If you include any Windows specific code (or a derivative thereof) from | |
40 * the apps directory (application code) you must include an acknowledgement: | |
41 * "This product includes software written by Tim Hudson (tjh@cryptsoft.com)" | |
42 * | |
43 * THIS SOFTWARE IS PROVIDED BY ERIC YOUNG ``AS IS'' AND | |
44 * ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE | |
45 * IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE | |
46 * ARE DISCLAIMED. IN NO EVENT SHALL THE AUTHOR OR CONTRIBUTORS BE LIABLE | |
47 * FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL | |
48 * DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS | |
49 * OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) | |
50 * HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT | |
51 * LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY | |
52 * OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF | |
53 * SUCH DAMAGE. | |
54 * | |
55 * The licence and distribution terms for any publically available version or | |
56 * derivative of this code cannot be changed. i.e. this code cannot simply be | |
57 * copied and put under another distribution licence | |
58 * [including the GNU Public Licence.] | |
59 */ | |
60 | |
61 /* most of this code has been pilfered from my libdes speed.c program */ | |
62 | |
63 #define BASENUM 5000 | |
64 #define NUM_START 0 | |
65 | |
66 | |
67 /* determine timings for modexp, modmul, modsqr, gcd, Kronecker symbol, | |
68 * modular inverse, or modular square roots */ | |
69 #define TEST_EXP | |
70 #undef TEST_MUL | |
71 #undef TEST_SQR | |
72 #undef TEST_GCD | |
73 #undef TEST_KRON | |
74 #undef TEST_INV | |
75 #undef TEST_SQRT | |
76 #define P_MOD_64 9 /* least significant 6 bits for prime to be used for BN_sqrt
timings */ | |
77 | |
78 #if defined(TEST_EXP) + defined(TEST_MUL) + defined(TEST_SQR) + defined(TEST_GCD
) + defined(TEST_KRON) + defined(TEST_INV) +defined(TEST_SQRT) != 1 | |
79 # error "choose one test" | |
80 #endif | |
81 | |
82 #if defined(TEST_INV) || defined(TEST_SQRT) | |
83 # define C_PRIME | |
84 static void genprime_cb(int p, int n, void *arg); | |
85 #endif | |
86 | |
87 | |
88 | |
89 #undef PROG | |
90 #define PROG bnspeed_main | |
91 | |
92 #include <stdio.h> | |
93 #include <stdlib.h> | |
94 #include <signal.h> | |
95 #include <string.h> | |
96 #include <openssl/crypto.h> | |
97 #include <openssl/err.h> | |
98 #include <openssl/rand.h> | |
99 | |
100 #if !defined(OPENSSL_SYS_MSDOS) && (!defined(OPENSSL_SYS_VMS) || defined(__DECC)
) && !defined(OPENSSL_SYS_MACOSX) | |
101 #define TIMES | |
102 #endif | |
103 | |
104 #ifndef _IRIX | |
105 #include <time.h> | |
106 #endif | |
107 #ifdef TIMES | |
108 #include <sys/types.h> | |
109 #include <sys/times.h> | |
110 #endif | |
111 | |
112 /* Depending on the VMS version, the tms structure is perhaps defined. | |
113 The __TMS macro will show if it was. If it wasn't defined, we should | |
114 undefine TIMES, since that tells the rest of the program how things | |
115 should be handled. -- Richard Levitte */ | |
116 #if defined(OPENSSL_SYS_VMS_DECC) && !defined(__TMS) | |
117 #undef TIMES | |
118 #endif | |
119 | |
120 #ifndef TIMES | |
121 #include <sys/timeb.h> | |
122 #endif | |
123 | |
124 #if defined(sun) || defined(__ultrix) | |
125 #define _POSIX_SOURCE | |
126 #include <limits.h> | |
127 #include <sys/param.h> | |
128 #endif | |
129 | |
130 #include <openssl/bn.h> | |
131 #include <openssl/x509.h> | |
132 | |
133 /* The following if from times(3) man page. It may need to be changed */ | |
134 #ifndef HZ | |
135 # ifndef CLK_TCK | |
136 # ifndef _BSD_CLK_TCK_ /* FreeBSD hack */ | |
137 # define HZ 100.0 | |
138 # else /* _BSD_CLK_TCK_ */ | |
139 # define HZ ((double)_BSD_CLK_TCK_) | |
140 # endif | |
141 # else /* CLK_TCK */ | |
142 # define HZ ((double)CLK_TCK) | |
143 # endif | |
144 #endif | |
145 | |
146 #undef BUFSIZE | |
147 #define BUFSIZE ((long)1024*8) | |
148 int run=0; | |
149 | |
150 static double Time_F(int s); | |
151 #define START 0 | |
152 #define STOP 1 | |
153 | |
154 static double Time_F(int s) | |
155 { | |
156 double ret; | |
157 #ifdef TIMES | |
158 static struct tms tstart,tend; | |
159 | |
160 if (s == START) | |
161 { | |
162 times(&tstart); | |
163 return(0); | |
164 } | |
165 else | |
166 { | |
167 times(&tend); | |
168 ret=((double)(tend.tms_utime-tstart.tms_utime))/HZ; | |
169 return((ret < 1e-3)?1e-3:ret); | |
170 } | |
171 #else /* !times() */ | |
172 static struct timeb tstart,tend; | |
173 long i; | |
174 | |
175 if (s == START) | |
176 { | |
177 ftime(&tstart); | |
178 return(0); | |
179 } | |
180 else | |
181 { | |
182 ftime(&tend); | |
183 i=(long)tend.millitm-(long)tstart.millitm; | |
184 ret=((double)(tend.time-tstart.time))+((double)i)/1000.0; | |
185 return((ret < 0.001)?0.001:ret); | |
186 } | |
187 #endif | |
188 } | |
189 | |
190 #define NUM_SIZES 7 | |
191 #if NUM_START > NUM_SIZES | |
192 # error "NUM_START > NUM_SIZES" | |
193 #endif | |
194 static int sizes[NUM_SIZES]={128,256,512,1024,2048,4096,8192}; | |
195 static int mul_c[NUM_SIZES]={8*8*8*8*8*8,8*8*8*8*8,8*8*8*8,8*8*8,8*8,8,1}; | |
196 /*static int sizes[NUM_SIZES]={59,179,299,419,539}; */ | |
197 | |
198 #define RAND_SEED(string) { const char str[] = string; RAND_seed(string, sizeof
str); } | |
199 | |
200 void do_mul_exp(BIGNUM *r,BIGNUM *a,BIGNUM *b,BIGNUM *c,BN_CTX *ctx); | |
201 | |
202 int main(int argc, char **argv) | |
203 { | |
204 BN_CTX *ctx; | |
205 BIGNUM *a,*b,*c,*r; | |
206 | |
207 #if 1 | |
208 if (!CRYPTO_set_mem_debug_functions(0,0,0,0,0)) | |
209 abort(); | |
210 #endif | |
211 | |
212 ctx=BN_CTX_new(); | |
213 a=BN_new(); | |
214 b=BN_new(); | |
215 c=BN_new(); | |
216 r=BN_new(); | |
217 | |
218 while (!RAND_status()) | |
219 /* not enough bits */ | |
220 RAND_SEED("I demand a manual recount!"); | |
221 | |
222 do_mul_exp(r,a,b,c,ctx); | |
223 return 0; | |
224 } | |
225 | |
226 void do_mul_exp(BIGNUM *r, BIGNUM *a, BIGNUM *b, BIGNUM *c, BN_CTX *ctx) | |
227 { | |
228 int i,k; | |
229 double tm; | |
230 long num; | |
231 | |
232 num=BASENUM; | |
233 for (i=NUM_START; i<NUM_SIZES; i++) | |
234 { | |
235 #ifdef C_PRIME | |
236 # ifdef TEST_SQRT | |
237 if (!BN_set_word(a, 64)) goto err; | |
238 if (!BN_set_word(b, P_MOD_64)) goto err; | |
239 # define ADD a | |
240 # define REM b | |
241 # else | |
242 # define ADD NULL | |
243 # define REM NULL | |
244 # endif | |
245 if (!BN_generate_prime(c,sizes[i],0,ADD,REM,genprime_cb,NULL)) g
oto err; | |
246 putc('\n', stderr); | |
247 fflush(stderr); | |
248 #endif | |
249 | |
250 for (k=0; k<num; k++) | |
251 { | |
252 if (k%50 == 0) /* Average over num/50 different choices
of random numbers. */ | |
253 { | |
254 if (!BN_pseudo_rand(a,sizes[i],1,0)) goto err; | |
255 | |
256 if (!BN_pseudo_rand(b,sizes[i],1,0)) goto err; | |
257 | |
258 #ifndef C_PRIME | |
259 if (!BN_pseudo_rand(c,sizes[i],1,1)) goto err; | |
260 #endif | |
261 | |
262 #ifdef TEST_SQRT | |
263 if (!BN_mod_sqr(a,a,c,ctx)) goto err; | |
264 if (!BN_mod_sqr(b,b,c,ctx)) goto err; | |
265 #else | |
266 if (!BN_nnmod(a,a,c,ctx)) goto err; | |
267 if (!BN_nnmod(b,b,c,ctx)) goto err; | |
268 #endif | |
269 | |
270 if (k == 0) | |
271 Time_F(START); | |
272 } | |
273 | |
274 #if defined(TEST_EXP) | |
275 if (!BN_mod_exp(r,a,b,c,ctx)) goto err; | |
276 #elif defined(TEST_MUL) | |
277 { | |
278 int i = 0; | |
279 for (i = 0; i < 50; i++) | |
280 if (!BN_mod_mul(r,a,b,c,ctx)) goto err; | |
281 } | |
282 #elif defined(TEST_SQR) | |
283 { | |
284 int i = 0; | |
285 for (i = 0; i < 50; i++) | |
286 { | |
287 if (!BN_mod_sqr(r,a,c,ctx)) goto err; | |
288 if (!BN_mod_sqr(r,b,c,ctx)) goto err; | |
289 } | |
290 } | |
291 #elif defined(TEST_GCD) | |
292 if (!BN_gcd(r,a,b,ctx)) goto err; | |
293 if (!BN_gcd(r,b,c,ctx)) goto err; | |
294 if (!BN_gcd(r,c,a,ctx)) goto err; | |
295 #elif defined(TEST_KRON) | |
296 if (-2 == BN_kronecker(a,b,ctx)) goto err; | |
297 if (-2 == BN_kronecker(b,c,ctx)) goto err; | |
298 if (-2 == BN_kronecker(c,a,ctx)) goto err; | |
299 #elif defined(TEST_INV) | |
300 if (!BN_mod_inverse(r,a,c,ctx)) goto err; | |
301 if (!BN_mod_inverse(r,b,c,ctx)) goto err; | |
302 #else /* TEST_SQRT */ | |
303 if (!BN_mod_sqrt(r,a,c,ctx)) goto err; | |
304 if (!BN_mod_sqrt(r,b,c,ctx)) goto err; | |
305 #endif | |
306 } | |
307 tm=Time_F(STOP); | |
308 printf( | |
309 #if defined(TEST_EXP) | |
310 "modexp %4d ^ %4d %% %4d" | |
311 #elif defined(TEST_MUL) | |
312 "50*modmul %4d %4d %4d" | |
313 #elif defined(TEST_SQR) | |
314 "100*modsqr %4d %4d %4d" | |
315 #elif defined(TEST_GCD) | |
316 "3*gcd %4d %4d %4d" | |
317 #elif defined(TEST_KRON) | |
318 "3*kronecker %4d %4d %4d" | |
319 #elif defined(TEST_INV) | |
320 "2*inv %4d %4d mod %4d" | |
321 #else /* TEST_SQRT */ | |
322 "2*sqrt [prime == %d (mod 64)] %4d %4d mod %4d" | |
323 #endif | |
324 " -> %8.6fms %5.1f (%ld)\n", | |
325 #ifdef TEST_SQRT | |
326 P_MOD_64, | |
327 #endif | |
328 sizes[i],sizes[i],sizes[i],tm*1000.0/num,tm*mul_c[i]/num
, num); | |
329 num/=7; | |
330 if (num <= 0) num=1; | |
331 } | |
332 return; | |
333 | |
334 err: | |
335 ERR_print_errors_fp(stderr); | |
336 } | |
337 | |
338 | |
339 #ifdef C_PRIME | |
340 static void genprime_cb(int p, int n, void *arg) | |
341 { | |
342 char c='*'; | |
343 | |
344 if (p == 0) c='.'; | |
345 if (p == 1) c='+'; | |
346 if (p == 2) c='*'; | |
347 if (p == 3) c='\n'; | |
348 putc(c, stderr); | |
349 fflush(stderr); | |
350 (void)n; | |
351 (void)arg; | |
352 } | |
353 #endif | |
OLD | NEW |