1 /* Copyright (C) 1995-1998 Eric Young (eay@cryptsoft.com)
2 * All rights reserved.
3 *
4 * This package is an SSL implementation written
5 * by Eric Young (eay@cryptsoft.com).
6 * The implementation was written so as to conform with Netscapes SSL.
7 *
8 * This library is free for commercial and non-commercial use as long as
9 * the following conditions are aheared to. The following conditions
10 * apply to all code found in this distribution, be it the RC4, RSA,
11 * lhash, DES, etc., code; not just the SSL code. The SSL documentation
12 * included with this distribution is covered by the same copyright terms
13 * except that the holder is Tim Hudson (tjh@cryptsoft.com).
14 *
15 * Copyright remains Eric Young's, and as such any Copyright notices in
16 * the code are not to be removed.
17 * If this package is used in a product, Eric Young should be given attribution
18 * as the author of the parts of the library used.
19 * This can be in the form of a textual message at program startup or
20 * in documentation (online or textual) provided with the package.
21 *
22 * Redistribution and use in source and binary forms, with or without
23 * modification, are permitted provided that the following conditions
24 * are met:
25 * 1. Redistributions of source code must retain the copyright
26 * notice, this list of conditions and the following disclaimer.
27 * 2. Redistributions in binary form must reproduce the above copyright
28 * notice, this list of conditions and the following disclaimer in the
29 * documentation and/or other materials provided with the distribution.
30 * 3. All advertising materials mentioning features or use of this software
31 * must display the following acknowledgement:
32 * "This product includes cryptographic software written by
33 * Eric Young (eay@cryptsoft.com)"
34 * The word 'cryptographic' can be left out if the rouines from the library
35 * being used are not cryptographic related :-).
36 * 4. If you include any Windows specific code (or a derivative thereof) from
37 * the apps directory (application code) you must include an acknowledgement:
38 * "This product includes software written by Tim Hudson (tjh@cryptsoft.com)"
39 *
40 * THIS SOFTWARE IS PROVIDED BY ERIC YOUNG ``AS IS'' AND
41 * ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
42 * IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
43 * ARE DISCLAIMED. IN NO EVENT SHALL THE AUTHOR OR CONTRIBUTORS BE LIABLE
44 * FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL
45 * DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS
46 * OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION)
47 * HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT
48 * LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY
49 * OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF
50 * SUCH DAMAGE.
51 *
52 * The licence and distribution terms for any publically available version or
53 * derivative of this code cannot be changed. i.e. this code cannot simply be
54 * copied and put under another distribution licence
55 * [including the GNU Public Licence.] */
56
57 #include <openssl/md5.h>
58
59 #include <string.h>
60
61 #include <openssl/mem.h>
62
63 #include "../../internal.h"
64 #include "../digest/md32_common.h"
65 #include "internal.h"
66
67
MD5(const uint8_t * data,size_t len,uint8_t out[MD5_DIGEST_LENGTH])68 uint8_t *MD5(const uint8_t *data, size_t len, uint8_t out[MD5_DIGEST_LENGTH]) {
69 MD5_CTX ctx;
70 MD5_Init(&ctx);
71 MD5_Update(&ctx, data, len);
72 MD5_Final(out, &ctx);
73
74 return out;
75 }
76
MD5_Init(MD5_CTX * md5)77 int MD5_Init(MD5_CTX *md5) {
78 OPENSSL_memset(md5, 0, sizeof(MD5_CTX));
79 md5->h[0] = 0x67452301UL;
80 md5->h[1] = 0xefcdab89UL;
81 md5->h[2] = 0x98badcfeUL;
82 md5->h[3] = 0x10325476UL;
83 return 1;
84 }
85
86 #if defined(MD5_ASM)
87 #define md5_block_data_order md5_block_asm_data_order
88 #else
89 static void md5_block_data_order(uint32_t *state, const uint8_t *data,
90 size_t num);
91 #endif
92
MD5_Transform(MD5_CTX * c,const uint8_t data[MD5_CBLOCK])93 void MD5_Transform(MD5_CTX *c, const uint8_t data[MD5_CBLOCK]) {
94 md5_block_data_order(c->h, data, 1);
95 }
96
MD5_Update(MD5_CTX * c,const void * data,size_t len)97 int MD5_Update(MD5_CTX *c, const void *data, size_t len) {
98 crypto_md32_update(&md5_block_data_order, c->h, c->data, MD5_CBLOCK, &c->num,
99 &c->Nh, &c->Nl, data, len);
100 return 1;
101 }
102
MD5_Final(uint8_t out[MD5_DIGEST_LENGTH],MD5_CTX * c)103 int MD5_Final(uint8_t out[MD5_DIGEST_LENGTH], MD5_CTX *c) {
104 crypto_md32_final(&md5_block_data_order, c->h, c->data, MD5_CBLOCK, &c->num,
105 c->Nh, c->Nl, /*is_big_endian=*/0);
106
107 CRYPTO_store_u32_le(out, c->h[0]);
108 CRYPTO_store_u32_le(out + 4, c->h[1]);
109 CRYPTO_store_u32_le(out + 8, c->h[2]);
110 CRYPTO_store_u32_le(out + 12, c->h[3]);
111 return 1;
112 }
113
114 // As pointed out by Wei Dai <weidai@eskimo.com>, the above can be
115 // simplified to the code below. Wei attributes these optimizations
116 // to Peter Gutmann's SHS code, and he attributes it to Rich Schroeppel.
117 #define F(b, c, d) ((((c) ^ (d)) & (b)) ^ (d))
118 #define G(b, c, d) ((((b) ^ (c)) & (d)) ^ (c))
119 #define H(b, c, d) ((b) ^ (c) ^ (d))
120 #define I(b, c, d) (((~(d)) | (b)) ^ (c))
121
122 #define ROTATE(a, n) (((a) << (n)) | ((a) >> (32 - (n))))
123
124 #define R0(a, b, c, d, k, s, t) \
125 do { \
126 (a) += ((k) + (t) + F((b), (c), (d))); \
127 (a) = ROTATE(a, s); \
128 (a) += (b); \
129 } while (0)
130
131 #define R1(a, b, c, d, k, s, t) \
132 do { \
133 (a) += ((k) + (t) + G((b), (c), (d))); \
134 (a) = ROTATE(a, s); \
135 (a) += (b); \
136 } while (0)
137
138 #define R2(a, b, c, d, k, s, t) \
139 do { \
140 (a) += ((k) + (t) + H((b), (c), (d))); \
141 (a) = ROTATE(a, s); \
142 (a) += (b); \
143 } while (0)
144
145 #define R3(a, b, c, d, k, s, t) \
146 do { \
147 (a) += ((k) + (t) + I((b), (c), (d))); \
148 (a) = ROTATE(a, s); \
149 (a) += (b); \
150 } while (0)
151
152 #ifndef MD5_ASM
153 #ifdef X
154 #undef X
155 #endif
md5_block_data_order(uint32_t * state,const uint8_t * data,size_t num)156 static void md5_block_data_order(uint32_t *state, const uint8_t *data,
157 size_t num) {
158 uint32_t A, B, C, D;
159 uint32_t XX0, XX1, XX2, XX3, XX4, XX5, XX6, XX7, XX8, XX9, XX10, XX11, XX12,
160 XX13, XX14, XX15;
161 #define X(i) XX##i
162
163 A = state[0];
164 B = state[1];
165 C = state[2];
166 D = state[3];
167
168 for (; num--;) {
169 X(0) = CRYPTO_load_u32_le(data);
170 data += 4;
171 X(1) = CRYPTO_load_u32_le(data);
172 data += 4;
173 // Round 0
174 R0(A, B, C, D, X(0), 7, 0xd76aa478L);
175 X(2) = CRYPTO_load_u32_le(data);
176 data += 4;
177 R0(D, A, B, C, X(1), 12, 0xe8c7b756L);
178 X(3) = CRYPTO_load_u32_le(data);
179 data += 4;
180 R0(C, D, A, B, X(2), 17, 0x242070dbL);
181 X(4) = CRYPTO_load_u32_le(data);
182 data += 4;
183 R0(B, C, D, A, X(3), 22, 0xc1bdceeeL);
184 X(5) = CRYPTO_load_u32_le(data);
185 data += 4;
186 R0(A, B, C, D, X(4), 7, 0xf57c0fafL);
187 X(6) = CRYPTO_load_u32_le(data);
188 data += 4;
189 R0(D, A, B, C, X(5), 12, 0x4787c62aL);
190 X(7) = CRYPTO_load_u32_le(data);
191 data += 4;
192 R0(C, D, A, B, X(6), 17, 0xa8304613L);
193 X(8) = CRYPTO_load_u32_le(data);
194 data += 4;
195 R0(B, C, D, A, X(7), 22, 0xfd469501L);
196 X(9) = CRYPTO_load_u32_le(data);
197 data += 4;
198 R0(A, B, C, D, X(8), 7, 0x698098d8L);
199 X(10) = CRYPTO_load_u32_le(data);
200 data += 4;
201 R0(D, A, B, C, X(9), 12, 0x8b44f7afL);
202 X(11) = CRYPTO_load_u32_le(data);
203 data += 4;
204 R0(C, D, A, B, X(10), 17, 0xffff5bb1L);
205 X(12) = CRYPTO_load_u32_le(data);
206 data += 4;
207 R0(B, C, D, A, X(11), 22, 0x895cd7beL);
208 X(13) = CRYPTO_load_u32_le(data);
209 data += 4;
210 R0(A, B, C, D, X(12), 7, 0x6b901122L);
211 X(14) = CRYPTO_load_u32_le(data);
212 data += 4;
213 R0(D, A, B, C, X(13), 12, 0xfd987193L);
214 X(15) = CRYPTO_load_u32_le(data);
215 data += 4;
216 R0(C, D, A, B, X(14), 17, 0xa679438eL);
217 R0(B, C, D, A, X(15), 22, 0x49b40821L);
218 // Round 1
219 R1(A, B, C, D, X(1), 5, 0xf61e2562L);
220 R1(D, A, B, C, X(6), 9, 0xc040b340L);
221 R1(C, D, A, B, X(11), 14, 0x265e5a51L);
222 R1(B, C, D, A, X(0), 20, 0xe9b6c7aaL);
223 R1(A, B, C, D, X(5), 5, 0xd62f105dL);
224 R1(D, A, B, C, X(10), 9, 0x02441453L);
225 R1(C, D, A, B, X(15), 14, 0xd8a1e681L);
226 R1(B, C, D, A, X(4), 20, 0xe7d3fbc8L);
227 R1(A, B, C, D, X(9), 5, 0x21e1cde6L);
228 R1(D, A, B, C, X(14), 9, 0xc33707d6L);
229 R1(C, D, A, B, X(3), 14, 0xf4d50d87L);
230 R1(B, C, D, A, X(8), 20, 0x455a14edL);
231 R1(A, B, C, D, X(13), 5, 0xa9e3e905L);
232 R1(D, A, B, C, X(2), 9, 0xfcefa3f8L);
233 R1(C, D, A, B, X(7), 14, 0x676f02d9L);
234 R1(B, C, D, A, X(12), 20, 0x8d2a4c8aL);
235 // Round 2
236 R2(A, B, C, D, X(5), 4, 0xfffa3942L);
237 R2(D, A, B, C, X(8), 11, 0x8771f681L);
238 R2(C, D, A, B, X(11), 16, 0x6d9d6122L);
239 R2(B, C, D, A, X(14), 23, 0xfde5380cL);
240 R2(A, B, C, D, X(1), 4, 0xa4beea44L);
241 R2(D, A, B, C, X(4), 11, 0x4bdecfa9L);
242 R2(C, D, A, B, X(7), 16, 0xf6bb4b60L);
243 R2(B, C, D, A, X(10), 23, 0xbebfbc70L);
244 R2(A, B, C, D, X(13), 4, 0x289b7ec6L);
245 R2(D, A, B, C, X(0), 11, 0xeaa127faL);
246 R2(C, D, A, B, X(3), 16, 0xd4ef3085L);
247 R2(B, C, D, A, X(6), 23, 0x04881d05L);
248 R2(A, B, C, D, X(9), 4, 0xd9d4d039L);
249 R2(D, A, B, C, X(12), 11, 0xe6db99e5L);
250 R2(C, D, A, B, X(15), 16, 0x1fa27cf8L);
251 R2(B, C, D, A, X(2), 23, 0xc4ac5665L);
252 // Round 3
253 R3(A, B, C, D, X(0), 6, 0xf4292244L);
254 R3(D, A, B, C, X(7), 10, 0x432aff97L);
255 R3(C, D, A, B, X(14), 15, 0xab9423a7L);
256 R3(B, C, D, A, X(5), 21, 0xfc93a039L);
257 R3(A, B, C, D, X(12), 6, 0x655b59c3L);
258 R3(D, A, B, C, X(3), 10, 0x8f0ccc92L);
259 R3(C, D, A, B, X(10), 15, 0xffeff47dL);
260 R3(B, C, D, A, X(1), 21, 0x85845dd1L);
261 R3(A, B, C, D, X(8), 6, 0x6fa87e4fL);
262 R3(D, A, B, C, X(15), 10, 0xfe2ce6e0L);
263 R3(C, D, A, B, X(6), 15, 0xa3014314L);
264 R3(B, C, D, A, X(13), 21, 0x4e0811a1L);
265 R3(A, B, C, D, X(4), 6, 0xf7537e82L);
266 R3(D, A, B, C, X(11), 10, 0xbd3af235L);
267 R3(C, D, A, B, X(2), 15, 0x2ad7d2bbL);
268 R3(B, C, D, A, X(9), 21, 0xeb86d391L);
269
270 A = state[0] += A;
271 B = state[1] += B;
272 C = state[2] += C;
273 D = state[3] += D;
274 }
275 }
276 #undef X
277 #endif
278
279 #undef F
280 #undef G
281 #undef H
282 #undef I
283 #undef ROTATE
284 #undef R0
285 #undef R1
286 #undef R2
287 #undef R3
288