src/md5.c

   1 #include "first.h"
   2
   3 /* MD5C.C - RSA Data Security, Inc., MD5 message-digest algorithm
   4  */
   5
   6 /* Copyright (C) 1991-2, RSA Data Security, Inc. Created 1991. All
   7 rights reserved.
   8
   9 License to copy and use this software is granted provided that it
  10 is identified as the "RSA Data Security, Inc. MD5 Message-Digest
  11 Algorithm" in all material mentioning or referencing this software
  12 or this function.
  13
  14 License is also granted to make and use derivative works provided
  15 that such works are identified as "derived from the RSA Data
  16 Security, Inc. MD5 Message-Digest Algorithm" in all material
  17 mentioning or referencing the derived work.
  18
  19 RSA Data Security, Inc. makes no representations concerning either
  20 the merchantability of this software or the suitability of this
  21 software for any particular purpose. It is provided "as is"
  22 without express or implied warranty of any kind.
  23
  24 These notices must be retained in any copies of any part of this
  25 documentation and/or software.
  26  */
  27
  28
  29 #include "md5.h"
  30
  31 #ifndef USE_OPENSSL
  32 #include <string.h>
  33
  34 /* Constants for MD5Transform routine.
  35  */
  36
  37 #define S11 7
  38 #define S12 12
  39 #define S13 17
  40 #define S14 22
  41 #define S21 5
  42 #define S22 9
  43 #define S23 14
  44 #define S24 20
  45 #define S31 4
  46 #define S32 11
  47 #define S33 16
  48 #define S34 23
  49 #define S41 6
  50 #define S42 10
  51 #define S43 15
  52 #define S44 21
  53
  54 static void li_MD5Transform (UINT4 [4], const unsigned char [64]);
  55 static void Encode (unsigned char *, UINT4 *, unsigned int);
  56 static void Decode (UINT4 *, const unsigned char *, unsigned int);
  57
  58 #ifdef HAVE_MEMCPY
  59 #define MD5_memcpy(output, input, len) memcpy((output), (input), (len))
  60 #else
  61 static void MD5_memcpy (POINTER, POINTER, unsigned int);
  62 #endif
  63 #ifdef HAVE_MEMSET
  64 #define MD5_memset(output, value, len) memset((output), (value), (len))
  65 #else
  66 static void MD5_memset (POINTER, int, unsigned int);
  67 #endif
  68
  69 static unsigned char PADDING[64] = {
  70   0x80, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
  71   0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
  72   0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0
  73 };
  74
  75 /* F, G, H and I are basic MD5 functions.
  76  */
  77 #define F(x, y, z) (((x) & (y)) | ((~x) & (z)))
  78 #define G(x, y, z) (((x) & (z)) | ((y) & (~z)))
  79 #define H(x, y, z) ((x) ^ (y) ^ (z))
  80 #define I(x, y, z) ((y) ^ ((x) | (~z)))
  81
  82 /* ROTATE_LEFT rotates x left n bits.
  83  */
  84 #define ROTATE_LEFT(x, n) (((x) << (n)) | ((x) >> (32-(n))))
  85
  86 /* FF, GG, HH, and II transformations for rounds 1, 2, 3, and 4.
  87 Rotation is separate from addition to prevent recomputation.
  88  */
  89 #define FF(a, b, c, d, x, s, ac) { \
  90  (a) += F ((b), (c), (d)) + (x) + (UINT4)(ac); \
  91  (a) = ROTATE_LEFT ((a), (s)); \
  92  (a) += (b); \
  93   }
  94 #define GG(a, b, c, d, x, s, ac) { \
  95  (a) += G ((b), (c), (d)) + (x) + (UINT4)(ac); \
  96  (a) = ROTATE_LEFT ((a), (s)); \
  97  (a) += (b); \
  98   }
  99 #define HH(a, b, c, d, x, s, ac) { \
 100  (a) += H ((b), (c), (d)) + (x) + (UINT4)(ac); \
 101  (a) = ROTATE_LEFT ((a), (s)); \
 102  (a) += (b); \
 103   }
 104 #define II(a, b, c, d, x, s, ac) { \
 105  (a) += I ((b), (c), (d)) + (x) + (UINT4)(ac); \
 106  (a) = ROTATE_LEFT ((a), (s)); \
 107  (a) += (b); \
 108   }
 109
 110 /* MD5 initialization. Begins an MD5 operation, writing a new context.
 111  */
 112 void li_MD5_Init (li_MD5_CTX *context)
 113 {
 114   context->count[0] = context->count[1] = 0;
 115   /* Load magic initialization constants.
 116 */
 117   context->state[0] = 0x67452301;
 118   context->state[1] = 0xefcdab89;
 119   context->state[2] = 0x98badcfe;
 120   context->state[3] = 0x10325476;
 121 }
 122
 123 /* MD5 block update operation. Continues an MD5 message-digest
 124   operation, processing another message block, and updating the
 125   context.
 126  */
 127 void li_MD5_Update (li_MD5_CTX *context, const void *_input, unsigned int inputLen)
 128 {
 129   unsigned int i, ndx, partLen;
 130   const unsigned char *input = (const unsigned char*) _input;
 131
 132   /* Compute number of bytes mod 64 */
 133   ndx = (unsigned int)((context->count[0] >> 3) & 0x3F);
 134
 135   /* Update number of bits */
 136   if ((context->count[0] += ((UINT4)inputLen << 3))
 137
 138    < ((UINT4)inputLen << 3))
 139  context->count[1]++;
 140   context->count[1] += ((UINT4)inputLen >> 29);
 141
 142   partLen = 64 - ndx;
 143
 144   /* Transform as many times as possible.
 145 */
 146   if (inputLen >= partLen) {
 147  MD5_memcpy
 148    ((POINTER)&context->buffer[ndx], (POINTER)input, partLen);
 149  li_MD5Transform (context->state, context->buffer);
 150
 151  for (i = partLen; i + 63 < inputLen; i += 64)
 152    li_MD5Transform (context->state, &input[i]);
 153
 154  ndx = 0;
 155   }
 156   else
 157  i = 0;
 158
 159   /* Buffer remaining input */
 160   MD5_memcpy
 161  ((POINTER)&context->buffer[ndx], (POINTER)&input[i],
 162   inputLen-i);
 163 }
 164
 165 /* MD5 finalization. Ends an MD5 message-digest operation, writing the
 166   the message digest and zeroizing the context.
 167  */
 168 void li_MD5_Final (unsigned char digest[16], li_MD5_CTX *context)
 169 {
 170   unsigned char bits[8];
 171   unsigned int ndx, padLen;
 172
 173   /* Save number of bits */
 174   Encode (bits, context->count, 8);
 175
 176   /* Pad out to 56 mod 64.
 177 */
 178   ndx = (unsigned int)((context->count[0] >> 3) & 0x3f);
 179   padLen = (ndx < 56) ? (56 - ndx) : (120 - ndx);
 180   li_MD5_Update (context, PADDING, padLen);
 181
 182   /* Append length (before padding) */
 183   li_MD5_Update (context, bits, 8);
 184
 185   /* Store state in digest */
 186   Encode (digest, context->state, 16);
 187
 188   /* Zeroize sensitive information.
 189 */
 190   MD5_memset ((POINTER)context, 0, sizeof (*context));
 191 }
 192
 193 /* MD5 basic transformation. Transforms state based on block.
 194  */
 195 static void li_MD5Transform (UINT4 state[4], const unsigned char block[64])
 196 {
 197   UINT4 a = state[0], b = state[1], c = state[2], d = state[3], x[16];
 198
 199   Decode (x, block, 64);
 200
 201   /* Round 1 */
 202   FF (a, b, c, d, x[ 0], S11, 0xd76aa478); /* 1 */
 203   FF (d, a, b, c, x[ 1], S12, 0xe8c7b756); /* 2 */
 204   FF (c, d, a, b, x[ 2], S13, 0x242070db); /* 3 */
 205   FF (b, c, d, a, x[ 3], S14, 0xc1bdceee); /* 4 */
 206   FF (a, b, c, d, x[ 4], S11, 0xf57c0faf); /* 5 */
 207   FF (d, a, b, c, x[ 5], S12, 0x4787c62a); /* 6 */
 208   FF (c, d, a, b, x[ 6], S13, 0xa8304613); /* 7 */
 209   FF (b, c, d, a, x[ 7], S14, 0xfd469501); /* 8 */
 210   FF (a, b, c, d, x[ 8], S11, 0x698098d8); /* 9 */
 211   FF (d, a, b, c, x[ 9], S12, 0x8b44f7af); /* 10 */
 212   FF (c, d, a, b, x[10], S13, 0xffff5bb1); /* 11 */
 213   FF (b, c, d, a, x[11], S14, 0x895cd7be); /* 12 */
 214   FF (a, b, c, d, x[12], S11, 0x6b901122); /* 13 */
 215   FF (d, a, b, c, x[13], S12, 0xfd987193); /* 14 */
 216   FF (c, d, a, b, x[14], S13, 0xa679438e); /* 15 */
 217   FF (b, c, d, a, x[15], S14, 0x49b40821); /* 16 */
 218
 219  /* Round 2 */
 220   GG (a, b, c, d, x[ 1], S21, 0xf61e2562); /* 17 */
 221   GG (d, a, b, c, x[ 6], S22, 0xc040b340); /* 18 */
 222   GG (c, d, a, b, x[11], S23, 0x265e5a51); /* 19 */
 223   GG (b, c, d, a, x[ 0], S24, 0xe9b6c7aa); /* 20 */
 224   GG (a, b, c, d, x[ 5], S21, 0xd62f105d); /* 21 */
 225   GG (d, a, b, c, x[10], S22,  0x2441453); /* 22 */
 226   GG (c, d, a, b, x[15], S23, 0xd8a1e681); /* 23 */
 227   GG (b, c, d, a, x[ 4], S24, 0xe7d3fbc8); /* 24 */
 228   GG (a, b, c, d, x[ 9], S21, 0x21e1cde6); /* 25 */
 229   GG (d, a, b, c, x[14], S22, 0xc33707d6); /* 26 */
 230   GG (c, d, a, b, x[ 3], S23, 0xf4d50d87); /* 27 */
 231
 232   GG (b, c, d, a, x[ 8], S24, 0x455a14ed); /* 28 */
 233   GG (a, b, c, d, x[13], S21, 0xa9e3e905); /* 29 */
 234   GG (d, a, b, c, x[ 2], S22, 0xfcefa3f8); /* 30 */
 235   GG (c, d, a, b, x[ 7], S23, 0x676f02d9); /* 31 */
 236   GG (b, c, d, a, x[12], S24, 0x8d2a4c8a); /* 32 */
 237
 238   /* Round 3 */
 239   HH (a, b, c, d, x[ 5], S31, 0xfffa3942); /* 33 */
 240   HH (d, a, b, c, x[ 8], S32, 0x8771f681); /* 34 */
 241   HH (c, d, a, b, x[11], S33, 0x6d9d6122); /* 35 */
 242   HH (b, c, d, a, x[14], S34, 0xfde5380c); /* 36 */
 243   HH (a, b, c, d, x[ 1], S31, 0xa4beea44); /* 37 */
 244   HH (d, a, b, c, x[ 4], S32, 0x4bdecfa9); /* 38 */
 245   HH (c, d, a, b, x[ 7], S33, 0xf6bb4b60); /* 39 */
 246   HH (b, c, d, a, x[10], S34, 0xbebfbc70); /* 40 */
 247   HH (a, b, c, d, x[13], S31, 0x289b7ec6); /* 41 */
 248   HH (d, a, b, c, x[ 0], S32, 0xeaa127fa); /* 42 */
 249   HH (c, d, a, b, x[ 3], S33, 0xd4ef3085); /* 43 */
 250   HH (b, c, d, a, x[ 6], S34,  0x4881d05); /* 44 */
 251   HH (a, b, c, d, x[ 9], S31, 0xd9d4d039); /* 45 */
 252   HH (d, a, b, c, x[12], S32, 0xe6db99e5); /* 46 */
 253   HH (c, d, a, b, x[15], S33, 0x1fa27cf8); /* 47 */
 254   HH (b, c, d, a, x[ 2], S34, 0xc4ac5665); /* 48 */
 255
 256   /* Round 4 */
 257   II (a, b, c, d, x[ 0], S41, 0xf4292244); /* 49 */
 258   II (d, a, b, c, x[ 7], S42, 0x432aff97); /* 50 */
 259   II (c, d, a, b, x[14], S43, 0xab9423a7); /* 51 */
 260   II (b, c, d, a, x[ 5], S44, 0xfc93a039); /* 52 */
 261   II (a, b, c, d, x[12], S41, 0x655b59c3); /* 53 */
 262   II (d, a, b, c, x[ 3], S42, 0x8f0ccc92); /* 54 */
 263   II (c, d, a, b, x[10], S43, 0xffeff47d); /* 55 */
 264   II (b, c, d, a, x[ 1], S44, 0x85845dd1); /* 56 */
 265   II (a, b, c, d, x[ 8], S41, 0x6fa87e4f); /* 57 */
 266   II (d, a, b, c, x[15], S42, 0xfe2ce6e0); /* 58 */
 267   II (c, d, a, b, x[ 6], S43, 0xa3014314); /* 59 */
 268   II (b, c, d, a, x[13], S44, 0x4e0811a1); /* 60 */
 269   II (a, b, c, d, x[ 4], S41, 0xf7537e82); /* 61 */
 270   II (d, a, b, c, x[11], S42, 0xbd3af235); /* 62 */
 271   II (c, d, a, b, x[ 2], S43, 0x2ad7d2bb); /* 63 */
 272   II (b, c, d, a, x[ 9], S44, 0xeb86d391); /* 64 */
 273
 274   state[0] += a;
 275   state[1] += b;
 276   state[2] += c;
 277   state[3] += d;
 278
 279   /* Zeroize sensitive information.
 280
 281 */
 282   MD5_memset ((POINTER)x, 0, sizeof (x));
 283 }
 284
 285 /* Encodes input (UINT4) into output (unsigned char). Assumes len is
 286   a multiple of 4.
 287  */
 288 static void Encode (unsigned char *output, UINT4 *input, unsigned int len)
 289 {
 290   unsigned int i, j;
 291
 292   for (i = 0, j = 0; j < len; i++, j += 4) {
 293  output[j] = (unsigned char)(input[i] & 0xff);
 294  output[j+1] = (unsigned char)((input[i] >> 8) & 0xff);
 295  output[j+2] = (unsigned char)((input[i] >> 16) & 0xff);
 296  output[j+3] = (unsigned char)((input[i] >> 24) & 0xff);
 297   }
 298 }
 299
 300 /* Decodes input (unsigned char) into output (UINT4). Assumes len is
 301   a multiple of 4.
 302  */
 303 static void Decode (UINT4 *output, const unsigned char *input, unsigned int len)
 304 {
 305   unsigned int i, j;
 306
 307   for (i = 0, j = 0; j < len; i++, j += 4)
 308  output[i] = ((UINT4)input[j]) | (((UINT4)input[j+1]) << 8) |
 309    (((UINT4)input[j+2]) << 16) | (((UINT4)input[j+3]) << 24);
 310 }
 311
 312 /* Note: Replace "for loop" with standard memcpy if possible.
 313  */
 314 #ifndef HAVE_MEMCPY
 315 static void MD5_memcpy (POINTER output, POINTER input, unsigned int len)
 316 {
 317   unsigned int i;
 318
 319   for (i = 0; i < len; i++)
 320
 321  output[i] = input[i];
 322 }
 323 #endif
 324 /* Note: Replace "for loop" with standard memset if possible.
 325  */
 326 #ifndef HAVE_MEMSET
 327 static void MD5_memset (POINTER output, int value, unsigned int len)
 328 {
 329   unsigned int i;
 330
 331   for (i = 0; i < len; i++)
 332  ((char *)output)[i] = (char)value;
 333 }
 334 #endif
 335 #endif