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