2 * Copyright (c) 2016 Thomas Pornin <pornin@bolet.org>
4 * Permission is hereby granted, free of charge, to any person obtaining
5 * a copy of this software and associated documentation files (the
6 * "Software"), to deal in the Software without restriction, including
7 * without limitation the rights to use, copy, modify, merge, publish,
8 * distribute, sublicense, and/or sell copies of the Software, and to
9 * permit persons to whom the Software is furnished to do so, subject to
10 * the following conditions:
12 * The above copyright notice and this permission notice shall be
13 * included in all copies or substantial portions of the Software.
15 * THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND,
16 * EXPRESS OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF
17 * MERCHANTABILITY, FITNESS FOR A PARTICULAR PURPOSE AND
18 * NONINFRINGEMENT. IN NO EVENT SHALL THE AUTHORS OR COPYRIGHT HOLDERS
19 * BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER LIABILITY, WHETHER IN AN
20 * ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM, OUT OF OR IN
21 * CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE
27 #define F(B, C, D) ((((C) ^ (D)) & (B)) ^ (D))
28 #define G(B, C, D) ((((C) ^ (B)) & (D)) ^ (C))
29 #define H(B, C, D) ((B) ^ (C) ^ (D))
30 #define I(B, C, D) ((C) ^ ((B) | ~(D)))
32 #define ROTL(x, n) (((x) << (n)) | ((x) >> (32 - (n))))
35 const uint32_t br_md5_IV
[4] = {
36 0x67452301, 0xEFCDAB89, 0x98BADCFE, 0x10325476
39 static const uint32_t K
[64] = {
40 0xD76AA478, 0xE8C7B756, 0x242070DB, 0xC1BDCEEE,
41 0xF57C0FAF, 0x4787C62A, 0xA8304613, 0xFD469501,
42 0x698098D8, 0x8B44F7AF, 0xFFFF5BB1, 0x895CD7BE,
43 0x6B901122, 0xFD987193, 0xA679438E, 0x49B40821,
45 0xF61E2562, 0xC040B340, 0x265E5A51, 0xE9B6C7AA,
46 0xD62F105D, 0x02441453, 0xD8A1E681, 0xE7D3FBC8,
47 0x21E1CDE6, 0xC33707D6, 0xF4D50D87, 0x455A14ED,
48 0xA9E3E905, 0xFCEFA3F8, 0x676F02D9, 0x8D2A4C8A,
50 0xFFFA3942, 0x8771F681, 0x6D9D6122, 0xFDE5380C,
51 0xA4BEEA44, 0x4BDECFA9, 0xF6BB4B60, 0xBEBFBC70,
52 0x289B7EC6, 0xEAA127FA, 0xD4EF3085, 0x04881D05,
53 0xD9D4D039, 0xE6DB99E5, 0x1FA27CF8, 0xC4AC5665,
55 0xF4292244, 0x432AFF97, 0xAB9423A7, 0xFC93A039,
56 0x655B59C3, 0x8F0CCC92, 0xFFEFF47D, 0x85845DD1,
57 0x6FA87E4F, 0xFE2CE6E0, 0xA3014314, 0x4E0811A1,
58 0xF7537E82, 0xBD3AF235, 0x2AD7D2BB, 0xEB86D391
61 static const unsigned char MP
[48] = {
62 1, 6, 11, 0, 5, 10, 15, 4, 9, 14, 3, 8, 13, 2, 7, 12,
63 5, 8, 11, 14, 1, 4, 7, 10, 13, 0, 3, 6, 9, 12, 15, 2,
64 0, 7, 14, 5, 12, 3, 10, 1, 8, 15, 6, 13, 4, 11, 2, 9
69 br_md5_round(const unsigned char *buf
, uint32_t *val
)
80 for (i = 0; i < 16; i ++) {
81 m[i] = br_dec32le(buf + (i << 2));
84 br_range_dec32le(m
, 16, buf
);
86 for (i
= 0; i
< 16; i
+= 4) {
87 a
= b
+ ROTL(a
+ F(b
, c
, d
) + m
[i
+ 0] + K
[i
+ 0], 7);
88 d
= a
+ ROTL(d
+ F(a
, b
, c
) + m
[i
+ 1] + K
[i
+ 1], 12);
89 c
= d
+ ROTL(c
+ F(d
, a
, b
) + m
[i
+ 2] + K
[i
+ 2], 17);
90 b
= c
+ ROTL(b
+ F(c
, d
, a
) + m
[i
+ 3] + K
[i
+ 3], 22);
92 for (i
= 16; i
< 32; i
+= 4) {
93 a
= b
+ ROTL(a
+ G(b
, c
, d
) + m
[MP
[i
- 16]] + K
[i
+ 0], 5);
94 d
= a
+ ROTL(d
+ G(a
, b
, c
) + m
[MP
[i
- 15]] + K
[i
+ 1], 9);
95 c
= d
+ ROTL(c
+ G(d
, a
, b
) + m
[MP
[i
- 14]] + K
[i
+ 2], 14);
96 b
= c
+ ROTL(b
+ G(c
, d
, a
) + m
[MP
[i
- 13]] + K
[i
+ 3], 20);
98 for (i
= 32; i
< 48; i
+= 4) {
99 a
= b
+ ROTL(a
+ H(b
, c
, d
) + m
[MP
[i
- 16]] + K
[i
+ 0], 4);
100 d
= a
+ ROTL(d
+ H(a
, b
, c
) + m
[MP
[i
- 15]] + K
[i
+ 1], 11);
101 c
= d
+ ROTL(c
+ H(d
, a
, b
) + m
[MP
[i
- 14]] + K
[i
+ 2], 16);
102 b
= c
+ ROTL(b
+ H(c
, d
, a
) + m
[MP
[i
- 13]] + K
[i
+ 3], 23);
104 for (i
= 48; i
< 64; i
+= 4) {
105 a
= b
+ ROTL(a
+ I(b
, c
, d
) + m
[MP
[i
- 16]] + K
[i
+ 0], 6);
106 d
= a
+ ROTL(d
+ I(a
, b
, c
) + m
[MP
[i
- 15]] + K
[i
+ 1], 10);
107 c
= d
+ ROTL(c
+ I(d
, a
, b
) + m
[MP
[i
- 14]] + K
[i
+ 2], 15);
108 b
= c
+ ROTL(b
+ I(c
, d
, a
) + m
[MP
[i
- 13]] + K
[i
+ 3], 21);
119 br_md5_init(br_md5_context
*cc
)
121 cc
->vtable
= &br_md5_vtable
;
122 memcpy(cc
->val
, br_md5_IV
, sizeof cc
->val
);
128 br_md5_update(br_md5_context
*cc
, const void *data
, size_t len
)
130 const unsigned char *buf
;
134 ptr
= (size_t)cc
->count
& 63;
142 memcpy(cc
->buf
+ ptr
, buf
, clen
);
146 cc
->count
+= (uint64_t)clen
;
148 br_md5_round(cc
->buf
, cc
->val
);
156 br_md5_out(const br_md5_context
*cc
, void *dst
)
158 unsigned char buf
[64];
162 ptr
= (size_t)cc
->count
& 63;
163 memcpy(buf
, cc
->buf
, ptr
);
164 memcpy(val
, cc
->val
, sizeof val
);
167 memset(buf
+ ptr
, 0, 64 - ptr
);
168 br_md5_round(buf
, val
);
171 memset(buf
+ ptr
, 0, 56 - ptr
);
173 br_enc64le(buf
+ 56, cc
->count
<< 3);
174 br_md5_round(buf
, val
);
175 br_range_enc32le(dst
, val
, 4);
180 br_md5_state(const br_md5_context
*cc
, void *dst
)
182 br_range_enc32le(dst
, cc
->val
, 4);
188 br_md5_set_state(br_md5_context
*cc
, const void *stb
, uint64_t count
)
190 br_range_dec32le(cc
->val
, 4, stb
);
195 const br_hash_class br_md5_vtable
= {
196 sizeof(br_md5_context
),
197 BR_HASHDESC_ID(br_md5_ID
)
198 | BR_HASHDESC_OUT(16)
199 | BR_HASHDESC_STATE(16)
200 | BR_HASHDESC_LBLEN(6)
201 | BR_HASHDESC_MD_PADDING
,
202 (void (*)(const br_hash_class
**))&br_md5_init
,
203 (void (*)(const br_hash_class
**, const void *, size_t))&br_md5_update
,
204 (void (*)(const br_hash_class
*const *, void *))&br_md5_out
,
205 (uint64_t (*)(const br_hash_class
*const *, void *))&br_md5_state
,
206 (void (*)(const br_hash_class
**, const void *, uint64_t))