See ChangeLog: Tue May 4 15:49:29 CEST 1999 Werner Koch
[gnupg.git] / cipher / rmd160.c
1 /* rmd160.c  -  RIPE-MD160
2  *      Copyright (C) 1998 Free Software Foundation, Inc.
3  *
4  * This file is part of GnuPG.
5  *
6  * GnuPG is free software; you can redistribute it and/or modify
7  * it under the terms of the GNU General Public License as published by
8  * the Free Software Foundation; either version 2 of the License, or
9  * (at your option) any later version.
10  *
11  * GnuPG is distributed in the hope that it will be useful,
12  * but WITHOUT ANY WARRANTY; without even the implied warranty of
13  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
14  * GNU General Public License for more details.
15  *
16  * You should have received a copy of the GNU General Public License
17  * along with this program; if not, write to the Free Software
18  * Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA 02111-1307, USA
19  */
20
21 #include <config.h>
22 #include <stdio.h>
23 #include <stdlib.h>
24 #include <string.h>
25 #include <assert.h>
26 #include "util.h"
27 #include "memory.h"
28 #include "rmd.h"
29 #include "cipher.h" /* only used for the rmd160_hash_buffer() prototype */
30 #include "dynload.h"
31
32 /*********************************
33  * RIPEMD-160 is not patented, see (as of 25.10.97)
34  *   http://www.esat.kuleuven.ac.be/~bosselae/ripemd160.html
35  * Note that the code uses Little Endian byteorder, which is good for
36  * 386 etc, but we must add some conversion when used on a big endian box.
37  *
38  *
39  * Pseudo-code for RIPEMD-160
40  *
41  * RIPEMD-160 is an iterative hash function that operates on 32-bit words.
42  * The round function takes as input a 5-word chaining variable and a 16-word
43  * message block and maps this to a new chaining variable. All operations are
44  * defined on 32-bit words. Padding is identical to that of MD4.
45  *
46  *
47  * RIPEMD-160: definitions
48  *
49  *
50  *   nonlinear functions at bit level: exor, mux, -, mux, -
51  *
52  *   f(j, x, y, z) = x XOR y XOR z                (0 <= j <= 15)
53  *   f(j, x, y, z) = (x AND y) OR (NOT(x) AND z)  (16 <= j <= 31)
54  *   f(j, x, y, z) = (x OR NOT(y)) XOR z          (32 <= j <= 47)
55  *   f(j, x, y, z) = (x AND z) OR (y AND NOT(z))  (48 <= j <= 63)
56  *   f(j, x, y, z) = x XOR (y OR NOT(z))          (64 <= j <= 79)
57  *
58  *
59  *   added constants (hexadecimal)
60  *
61  *   K(j) = 0x00000000      (0 <= j <= 15)
62  *   K(j) = 0x5A827999     (16 <= j <= 31)      int(2**30 x sqrt(2))
63  *   K(j) = 0x6ED9EBA1     (32 <= j <= 47)      int(2**30 x sqrt(3))
64  *   K(j) = 0x8F1BBCDC     (48 <= j <= 63)      int(2**30 x sqrt(5))
65  *   K(j) = 0xA953FD4E     (64 <= j <= 79)      int(2**30 x sqrt(7))
66  *   K'(j) = 0x50A28BE6     (0 <= j <= 15)      int(2**30 x cbrt(2))
67  *   K'(j) = 0x5C4DD124    (16 <= j <= 31)      int(2**30 x cbrt(3))
68  *   K'(j) = 0x6D703EF3    (32 <= j <= 47)      int(2**30 x cbrt(5))
69  *   K'(j) = 0x7A6D76E9    (48 <= j <= 63)      int(2**30 x cbrt(7))
70  *   K'(j) = 0x00000000    (64 <= j <= 79)
71  *
72  *
73  *   selection of message word
74  *
75  *   r(j)      = j                    (0 <= j <= 15)
76  *   r(16..31) = 7, 4, 13, 1, 10, 6, 15, 3, 12, 0, 9, 5, 2, 14, 11, 8
77  *   r(32..47) = 3, 10, 14, 4, 9, 15, 8, 1, 2, 7, 0, 6, 13, 11, 5, 12
78  *   r(48..63) = 1, 9, 11, 10, 0, 8, 12, 4, 13, 3, 7, 15, 14, 5, 6, 2
79  *   r(64..79) = 4, 0, 5, 9, 7, 12, 2, 10, 14, 1, 3, 8, 11, 6, 15, 13
80  *   r0(0..15) = 5, 14, 7, 0, 9, 2, 11, 4, 13, 6, 15, 8, 1, 10, 3, 12
81  *   r0(16..31)= 6, 11, 3, 7, 0, 13, 5, 10, 14, 15, 8, 12, 4, 9, 1, 2
82  *   r0(32..47)= 15, 5, 1, 3, 7, 14, 6, 9, 11, 8, 12, 2, 10, 0, 4, 13
83  *   r0(48..63)= 8, 6, 4, 1, 3, 11, 15, 0, 5, 12, 2, 13, 9, 7, 10, 14
84  *   r0(64..79)= 12, 15, 10, 4, 1, 5, 8, 7, 6, 2, 13, 14, 0, 3, 9, 11
85  *
86  *
87  *   amount for rotate left (rol)
88  *
89  *   s(0..15)  = 11, 14, 15, 12, 5, 8, 7, 9, 11, 13, 14, 15, 6, 7, 9, 8
90  *   s(16..31) = 7, 6, 8, 13, 11, 9, 7, 15, 7, 12, 15, 9, 11, 7, 13, 12
91  *   s(32..47) = 11, 13, 6, 7, 14, 9, 13, 15, 14, 8, 13, 6, 5, 12, 7, 5
92  *   s(48..63) = 11, 12, 14, 15, 14, 15, 9, 8, 9, 14, 5, 6, 8, 6, 5, 12
93  *   s(64..79) = 9, 15, 5, 11, 6, 8, 13, 12, 5, 12, 13, 14, 11, 8, 5, 6
94  *   s'(0..15) = 8, 9, 9, 11, 13, 15, 15, 5, 7, 7, 8, 11, 14, 14, 12, 6
95  *   s'(16..31)= 9, 13, 15, 7, 12, 8, 9, 11, 7, 7, 12, 7, 6, 15, 13, 11
96  *   s'(32..47)= 9, 7, 15, 11, 8, 6, 6, 14, 12, 13, 5, 14, 13, 13, 7, 5
97  *   s'(48..63)= 15, 5, 8, 11, 14, 14, 6, 14, 6, 9, 12, 9, 12, 5, 15, 8
98  *   s'(64..79)= 8, 5, 12, 9, 12, 5, 14, 6, 8, 13, 6, 5, 15, 13, 11, 11
99  *
100  *
101  *   initial value (hexadecimal)
102  *
103  *   h0 = 0x67452301; h1 = 0xEFCDAB89; h2 = 0x98BADCFE; h3 = 0x10325476;
104  *                                                      h4 = 0xC3D2E1F0;
105  *
106  *
107  * RIPEMD-160: pseudo-code
108  *
109  *   It is assumed that the message after padding consists of t 16-word blocks
110  *   that will be denoted with X[i][j], with 0 <= i <= t-1 and 0 <= j <= 15.
111  *   The symbol [+] denotes addition modulo 2**32 and rol_s denotes cyclic left
112  *   shift (rotate) over s positions.
113  *
114  *
115  *   for i := 0 to t-1 {
116  *       A := h0; B := h1; C := h2; D = h3; E = h4;
117  *       A' := h0; B' := h1; C' := h2; D' = h3; E' = h4;
118  *       for j := 0 to 79 {
119  *           T := rol_s(j)(A [+] f(j, B, C, D) [+] X[i][r(j)] [+] K(j)) [+] E;
120  *           A := E; E := D; D := rol_10(C); C := B; B := T;
121  *           T := rol_s'(j)(A' [+] f(79-j, B', C', D') [+] X[i][r'(j)]
122                                                        [+] K'(j)) [+] E';
123  *           A' := E'; E' := D'; D' := rol_10(C'); C' := B'; B' := T;
124  *       }
125  *       T := h1 [+] C [+] D'; h1 := h2 [+] D [+] E'; h2 := h3 [+] E [+] A';
126  *       h3 := h4 [+] A [+] B'; h4 := h0 [+] B [+] C'; h0 := T;
127  *   }
128  */
129
130 /* Some examples:
131  * ""                    9c1185a5c5e9fc54612808977ee8f548b2258d31
132  * "a"                   0bdc9d2d256b3ee9daae347be6f4dc835a467ffe
133  * "abc"                 8eb208f7e05d987a9b044a8e98c6b087f15a0bfc
134  * "message digest"      5d0689ef49d2fae572b881b123a85ffa21595f36
135  * "a...z"               f71c27109c692c1b56bbdceb5b9d2865b3708dbc
136  * "abcdbcde...nopq"     12a053384a9c0c88e405a06c27dcf49ada62eb2b
137  * "A...Za...z0...9"     b0e20b6e3116640286ed3a87a5713079b21f5189
138  * 8 times "1234567890"  9b752e45573d4b39f4dbd3323cab82bf63326bfb
139  * 1 million times "a"   52783243c1697bdbe16d37f97f68f08325dc1528
140  */
141
142
143 void
144 rmd160_init( RMD160_CONTEXT *hd )
145 {
146     hd->h0 = 0x67452301;
147     hd->h1 = 0xEFCDAB89;
148     hd->h2 = 0x98BADCFE;
149     hd->h3 = 0x10325476;
150     hd->h4 = 0xC3D2E1F0;
151     hd->nblocks = 0;
152     hd->count = 0;
153 }
154
155
156 #if defined(__GNUC__) && defined(__i386__)
157 static inline u32
158 rol(int n, u32 x)
159 {
160         __asm__("roll %%cl,%0"
161                 :"=r" (x)
162                 :"0" (x),"c" (n));
163         return x;
164 }
165 #else
166   #define rol(n,x) ( ((x) << (n)) | ((x) >> (32-(n))) )
167 #endif
168
169
170 /****************
171  * Transform the message X which consists of 16 32-bit-words
172  */
173 static void
174 transform( RMD160_CONTEXT *hd, byte *data )
175 {
176     u32 a,b,c,d,e,aa,bb,cc,dd,ee,t;
177   #ifdef BIG_ENDIAN_HOST
178     u32 x[16];
179     { int i;
180       byte *p2, *p1;
181       for(i=0, p1=data, p2=(byte*)x; i < 16; i++, p2 += 4 ) {
182         p2[3] = *p1++;
183         p2[2] = *p1++;
184         p2[1] = *p1++;
185         p2[0] = *p1++;
186       }
187     }
188   #else
189    #if 0
190     u32 *x =(u32*)data;
191    #else
192     /* this version is better because it is always aligned;
193      * The performance penalty on a 586-100 is about 6% which
194      * is acceptable - because the data is more local it might
195      * also be possible that this is faster on some machines.
196      * This function (when compiled with -02 on gcc 2.7.2)
197      * executes on a 586-100 (39.73 bogomips) at about 1900kb/sec;
198      * [measured with a 4MB data and "gpgm --print-md rmd160"] */
199     u32 x[16];
200     memcpy( x, data, 64 );
201    #endif
202   #endif
203
204
205 #define K0  0x00000000
206 #define K1  0x5A827999
207 #define K2  0x6ED9EBA1
208 #define K3  0x8F1BBCDC
209 #define K4  0xA953FD4E
210 #define KK0 0x50A28BE6
211 #define KK1 0x5C4DD124
212 #define KK2 0x6D703EF3
213 #define KK3 0x7A6D76E9
214 #define KK4 0x00000000
215 #define F0(x,y,z)   ( (x) ^ (y) ^ (z) )
216 #define F1(x,y,z)   ( ((x) & (y)) | (~(x) & (z)) )
217 #define F2(x,y,z)   ( ((x) | ~(y)) ^ (z) )
218 #define F3(x,y,z)   ( ((x) & (z)) | ((y) & ~(z)) )
219 #define F4(x,y,z)   ( (x) ^ ((y) | ~(z)) )
220 #define R(a,b,c,d,e,f,k,r,s) do { t = a + f(b,c,d) + k + x[r]; \
221                                   a = rol(s,t) + e;            \
222                                   c = rol(10,c);               \
223                                 } while(0)
224
225     /* left lane */
226     a = hd->h0;
227     b = hd->h1;
228     c = hd->h2;
229     d = hd->h3;
230     e = hd->h4;
231     R( a, b, c, d, e, F0, K0,  0, 11 );
232     R( e, a, b, c, d, F0, K0,  1, 14 );
233     R( d, e, a, b, c, F0, K0,  2, 15 );
234     R( c, d, e, a, b, F0, K0,  3, 12 );
235     R( b, c, d, e, a, F0, K0,  4,  5 );
236     R( a, b, c, d, e, F0, K0,  5,  8 );
237     R( e, a, b, c, d, F0, K0,  6,  7 );
238     R( d, e, a, b, c, F0, K0,  7,  9 );
239     R( c, d, e, a, b, F0, K0,  8, 11 );
240     R( b, c, d, e, a, F0, K0,  9, 13 );
241     R( a, b, c, d, e, F0, K0, 10, 14 );
242     R( e, a, b, c, d, F0, K0, 11, 15 );
243     R( d, e, a, b, c, F0, K0, 12,  6 );
244     R( c, d, e, a, b, F0, K0, 13,  7 );
245     R( b, c, d, e, a, F0, K0, 14,  9 );
246     R( a, b, c, d, e, F0, K0, 15,  8 );
247     R( e, a, b, c, d, F1, K1,  7,  7 );
248     R( d, e, a, b, c, F1, K1,  4,  6 );
249     R( c, d, e, a, b, F1, K1, 13,  8 );
250     R( b, c, d, e, a, F1, K1,  1, 13 );
251     R( a, b, c, d, e, F1, K1, 10, 11 );
252     R( e, a, b, c, d, F1, K1,  6,  9 );
253     R( d, e, a, b, c, F1, K1, 15,  7 );
254     R( c, d, e, a, b, F1, K1,  3, 15 );
255     R( b, c, d, e, a, F1, K1, 12,  7 );
256     R( a, b, c, d, e, F1, K1,  0, 12 );
257     R( e, a, b, c, d, F1, K1,  9, 15 );
258     R( d, e, a, b, c, F1, K1,  5,  9 );
259     R( c, d, e, a, b, F1, K1,  2, 11 );
260     R( b, c, d, e, a, F1, K1, 14,  7 );
261     R( a, b, c, d, e, F1, K1, 11, 13 );
262     R( e, a, b, c, d, F1, K1,  8, 12 );
263     R( d, e, a, b, c, F2, K2,  3, 11 );
264     R( c, d, e, a, b, F2, K2, 10, 13 );
265     R( b, c, d, e, a, F2, K2, 14,  6 );
266     R( a, b, c, d, e, F2, K2,  4,  7 );
267     R( e, a, b, c, d, F2, K2,  9, 14 );
268     R( d, e, a, b, c, F2, K2, 15,  9 );
269     R( c, d, e, a, b, F2, K2,  8, 13 );
270     R( b, c, d, e, a, F2, K2,  1, 15 );
271     R( a, b, c, d, e, F2, K2,  2, 14 );
272     R( e, a, b, c, d, F2, K2,  7,  8 );
273     R( d, e, a, b, c, F2, K2,  0, 13 );
274     R( c, d, e, a, b, F2, K2,  6,  6 );
275     R( b, c, d, e, a, F2, K2, 13,  5 );
276     R( a, b, c, d, e, F2, K2, 11, 12 );
277     R( e, a, b, c, d, F2, K2,  5,  7 );
278     R( d, e, a, b, c, F2, K2, 12,  5 );
279     R( c, d, e, a, b, F3, K3,  1, 11 );
280     R( b, c, d, e, a, F3, K3,  9, 12 );
281     R( a, b, c, d, e, F3, K3, 11, 14 );
282     R( e, a, b, c, d, F3, K3, 10, 15 );
283     R( d, e, a, b, c, F3, K3,  0, 14 );
284     R( c, d, e, a, b, F3, K3,  8, 15 );
285     R( b, c, d, e, a, F3, K3, 12,  9 );
286     R( a, b, c, d, e, F3, K3,  4,  8 );
287     R( e, a, b, c, d, F3, K3, 13,  9 );
288     R( d, e, a, b, c, F3, K3,  3, 14 );
289     R( c, d, e, a, b, F3, K3,  7,  5 );
290     R( b, c, d, e, a, F3, K3, 15,  6 );
291     R( a, b, c, d, e, F3, K3, 14,  8 );
292     R( e, a, b, c, d, F3, K3,  5,  6 );
293     R( d, e, a, b, c, F3, K3,  6,  5 );
294     R( c, d, e, a, b, F3, K3,  2, 12 );
295     R( b, c, d, e, a, F4, K4,  4,  9 );
296     R( a, b, c, d, e, F4, K4,  0, 15 );
297     R( e, a, b, c, d, F4, K4,  5,  5 );
298     R( d, e, a, b, c, F4, K4,  9, 11 );
299     R( c, d, e, a, b, F4, K4,  7,  6 );
300     R( b, c, d, e, a, F4, K4, 12,  8 );
301     R( a, b, c, d, e, F4, K4,  2, 13 );
302     R( e, a, b, c, d, F4, K4, 10, 12 );
303     R( d, e, a, b, c, F4, K4, 14,  5 );
304     R( c, d, e, a, b, F4, K4,  1, 12 );
305     R( b, c, d, e, a, F4, K4,  3, 13 );
306     R( a, b, c, d, e, F4, K4,  8, 14 );
307     R( e, a, b, c, d, F4, K4, 11, 11 );
308     R( d, e, a, b, c, F4, K4,  6,  8 );
309     R( c, d, e, a, b, F4, K4, 15,  5 );
310     R( b, c, d, e, a, F4, K4, 13,  6 );
311
312     aa = a; bb = b; cc = c; dd = d; ee = e;
313
314     /* right lane */
315     a = hd->h0;
316     b = hd->h1;
317     c = hd->h2;
318     d = hd->h3;
319     e = hd->h4;
320     R( a, b, c, d, e, F4, KK0,  5,  8);
321     R( e, a, b, c, d, F4, KK0, 14,  9);
322     R( d, e, a, b, c, F4, KK0,  7,  9);
323     R( c, d, e, a, b, F4, KK0,  0, 11);
324     R( b, c, d, e, a, F4, KK0,  9, 13);
325     R( a, b, c, d, e, F4, KK0,  2, 15);
326     R( e, a, b, c, d, F4, KK0, 11, 15);
327     R( d, e, a, b, c, F4, KK0,  4,  5);
328     R( c, d, e, a, b, F4, KK0, 13,  7);
329     R( b, c, d, e, a, F4, KK0,  6,  7);
330     R( a, b, c, d, e, F4, KK0, 15,  8);
331     R( e, a, b, c, d, F4, KK0,  8, 11);
332     R( d, e, a, b, c, F4, KK0,  1, 14);
333     R( c, d, e, a, b, F4, KK0, 10, 14);
334     R( b, c, d, e, a, F4, KK0,  3, 12);
335     R( a, b, c, d, e, F4, KK0, 12,  6);
336     R( e, a, b, c, d, F3, KK1,  6,  9);
337     R( d, e, a, b, c, F3, KK1, 11, 13);
338     R( c, d, e, a, b, F3, KK1,  3, 15);
339     R( b, c, d, e, a, F3, KK1,  7,  7);
340     R( a, b, c, d, e, F3, KK1,  0, 12);
341     R( e, a, b, c, d, F3, KK1, 13,  8);
342     R( d, e, a, b, c, F3, KK1,  5,  9);
343     R( c, d, e, a, b, F3, KK1, 10, 11);
344     R( b, c, d, e, a, F3, KK1, 14,  7);
345     R( a, b, c, d, e, F3, KK1, 15,  7);
346     R( e, a, b, c, d, F3, KK1,  8, 12);
347     R( d, e, a, b, c, F3, KK1, 12,  7);
348     R( c, d, e, a, b, F3, KK1,  4,  6);
349     R( b, c, d, e, a, F3, KK1,  9, 15);
350     R( a, b, c, d, e, F3, KK1,  1, 13);
351     R( e, a, b, c, d, F3, KK1,  2, 11);
352     R( d, e, a, b, c, F2, KK2, 15,  9);
353     R( c, d, e, a, b, F2, KK2,  5,  7);
354     R( b, c, d, e, a, F2, KK2,  1, 15);
355     R( a, b, c, d, e, F2, KK2,  3, 11);
356     R( e, a, b, c, d, F2, KK2,  7,  8);
357     R( d, e, a, b, c, F2, KK2, 14,  6);
358     R( c, d, e, a, b, F2, KK2,  6,  6);
359     R( b, c, d, e, a, F2, KK2,  9, 14);
360     R( a, b, c, d, e, F2, KK2, 11, 12);
361     R( e, a, b, c, d, F2, KK2,  8, 13);
362     R( d, e, a, b, c, F2, KK2, 12,  5);
363     R( c, d, e, a, b, F2, KK2,  2, 14);
364     R( b, c, d, e, a, F2, KK2, 10, 13);
365     R( a, b, c, d, e, F2, KK2,  0, 13);
366     R( e, a, b, c, d, F2, KK2,  4,  7);
367     R( d, e, a, b, c, F2, KK2, 13,  5);
368     R( c, d, e, a, b, F1, KK3,  8, 15);
369     R( b, c, d, e, a, F1, KK3,  6,  5);
370     R( a, b, c, d, e, F1, KK3,  4,  8);
371     R( e, a, b, c, d, F1, KK3,  1, 11);
372     R( d, e, a, b, c, F1, KK3,  3, 14);
373     R( c, d, e, a, b, F1, KK3, 11, 14);
374     R( b, c, d, e, a, F1, KK3, 15,  6);
375     R( a, b, c, d, e, F1, KK3,  0, 14);
376     R( e, a, b, c, d, F1, KK3,  5,  6);
377     R( d, e, a, b, c, F1, KK3, 12,  9);
378     R( c, d, e, a, b, F1, KK3,  2, 12);
379     R( b, c, d, e, a, F1, KK3, 13,  9);
380     R( a, b, c, d, e, F1, KK3,  9, 12);
381     R( e, a, b, c, d, F1, KK3,  7,  5);
382     R( d, e, a, b, c, F1, KK3, 10, 15);
383     R( c, d, e, a, b, F1, KK3, 14,  8);
384     R( b, c, d, e, a, F0, KK4, 12,  8);
385     R( a, b, c, d, e, F0, KK4, 15,  5);
386     R( e, a, b, c, d, F0, KK4, 10, 12);
387     R( d, e, a, b, c, F0, KK4,  4,  9);
388     R( c, d, e, a, b, F0, KK4,  1, 12);
389     R( b, c, d, e, a, F0, KK4,  5,  5);
390     R( a, b, c, d, e, F0, KK4,  8, 14);
391     R( e, a, b, c, d, F0, KK4,  7,  6);
392     R( d, e, a, b, c, F0, KK4,  6,  8);
393     R( c, d, e, a, b, F0, KK4,  2, 13);
394     R( b, c, d, e, a, F0, KK4, 13,  6);
395     R( a, b, c, d, e, F0, KK4, 14,  5);
396     R( e, a, b, c, d, F0, KK4,  0, 15);
397     R( d, e, a, b, c, F0, KK4,  3, 13);
398     R( c, d, e, a, b, F0, KK4,  9, 11);
399     R( b, c, d, e, a, F0, KK4, 11, 11);
400
401
402     t      = hd->h1 + d + cc;
403     hd->h1 = hd->h2 + e + dd;
404     hd->h2 = hd->h3 + a + ee;
405     hd->h3 = hd->h4 + b + aa;
406     hd->h4 = hd->h0 + c + bb;
407     hd->h0 = t;
408 }
409
410
411 /* Update the message digest with the contents
412  * of INBUF with length INLEN.
413  */
414 static void
415 rmd160_write( RMD160_CONTEXT *hd, byte *inbuf, size_t inlen)
416 {
417     if( hd->count == 64 ) { /* flush the buffer */
418         transform( hd, hd->buf );
419         hd->count = 0;
420         hd->nblocks++;
421     }
422     if( !inbuf )
423         return;
424     if( hd->count ) {
425         for( ; inlen && hd->count < 64; inlen-- )
426             hd->buf[hd->count++] = *inbuf++;
427         rmd160_write( hd, NULL, 0 );
428         if( !inlen )
429             return;
430     }
431
432     while( inlen >= 64 ) {
433         transform( hd, inbuf );
434         hd->count = 0;
435         hd->nblocks++;
436         inlen -= 64;
437         inbuf += 64;
438     }
439     for( ; inlen && hd->count < 64; inlen-- )
440         hd->buf[hd->count++] = *inbuf++;
441 }
442
443 /****************
444  * Apply the rmd160 transform function on the buffer which must have
445  * a length 64 bytes. Do not use this function together with the
446  * other functions, use rmd160_init to initialize internal variables.
447  * Returns: 16 bytes in buffer with the mixed contentes of buffer.
448  */
449 void
450 rmd160_mixblock( RMD160_CONTEXT *hd, char *buffer )
451 {
452     char *p = buffer;
453     transform( hd, buffer );
454   #define X(a) do { *(u32*)p = hd->h##a ; p += 4; } while(0)
455     X(0);
456     X(1);
457     X(2);
458     X(3);
459     X(4);
460   #undef X
461 }
462
463
464 /* The routine terminates the computation
465  */
466
467 static void
468 rmd160_final( RMD160_CONTEXT *hd )
469 {
470     u32 t, msb, lsb;
471     byte *p;
472
473     rmd160_write(hd, NULL, 0); /* flush */;
474
475     msb = 0;
476     t = hd->nblocks;
477     if( (lsb = t << 6) < t ) /* multiply by 64 to make a byte count */
478         msb++;
479     msb += t >> 26;
480     t = lsb;
481     if( (lsb = t + hd->count) < t ) /* add the count */
482         msb++;
483     t = lsb;
484     if( (lsb = t << 3) < t ) /* multiply by 8 to make a bit count */
485         msb++;
486     msb += t >> 29;
487
488     if( hd->count < 56 ) { /* enough room */
489         hd->buf[hd->count++] = 0x80; /* pad */
490         while( hd->count < 56 )
491             hd->buf[hd->count++] = 0;  /* pad */
492     }
493     else { /* need one extra block */
494         hd->buf[hd->count++] = 0x80; /* pad character */
495         while( hd->count < 64 )
496             hd->buf[hd->count++] = 0;
497         rmd160_write(hd, NULL, 0);  /* flush */;
498         memset(hd->buf, 0, 56 ); /* fill next block with zeroes */
499     }
500     /* append the 64 bit count */
501     hd->buf[56] = lsb      ;
502     hd->buf[57] = lsb >>  8;
503     hd->buf[58] = lsb >> 16;
504     hd->buf[59] = lsb >> 24;
505     hd->buf[60] = msb      ;
506     hd->buf[61] = msb >>  8;
507     hd->buf[62] = msb >> 16;
508     hd->buf[63] = msb >> 24;
509     transform( hd, hd->buf );
510
511     p = hd->buf;
512   #ifdef BIG_ENDIAN_HOST
513     #define X(a) do { *p++ = hd->h##a      ; *p++ = hd->h##a >> 8;      \
514                       *p++ = hd->h##a >> 16; *p++ = hd->h##a >> 24; } while(0)
515   #else /* little endian */
516     #define X(a) do { *(u32*)p = hd->h##a ; p += 4; } while(0)
517   #endif
518     X(0);
519     X(1);
520     X(2);
521     X(3);
522     X(4);
523   #undef X
524 }
525
526 static byte *
527 rmd160_read( RMD160_CONTEXT *hd )
528 {
529     return hd->buf;
530 }
531
532
533
534 /****************
535  * Shortcut functions which puts the hash value of the supplied buffer
536  * into outbuf which must have a size of 20 bytes.
537  */
538 void
539 rmd160_hash_buffer( char *outbuf, const char *buffer, size_t length )
540 {
541     RMD160_CONTEXT hd;
542
543     rmd160_init( &hd );
544     rmd160_write( &hd, (byte*)buffer, length );
545     rmd160_final( &hd );
546     memcpy( outbuf, hd.buf, 20 );
547 }
548
549
550 /****************
551  * Return some information about the algorithm.  We need algo here to
552  * distinguish different flavors of the algorithm.
553  * Returns: A pointer to string describing the algorithm or NULL if
554  *          the ALGO is invalid.
555  */
556 static const char *
557 rmd160_get_info( int algo, size_t *contextsize,
558                byte **r_asnoid, int *r_asnlen, int *r_mdlen,
559                void (**r_init)( void *c ),
560                void (**r_write)( void *c, byte *buf, size_t nbytes ),
561                void (**r_final)( void *c ),
562                byte *(**r_read)( void *c )
563              )
564 {
565     static byte asn[15] = /* Object ID is 1.3.36.3.2.1 */
566           { 0x30, 0x21, 0x30, 0x09, 0x06, 0x05, 0x2b, 0x24, 0x03,
567             0x02, 0x01, 0x05, 0x00, 0x04, 0x14 };
568
569     if( algo != 3 )
570         return NULL;
571
572     *contextsize = sizeof(RMD160_CONTEXT);
573     *r_asnoid = asn;
574     *r_asnlen = DIM(asn);
575     *r_mdlen = 20;
576     *r_init  = (void (*)(void *))rmd160_init;
577     *r_write = (void (*)(void *, byte*, size_t))rmd160_write;
578     *r_final = (void (*)(void *))rmd160_final;
579     *r_read  = (byte *(*)(void *))rmd160_read;
580
581     return "RIPEMD160";
582 }
583
584
585 #ifndef IS_MODULE
586 static
587 #endif
588 const char * const gnupgext_version = "RMD160 ($Revision$)";
589
590 static struct {
591     int class;
592     int version;
593     int  value;
594     void (*func)(void);
595 } func_table[] = {
596     { 10, 1, 0, (void(*)(void))rmd160_get_info },
597     { 11, 1, 3 },
598 };
599
600
601 #ifndef IS_MODULE
602 static
603 #endif
604 void *
605 gnupgext_enum_func( int what, int *sequence, int *class, int *vers )
606 {
607     void *ret;
608     int i = *sequence;
609
610     do {
611         if( i >= DIM(func_table) || i < 0 ) {
612             return NULL;
613         }
614         *class = func_table[i].class;
615         *vers  = func_table[i].version;
616         switch( *class ) {
617           case 11:
618           case 21:
619           case 31:
620             ret = &func_table[i].value;
621             break;
622           default:
623             ret = func_table[i].func;
624             break;
625         }
626         i++;
627     } while( what && what != *class );
628
629     *sequence = i;
630     return ret;
631 }
632
633
634
635
636 #ifndef IS_MODULE
637 void
638 rmd160_constructor(void)
639 {
640     register_internal_cipher_extension( gnupgext_version, gnupgext_enum_func );
641 }
642 #endif
643