diff options
author | miod <> | 2014-04-22 20:31:38 +0000 |
---|---|---|
committer | miod <> | 2014-04-22 20:31:38 +0000 |
commit | 1c26126a31b9d310c9fb3e33570fcbbe14676b11 (patch) | |
tree | 4598ce10bcfadfefbee541f6f4d2d7f25b1762a7 /src/lib/libcrypto/bn | |
parent | d9e10f75814922e77d3c507d0bed87de0e2095c6 (diff) | |
download | openbsd-1c26126a31b9d310c9fb3e33570fcbbe14676b11.tar.gz openbsd-1c26126a31b9d310c9fb3e33570fcbbe14676b11.tar.bz2 openbsd-1c26126a31b9d310c9fb3e33570fcbbe14676b11.zip |
Remove meat which either duplicates code found in apps/, or is only of value
for 20th century historians, and can be put in the Attic.
Diffstat (limited to 'src/lib/libcrypto/bn')
-rw-r--r-- | src/lib/libcrypto/bn/bn.mul | 19 | ||||
-rw-r--r-- | src/lib/libcrypto/bn/bnspeed.c | 233 | ||||
-rw-r--r-- | src/lib/libcrypto/bn/exp.c | 62 | ||||
-rw-r--r-- | src/lib/libcrypto/bn/expspeed.c | 347 |
4 files changed, 0 insertions, 661 deletions
diff --git a/src/lib/libcrypto/bn/bn.mul b/src/lib/libcrypto/bn/bn.mul deleted file mode 100644 index 9728870d38..0000000000 --- a/src/lib/libcrypto/bn/bn.mul +++ /dev/null | |||
@@ -1,19 +0,0 @@ | |||
1 | We need | ||
2 | |||
3 | * bn_mul_comba8 | ||
4 | * bn_mul_comba4 | ||
5 | * bn_mul_normal | ||
6 | * bn_mul_recursive | ||
7 | |||
8 | * bn_sqr_comba8 | ||
9 | * bn_sqr_comba4 | ||
10 | bn_sqr_normal -> BN_sqr | ||
11 | * bn_sqr_recursive | ||
12 | |||
13 | * bn_mul_low_recursive | ||
14 | * bn_mul_low_normal | ||
15 | * bn_mul_high | ||
16 | |||
17 | * bn_mul_part_recursive # symetric but not power of 2 | ||
18 | |||
19 | bn_mul_asymetric_recursive # uneven, but do the chop up. | ||
diff --git a/src/lib/libcrypto/bn/bnspeed.c b/src/lib/libcrypto/bn/bnspeed.c deleted file mode 100644 index b554ac8cf8..0000000000 --- a/src/lib/libcrypto/bn/bnspeed.c +++ /dev/null | |||
@@ -1,233 +0,0 @@ | |||
1 | /* unused */ | ||
2 | |||
3 | /* crypto/bn/bnspeed.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 1000000 | ||
64 | #undef PROG | ||
65 | #define PROG bnspeed_main | ||
66 | |||
67 | #include <stdio.h> | ||
68 | #include <stdlib.h> | ||
69 | #include <signal.h> | ||
70 | #include <string.h> | ||
71 | #include <openssl/crypto.h> | ||
72 | #include <openssl/err.h> | ||
73 | |||
74 | #if !defined(OPENSSL_SYS_MSDOS) && (!defined(OPENSSL_SYS_VMS) || defined(__DECC)) && !defined(OPENSSL_SYS_MACOSX) | ||
75 | #define TIMES | ||
76 | #endif | ||
77 | |||
78 | #ifndef _IRIX | ||
79 | #include <time.h> | ||
80 | #endif | ||
81 | #ifdef TIMES | ||
82 | #include <sys/types.h> | ||
83 | #include <sys/times.h> | ||
84 | #endif | ||
85 | |||
86 | /* Depending on the VMS version, the tms structure is perhaps defined. | ||
87 | The __TMS macro will show if it was. If it wasn't defined, we should | ||
88 | undefine TIMES, since that tells the rest of the program how things | ||
89 | should be handled. -- Richard Levitte */ | ||
90 | #if defined(OPENSSL_SYS_VMS_DECC) && !defined(__TMS) | ||
91 | #undef TIMES | ||
92 | #endif | ||
93 | |||
94 | #ifndef TIMES | ||
95 | #include <sys/timeb.h> | ||
96 | #endif | ||
97 | |||
98 | #if defined(sun) || defined(__ultrix) | ||
99 | #define _POSIX_SOURCE | ||
100 | #include <limits.h> | ||
101 | #include <sys/param.h> | ||
102 | #endif | ||
103 | |||
104 | #include <openssl/bn.h> | ||
105 | #include <openssl/x509.h> | ||
106 | |||
107 | /* The following if from times(3) man page. It may need to be changed */ | ||
108 | #ifndef HZ | ||
109 | # ifndef CLK_TCK | ||
110 | # ifndef _BSD_CLK_TCK_ /* FreeBSD hack */ | ||
111 | # define HZ 100.0 | ||
112 | # else /* _BSD_CLK_TCK_ */ | ||
113 | # define HZ ((double)_BSD_CLK_TCK_) | ||
114 | # endif | ||
115 | # else /* CLK_TCK */ | ||
116 | # define HZ ((double)CLK_TCK) | ||
117 | # endif | ||
118 | #endif | ||
119 | |||
120 | #undef BUFSIZE | ||
121 | #define BUFSIZE ((long)1024*8) | ||
122 | int run=0; | ||
123 | |||
124 | static double Time_F(int s); | ||
125 | #define START 0 | ||
126 | #define STOP 1 | ||
127 | |||
128 | static double Time_F(int s) | ||
129 | { | ||
130 | double ret; | ||
131 | #ifdef TIMES | ||
132 | static struct tms tstart,tend; | ||
133 | |||
134 | if (s == START) | ||
135 | { | ||
136 | times(&tstart); | ||
137 | return(0); | ||
138 | } | ||
139 | else | ||
140 | { | ||
141 | times(&tend); | ||
142 | ret=((double)(tend.tms_utime-tstart.tms_utime))/HZ; | ||
143 | return((ret < 1e-3)?1e-3:ret); | ||
144 | } | ||
145 | #else /* !times() */ | ||
146 | static struct timeb tstart,tend; | ||
147 | long i; | ||
148 | |||
149 | if (s == START) | ||
150 | { | ||
151 | ftime(&tstart); | ||
152 | return(0); | ||
153 | } | ||
154 | else | ||
155 | { | ||
156 | ftime(&tend); | ||
157 | i=(long)tend.millitm-(long)tstart.millitm; | ||
158 | ret=((double)(tend.time-tstart.time))+((double)i)/1000.0; | ||
159 | return((ret < 0.001)?0.001:ret); | ||
160 | } | ||
161 | #endif | ||
162 | } | ||
163 | |||
164 | #define NUM_SIZES 5 | ||
165 | static int sizes[NUM_SIZES]={128,256,512,1024,2048}; | ||
166 | /*static int sizes[NUM_SIZES]={59,179,299,419,539}; */ | ||
167 | |||
168 | void do_mul(BIGNUM *r,BIGNUM *a,BIGNUM *b,BN_CTX *ctx); | ||
169 | |||
170 | int main(int argc, char **argv) | ||
171 | { | ||
172 | BN_CTX *ctx; | ||
173 | BIGNUM a,b,c; | ||
174 | |||
175 | ctx=BN_CTX_new(); | ||
176 | BN_init(&a); | ||
177 | BN_init(&b); | ||
178 | BN_init(&c); | ||
179 | |||
180 | do_mul(&a,&b,&c,ctx); | ||
181 | } | ||
182 | |||
183 | void do_mul(BIGNUM *r, BIGNUM *a, BIGNUM *b, BN_CTX *ctx) | ||
184 | { | ||
185 | int i,j,k; | ||
186 | double tm; | ||
187 | long num; | ||
188 | |||
189 | for (i=0; i<NUM_SIZES; i++) | ||
190 | { | ||
191 | num=BASENUM; | ||
192 | if (i) num/=(i*3); | ||
193 | BN_rand(a,sizes[i],1,0); | ||
194 | for (j=i; j<NUM_SIZES; j++) | ||
195 | { | ||
196 | BN_rand(b,sizes[j],1,0); | ||
197 | Time_F(START); | ||
198 | for (k=0; k<num; k++) | ||
199 | BN_mul(r,b,a,ctx); | ||
200 | tm=Time_F(STOP); | ||
201 | printf("mul %4d x %4d -> %8.3fms\n",sizes[i],sizes[j],tm*1000.0/num); | ||
202 | } | ||
203 | } | ||
204 | |||
205 | for (i=0; i<NUM_SIZES; i++) | ||
206 | { | ||
207 | num=BASENUM; | ||
208 | if (i) num/=(i*3); | ||
209 | BN_rand(a,sizes[i],1,0); | ||
210 | Time_F(START); | ||
211 | for (k=0; k<num; k++) | ||
212 | BN_sqr(r,a,ctx); | ||
213 | tm=Time_F(STOP); | ||
214 | printf("sqr %4d x %4d -> %8.3fms\n",sizes[i],sizes[i],tm*1000.0/num); | ||
215 | } | ||
216 | |||
217 | for (i=0; i<NUM_SIZES; i++) | ||
218 | { | ||
219 | num=BASENUM/10; | ||
220 | if (i) num/=(i*3); | ||
221 | BN_rand(a,sizes[i]-1,1,0); | ||
222 | for (j=i; j<NUM_SIZES; j++) | ||
223 | { | ||
224 | BN_rand(b,sizes[j],1,0); | ||
225 | Time_F(START); | ||
226 | for (k=0; k<100000; k++) | ||
227 | BN_div(r, NULL, b, a,ctx); | ||
228 | tm=Time_F(STOP); | ||
229 | printf("div %4d / %4d -> %8.3fms\n",sizes[j],sizes[i]-1,tm*1000.0/num); | ||
230 | } | ||
231 | } | ||
232 | } | ||
233 | |||
diff --git a/src/lib/libcrypto/bn/exp.c b/src/lib/libcrypto/bn/exp.c deleted file mode 100644 index 4865b0ef74..0000000000 --- a/src/lib/libcrypto/bn/exp.c +++ /dev/null | |||
@@ -1,62 +0,0 @@ | |||
1 | /* unused */ | ||
2 | |||
3 | #include <stdio.h> | ||
4 | #include <openssl/tmdiff.h> | ||
5 | #include "bn_lcl.h" | ||
6 | |||
7 | #define SIZE 256 | ||
8 | #define NUM (8*8*8) | ||
9 | #define MOD (8*8*8*8*8) | ||
10 | |||
11 | main(argc,argv) | ||
12 | int argc; | ||
13 | char *argv[]; | ||
14 | { | ||
15 | BN_CTX ctx; | ||
16 | BIGNUM a,b,c,r,rr,t,l; | ||
17 | int j,i,size=SIZE,num=NUM,mod=MOD; | ||
18 | char *start,*end; | ||
19 | BN_MONT_CTX mont; | ||
20 | double d,md; | ||
21 | |||
22 | BN_MONT_CTX_init(&mont); | ||
23 | BN_CTX_init(&ctx); | ||
24 | BN_init(&a); | ||
25 | BN_init(&b); | ||
26 | BN_init(&c); | ||
27 | BN_init(&r); | ||
28 | |||
29 | start=ms_time_new(); | ||
30 | end=ms_time_new(); | ||
31 | while (size <= 1024*8) | ||
32 | { | ||
33 | BN_rand(&a,size,0,0); | ||
34 | BN_rand(&b,size,1,0); | ||
35 | BN_rand(&c,size,0,1); | ||
36 | |||
37 | BN_mod(&a,&a,&c,&ctx); | ||
38 | |||
39 | ms_time_get(start); | ||
40 | for (i=0; i<10; i++) | ||
41 | BN_MONT_CTX_set(&mont,&c,&ctx); | ||
42 | ms_time_get(end); | ||
43 | md=ms_time_diff(start,end); | ||
44 | |||
45 | ms_time_get(start); | ||
46 | for (i=0; i<num; i++) | ||
47 | { | ||
48 | /* bn_mull(&r,&a,&b,&ctx); */ | ||
49 | /* BN_sqr(&r,&a,&ctx); */ | ||
50 | BN_mod_exp_mont(&r,&a,&b,&c,&ctx,&mont); | ||
51 | } | ||
52 | ms_time_get(end); | ||
53 | d=ms_time_diff(start,end)/* *50/33 */; | ||
54 | printf("%5d bit:%6.2f %6d %6.4f %4d m_set(%5.4f)\n",size, | ||
55 | d,num,d/num,(int)((d/num)*mod),md/10.0); | ||
56 | num/=8; | ||
57 | mod/=8; | ||
58 | if (num <= 0) num=1; | ||
59 | size*=2; | ||
60 | } | ||
61 | |||
62 | } | ||
diff --git a/src/lib/libcrypto/bn/expspeed.c b/src/lib/libcrypto/bn/expspeed.c deleted file mode 100644 index 1b434db1ab..0000000000 --- a/src/lib/libcrypto/bn/expspeed.c +++ /dev/null | |||
@@ -1,347 +0,0 @@ | |||
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 | void do_mul_exp(BIGNUM *r,BIGNUM *a,BIGNUM *b,BIGNUM *c,BN_CTX *ctx); | ||
199 | |||
200 | int main(int argc, char **argv) | ||
201 | { | ||
202 | BN_CTX *ctx; | ||
203 | BIGNUM *a,*b,*c,*r; | ||
204 | |||
205 | #if 1 | ||
206 | if (!CRYPTO_set_mem_debug_functions(0,0,0,0,0)) | ||
207 | abort(); | ||
208 | #endif | ||
209 | |||
210 | ctx=BN_CTX_new(); | ||
211 | a=BN_new(); | ||
212 | b=BN_new(); | ||
213 | c=BN_new(); | ||
214 | r=BN_new(); | ||
215 | |||
216 | do_mul_exp(r,a,b,c,ctx); | ||
217 | return 0; | ||
218 | } | ||
219 | |||
220 | void do_mul_exp(BIGNUM *r, BIGNUM *a, BIGNUM *b, BIGNUM *c, BN_CTX *ctx) | ||
221 | { | ||
222 | int i,k; | ||
223 | double tm; | ||
224 | long num; | ||
225 | |||
226 | num=BASENUM; | ||
227 | for (i=NUM_START; i<NUM_SIZES; i++) | ||
228 | { | ||
229 | #ifdef C_PRIME | ||
230 | # ifdef TEST_SQRT | ||
231 | if (!BN_set_word(a, 64)) goto err; | ||
232 | if (!BN_set_word(b, P_MOD_64)) goto err; | ||
233 | # define ADD a | ||
234 | # define REM b | ||
235 | # else | ||
236 | # define ADD NULL | ||
237 | # define REM NULL | ||
238 | # endif | ||
239 | if (!BN_generate_prime(c,sizes[i],0,ADD,REM,genprime_cb,NULL)) goto err; | ||
240 | putc('\n', stderr); | ||
241 | fflush(stderr); | ||
242 | #endif | ||
243 | |||
244 | for (k=0; k<num; k++) | ||
245 | { | ||
246 | if (k%50 == 0) /* Average over num/50 different choices of random numbers. */ | ||
247 | { | ||
248 | if (!BN_pseudo_rand(a,sizes[i],1,0)) goto err; | ||
249 | |||
250 | if (!BN_pseudo_rand(b,sizes[i],1,0)) goto err; | ||
251 | |||
252 | #ifndef C_PRIME | ||
253 | if (!BN_pseudo_rand(c,sizes[i],1,1)) goto err; | ||
254 | #endif | ||
255 | |||
256 | #ifdef TEST_SQRT | ||
257 | if (!BN_mod_sqr(a,a,c,ctx)) goto err; | ||
258 | if (!BN_mod_sqr(b,b,c,ctx)) goto err; | ||
259 | #else | ||
260 | if (!BN_nnmod(a,a,c,ctx)) goto err; | ||
261 | if (!BN_nnmod(b,b,c,ctx)) goto err; | ||
262 | #endif | ||
263 | |||
264 | if (k == 0) | ||
265 | Time_F(START); | ||
266 | } | ||
267 | |||
268 | #if defined(TEST_EXP) | ||
269 | if (!BN_mod_exp(r,a,b,c,ctx)) goto err; | ||
270 | #elif defined(TEST_MUL) | ||
271 | { | ||
272 | int i = 0; | ||
273 | for (i = 0; i < 50; i++) | ||
274 | if (!BN_mod_mul(r,a,b,c,ctx)) goto err; | ||
275 | } | ||
276 | #elif defined(TEST_SQR) | ||
277 | { | ||
278 | int i = 0; | ||
279 | for (i = 0; i < 50; i++) | ||
280 | { | ||
281 | if (!BN_mod_sqr(r,a,c,ctx)) goto err; | ||
282 | if (!BN_mod_sqr(r,b,c,ctx)) goto err; | ||
283 | } | ||
284 | } | ||
285 | #elif defined(TEST_GCD) | ||
286 | if (!BN_gcd(r,a,b,ctx)) goto err; | ||
287 | if (!BN_gcd(r,b,c,ctx)) goto err; | ||
288 | if (!BN_gcd(r,c,a,ctx)) goto err; | ||
289 | #elif defined(TEST_KRON) | ||
290 | if (-2 == BN_kronecker(a,b,ctx)) goto err; | ||
291 | if (-2 == BN_kronecker(b,c,ctx)) goto err; | ||
292 | if (-2 == BN_kronecker(c,a,ctx)) goto err; | ||
293 | #elif defined(TEST_INV) | ||
294 | if (!BN_mod_inverse(r,a,c,ctx)) goto err; | ||
295 | if (!BN_mod_inverse(r,b,c,ctx)) goto err; | ||
296 | #else /* TEST_SQRT */ | ||
297 | if (!BN_mod_sqrt(r,a,c,ctx)) goto err; | ||
298 | if (!BN_mod_sqrt(r,b,c,ctx)) goto err; | ||
299 | #endif | ||
300 | } | ||
301 | tm=Time_F(STOP); | ||
302 | printf( | ||
303 | #if defined(TEST_EXP) | ||
304 | "modexp %4d ^ %4d %% %4d" | ||
305 | #elif defined(TEST_MUL) | ||
306 | "50*modmul %4d %4d %4d" | ||
307 | #elif defined(TEST_SQR) | ||
308 | "100*modsqr %4d %4d %4d" | ||
309 | #elif defined(TEST_GCD) | ||
310 | "3*gcd %4d %4d %4d" | ||
311 | #elif defined(TEST_KRON) | ||
312 | "3*kronecker %4d %4d %4d" | ||
313 | #elif defined(TEST_INV) | ||
314 | "2*inv %4d %4d mod %4d" | ||
315 | #else /* TEST_SQRT */ | ||
316 | "2*sqrt [prime == %d (mod 64)] %4d %4d mod %4d" | ||
317 | #endif | ||
318 | " -> %8.6fms %5.1f (%ld)\n", | ||
319 | #ifdef TEST_SQRT | ||
320 | P_MOD_64, | ||
321 | #endif | ||
322 | sizes[i],sizes[i],sizes[i],tm*1000.0/num,tm*mul_c[i]/num, num); | ||
323 | num/=7; | ||
324 | if (num <= 0) num=1; | ||
325 | } | ||
326 | return; | ||
327 | |||
328 | err: | ||
329 | ERR_print_errors_fp(stderr); | ||
330 | } | ||
331 | |||
332 | |||
333 | #ifdef C_PRIME | ||
334 | static void genprime_cb(int p, int n, void *arg) | ||
335 | { | ||
336 | char c='*'; | ||
337 | |||
338 | if (p == 0) c='.'; | ||
339 | if (p == 1) c='+'; | ||
340 | if (p == 2) c='*'; | ||
341 | if (p == 3) c='\n'; | ||
342 | putc(c, stderr); | ||
343 | fflush(stderr); | ||
344 | (void)n; | ||
345 | (void)arg; | ||
346 | } | ||
347 | #endif | ||