1 /* crypto/sha/sha512.c */
2 /* ====================================================================
3 * Copyright (c) 2004 The OpenSSL Project. All rights reserved
4 * according to the OpenSSL license [found in ../../LICENSE].
5 * ====================================================================
7 #include <openssl/opensslconf.h>
8 #if !defined(OPENSSL_NO_SHA) && !defined(OPENSSL_NO_SHA512)
10 * IMPLEMENTATION NOTES.
12 * As you might have noticed 32-bit hash algorithms:
14 * - permit SHA_LONG to be wider than 32-bit (case on CRAY);
15 * - optimized versions implement two transform functions: one operating
16 * on [aligned] data in host byte order and one - on data in input
18 * - share common byte-order neutral collector and padding function
19 * implementations, ../md32_common.h;
21 * Neither of the above applies to this SHA-512 implementations. Reasons
22 * [in reverse order] are:
24 * - it's the only 64-bit hash algorithm for the moment of this writing,
25 * there is no need for common collector/padding implementation [yet];
26 * - by supporting only one transform function [which operates on
27 * *aligned* data in input stream byte order, big-endian in this case]
28 * we minimize burden of maintenance in two ways: a) collector/padding
29 * function is simpler; b) only one transform function to stare at;
30 * - SHA_LONG64 is required to be exactly 64-bit in order to be able to
31 * apply a number of optimizations to mitigate potential performance
32 * penalties caused by previous design decision;
36 * Implementation relies on the fact that "long long" is 64-bit on
37 * both 32- and 64-bit platforms. If some compiler vendor comes up
38 * with 128-bit long long, adjustment to sha.h would be required.
39 * As this implementation relies on 64-bit integer type, it's totally
40 * inappropriate for platforms which don't support it, most notably
42 * <appro@fy.chalmers.se>
47 # include <openssl/crypto.h>
48 # include <openssl/sha.h>
49 # include <openssl/opensslv.h>
51 # include "cryptlib.h"
53 const char SHA512_version[] = "SHA-512" OPENSSL_VERSION_PTEXT;
55 # if defined(__i386) || defined(__i386__) || defined(_M_IX86) || \
56 defined(__x86_64) || defined(_M_AMD64) || defined(_M_X64) || \
57 defined(__s390__) || defined(__s390x__) || \
59 # define SHA512_BLOCK_CAN_MANAGE_UNALIGNED_DATA
62 fips_md_init_ctx(SHA384, SHA512)
64 c->h[0] = U64(0xcbbb9d5dc1059ed8);
65 c->h[1] = U64(0x629a292a367cd507);
66 c->h[2] = U64(0x9159015a3070dd17);
67 c->h[3] = U64(0x152fecd8f70e5939);
68 c->h[4] = U64(0x67332667ffc00b31);
69 c->h[5] = U64(0x8eb44a8768581511);
70 c->h[6] = U64(0xdb0c2e0d64f98fa7);
71 c->h[7] = U64(0x47b5481dbefa4fa4);
76 c->md_len = SHA384_DIGEST_LENGTH;
82 c->h[0] = U64(0x6a09e667f3bcc908);
83 c->h[1] = U64(0xbb67ae8584caa73b);
84 c->h[2] = U64(0x3c6ef372fe94f82b);
85 c->h[3] = U64(0xa54ff53a5f1d36f1);
86 c->h[4] = U64(0x510e527fade682d1);
87 c->h[5] = U64(0x9b05688c2b3e6c1f);
88 c->h[6] = U64(0x1f83d9abfb41bd6b);
89 c->h[7] = U64(0x5be0cd19137e2179);
94 c->md_len = SHA512_DIGEST_LENGTH;
101 void sha512_block_data_order(SHA512_CTX *ctx, const void *in, size_t num);
103 int SHA512_Final(unsigned char *md, SHA512_CTX *c)
105 unsigned char *p = (unsigned char *)c->u.p;
108 p[n] = 0x80; /* There always is a room for one */
110 if (n > (sizeof(c->u) - 16))
111 memset(p + n, 0, sizeof(c->u) - n), n = 0,
112 sha512_block_data_order(c, p, 1);
114 memset(p + n, 0, sizeof(c->u) - 16 - n);
116 c->u.d[SHA_LBLOCK - 2] = c->Nh;
117 c->u.d[SHA_LBLOCK - 1] = c->Nl;
119 p[sizeof(c->u) - 1] = (unsigned char)(c->Nl);
120 p[sizeof(c->u) - 2] = (unsigned char)(c->Nl >> 8);
121 p[sizeof(c->u) - 3] = (unsigned char)(c->Nl >> 16);
122 p[sizeof(c->u) - 4] = (unsigned char)(c->Nl >> 24);
123 p[sizeof(c->u) - 5] = (unsigned char)(c->Nl >> 32);
124 p[sizeof(c->u) - 6] = (unsigned char)(c->Nl >> 40);
125 p[sizeof(c->u) - 7] = (unsigned char)(c->Nl >> 48);
126 p[sizeof(c->u) - 8] = (unsigned char)(c->Nl >> 56);
127 p[sizeof(c->u) - 9] = (unsigned char)(c->Nh);
128 p[sizeof(c->u) - 10] = (unsigned char)(c->Nh >> 8);
129 p[sizeof(c->u) - 11] = (unsigned char)(c->Nh >> 16);
130 p[sizeof(c->u) - 12] = (unsigned char)(c->Nh >> 24);
131 p[sizeof(c->u) - 13] = (unsigned char)(c->Nh >> 32);
132 p[sizeof(c->u) - 14] = (unsigned char)(c->Nh >> 40);
133 p[sizeof(c->u) - 15] = (unsigned char)(c->Nh >> 48);
134 p[sizeof(c->u) - 16] = (unsigned char)(c->Nh >> 56);
137 sha512_block_data_order(c, p, 1);
143 /* Let compiler decide if it's appropriate to unroll... */
144 case SHA384_DIGEST_LENGTH:
145 for (n = 0; n < SHA384_DIGEST_LENGTH / 8; n++) {
146 SHA_LONG64 t = c->h[n];
148 *(md++) = (unsigned char)(t >> 56);
149 *(md++) = (unsigned char)(t >> 48);
150 *(md++) = (unsigned char)(t >> 40);
151 *(md++) = (unsigned char)(t >> 32);
152 *(md++) = (unsigned char)(t >> 24);
153 *(md++) = (unsigned char)(t >> 16);
154 *(md++) = (unsigned char)(t >> 8);
155 *(md++) = (unsigned char)(t);
158 case SHA512_DIGEST_LENGTH:
159 for (n = 0; n < SHA512_DIGEST_LENGTH / 8; n++) {
160 SHA_LONG64 t = c->h[n];
162 *(md++) = (unsigned char)(t >> 56);
163 *(md++) = (unsigned char)(t >> 48);
164 *(md++) = (unsigned char)(t >> 40);
165 *(md++) = (unsigned char)(t >> 32);
166 *(md++) = (unsigned char)(t >> 24);
167 *(md++) = (unsigned char)(t >> 16);
168 *(md++) = (unsigned char)(t >> 8);
169 *(md++) = (unsigned char)(t);
172 /* ... as well as make sure md_len is not abused. */
180 int SHA384_Final(unsigned char *md, SHA512_CTX *c)
182 return SHA512_Final(md, c);
185 int SHA512_Update(SHA512_CTX *c, const void *_data, size_t len)
188 unsigned char *p = c->u.p;
189 const unsigned char *data = (const unsigned char *)_data;
194 l = (c->Nl + (((SHA_LONG64) len) << 3)) & U64(0xffffffffffffffff);
197 if (sizeof(len) >= 8)
198 c->Nh += (((SHA_LONG64) len) >> 61);
202 size_t n = sizeof(c->u) - c->num;
205 memcpy(p + c->num, data, len), c->num += (unsigned int)len;
208 memcpy(p + c->num, data, n), c->num = 0;
210 sha512_block_data_order(c, p, 1);
214 if (len >= sizeof(c->u)) {
215 # ifndef SHA512_BLOCK_CAN_MANAGE_UNALIGNED_DATA
216 if ((size_t)data % sizeof(c->u.d[0]) != 0)
217 while (len >= sizeof(c->u))
218 memcpy(p, data, sizeof(c->u)),
219 sha512_block_data_order(c, p, 1),
220 len -= sizeof(c->u), data += sizeof(c->u);
223 sha512_block_data_order(c, data, len / sizeof(c->u)),
224 data += len, len %= sizeof(c->u), data -= len;
228 memcpy(p, data, len), c->num = (int)len;
233 int SHA384_Update(SHA512_CTX *c, const void *data, size_t len)
235 return SHA512_Update(c, data, len);
238 void SHA512_Transform(SHA512_CTX *c, const unsigned char *data)
240 # ifndef SHA512_BLOCK_CAN_MANAGE_UNALIGNED_DATA
241 if ((size_t)data % sizeof(c->u.d[0]) != 0)
242 memcpy(c->u.p, data, sizeof(c->u.p)), data = c->u.p;
244 sha512_block_data_order(c, data, 1);
247 unsigned char *SHA384(const unsigned char *d, size_t n, unsigned char *md)
250 static unsigned char m[SHA384_DIGEST_LENGTH];
255 SHA512_Update(&c, d, n);
256 SHA512_Final(md, &c);
257 OPENSSL_cleanse(&c, sizeof(c));
261 unsigned char *SHA512(const unsigned char *d, size_t n, unsigned char *md)
264 static unsigned char m[SHA512_DIGEST_LENGTH];
269 SHA512_Update(&c, d, n);
270 SHA512_Final(md, &c);
271 OPENSSL_cleanse(&c, sizeof(c));
276 static const SHA_LONG64 K512[80] = {
277 U64(0x428a2f98d728ae22), U64(0x7137449123ef65cd),
278 U64(0xb5c0fbcfec4d3b2f), U64(0xe9b5dba58189dbbc),
279 U64(0x3956c25bf348b538), U64(0x59f111f1b605d019),
280 U64(0x923f82a4af194f9b), U64(0xab1c5ed5da6d8118),
281 U64(0xd807aa98a3030242), U64(0x12835b0145706fbe),
282 U64(0x243185be4ee4b28c), U64(0x550c7dc3d5ffb4e2),
283 U64(0x72be5d74f27b896f), U64(0x80deb1fe3b1696b1),
284 U64(0x9bdc06a725c71235), U64(0xc19bf174cf692694),
285 U64(0xe49b69c19ef14ad2), U64(0xefbe4786384f25e3),
286 U64(0x0fc19dc68b8cd5b5), U64(0x240ca1cc77ac9c65),
287 U64(0x2de92c6f592b0275), U64(0x4a7484aa6ea6e483),
288 U64(0x5cb0a9dcbd41fbd4), U64(0x76f988da831153b5),
289 U64(0x983e5152ee66dfab), U64(0xa831c66d2db43210),
290 U64(0xb00327c898fb213f), U64(0xbf597fc7beef0ee4),
291 U64(0xc6e00bf33da88fc2), U64(0xd5a79147930aa725),
292 U64(0x06ca6351e003826f), U64(0x142929670a0e6e70),
293 U64(0x27b70a8546d22ffc), U64(0x2e1b21385c26c926),
294 U64(0x4d2c6dfc5ac42aed), U64(0x53380d139d95b3df),
295 U64(0x650a73548baf63de), U64(0x766a0abb3c77b2a8),
296 U64(0x81c2c92e47edaee6), U64(0x92722c851482353b),
297 U64(0xa2bfe8a14cf10364), U64(0xa81a664bbc423001),
298 U64(0xc24b8b70d0f89791), U64(0xc76c51a30654be30),
299 U64(0xd192e819d6ef5218), U64(0xd69906245565a910),
300 U64(0xf40e35855771202a), U64(0x106aa07032bbd1b8),
301 U64(0x19a4c116b8d2d0c8), U64(0x1e376c085141ab53),
302 U64(0x2748774cdf8eeb99), U64(0x34b0bcb5e19b48a8),
303 U64(0x391c0cb3c5c95a63), U64(0x4ed8aa4ae3418acb),
304 U64(0x5b9cca4f7763e373), U64(0x682e6ff3d6b2b8a3),
305 U64(0x748f82ee5defb2fc), U64(0x78a5636f43172f60),
306 U64(0x84c87814a1f0ab72), U64(0x8cc702081a6439ec),
307 U64(0x90befffa23631e28), U64(0xa4506cebde82bde9),
308 U64(0xbef9a3f7b2c67915), U64(0xc67178f2e372532b),
309 U64(0xca273eceea26619c), U64(0xd186b8c721c0c207),
310 U64(0xeada7dd6cde0eb1e), U64(0xf57d4f7fee6ed178),
311 U64(0x06f067aa72176fba), U64(0x0a637dc5a2c898a6),
312 U64(0x113f9804bef90dae), U64(0x1b710b35131c471b),
313 U64(0x28db77f523047d84), U64(0x32caab7b40c72493),
314 U64(0x3c9ebe0a15c9bebc), U64(0x431d67c49c100d4c),
315 U64(0x4cc5d4becb3e42b6), U64(0x597f299cfc657e2a),
316 U64(0x5fcb6fab3ad6faec), U64(0x6c44198c4a475817)
320 # if defined(__GNUC__) && __GNUC__>=2 && !defined(OPENSSL_NO_ASM) && !defined(OPENSSL_NO_INLINE_ASM)
321 # if defined(__x86_64) || defined(__x86_64__)
322 # define ROTR(a,n) ({ SHA_LONG64 ret; \
327 # if !defined(B_ENDIAN)
328 # define PULL64(x) ({ SHA_LONG64 ret=*((const SHA_LONG64 *)(&(x))); \
333 # elif (defined(__i386) || defined(__i386__)) && !defined(B_ENDIAN)
334 # if defined(I386_ONLY)
335 # define PULL64(x) ({ const unsigned int *p=(const unsigned int *)(&(x));\
336 unsigned int hi=p[0],lo=p[1]; \
337 asm("xchgb %%ah,%%al;xchgb %%dh,%%dl;"\
338 "roll $16,%%eax; roll $16,%%edx; "\
339 "xchgb %%ah,%%al;xchgb %%dh,%%dl;" \
340 : "=a"(lo),"=d"(hi) \
341 : "0"(lo),"1"(hi) : "cc"); \
342 ((SHA_LONG64)hi)<<32|lo; })
344 # define PULL64(x) ({ const unsigned int *p=(const unsigned int *)(&(x));\
345 unsigned int hi=p[0],lo=p[1]; \
346 asm ("bswapl %0; bswapl %1;" \
347 : "=r"(lo),"=r"(hi) \
348 : "0"(lo),"1"(hi)); \
349 ((SHA_LONG64)hi)<<32|lo; })
351 # elif (defined(_ARCH_PPC) && defined(__64BIT__)) || defined(_ARCH_PPC64)
352 # define ROTR(a,n) ({ SHA_LONG64 ret; \
353 asm ("rotrdi %0,%1,%2" \
355 : "r"(a),"K"(n)); ret; })
357 # elif defined(_MSC_VER)
358 # if defined(_WIN64) /* applies to both IA-64 and AMD64 */
359 # pragma intrinsic(_rotr64)
360 # define ROTR(a,n) _rotr64((a),n)
362 # if defined(_M_IX86) && !defined(OPENSSL_NO_ASM) && !defined(OPENSSL_NO_INLINE_ASM)
363 # if defined(I386_ONLY)
364 static SHA_LONG64 __fastcall __pull64be(const void *x)
366 _asm mov edx,[ecx + 0]
367 _asm mov eax,[ecx + 4]
370 _asm rol edx, 16 _asm rol eax, 16 _asm xchg dh, dl _asm xchg ah, al}
372 static SHA_LONG64 __fastcall __pull64be(const void *x)
374 _asm mov edx,[ecx + 0]
375 _asm mov eax,[ecx + 4]
376 _asm bswap edx _asm bswap eax}
378 # define PULL64(x) __pull64be(&(x))
380 # pragma inline_depth(0)
386 # define B(x,j) (((SHA_LONG64)(*(((const unsigned char *)(&x))+j)))<<((7-j)*8))
387 # define PULL64(x) (B(x,0)|B(x,1)|B(x,2)|B(x,3)|B(x,4)|B(x,5)|B(x,6)|B(x,7))
390 # define ROTR(x,s) (((x)>>s) | (x)<<(64-s))
392 # define Sigma0(x) (ROTR((x),28) ^ ROTR((x),34) ^ ROTR((x),39))
393 # define Sigma1(x) (ROTR((x),14) ^ ROTR((x),18) ^ ROTR((x),41))
394 # define sigma0(x) (ROTR((x),1) ^ ROTR((x),8) ^ ((x)>>7))
395 # define sigma1(x) (ROTR((x),19) ^ ROTR((x),61) ^ ((x)>>6))
396 # define Ch(x,y,z) (((x) & (y)) ^ ((~(x)) & (z)))
397 # define Maj(x,y,z) (((x) & (y)) ^ ((x) & (z)) ^ ((y) & (z)))
398 # if defined(__i386) || defined(__i386__) || defined(_M_IX86)
400 * This code should give better results on 32-bit CPU with less than
401 * ~24 registers, both size and performance wise...
402 */ static void sha512_block_data_order(SHA512_CTX *ctx, const void *in,
405 const SHA_LONG64 *W = in;
407 SHA_LONG64 X[9 + 80], *F;
422 for (i = 0; i < 16; i++, F--) {
431 T += F[7] + Sigma1(E) + Ch(E, F[5], F[6]) + K512[i];
433 A = T + Sigma0(A) + Maj(A, F[1], F[2]);
436 for (; i < 80; i++, F--) {
437 T = sigma0(F[8 + 16 - 1]);
438 T += sigma1(F[8 + 16 - 14]);
439 T += F[8 + 16] + F[8 + 16 - 9];
444 T += F[7] + Sigma1(E) + Ch(E, F[5], F[6]) + K512[i];
446 A = T + Sigma0(A) + Maj(A, F[1], F[2]);
462 # elif defined(OPENSSL_SMALL_FOOTPRINT)
463 static void sha512_block_data_order(SHA512_CTX *ctx, const void *in,
466 const SHA_LONG64 *W = in;
467 SHA_LONG64 a, b, c, d, e, f, g, h, s0, s1, T1, T2;
482 for (i = 0; i < 16; i++) {
486 T1 = X[i] = PULL64(W[i]);
488 T1 += h + Sigma1(e) + Ch(e, f, g) + K512[i];
489 T2 = Sigma0(a) + Maj(a, b, c);
500 for (; i < 80; i++) {
501 s0 = X[(i + 1) & 0x0f];
503 s1 = X[(i + 14) & 0x0f];
506 T1 = X[i & 0xf] += s0 + s1 + X[(i + 9) & 0xf];
507 T1 += h + Sigma1(e) + Ch(e, f, g) + K512[i];
508 T2 = Sigma0(a) + Maj(a, b, c);
533 # define ROUND_00_15(i,a,b,c,d,e,f,g,h) do { \
534 T1 += h + Sigma1(e) + Ch(e,f,g) + K512[i]; \
535 h = Sigma0(a) + Maj(a,b,c); \
536 d += T1; h += T1; } while (0)
537 # define ROUND_16_80(i,j,a,b,c,d,e,f,g,h,X) do { \
538 s0 = X[(j+1)&0x0f]; s0 = sigma0(s0); \
539 s1 = X[(j+14)&0x0f]; s1 = sigma1(s1); \
540 T1 = X[(j)&0x0f] += s0 + s1 + X[(j+9)&0x0f]; \
541 ROUND_00_15(i+j,a,b,c,d,e,f,g,h); } while (0)
542 static void sha512_block_data_order(SHA512_CTX *ctx, const void *in,
545 const SHA_LONG64 *W = in;
546 SHA_LONG64 a, b, c, d, e, f, g, h, s0, s1, T1;
563 ROUND_00_15(0, a, b, c, d, e, f, g, h);
565 ROUND_00_15(1, h, a, b, c, d, e, f, g);
567 ROUND_00_15(2, g, h, a, b, c, d, e, f);
569 ROUND_00_15(3, f, g, h, a, b, c, d, e);
571 ROUND_00_15(4, e, f, g, h, a, b, c, d);
573 ROUND_00_15(5, d, e, f, g, h, a, b, c);
575 ROUND_00_15(6, c, d, e, f, g, h, a, b);
577 ROUND_00_15(7, b, c, d, e, f, g, h, a);
579 ROUND_00_15(8, a, b, c, d, e, f, g, h);
581 ROUND_00_15(9, h, a, b, c, d, e, f, g);
583 ROUND_00_15(10, g, h, a, b, c, d, e, f);
585 ROUND_00_15(11, f, g, h, a, b, c, d, e);
587 ROUND_00_15(12, e, f, g, h, a, b, c, d);
589 ROUND_00_15(13, d, e, f, g, h, a, b, c);
591 ROUND_00_15(14, c, d, e, f, g, h, a, b);
593 ROUND_00_15(15, b, c, d, e, f, g, h, a);
595 T1 = X[0] = PULL64(W[0]);
596 ROUND_00_15(0, a, b, c, d, e, f, g, h);
597 T1 = X[1] = PULL64(W[1]);
598 ROUND_00_15(1, h, a, b, c, d, e, f, g);
599 T1 = X[2] = PULL64(W[2]);
600 ROUND_00_15(2, g, h, a, b, c, d, e, f);
601 T1 = X[3] = PULL64(W[3]);
602 ROUND_00_15(3, f, g, h, a, b, c, d, e);
603 T1 = X[4] = PULL64(W[4]);
604 ROUND_00_15(4, e, f, g, h, a, b, c, d);
605 T1 = X[5] = PULL64(W[5]);
606 ROUND_00_15(5, d, e, f, g, h, a, b, c);
607 T1 = X[6] = PULL64(W[6]);
608 ROUND_00_15(6, c, d, e, f, g, h, a, b);
609 T1 = X[7] = PULL64(W[7]);
610 ROUND_00_15(7, b, c, d, e, f, g, h, a);
611 T1 = X[8] = PULL64(W[8]);
612 ROUND_00_15(8, a, b, c, d, e, f, g, h);
613 T1 = X[9] = PULL64(W[9]);
614 ROUND_00_15(9, h, a, b, c, d, e, f, g);
615 T1 = X[10] = PULL64(W[10]);
616 ROUND_00_15(10, g, h, a, b, c, d, e, f);
617 T1 = X[11] = PULL64(W[11]);
618 ROUND_00_15(11, f, g, h, a, b, c, d, e);
619 T1 = X[12] = PULL64(W[12]);
620 ROUND_00_15(12, e, f, g, h, a, b, c, d);
621 T1 = X[13] = PULL64(W[13]);
622 ROUND_00_15(13, d, e, f, g, h, a, b, c);
623 T1 = X[14] = PULL64(W[14]);
624 ROUND_00_15(14, c, d, e, f, g, h, a, b);
625 T1 = X[15] = PULL64(W[15]);
626 ROUND_00_15(15, b, c, d, e, f, g, h, a);
629 for (i = 16; i < 80; i += 16) {
630 ROUND_16_80(i, 0, a, b, c, d, e, f, g, h, X);
631 ROUND_16_80(i, 1, h, a, b, c, d, e, f, g, X);
632 ROUND_16_80(i, 2, g, h, a, b, c, d, e, f, X);
633 ROUND_16_80(i, 3, f, g, h, a, b, c, d, e, X);
634 ROUND_16_80(i, 4, e, f, g, h, a, b, c, d, X);
635 ROUND_16_80(i, 5, d, e, f, g, h, a, b, c, X);
636 ROUND_16_80(i, 6, c, d, e, f, g, h, a, b, X);
637 ROUND_16_80(i, 7, b, c, d, e, f, g, h, a, X);
638 ROUND_16_80(i, 8, a, b, c, d, e, f, g, h, X);
639 ROUND_16_80(i, 9, h, a, b, c, d, e, f, g, X);
640 ROUND_16_80(i, 10, g, h, a, b, c, d, e, f, X);
641 ROUND_16_80(i, 11, f, g, h, a, b, c, d, e, X);
642 ROUND_16_80(i, 12, e, f, g, h, a, b, c, d, X);
643 ROUND_16_80(i, 13, d, e, f, g, h, a, b, c, X);
644 ROUND_16_80(i, 14, c, d, e, f, g, h, a, b, X);
645 ROUND_16_80(i, 15, b, c, d, e, f, g, h, a, X);
663 # endif /* SHA512_ASM */
665 #else /* !OPENSSL_NO_SHA512 */
667 # if defined(PEDANTIC) || defined(__DECC) || defined(OPENSSL_SYS_MACOSX)
668 static void *dummy = &dummy;
671 #endif /* !OPENSSL_NO_SHA512 */