sha512.c 23 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190191192193194195196197198199200201202203204205206207208209210211212213214215216217218219220221222223224225226227228229230231232233234235236237238239240241242243244245246247248249250251252253254255256257258259260261262263264265266267268269270271272273274275276277278279280281282283284285286287288289290291292293294295296297298299300301302303304305306307308309310311312313314315316317318319320321322323324325326327328329330331332333334335336337338339340341342343344345346347348349350351352353354355356357358359360361362363364365366367368369370371372373374375376377378379380381382383384385386387388389390391392393394395396397398399400401402403404405406407408409410411412413414415416417418419420421422423424425426427428429430431432433434435436437438439440441442443444445446447448449450451452453454455456457458459460461462463464465466467468469470471472473474475476477478479480481482483484485486487488489490491492493494495496497498499500501502503504505506507508509510511512513514515516517518519520521522523524525526527528529530531532533534535536537538539540541542543544545546547548549550551552553554555556557558559560561562563564565566567568569570571572573574575576577578579580581582583584585586587588589590591592593594595596597598599600601602603604605606607608609610611612613614615616617618619620621622623624625626627628629630631632633634635636637638639640641642643644645646647648649650651652653654655656657658659660661662663664665666667668669670671672673674675676677678679680681682683684685
  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. * ====================================================================
  6. */
  7. #include <openssl/opensslconf.h>
  8. #if !defined(OPENSSL_NO_SHA) && !defined(OPENSSL_NO_SHA512)
  9. /*-
  10. * IMPLEMENTATION NOTES.
  11. *
  12. * As you might have noticed 32-bit hash algorithms:
  13. *
  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
  17. * stream byte order;
  18. * - share common byte-order neutral collector and padding function
  19. * implementations, ../md32_common.h;
  20. *
  21. * Neither of the above applies to this SHA-512 implementations. Reasons
  22. * [in reverse order] are:
  23. *
  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;
  33. *
  34. * Caveat lector.
  35. *
  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
  41. * 16-bit platforms.
  42. * <appro@fy.chalmers.se>
  43. */
  44. # include <stdlib.h>
  45. # include <string.h>
  46. # include <openssl/crypto.h>
  47. # include <openssl/sha.h>
  48. # include <openssl/opensslv.h>
  49. # include "cryptlib.h"
  50. const char SHA512_version[] = "SHA-512" OPENSSL_VERSION_PTEXT;
  51. # if defined(__i386) || defined(__i386__) || defined(_M_IX86) || \
  52. defined(__x86_64) || defined(_M_AMD64) || defined(_M_X64) || \
  53. defined(__s390__) || defined(__s390x__) || \
  54. defined(__aarch64__) || \
  55. defined(SHA512_ASM)
  56. # define SHA512_BLOCK_CAN_MANAGE_UNALIGNED_DATA
  57. # endif
  58. fips_md_init_ctx(SHA384, SHA512)
  59. {
  60. c->h[0] = U64(0xcbbb9d5dc1059ed8);
  61. c->h[1] = U64(0x629a292a367cd507);
  62. c->h[2] = U64(0x9159015a3070dd17);
  63. c->h[3] = U64(0x152fecd8f70e5939);
  64. c->h[4] = U64(0x67332667ffc00b31);
  65. c->h[5] = U64(0x8eb44a8768581511);
  66. c->h[6] = U64(0xdb0c2e0d64f98fa7);
  67. c->h[7] = U64(0x47b5481dbefa4fa4);
  68. c->Nl = 0;
  69. c->Nh = 0;
  70. c->num = 0;
  71. c->md_len = SHA384_DIGEST_LENGTH;
  72. return 1;
  73. }
  74. fips_md_init(SHA512)
  75. {
  76. c->h[0] = U64(0x6a09e667f3bcc908);
  77. c->h[1] = U64(0xbb67ae8584caa73b);
  78. c->h[2] = U64(0x3c6ef372fe94f82b);
  79. c->h[3] = U64(0xa54ff53a5f1d36f1);
  80. c->h[4] = U64(0x510e527fade682d1);
  81. c->h[5] = U64(0x9b05688c2b3e6c1f);
  82. c->h[6] = U64(0x1f83d9abfb41bd6b);
  83. c->h[7] = U64(0x5be0cd19137e2179);
  84. c->Nl = 0;
  85. c->Nh = 0;
  86. c->num = 0;
  87. c->md_len = SHA512_DIGEST_LENGTH;
  88. return 1;
  89. }
  90. # ifndef SHA512_ASM
  91. static
  92. # endif
  93. void sha512_block_data_order(SHA512_CTX *ctx, const void *in, size_t num);
  94. int SHA512_Final(unsigned char *md, SHA512_CTX *c)
  95. {
  96. unsigned char *p = (unsigned char *)c->u.p;
  97. size_t n = c->num;
  98. p[n] = 0x80; /* There always is a room for one */
  99. n++;
  100. if (n > (sizeof(c->u) - 16))
  101. memset(p + n, 0, sizeof(c->u) - n), n = 0,
  102. sha512_block_data_order(c, p, 1);
  103. memset(p + n, 0, sizeof(c->u) - 16 - n);
  104. # ifdef B_ENDIAN
  105. c->u.d[SHA_LBLOCK - 2] = c->Nh;
  106. c->u.d[SHA_LBLOCK - 1] = c->Nl;
  107. # else
  108. p[sizeof(c->u) - 1] = (unsigned char)(c->Nl);
  109. p[sizeof(c->u) - 2] = (unsigned char)(c->Nl >> 8);
  110. p[sizeof(c->u) - 3] = (unsigned char)(c->Nl >> 16);
  111. p[sizeof(c->u) - 4] = (unsigned char)(c->Nl >> 24);
  112. p[sizeof(c->u) - 5] = (unsigned char)(c->Nl >> 32);
  113. p[sizeof(c->u) - 6] = (unsigned char)(c->Nl >> 40);
  114. p[sizeof(c->u) - 7] = (unsigned char)(c->Nl >> 48);
  115. p[sizeof(c->u) - 8] = (unsigned char)(c->Nl >> 56);
  116. p[sizeof(c->u) - 9] = (unsigned char)(c->Nh);
  117. p[sizeof(c->u) - 10] = (unsigned char)(c->Nh >> 8);
  118. p[sizeof(c->u) - 11] = (unsigned char)(c->Nh >> 16);
  119. p[sizeof(c->u) - 12] = (unsigned char)(c->Nh >> 24);
  120. p[sizeof(c->u) - 13] = (unsigned char)(c->Nh >> 32);
  121. p[sizeof(c->u) - 14] = (unsigned char)(c->Nh >> 40);
  122. p[sizeof(c->u) - 15] = (unsigned char)(c->Nh >> 48);
  123. p[sizeof(c->u) - 16] = (unsigned char)(c->Nh >> 56);
  124. # endif
  125. sha512_block_data_order(c, p, 1);
  126. if (md == 0)
  127. return 0;
  128. switch (c->md_len) {
  129. /* Let compiler decide if it's appropriate to unroll... */
  130. case SHA384_DIGEST_LENGTH:
  131. for (n = 0; n < SHA384_DIGEST_LENGTH / 8; n++) {
  132. SHA_LONG64 t = c->h[n];
  133. *(md++) = (unsigned char)(t >> 56);
  134. *(md++) = (unsigned char)(t >> 48);
  135. *(md++) = (unsigned char)(t >> 40);
  136. *(md++) = (unsigned char)(t >> 32);
  137. *(md++) = (unsigned char)(t >> 24);
  138. *(md++) = (unsigned char)(t >> 16);
  139. *(md++) = (unsigned char)(t >> 8);
  140. *(md++) = (unsigned char)(t);
  141. }
  142. break;
  143. case SHA512_DIGEST_LENGTH:
  144. for (n = 0; n < SHA512_DIGEST_LENGTH / 8; n++) {
  145. SHA_LONG64 t = c->h[n];
  146. *(md++) = (unsigned char)(t >> 56);
  147. *(md++) = (unsigned char)(t >> 48);
  148. *(md++) = (unsigned char)(t >> 40);
  149. *(md++) = (unsigned char)(t >> 32);
  150. *(md++) = (unsigned char)(t >> 24);
  151. *(md++) = (unsigned char)(t >> 16);
  152. *(md++) = (unsigned char)(t >> 8);
  153. *(md++) = (unsigned char)(t);
  154. }
  155. break;
  156. /* ... as well as make sure md_len is not abused. */
  157. default:
  158. return 0;
  159. }
  160. return 1;
  161. }
  162. int SHA384_Final(unsigned char *md, SHA512_CTX *c)
  163. {
  164. return SHA512_Final(md, c);
  165. }
  166. int SHA512_Update(SHA512_CTX *c, const void *_data, size_t len)
  167. {
  168. SHA_LONG64 l;
  169. unsigned char *p = c->u.p;
  170. const unsigned char *data = (const unsigned char *)_data;
  171. if (len == 0)
  172. return 1;
  173. l = (c->Nl + (((SHA_LONG64) len) << 3)) & U64(0xffffffffffffffff);
  174. if (l < c->Nl)
  175. c->Nh++;
  176. if (sizeof(len) >= 8)
  177. c->Nh += (((SHA_LONG64) len) >> 61);
  178. c->Nl = l;
  179. if (c->num != 0) {
  180. size_t n = sizeof(c->u) - c->num;
  181. if (len < n) {
  182. memcpy(p + c->num, data, len), c->num += (unsigned int)len;
  183. return 1;
  184. } else {
  185. memcpy(p + c->num, data, n), c->num = 0;
  186. len -= n, data += n;
  187. sha512_block_data_order(c, p, 1);
  188. }
  189. }
  190. if (len >= sizeof(c->u)) {
  191. # ifndef SHA512_BLOCK_CAN_MANAGE_UNALIGNED_DATA
  192. if ((size_t)data % sizeof(c->u.d[0]) != 0)
  193. while (len >= sizeof(c->u))
  194. memcpy(p, data, sizeof(c->u)),
  195. sha512_block_data_order(c, p, 1),
  196. len -= sizeof(c->u), data += sizeof(c->u);
  197. else
  198. # endif
  199. sha512_block_data_order(c, data, len / sizeof(c->u)),
  200. data += len, len %= sizeof(c->u), data -= len;
  201. }
  202. if (len != 0)
  203. memcpy(p, data, len), c->num = (int)len;
  204. return 1;
  205. }
  206. int SHA384_Update(SHA512_CTX *c, const void *data, size_t len)
  207. {
  208. return SHA512_Update(c, data, len);
  209. }
  210. void SHA512_Transform(SHA512_CTX *c, const unsigned char *data)
  211. {
  212. # ifndef SHA512_BLOCK_CAN_MANAGE_UNALIGNED_DATA
  213. if ((size_t)data % sizeof(c->u.d[0]) != 0)
  214. memcpy(c->u.p, data, sizeof(c->u.p)), data = c->u.p;
  215. # endif
  216. sha512_block_data_order(c, data, 1);
  217. }
  218. unsigned char *SHA384(const unsigned char *d, size_t n, unsigned char *md)
  219. {
  220. SHA512_CTX c;
  221. static unsigned char m[SHA384_DIGEST_LENGTH];
  222. if (md == NULL)
  223. md = m;
  224. SHA384_Init(&c);
  225. SHA512_Update(&c, d, n);
  226. SHA512_Final(md, &c);
  227. OPENSSL_cleanse(&c, sizeof(c));
  228. return (md);
  229. }
  230. unsigned char *SHA512(const unsigned char *d, size_t n, unsigned char *md)
  231. {
  232. SHA512_CTX c;
  233. static unsigned char m[SHA512_DIGEST_LENGTH];
  234. if (md == NULL)
  235. md = m;
  236. SHA512_Init(&c);
  237. SHA512_Update(&c, d, n);
  238. SHA512_Final(md, &c);
  239. OPENSSL_cleanse(&c, sizeof(c));
  240. return (md);
  241. }
  242. # ifndef SHA512_ASM
  243. static const SHA_LONG64 K512[80] = {
  244. U64(0x428a2f98d728ae22), U64(0x7137449123ef65cd),
  245. U64(0xb5c0fbcfec4d3b2f), U64(0xe9b5dba58189dbbc),
  246. U64(0x3956c25bf348b538), U64(0x59f111f1b605d019),
  247. U64(0x923f82a4af194f9b), U64(0xab1c5ed5da6d8118),
  248. U64(0xd807aa98a3030242), U64(0x12835b0145706fbe),
  249. U64(0x243185be4ee4b28c), U64(0x550c7dc3d5ffb4e2),
  250. U64(0x72be5d74f27b896f), U64(0x80deb1fe3b1696b1),
  251. U64(0x9bdc06a725c71235), U64(0xc19bf174cf692694),
  252. U64(0xe49b69c19ef14ad2), U64(0xefbe4786384f25e3),
  253. U64(0x0fc19dc68b8cd5b5), U64(0x240ca1cc77ac9c65),
  254. U64(0x2de92c6f592b0275), U64(0x4a7484aa6ea6e483),
  255. U64(0x5cb0a9dcbd41fbd4), U64(0x76f988da831153b5),
  256. U64(0x983e5152ee66dfab), U64(0xa831c66d2db43210),
  257. U64(0xb00327c898fb213f), U64(0xbf597fc7beef0ee4),
  258. U64(0xc6e00bf33da88fc2), U64(0xd5a79147930aa725),
  259. U64(0x06ca6351e003826f), U64(0x142929670a0e6e70),
  260. U64(0x27b70a8546d22ffc), U64(0x2e1b21385c26c926),
  261. U64(0x4d2c6dfc5ac42aed), U64(0x53380d139d95b3df),
  262. U64(0x650a73548baf63de), U64(0x766a0abb3c77b2a8),
  263. U64(0x81c2c92e47edaee6), U64(0x92722c851482353b),
  264. U64(0xa2bfe8a14cf10364), U64(0xa81a664bbc423001),
  265. U64(0xc24b8b70d0f89791), U64(0xc76c51a30654be30),
  266. U64(0xd192e819d6ef5218), U64(0xd69906245565a910),
  267. U64(0xf40e35855771202a), U64(0x106aa07032bbd1b8),
  268. U64(0x19a4c116b8d2d0c8), U64(0x1e376c085141ab53),
  269. U64(0x2748774cdf8eeb99), U64(0x34b0bcb5e19b48a8),
  270. U64(0x391c0cb3c5c95a63), U64(0x4ed8aa4ae3418acb),
  271. U64(0x5b9cca4f7763e373), U64(0x682e6ff3d6b2b8a3),
  272. U64(0x748f82ee5defb2fc), U64(0x78a5636f43172f60),
  273. U64(0x84c87814a1f0ab72), U64(0x8cc702081a6439ec),
  274. U64(0x90befffa23631e28), U64(0xa4506cebde82bde9),
  275. U64(0xbef9a3f7b2c67915), U64(0xc67178f2e372532b),
  276. U64(0xca273eceea26619c), U64(0xd186b8c721c0c207),
  277. U64(0xeada7dd6cde0eb1e), U64(0xf57d4f7fee6ed178),
  278. U64(0x06f067aa72176fba), U64(0x0a637dc5a2c898a6),
  279. U64(0x113f9804bef90dae), U64(0x1b710b35131c471b),
  280. U64(0x28db77f523047d84), U64(0x32caab7b40c72493),
  281. U64(0x3c9ebe0a15c9bebc), U64(0x431d67c49c100d4c),
  282. U64(0x4cc5d4becb3e42b6), U64(0x597f299cfc657e2a),
  283. U64(0x5fcb6fab3ad6faec), U64(0x6c44198c4a475817)
  284. };
  285. # ifndef PEDANTIC
  286. # if defined(__GNUC__) && __GNUC__>=2 && !defined(OPENSSL_NO_ASM) && !defined(OPENSSL_NO_INLINE_ASM)
  287. # if defined(__x86_64) || defined(__x86_64__)
  288. # define ROTR(a,n) ({ SHA_LONG64 ret; \
  289. asm ("rorq %1,%0" \
  290. : "=r"(ret) \
  291. : "J"(n),"0"(a) \
  292. : "cc"); ret; })
  293. # if !defined(B_ENDIAN)
  294. # define PULL64(x) ({ SHA_LONG64 ret=*((const SHA_LONG64 *)(&(x))); \
  295. asm ("bswapq %0" \
  296. : "=r"(ret) \
  297. : "0"(ret)); ret; })
  298. # endif
  299. # elif (defined(__i386) || defined(__i386__)) && !defined(B_ENDIAN)
  300. # if defined(I386_ONLY)
  301. # define PULL64(x) ({ const unsigned int *p=(const unsigned int *)(&(x));\
  302. unsigned int hi=p[0],lo=p[1]; \
  303. asm("xchgb %%ah,%%al;xchgb %%dh,%%dl;"\
  304. "roll $16,%%eax; roll $16,%%edx; "\
  305. "xchgb %%ah,%%al;xchgb %%dh,%%dl;" \
  306. : "=a"(lo),"=d"(hi) \
  307. : "0"(lo),"1"(hi) : "cc"); \
  308. ((SHA_LONG64)hi)<<32|lo; })
  309. # else
  310. # define PULL64(x) ({ const unsigned int *p=(const unsigned int *)(&(x));\
  311. unsigned int hi=p[0],lo=p[1]; \
  312. asm ("bswapl %0; bswapl %1;" \
  313. : "=r"(lo),"=r"(hi) \
  314. : "0"(lo),"1"(hi)); \
  315. ((SHA_LONG64)hi)<<32|lo; })
  316. # endif
  317. # elif (defined(_ARCH_PPC) && defined(__64BIT__)) || defined(_ARCH_PPC64)
  318. # define ROTR(a,n) ({ SHA_LONG64 ret; \
  319. asm ("rotrdi %0,%1,%2" \
  320. : "=r"(ret) \
  321. : "r"(a),"K"(n)); ret; })
  322. # elif defined(__aarch64__)
  323. # define ROTR(a,n) ({ SHA_LONG64 ret; \
  324. asm ("ror %0,%1,%2" \
  325. : "=r"(ret) \
  326. : "r"(a),"I"(n)); ret; })
  327. # if defined(__BYTE_ORDER__) && defined(__ORDER_LITTLE_ENDIAN__) && \
  328. __BYTE_ORDER__==__ORDER_LITTLE_ENDIAN__
  329. # define PULL64(x) ({ SHA_LONG64 ret; \
  330. asm ("rev %0,%1" \
  331. : "=r"(ret) \
  332. : "r"(*((const SHA_LONG64 *)(&(x))))); ret; })
  333. # endif
  334. # endif
  335. # elif defined(_MSC_VER)
  336. # if defined(_WIN64) /* applies to both IA-64 and AMD64 */
  337. # pragma intrinsic(_rotr64)
  338. # define ROTR(a,n) _rotr64((a),n)
  339. # endif
  340. # if defined(_M_IX86) && !defined(OPENSSL_NO_ASM) && !defined(OPENSSL_NO_INLINE_ASM)
  341. # if defined(I386_ONLY)
  342. static SHA_LONG64 __fastcall __pull64be(const void *x)
  343. {
  344. _asm mov edx,[ecx + 0]
  345. _asm mov eax,[ecx + 4]
  346. _asm xchg dh, dl
  347. _asm xchg ah, al
  348. _asm rol edx, 16 _asm rol eax, 16 _asm xchg dh, dl _asm xchg ah, al}
  349. # else
  350. static SHA_LONG64 __fastcall __pull64be(const void *x)
  351. {
  352. _asm mov edx,[ecx + 0]
  353. _asm mov eax,[ecx + 4]
  354. _asm bswap edx _asm bswap eax}
  355. # endif
  356. # define PULL64(x) __pull64be(&(x))
  357. # if _MSC_VER<=1200
  358. # pragma inline_depth(0)
  359. # endif
  360. # endif
  361. # endif
  362. # endif
  363. # ifndef PULL64
  364. # define B(x,j) (((SHA_LONG64)(*(((const unsigned char *)(&x))+j)))<<((7-j)*8))
  365. # 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))
  366. # endif
  367. # ifndef ROTR
  368. # define ROTR(x,s) (((x)>>s) | (x)<<(64-s))
  369. # endif
  370. # define Sigma0(x) (ROTR((x),28) ^ ROTR((x),34) ^ ROTR((x),39))
  371. # define Sigma1(x) (ROTR((x),14) ^ ROTR((x),18) ^ ROTR((x),41))
  372. # define sigma0(x) (ROTR((x),1) ^ ROTR((x),8) ^ ((x)>>7))
  373. # define sigma1(x) (ROTR((x),19) ^ ROTR((x),61) ^ ((x)>>6))
  374. # define Ch(x,y,z) (((x) & (y)) ^ ((~(x)) & (z)))
  375. # define Maj(x,y,z) (((x) & (y)) ^ ((x) & (z)) ^ ((y) & (z)))
  376. # if defined(__i386) || defined(__i386__) || defined(_M_IX86)
  377. /*
  378. * This code should give better results on 32-bit CPU with less than
  379. * ~24 registers, both size and performance wise...
  380. */ static void sha512_block_data_order(SHA512_CTX *ctx, const void *in,
  381. size_t num)
  382. {
  383. const SHA_LONG64 *W = in;
  384. SHA_LONG64 A, E, T;
  385. SHA_LONG64 X[9 + 80], *F;
  386. int i;
  387. while (num--) {
  388. F = X + 80;
  389. A = ctx->h[0];
  390. F[1] = ctx->h[1];
  391. F[2] = ctx->h[2];
  392. F[3] = ctx->h[3];
  393. E = ctx->h[4];
  394. F[5] = ctx->h[5];
  395. F[6] = ctx->h[6];
  396. F[7] = ctx->h[7];
  397. for (i = 0; i < 16; i++, F--) {
  398. # ifdef B_ENDIAN
  399. T = W[i];
  400. # else
  401. T = PULL64(W[i]);
  402. # endif
  403. F[0] = A;
  404. F[4] = E;
  405. F[8] = T;
  406. T += F[7] + Sigma1(E) + Ch(E, F[5], F[6]) + K512[i];
  407. E = F[3] + T;
  408. A = T + Sigma0(A) + Maj(A, F[1], F[2]);
  409. }
  410. for (; i < 80; i++, F--) {
  411. T = sigma0(F[8 + 16 - 1]);
  412. T += sigma1(F[8 + 16 - 14]);
  413. T += F[8 + 16] + F[8 + 16 - 9];
  414. F[0] = A;
  415. F[4] = E;
  416. F[8] = T;
  417. T += F[7] + Sigma1(E) + Ch(E, F[5], F[6]) + K512[i];
  418. E = F[3] + T;
  419. A = T + Sigma0(A) + Maj(A, F[1], F[2]);
  420. }
  421. ctx->h[0] += A;
  422. ctx->h[1] += F[1];
  423. ctx->h[2] += F[2];
  424. ctx->h[3] += F[3];
  425. ctx->h[4] += E;
  426. ctx->h[5] += F[5];
  427. ctx->h[6] += F[6];
  428. ctx->h[7] += F[7];
  429. W += SHA_LBLOCK;
  430. }
  431. }
  432. # elif defined(OPENSSL_SMALL_FOOTPRINT)
  433. static void sha512_block_data_order(SHA512_CTX *ctx, const void *in,
  434. size_t num)
  435. {
  436. const SHA_LONG64 *W = in;
  437. SHA_LONG64 a, b, c, d, e, f, g, h, s0, s1, T1, T2;
  438. SHA_LONG64 X[16];
  439. int i;
  440. while (num--) {
  441. a = ctx->h[0];
  442. b = ctx->h[1];
  443. c = ctx->h[2];
  444. d = ctx->h[3];
  445. e = ctx->h[4];
  446. f = ctx->h[5];
  447. g = ctx->h[6];
  448. h = ctx->h[7];
  449. for (i = 0; i < 16; i++) {
  450. # ifdef B_ENDIAN
  451. T1 = X[i] = W[i];
  452. # else
  453. T1 = X[i] = PULL64(W[i]);
  454. # endif
  455. T1 += h + Sigma1(e) + Ch(e, f, g) + K512[i];
  456. T2 = Sigma0(a) + Maj(a, b, c);
  457. h = g;
  458. g = f;
  459. f = e;
  460. e = d + T1;
  461. d = c;
  462. c = b;
  463. b = a;
  464. a = T1 + T2;
  465. }
  466. for (; i < 80; i++) {
  467. s0 = X[(i + 1) & 0x0f];
  468. s0 = sigma0(s0);
  469. s1 = X[(i + 14) & 0x0f];
  470. s1 = sigma1(s1);
  471. T1 = X[i & 0xf] += s0 + s1 + X[(i + 9) & 0xf];
  472. T1 += h + Sigma1(e) + Ch(e, f, g) + K512[i];
  473. T2 = Sigma0(a) + Maj(a, b, c);
  474. h = g;
  475. g = f;
  476. f = e;
  477. e = d + T1;
  478. d = c;
  479. c = b;
  480. b = a;
  481. a = T1 + T2;
  482. }
  483. ctx->h[0] += a;
  484. ctx->h[1] += b;
  485. ctx->h[2] += c;
  486. ctx->h[3] += d;
  487. ctx->h[4] += e;
  488. ctx->h[5] += f;
  489. ctx->h[6] += g;
  490. ctx->h[7] += h;
  491. W += SHA_LBLOCK;
  492. }
  493. }
  494. # else
  495. # define ROUND_00_15(i,a,b,c,d,e,f,g,h) do { \
  496. T1 += h + Sigma1(e) + Ch(e,f,g) + K512[i]; \
  497. h = Sigma0(a) + Maj(a,b,c); \
  498. d += T1; h += T1; } while (0)
  499. # define ROUND_16_80(i,j,a,b,c,d,e,f,g,h,X) do { \
  500. s0 = X[(j+1)&0x0f]; s0 = sigma0(s0); \
  501. s1 = X[(j+14)&0x0f]; s1 = sigma1(s1); \
  502. T1 = X[(j)&0x0f] += s0 + s1 + X[(j+9)&0x0f]; \
  503. ROUND_00_15(i+j,a,b,c,d,e,f,g,h); } while (0)
  504. static void sha512_block_data_order(SHA512_CTX *ctx, const void *in,
  505. size_t num)
  506. {
  507. const SHA_LONG64 *W = in;
  508. SHA_LONG64 a, b, c, d, e, f, g, h, s0, s1, T1;
  509. SHA_LONG64 X[16];
  510. int i;
  511. while (num--) {
  512. a = ctx->h[0];
  513. b = ctx->h[1];
  514. c = ctx->h[2];
  515. d = ctx->h[3];
  516. e = ctx->h[4];
  517. f = ctx->h[5];
  518. g = ctx->h[6];
  519. h = ctx->h[7];
  520. # ifdef B_ENDIAN
  521. T1 = X[0] = W[0];
  522. ROUND_00_15(0, a, b, c, d, e, f, g, h);
  523. T1 = X[1] = W[1];
  524. ROUND_00_15(1, h, a, b, c, d, e, f, g);
  525. T1 = X[2] = W[2];
  526. ROUND_00_15(2, g, h, a, b, c, d, e, f);
  527. T1 = X[3] = W[3];
  528. ROUND_00_15(3, f, g, h, a, b, c, d, e);
  529. T1 = X[4] = W[4];
  530. ROUND_00_15(4, e, f, g, h, a, b, c, d);
  531. T1 = X[5] = W[5];
  532. ROUND_00_15(5, d, e, f, g, h, a, b, c);
  533. T1 = X[6] = W[6];
  534. ROUND_00_15(6, c, d, e, f, g, h, a, b);
  535. T1 = X[7] = W[7];
  536. ROUND_00_15(7, b, c, d, e, f, g, h, a);
  537. T1 = X[8] = W[8];
  538. ROUND_00_15(8, a, b, c, d, e, f, g, h);
  539. T1 = X[9] = W[9];
  540. ROUND_00_15(9, h, a, b, c, d, e, f, g);
  541. T1 = X[10] = W[10];
  542. ROUND_00_15(10, g, h, a, b, c, d, e, f);
  543. T1 = X[11] = W[11];
  544. ROUND_00_15(11, f, g, h, a, b, c, d, e);
  545. T1 = X[12] = W[12];
  546. ROUND_00_15(12, e, f, g, h, a, b, c, d);
  547. T1 = X[13] = W[13];
  548. ROUND_00_15(13, d, e, f, g, h, a, b, c);
  549. T1 = X[14] = W[14];
  550. ROUND_00_15(14, c, d, e, f, g, h, a, b);
  551. T1 = X[15] = W[15];
  552. ROUND_00_15(15, b, c, d, e, f, g, h, a);
  553. # else
  554. T1 = X[0] = PULL64(W[0]);
  555. ROUND_00_15(0, a, b, c, d, e, f, g, h);
  556. T1 = X[1] = PULL64(W[1]);
  557. ROUND_00_15(1, h, a, b, c, d, e, f, g);
  558. T1 = X[2] = PULL64(W[2]);
  559. ROUND_00_15(2, g, h, a, b, c, d, e, f);
  560. T1 = X[3] = PULL64(W[3]);
  561. ROUND_00_15(3, f, g, h, a, b, c, d, e);
  562. T1 = X[4] = PULL64(W[4]);
  563. ROUND_00_15(4, e, f, g, h, a, b, c, d);
  564. T1 = X[5] = PULL64(W[5]);
  565. ROUND_00_15(5, d, e, f, g, h, a, b, c);
  566. T1 = X[6] = PULL64(W[6]);
  567. ROUND_00_15(6, c, d, e, f, g, h, a, b);
  568. T1 = X[7] = PULL64(W[7]);
  569. ROUND_00_15(7, b, c, d, e, f, g, h, a);
  570. T1 = X[8] = PULL64(W[8]);
  571. ROUND_00_15(8, a, b, c, d, e, f, g, h);
  572. T1 = X[9] = PULL64(W[9]);
  573. ROUND_00_15(9, h, a, b, c, d, e, f, g);
  574. T1 = X[10] = PULL64(W[10]);
  575. ROUND_00_15(10, g, h, a, b, c, d, e, f);
  576. T1 = X[11] = PULL64(W[11]);
  577. ROUND_00_15(11, f, g, h, a, b, c, d, e);
  578. T1 = X[12] = PULL64(W[12]);
  579. ROUND_00_15(12, e, f, g, h, a, b, c, d);
  580. T1 = X[13] = PULL64(W[13]);
  581. ROUND_00_15(13, d, e, f, g, h, a, b, c);
  582. T1 = X[14] = PULL64(W[14]);
  583. ROUND_00_15(14, c, d, e, f, g, h, a, b);
  584. T1 = X[15] = PULL64(W[15]);
  585. ROUND_00_15(15, b, c, d, e, f, g, h, a);
  586. # endif
  587. for (i = 16; i < 80; i += 16) {
  588. ROUND_16_80(i, 0, a, b, c, d, e, f, g, h, X);
  589. ROUND_16_80(i, 1, h, a, b, c, d, e, f, g, X);
  590. ROUND_16_80(i, 2, g, h, a, b, c, d, e, f, X);
  591. ROUND_16_80(i, 3, f, g, h, a, b, c, d, e, X);
  592. ROUND_16_80(i, 4, e, f, g, h, a, b, c, d, X);
  593. ROUND_16_80(i, 5, d, e, f, g, h, a, b, c, X);
  594. ROUND_16_80(i, 6, c, d, e, f, g, h, a, b, X);
  595. ROUND_16_80(i, 7, b, c, d, e, f, g, h, a, X);
  596. ROUND_16_80(i, 8, a, b, c, d, e, f, g, h, X);
  597. ROUND_16_80(i, 9, h, a, b, c, d, e, f, g, X);
  598. ROUND_16_80(i, 10, g, h, a, b, c, d, e, f, X);
  599. ROUND_16_80(i, 11, f, g, h, a, b, c, d, e, X);
  600. ROUND_16_80(i, 12, e, f, g, h, a, b, c, d, X);
  601. ROUND_16_80(i, 13, d, e, f, g, h, a, b, c, X);
  602. ROUND_16_80(i, 14, c, d, e, f, g, h, a, b, X);
  603. ROUND_16_80(i, 15, b, c, d, e, f, g, h, a, X);
  604. }
  605. ctx->h[0] += a;
  606. ctx->h[1] += b;
  607. ctx->h[2] += c;
  608. ctx->h[3] += d;
  609. ctx->h[4] += e;
  610. ctx->h[5] += f;
  611. ctx->h[6] += g;
  612. ctx->h[7] += h;
  613. W += SHA_LBLOCK;
  614. }
  615. }
  616. # endif
  617. # endif /* SHA512_ASM */
  618. #else /* !OPENSSL_NO_SHA512 */
  619. # if defined(PEDANTIC) || defined(__DECC) || defined(OPENSSL_SYS_MACOSX)
  620. static void *dummy = &dummy;
  621. # endif
  622. #endif /* !OPENSSL_NO_SHA512 */