1 /* unused */ 2 3 /* crypto/bn/expspeed.c */ 4 /* Copyright (C) 1995-1998 Eric Young (eay (at) cryptsoft.com) 5 * All rights reserved. 6 * 7 * This package is an SSL implementation written 8 * by Eric Young (eay (at) 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 (at) 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 (at) 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 (at) 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)) goto 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 354
