1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21
22
23
24
25 #ifdef HAVE_CONFIG_H
26 # include <config.h>
27 #endif
28
29 #include <sys/types.h>
30
31 #if STDC_HEADERS || defined _LIBC
32 # include <stdlib.h>
33 # include <string.h>
34 #else
35 # ifndef HAVE_MEMCPY
36 # define memcpy(d, s, n) bcopy ((s), (d), (n))
37 # endif
38 #endif
39
40 #ifdef _LIBC
41 # include <endian.h>
42 # if __BYTE_ORDER == __BIG_ENDIAN
43 # define WORDS_BIG_ENDIAN 1
44 # endif
45 46
47 # define md5_init_ctx __md5_init_ctx
48 # define md5_process_block __md5_process_block
49 # define md5_process_bytes __md5_process_bytes
50 # define md5_finish_ctx __md5_finish_ctx
51 # define md5_read_ctx __md5_read_ctx
52 # define md5_stream __md5_stream
53 # define md5_buffer __md5_buffer
54 #endif
55
56 #include "md5.h"
57
58 #ifdef WORDS_BIG_ENDIAN
59 # define SWAP(n) \
60 (((n) << 24) | (((n) & 0xff00) << 8) | (((n) >> 8) & 0xff00) | ((n) >> 24))
61 #else
62 # define SWAP(n) (n)
63 #endif
64
65
66 67
68 static const unsigned char fillbuf[64] = { 0x80, 0 };
69
70
71 72
73 void
74 md5_init_ctx (ctx)
75 struct md5_ctx *ctx;
76 {
77 ctx->A = 0x67452301;
78 ctx->B = 0xefcdab89;
79 ctx->C = 0x98badcfe;
80 ctx->D = 0x10325476;
81
82 ctx->total[0] = ctx->total[1] = 0;
83 ctx->buflen = 0;
84 }
85
86 87 88 89 90
91 void *
92 md5_read_ctx (ctx, resbuf)
93 const struct md5_ctx *ctx;
94 void *resbuf;
95 {
96 ((md5_uint32 *) resbuf)[0] = SWAP (ctx->A);
97 ((md5_uint32 *) resbuf)[1] = SWAP (ctx->B);
98 ((md5_uint32 *) resbuf)[2] = SWAP (ctx->C);
99 ((md5_uint32 *) resbuf)[3] = SWAP (ctx->D);
100
101 return resbuf;
102 }
103
104 105 106 107 108
109 void *
110 md5_finish_ctx (ctx, resbuf)
111 struct md5_ctx *ctx;
112 void *resbuf;
113 {
114
115 md5_uint32 bytes = ctx->buflen;
116 size_t pad;
117
118
119 ctx->total[0] += bytes;
120 if (ctx->total[0] < bytes)
121 ++ctx->total[1];
122
123 pad = bytes >= 56 ? 64 + 56 - bytes : 56 - bytes;
124 memcpy (&ctx->buffer[bytes], fillbuf, pad);
125
126
127 *(md5_uint32 *) &ctx->buffer[bytes + pad] = SWAP (ctx->total[0] << 3);
128 *(md5_uint32 *) &ctx->buffer[bytes + pad + 4] = SWAP ((ctx->total[1] << 3) |
129 (ctx->total[0] >> 29));
130
131
132 md5_process_block (ctx->buffer, bytes + pad + 8, ctx);
133
134 return md5_read_ctx (ctx, resbuf);
135 }
136
137 138 139
140 int
141 md5_stream (stream, resblock)
142 FILE *stream;
143 void *resblock;
144 {
145
146 #define BLOCKSIZE 4096
147 struct md5_ctx ctx;
148 char buffer[BLOCKSIZE + 72];
149 size_t sum;
150
151
152 md5_init_ctx (&ctx);
153
154
155 while (1)
156 {
157 158 159
160 size_t n;
161 sum = 0;
162
163
164 do
165 {
166 n = fread (buffer + sum, 1, BLOCKSIZE - sum, stream);
167
168 sum += n;
169 }
170 while (sum < BLOCKSIZE && n != 0);
171 if (n == 0 && ferror (stream))
172 return 1;
173
174
175 if (n == 0)
176 break;
177
178 179 180
181 md5_process_block (buffer, BLOCKSIZE, &ctx);
182 }
183
184
185 if (sum > 0)
186 md5_process_bytes (buffer, sum, &ctx);
187
188
189 md5_finish_ctx (&ctx, resblock);
190 return 0;
191 }
192
193 194 195 196
197 void *
198 md5_buffer (buffer, len, resblock)
199 const char *buffer;
200 size_t len;
201 void *resblock;
202 {
203 struct md5_ctx ctx;
204
205
206 md5_init_ctx (&ctx);
207
208
209 md5_process_bytes (buffer, len, &ctx);
210
211
212 return md5_finish_ctx (&ctx, resblock);
213 }
214
215
216 void
217 md5_process_bytes (buffer, len, ctx)
218 const void *buffer;
219 size_t len;
220 struct md5_ctx *ctx;
221 {
222
223
224 225
226 if (ctx->buflen != 0)
227 {
228 size_t left_over = ctx->buflen;
229 size_t add = 128 - left_over > len ? len : 128 - left_over;
230
231
232 add -= add % __alignof__ (md5_uint32);
233
234 memcpy (&ctx->buffer[left_over], buffer, add);
235 ctx->buflen += add;
236
237 if (ctx->buflen > 64)
238 {
239 md5_process_block (ctx->buffer, ctx->buflen & ~63, ctx);
240
241 ctx->buflen &= 63;
242
243 memcpy (ctx->buffer, &ctx->buffer[(left_over + add) & ~63],
244 ctx->buflen);
245 }
246
247 buffer = (const char *) buffer + add;
248 len -= add;
249 }
250
251
252 if (len > 64)
253 {
254 md5_process_block (buffer, len & ~63, ctx);
255 buffer = (const char *) buffer + (len & ~63);
256 len &= 63;
257 }
258
259
260 if (len > 0)
261 {
262 size_t left_over = ctx->buflen;
263
264 memcpy (&ctx->buffer[left_over], buffer, len);
265 left_over += len;
266 if (left_over >= 64)
267 {
268 md5_process_block (ctx->buffer, 64, ctx);
269 left_over -= 64;
270 memcpy (ctx->buffer, &ctx->buffer[64], left_over);
271 }
272 ctx->buflen = left_over;
273 }
274 }
275
276
277 278 279
280
281 #define FF(b, c, d) (d ^ (b & (c ^ d)))
282 #define FG(b, c, d) FF (d, b, c)
283 #define FH(b, c, d) (b ^ c ^ d)
284 #define FI(b, c, d) (c ^ (b | ~d))
285
286 287
288
289 void
290 md5_process_block (buffer, len, ctx)
291 const void *buffer;
292 size_t len;
293 struct md5_ctx *ctx;
294 {
295 md5_uint32 correct_words[16];
296 const md5_uint32 *words = buffer;
297 size_t nwords = len / sizeof (md5_uint32);
298 const md5_uint32 *endp = words + nwords;
299 md5_uint32 A = ctx->A;
300 md5_uint32 B = ctx->B;
301 md5_uint32 C = ctx->C;
302 md5_uint32 D = ctx->D;
303
304 305 306
307 ctx->total[0] += len;
308 if (ctx->total[0] < len)
309 ++ctx->total[1];
310
311 312
313 while (words < endp)
314 {
315 md5_uint32 *cwp = correct_words;
316 md5_uint32 A_save = A;
317 md5_uint32 B_save = B;
318 md5_uint32 C_save = C;
319 md5_uint32 D_save = D;
320
321 322 323 324 325 326
327
328 #define OP(a, b, c, d, s, T) \
329 do \
330 { \
331 a += FF (b, c, d) + (*cwp++ = SWAP (*words)) + T; \
332 ++words; \
333 CYCLIC (a, s); \
334 a += b; \
335 } \
336 while (0)
337
338 339
340 #define CYCLIC(w, s) (w = (w << s) | (w >> (32 - s)))
341
342 343 344 345 346
347
348
349 OP (A, B, C, D, 7, 0xd76aa478);
350 OP (D, A, B, C, 12, 0xe8c7b756);
351 OP (C, D, A, B, 17, 0x242070db);
352 OP (B, C, D, A, 22, 0xc1bdceee);
353 OP (A, B, C, D, 7, 0xf57c0faf);
354 OP (D, A, B, C, 12, 0x4787c62a);
355 OP (C, D, A, B, 17, 0xa8304613);
356 OP (B, C, D, A, 22, 0xfd469501);
357 OP (A, B, C, D, 7, 0x698098d8);
358 OP (D, A, B, C, 12, 0x8b44f7af);
359 OP (C, D, A, B, 17, 0xffff5bb1);
360 OP (B, C, D, A, 22, 0x895cd7be);
361 OP (A, B, C, D, 7, 0x6b901122);
362 OP (D, A, B, C, 12, 0xfd987193);
363 OP (C, D, A, B, 17, 0xa679438e);
364 OP (B, C, D, A, 22, 0x49b40821);
365
366 367 368
369 #undef OP
370 #define OP(f, a, b, c, d, k, s, T) \
371 do \
372 { \
373 a += f (b, c, d) + correct_words[k] + T; \
374 CYCLIC (a, s); \
375 a += b; \
376 } \
377 while (0)
378
379
380 OP (FG, A, B, C, D, 1, 5, 0xf61e2562);
381 OP (FG, D, A, B, C, 6, 9, 0xc040b340);
382 OP (FG, C, D, A, B, 11, 14, 0x265e5a51);
383 OP (FG, B, C, D, A, 0, 20, 0xe9b6c7aa);
384 OP (FG, A, B, C, D, 5, 5, 0xd62f105d);
385 OP (FG, D, A, B, C, 10, 9, 0x02441453);
386 OP (FG, C, D, A, B, 15, 14, 0xd8a1e681);
387 OP (FG, B, C, D, A, 4, 20, 0xe7d3fbc8);
388 OP (FG, A, B, C, D, 9, 5, 0x21e1cde6);
389 OP (FG, D, A, B, C, 14, 9, 0xc33707d6);
390 OP (FG, C, D, A, B, 3, 14, 0xf4d50d87);
391 OP (FG, B, C, D, A, 8, 20, 0x455a14ed);
392 OP (FG, A, B, C, D, 13, 5, 0xa9e3e905);
393 OP (FG, D, A, B, C, 2, 9, 0xfcefa3f8);
394 OP (FG, C, D, A, B, 7, 14, 0x676f02d9);
395 OP (FG, B, C, D, A, 12, 20, 0x8d2a4c8a);
396
397
398 OP (FH, A, B, C, D, 5, 4, 0xfffa3942);
399 OP (FH, D, A, B, C, 8, 11, 0x8771f681);
400 OP (FH, C, D, A, B, 11, 16, 0x6d9d6122);
401 OP (FH, B, C, D, A, 14, 23, 0xfde5380c);
402 OP (FH, A, B, C, D, 1, 4, 0xa4beea44);
403 OP (FH, D, A, B, C, 4, 11, 0x4bdecfa9);
404 OP (FH, C, D, A, B, 7, 16, 0xf6bb4b60);
405 OP (FH, B, C, D, A, 10, 23, 0xbebfbc70);
406 OP (FH, A, B, C, D, 13, 4, 0x289b7ec6);
407 OP (FH, D, A, B, C, 0, 11, 0xeaa127fa);
408 OP (FH, C, D, A, B, 3, 16, 0xd4ef3085);
409 OP (FH, B, C, D, A, 6, 23, 0x04881d05);
410 OP (FH, A, B, C, D, 9, 4, 0xd9d4d039);
411 OP (FH, D, A, B, C, 12, 11, 0xe6db99e5);
412 OP (FH, C, D, A, B, 15, 16, 0x1fa27cf8);
413 OP (FH, B, C, D, A, 2, 23, 0xc4ac5665);
414
415
416 OP (FI, A, B, C, D, 0, 6, 0xf4292244);
417 OP (FI, D, A, B, C, 7, 10, 0x432aff97);
418 OP (FI, C, D, A, B, 14, 15, 0xab9423a7);
419 OP (FI, B, C, D, A, 5, 21, 0xfc93a039);
420 OP (FI, A, B, C, D, 12, 6, 0x655b59c3);
421 OP (FI, D, A, B, C, 3, 10, 0x8f0ccc92);
422 OP (FI, C, D, A, B, 10, 15, 0xffeff47d);
423 OP (FI, B, C, D, A, 1, 21, 0x85845dd1);
424 OP (FI, A, B, C, D, 8, 6, 0x6fa87e4f);
425 OP (FI, D, A, B, C, 15, 10, 0xfe2ce6e0);
426 OP (FI, C, D, A, B, 6, 15, 0xa3014314);
427 OP (FI, B, C, D, A, 13, 21, 0x4e0811a1);
428 OP (FI, A, B, C, D, 4, 6, 0xf7537e82);
429 OP (FI, D, A, B, C, 11, 10, 0xbd3af235);
430 OP (FI, C, D, A, B, 2, 15, 0x2ad7d2bb);
431 OP (FI, B, C, D, A, 9, 21, 0xeb86d391);
432
433
434 A += A_save;
435 B += B_save;
436 C += C_save;
437 D += D_save;
438 }
439
440
441 ctx->A = A;
442 ctx->B = B;
443 ctx->C = C;
444 ctx->D = D;
445 }
446
447 448