1 db6608bd 2003-11-23 devnull #include "os.h"
2 db6608bd 2003-11-23 devnull #include <libsec.h>
5 db6608bd 2003-11-23 devnull * rfc1321 requires that I include this. The code is new. The constants
6 db6608bd 2003-11-23 devnull * all come from the rfc (hence the copyright). We trade a table for the
7 db6608bd 2003-11-23 devnull * macros in rfc. The total size is a lot less. -- presotto
9 db6608bd 2003-11-23 devnull * Copyright (C) 1991-2, RSA Data Security, Inc. Created 1991. All
10 db6608bd 2003-11-23 devnull * rights reserved.
12 db6608bd 2003-11-23 devnull * License to copy and use this software is granted provided that it
13 db6608bd 2003-11-23 devnull * is identified as the "RSA Data Security, Inc. MD5 Message-Digest
14 db6608bd 2003-11-23 devnull * Algorithm" in all material mentioning or referencing this software
15 db6608bd 2003-11-23 devnull * or this function.
17 db6608bd 2003-11-23 devnull * License is also granted to make and use derivative works provided
18 db6608bd 2003-11-23 devnull * that such works are identified as "derived from the RSA Data
19 db6608bd 2003-11-23 devnull * Security, Inc. MD5 Message-Digest Algorithm" in all material
20 db6608bd 2003-11-23 devnull * mentioning or referencing the derived work.
22 db6608bd 2003-11-23 devnull * RSA Data Security, Inc. makes no representations concerning either
23 db6608bd 2003-11-23 devnull * the merchantability of this software or the suitability of this
24 db6608bd 2003-11-23 devnull * software forany particular purpose. It is provided "as is"
25 db6608bd 2003-11-23 devnull * without express or implied warranty of any kind.
26 db6608bd 2003-11-23 devnull * These notices must be retained in any copies of any part of this
27 db6608bd 2003-11-23 devnull * documentation and/or software.
31 db6608bd 2003-11-23 devnull * Rotate ammounts used in the algorithm
56 db6608bd 2003-11-23 devnull static u32int md5tab[] =
58 db6608bd 2003-11-23 devnull /* round 1 */
59 db6608bd 2003-11-23 devnull /*[0]*/ 0xd76aa478,
60 db6608bd 2003-11-23 devnull 0xe8c7b756,
61 db6608bd 2003-11-23 devnull 0x242070db,
62 db6608bd 2003-11-23 devnull 0xc1bdceee,
63 db6608bd 2003-11-23 devnull 0xf57c0faf,
64 db6608bd 2003-11-23 devnull 0x4787c62a,
65 db6608bd 2003-11-23 devnull 0xa8304613,
66 db6608bd 2003-11-23 devnull 0xfd469501,
67 db6608bd 2003-11-23 devnull 0x698098d8,
68 db6608bd 2003-11-23 devnull 0x8b44f7af,
69 db6608bd 2003-11-23 devnull 0xffff5bb1,
70 db6608bd 2003-11-23 devnull 0x895cd7be,
71 db6608bd 2003-11-23 devnull 0x6b901122,
72 db6608bd 2003-11-23 devnull 0xfd987193,
73 db6608bd 2003-11-23 devnull 0xa679438e,
74 db6608bd 2003-11-23 devnull 0x49b40821,
76 db6608bd 2003-11-23 devnull /* round 2 */
77 db6608bd 2003-11-23 devnull /*[16]*/0xf61e2562,
78 db6608bd 2003-11-23 devnull 0xc040b340,
79 db6608bd 2003-11-23 devnull 0x265e5a51,
80 db6608bd 2003-11-23 devnull 0xe9b6c7aa,
81 db6608bd 2003-11-23 devnull 0xd62f105d,
82 db6608bd 2003-11-23 devnull 0x2441453,
83 db6608bd 2003-11-23 devnull 0xd8a1e681,
84 db6608bd 2003-11-23 devnull 0xe7d3fbc8,
85 db6608bd 2003-11-23 devnull 0x21e1cde6,
86 db6608bd 2003-11-23 devnull 0xc33707d6,
87 db6608bd 2003-11-23 devnull 0xf4d50d87,
88 db6608bd 2003-11-23 devnull 0x455a14ed,
89 db6608bd 2003-11-23 devnull 0xa9e3e905,
90 db6608bd 2003-11-23 devnull 0xfcefa3f8,
91 db6608bd 2003-11-23 devnull 0x676f02d9,
92 db6608bd 2003-11-23 devnull 0x8d2a4c8a,
94 db6608bd 2003-11-23 devnull /* round 3 */
95 db6608bd 2003-11-23 devnull /*[32]*/0xfffa3942,
96 db6608bd 2003-11-23 devnull 0x8771f681,
97 db6608bd 2003-11-23 devnull 0x6d9d6122,
98 db6608bd 2003-11-23 devnull 0xfde5380c,
99 db6608bd 2003-11-23 devnull 0xa4beea44,
100 db6608bd 2003-11-23 devnull 0x4bdecfa9,
101 db6608bd 2003-11-23 devnull 0xf6bb4b60,
102 db6608bd 2003-11-23 devnull 0xbebfbc70,
103 db6608bd 2003-11-23 devnull 0x289b7ec6,
104 db6608bd 2003-11-23 devnull 0xeaa127fa,
105 db6608bd 2003-11-23 devnull 0xd4ef3085,
106 db6608bd 2003-11-23 devnull 0x4881d05,
107 db6608bd 2003-11-23 devnull 0xd9d4d039,
108 db6608bd 2003-11-23 devnull 0xe6db99e5,
109 db6608bd 2003-11-23 devnull 0x1fa27cf8,
110 db6608bd 2003-11-23 devnull 0xc4ac5665,
112 db6608bd 2003-11-23 devnull /* round 4 */
113 db6608bd 2003-11-23 devnull /*[48]*/0xf4292244,
114 db6608bd 2003-11-23 devnull 0x432aff97,
115 db6608bd 2003-11-23 devnull 0xab9423a7,
116 db6608bd 2003-11-23 devnull 0xfc93a039,
117 db6608bd 2003-11-23 devnull 0x655b59c3,
118 db6608bd 2003-11-23 devnull 0x8f0ccc92,
119 db6608bd 2003-11-23 devnull 0xffeff47d,
120 db6608bd 2003-11-23 devnull 0x85845dd1,
121 db6608bd 2003-11-23 devnull 0x6fa87e4f,
122 db6608bd 2003-11-23 devnull 0xfe2ce6e0,
123 db6608bd 2003-11-23 devnull 0xa3014314,
124 db6608bd 2003-11-23 devnull 0x4e0811a1,
125 db6608bd 2003-11-23 devnull 0xf7537e82,
126 db6608bd 2003-11-23 devnull 0xbd3af235,
127 db6608bd 2003-11-23 devnull 0x2ad7d2bb,
128 db6608bd 2003-11-23 devnull 0xeb86d391,
131 db6608bd 2003-11-23 devnull static void decode(u32int*, uchar*, ulong);
132 db6608bd 2003-11-23 devnull extern void _md5block(uchar *p, ulong len, u32int *s);
135 db6608bd 2003-11-23 devnull _md5block(uchar *p, ulong len, u32int *s)
137 db6608bd 2003-11-23 devnull u32int a, b, c, d, sh;
138 db6608bd 2003-11-23 devnull u32int *t;
139 db6608bd 2003-11-23 devnull uchar *end;
140 db6608bd 2003-11-23 devnull u32int x[16];
142 db6608bd 2003-11-23 devnull for(end = p+len; p < end; p += 64){
143 db6608bd 2003-11-23 devnull a = s[0];
144 db6608bd 2003-11-23 devnull b = s[1];
145 db6608bd 2003-11-23 devnull c = s[2];
146 db6608bd 2003-11-23 devnull d = s[3];
148 db6608bd 2003-11-23 devnull decode(x, p, 64);
150 db6608bd 2003-11-23 devnull t = md5tab;
152 db6608bd 2003-11-23 devnull for(; sh != 16; t += 4){
153 db6608bd 2003-11-23 devnull a += ((c ^ d) & b) ^ d;
154 db6608bd 2003-11-23 devnull a += x[sh] + t[0];
155 db6608bd 2003-11-23 devnull a = (a << S11) | (a >> (32 - S11));
158 db6608bd 2003-11-23 devnull d += ((b ^ c) & a) ^ c;
159 db6608bd 2003-11-23 devnull d += x[sh + 1] + t[1];
160 db6608bd 2003-11-23 devnull d = (d << S12) | (d >> (32 - S12));
163 db6608bd 2003-11-23 devnull c += ((a ^ b) & d) ^ b;
164 db6608bd 2003-11-23 devnull c += x[sh + 2] + t[2];
165 db6608bd 2003-11-23 devnull c = (c << S13) | (c >> (32 - S13));
168 db6608bd 2003-11-23 devnull b += ((d ^ a) & c) ^ a;
169 db6608bd 2003-11-23 devnull b += x[sh + 3] + t[3];
170 db6608bd 2003-11-23 devnull b = (b << S14) | (b >> (32 - S14));
173 db6608bd 2003-11-23 devnull sh += 4;
176 db6608bd 2003-11-23 devnull for(; sh != 1+20*4; t += 4){
177 db6608bd 2003-11-23 devnull a += ((b ^ c) & d) ^ c;
178 db6608bd 2003-11-23 devnull a += x[sh & 0xf] + t[0];
179 db6608bd 2003-11-23 devnull a = (a << S21) | (a >> (32 - S21));
182 db6608bd 2003-11-23 devnull d += ((a ^ b) & c) ^ b;
183 db6608bd 2003-11-23 devnull d += x[(sh + 5) & 0xf] + t[1];
184 db6608bd 2003-11-23 devnull d = (d << S22) | (d >> (32 - S22));
187 db6608bd 2003-11-23 devnull c += ((d ^ a) & b) ^ a;
188 db6608bd 2003-11-23 devnull c += x[(sh + 10) & 0xf] + t[2];
189 db6608bd 2003-11-23 devnull c = (c << S23) | (c >> (32 - S23));
192 db6608bd 2003-11-23 devnull b += ((c ^ d) & a) ^ d;
193 db6608bd 2003-11-23 devnull b += x[(sh + 15) & 0xf] + t[3];
194 db6608bd 2003-11-23 devnull b = (b << S24) | (b >> (32 - S24));
197 db6608bd 2003-11-23 devnull sh += 20;
200 db6608bd 2003-11-23 devnull for(; sh != 5+12*4; t += 4){
201 db6608bd 2003-11-23 devnull a += b ^ c ^ d;
202 db6608bd 2003-11-23 devnull a += x[sh & 0xf] + t[0];
203 db6608bd 2003-11-23 devnull a = (a << S31) | (a >> (32 - S31));
206 db6608bd 2003-11-23 devnull d += a ^ b ^ c;
207 db6608bd 2003-11-23 devnull d += x[(sh + 3) & 0xf] + t[1];
208 db6608bd 2003-11-23 devnull d = (d << S32) | (d >> (32 - S32));
211 db6608bd 2003-11-23 devnull c += d ^ a ^ b;
212 db6608bd 2003-11-23 devnull c += x[(sh + 6) & 0xf] + t[2];
213 db6608bd 2003-11-23 devnull c = (c << S33) | (c >> (32 - S33));
216 db6608bd 2003-11-23 devnull b += c ^ d ^ a;
217 db6608bd 2003-11-23 devnull b += x[(sh + 9) & 0xf] + t[3];
218 db6608bd 2003-11-23 devnull b = (b << S34) | (b >> (32 - S34));
221 db6608bd 2003-11-23 devnull sh += 12;
224 db6608bd 2003-11-23 devnull for(; sh != 28*4; t += 4){
225 db6608bd 2003-11-23 devnull a += c ^ (b | ~d);
226 db6608bd 2003-11-23 devnull a += x[sh & 0xf] + t[0];
227 db6608bd 2003-11-23 devnull a = (a << S41) | (a >> (32 - S41));
230 db6608bd 2003-11-23 devnull d += b ^ (a | ~c);
231 db6608bd 2003-11-23 devnull d += x[(sh + 7) & 0xf] + t[1];
232 db6608bd 2003-11-23 devnull d = (d << S42) | (d >> (32 - S42));
235 db6608bd 2003-11-23 devnull c += a ^ (d | ~b);
236 db6608bd 2003-11-23 devnull c += x[(sh + 14) & 0xf] + t[2];
237 db6608bd 2003-11-23 devnull c = (c << S43) | (c >> (32 - S43));
240 db6608bd 2003-11-23 devnull b += d ^ (c | ~a);
241 db6608bd 2003-11-23 devnull b += x[(sh + 21) & 0xf] + t[3];
242 db6608bd 2003-11-23 devnull b = (b << S44) | (b >> (32 - S44));
245 db6608bd 2003-11-23 devnull sh += 28;
248 db6608bd 2003-11-23 devnull s[0] += a;
249 db6608bd 2003-11-23 devnull s[1] += b;
250 db6608bd 2003-11-23 devnull s[2] += c;
251 db6608bd 2003-11-23 devnull s[3] += d;
256 db6608bd 2003-11-23 devnull * decodes input (uchar) into output (u32int). Assumes len is
257 db6608bd 2003-11-23 devnull * a multiple of 4.
259 db6608bd 2003-11-23 devnull static void
260 db6608bd 2003-11-23 devnull decode(u32int *output, uchar *input, ulong len)
262 db6608bd 2003-11-23 devnull uchar *e;
264 db6608bd 2003-11-23 devnull for(e = input+len; input < e; input += 4)
265 db6608bd 2003-11-23 devnull *output++ = input[0] | (input[1] << 8) |
266 db6608bd 2003-11-23 devnull (input[2] << 16) | (input[3] << 24);
