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