SHA512.cpp 8.4 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190191192193194195196197198199200201202203204205206207208209210211212213214
  1. // Code taken from NaCl by D. J. Bernstein and others
  2. // Public domain
  3. /*
  4. 20080913
  5. D. J. Bernstein
  6. Public domain.
  7. */
  8. #include <stdint.h>
  9. #include <stdlib.h>
  10. #include <string.h>
  11. #include "SHA512.hpp"
  12. #include "Utils.hpp"
  13. #ifndef ZT_HAVE_NATIVE_SHA512
  14. namespace ZeroTier {
  15. namespace {
  16. static inline void sha512_encode(uint64_t input, uint8_t *output, uint32_t idx)
  17. {
  18. output[idx + 0] = (uint8_t)(input >> 56);
  19. output[idx + 1] = (uint8_t)(input >> 48);
  20. output[idx + 2] = (uint8_t)(input >> 40);
  21. output[idx + 3] = (uint8_t)(input >> 32);
  22. output[idx + 4] = (uint8_t)(input >> 24);
  23. output[idx + 5] = (uint8_t)(input >> 16);
  24. output[idx + 6] = (uint8_t)(input >> 8);
  25. output[idx + 7] = (uint8_t)(input >> 0);
  26. }
  27. static inline void sha512_decode(uint64_t *output, uint8_t *input, uint32_t idx)
  28. {
  29. *output = ((uint64_t)input[idx + 0] << 56)
  30. | ((uint64_t)input[idx + 1] << 48)
  31. | ((uint64_t)input[idx + 2] << 40)
  32. | ((uint64_t)input[idx + 3] << 32)
  33. | ((uint64_t)input[idx + 4] << 24)
  34. | ((uint64_t)input[idx + 5] << 16)
  35. | ((uint64_t)input[idx + 6] << 8)
  36. | ((uint64_t)input[idx + 7] << 0);
  37. }
  38. typedef struct sha512_ctx_tag {
  39. uint32_t is_sha384;
  40. uint8_t block[128];
  41. uint64_t len[2];
  42. uint64_t val[8];
  43. uint8_t *payload_addr;
  44. uint64_t payload_len;
  45. } sha512_ctx_t;
  46. #define LSR(x,n) (x >> n)
  47. #define ROR(x,n) (LSR(x,n) | (x << (64 - n)))
  48. #define MA(x,y,z) ((x & y) | (z & (x | y)))
  49. #define CH(x,y,z) (z ^ (x & (y ^ z)))
  50. #define GAMMA0(x) (ROR(x, 1) ^ ROR(x, 8) ^ LSR(x, 7))
  51. #define GAMMA1(x) (ROR(x,19) ^ ROR(x,61) ^ LSR(x, 6))
  52. #define SIGMA0(x) (ROR(x,28) ^ ROR(x,34) ^ ROR(x,39))
  53. #define SIGMA1(x) (ROR(x,14) ^ ROR(x,18) ^ ROR(x,41))
  54. #define INIT_COMPRESSOR() uint64_t tmp0 = 0, tmp1 = 0
  55. #define COMPRESS( a, b, c, d, e, f, g, h, x, k) \
  56. tmp0 = h + SIGMA1(e) + CH(e,f,g) + k + x; \
  57. tmp1 = SIGMA0(a) + MA(a,b,c); d += tmp0; h = tmp0 + tmp1;
  58. static const uint8_t sha512_padding[128] = { 0x80, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 };
  59. static const uint64_t K[80] = {
  60. 0x428A2F98D728AE22ULL, 0x7137449123EF65CDULL, 0xB5C0FBCFEC4D3B2FULL, 0xE9B5DBA58189DBBCULL,
  61. 0x3956C25BF348B538ULL, 0x59F111F1B605D019ULL, 0x923F82A4AF194F9BULL, 0xAB1C5ED5DA6D8118ULL,
  62. 0xD807AA98A3030242ULL, 0x12835B0145706FBEULL, 0x243185BE4EE4B28CULL, 0x550C7DC3D5FFB4E2ULL,
  63. 0x72BE5D74F27B896FULL, 0x80DEB1FE3B1696B1ULL, 0x9BDC06A725C71235ULL, 0xC19BF174CF692694ULL,
  64. 0xE49B69C19EF14AD2ULL, 0xEFBE4786384F25E3ULL, 0x0FC19DC68B8CD5B5ULL, 0x240CA1CC77AC9C65ULL,
  65. 0x2DE92C6F592B0275ULL, 0x4A7484AA6EA6E483ULL, 0x5CB0A9DCBD41FBD4ULL, 0x76F988DA831153B5ULL,
  66. 0x983E5152EE66DFABULL, 0xA831C66D2DB43210ULL, 0xB00327C898FB213FULL, 0xBF597FC7BEEF0EE4ULL,
  67. 0xC6E00BF33DA88FC2ULL, 0xD5A79147930AA725ULL, 0x06CA6351E003826FULL, 0x142929670A0E6E70ULL,
  68. 0x27B70A8546D22FFCULL, 0x2E1B21385C26C926ULL, 0x4D2C6DFC5AC42AEDULL, 0x53380D139D95B3DFULL,
  69. 0x650A73548BAF63DEULL, 0x766A0ABB3C77B2A8ULL, 0x81C2C92E47EDAEE6ULL, 0x92722C851482353BULL,
  70. 0xA2BFE8A14CF10364ULL, 0xA81A664BBC423001ULL, 0xC24B8B70D0F89791ULL, 0xC76C51A30654BE30ULL,
  71. 0xD192E819D6EF5218ULL, 0xD69906245565A910ULL, 0xF40E35855771202AULL, 0x106AA07032BBD1B8ULL,
  72. 0x19A4C116B8D2D0C8ULL, 0x1E376C085141AB53ULL, 0x2748774CDF8EEB99ULL, 0x34B0BCB5E19B48A8ULL,
  73. 0x391C0CB3C5C95A63ULL, 0x4ED8AA4AE3418ACBULL, 0x5B9CCA4F7763E373ULL, 0x682E6FF3D6B2B8A3ULL,
  74. 0x748F82EE5DEFB2FCULL, 0x78A5636F43172F60ULL, 0x84C87814A1F0AB72ULL, 0x8CC702081A6439ECULL,
  75. 0x90BEFFFA23631E28ULL, 0xA4506CEBDE82BDE9ULL, 0xBEF9A3F7B2C67915ULL, 0xC67178F2E372532BULL,
  76. 0xCA273ECEEA26619CULL, 0xD186B8C721C0C207ULL, 0xEADA7DD6CDE0EB1EULL, 0xF57D4F7FEE6ED178ULL,
  77. 0x06F067AA72176FBAULL, 0x0A637DC5A2C898A6ULL, 0x113F9804BEF90DAEULL, 0x1B710B35131C471BULL,
  78. 0x28DB77F523047D84ULL, 0x32CAAB7B40C72493ULL, 0x3C9EBE0A15C9BEBCULL, 0x431D67C49C100D4CULL,
  79. 0x4CC5D4BECB3E42B6ULL, 0x597F299CFC657E2AULL, 0x5FCB6FAB3AD6FAECULL, 0x6C44198C4A475817ULL
  80. };
  81. static inline void sha512_memcpy(uint8_t *src, uint8_t *dst, uint32_t size)
  82. {
  83. uint32_t i = 0;
  84. for (;i < size;i++) { *dst++ = *src++; }
  85. }
  86. static inline void sha512_memclr(uint8_t *dst, uint32_t size)
  87. {
  88. uint32_t i = 0;
  89. for (;i < size;i++) { *dst++ = 0; }
  90. }
  91. static inline void sha512_init_512(sha512_ctx_t *sha512_ctx, uint8_t *payload_addr, uint64_t payload_len)
  92. {
  93. sha512_memclr((uint8_t *)sha512_ctx,sizeof(sha512_ctx_t));
  94. sha512_ctx->val[0] = 0x6A09E667F3BCC908ULL;
  95. sha512_ctx->val[1] = 0xBB67AE8584CAA73BULL;
  96. sha512_ctx->val[2] = 0x3C6EF372FE94F82BULL;
  97. sha512_ctx->val[3] = 0xA54FF53A5F1D36F1ULL;
  98. sha512_ctx->val[4] = 0x510E527FADE682D1ULL;
  99. sha512_ctx->val[5] = 0x9B05688C2B3E6C1FULL;
  100. sha512_ctx->val[6] = 0x1F83D9ABFB41BD6BULL;
  101. sha512_ctx->val[7] = 0x5BE0CD19137E2179ULL;
  102. sha512_ctx->is_sha384 = 0;
  103. sha512_ctx->payload_addr = payload_addr;
  104. sha512_ctx->payload_len = (uint64_t)payload_len;
  105. sha512_ctx->len[0] = payload_len << 3;
  106. sha512_ctx->len[1] = payload_len >> 61;
  107. }
  108. static inline void sha512_init_384(sha512_ctx_t *sha512_ctx, uint8_t *payload_addr, uint64_t payload_len)
  109. {
  110. sha512_memclr((uint8_t *)sha512_ctx,sizeof(sha512_ctx_t));
  111. sha512_ctx->val[0] = 0xCBBB9D5DC1059ED8ULL;
  112. sha512_ctx->val[1] = 0x629A292A367CD507ULL;
  113. sha512_ctx->val[2] = 0x9159015A3070DD17ULL;
  114. sha512_ctx->val[3] = 0x152FECD8F70E5939ULL;
  115. sha512_ctx->val[4] = 0x67332667FFC00B31ULL;
  116. sha512_ctx->val[5] = 0x8EB44A8768581511ULL;
  117. sha512_ctx->val[6] = 0xDB0C2E0D64F98FA7ULL;
  118. sha512_ctx->val[7] = 0x47B5481DBEFA4FA4ULL;
  119. sha512_ctx->is_sha384 = 1;
  120. sha512_ctx->payload_addr = payload_addr;
  121. sha512_ctx->payload_len = (uint64_t)payload_len;
  122. sha512_ctx->len[0] = payload_len << 3;
  123. sha512_ctx->len[1] = payload_len >> 61;
  124. }
  125. static inline void sha512_hash_factory(sha512_ctx_t *ctx, uint8_t data[128])
  126. {
  127. uint32_t i = 0;
  128. uint64_t W[80];
  129. uint64_t v[8];
  130. INIT_COMPRESSOR();
  131. for(i = 0; i < 16; i++) { sha512_decode(&W[i], data, i << 3 ); }
  132. for(; i < 80; i++) { W[i] = GAMMA1(W[i - 2]) + W[i - 7] + GAMMA0(W[i - 15]) + W[i - 16]; }
  133. for (i = 0;i < 8; i++) { v[i] = ctx->val[i]; }
  134. for(i = 0; i < 80;) {
  135. COMPRESS(v[0], v[1], v[2], v[3], v[4], v[5], v[6], v[7], W[i], K[i] ); i++;
  136. COMPRESS(v[7], v[0], v[1], v[2], v[3], v[4], v[5], v[6], W[i], K[i] ); i++;
  137. COMPRESS(v[6], v[7], v[0], v[1], v[2], v[3], v[4], v[5], W[i], K[i] ); i++;
  138. COMPRESS(v[5], v[6], v[7], v[0], v[1], v[2], v[3], v[4], W[i], K[i] ); i++;
  139. COMPRESS(v[4], v[5], v[6], v[7], v[0], v[1], v[2], v[3], W[i], K[i] ); i++;
  140. COMPRESS(v[3], v[4], v[5], v[6], v[7], v[0], v[1], v[2], W[i], K[i] ); i++;
  141. COMPRESS(v[2], v[3], v[4], v[5], v[6], v[7], v[0], v[1], W[i], K[i] ); i++;
  142. COMPRESS(v[1], v[2], v[3], v[4], v[5], v[6], v[7], v[0], W[i], K[i] ); i++;
  143. }
  144. for (i = 0; i < 8; i++) { ctx->val[i] += v[i]; }
  145. }
  146. static inline void sha512_stage1(sha512_ctx_t *sha512_ctx)
  147. {
  148. while (sha512_ctx->payload_len >= 128) {
  149. sha512_hash_factory(sha512_ctx, sha512_ctx->payload_addr);
  150. sha512_ctx->payload_addr += 128;
  151. sha512_ctx->payload_len -= 128;
  152. }
  153. }
  154. static inline void sha512_stage2(sha512_ctx_t *sha512_ctx, uint8_t output[64])
  155. {
  156. uint32_t block_pos = sha512_ctx->payload_len;
  157. uint32_t padding_bytes = 0;
  158. uint8_t temp_data[128] = {0};
  159. uint8_t *temp_data_p = (uint8_t *)&temp_data[0];
  160. uint8_t len_be[16] = {0};
  161. uint8_t i = 0;
  162. sha512_memcpy(sha512_ctx->payload_addr, temp_data_p, sha512_ctx->payload_len);
  163. padding_bytes = 112 - block_pos;
  164. temp_data_p += block_pos;
  165. sha512_memcpy((uint8_t *)sha512_padding, temp_data_p, padding_bytes);
  166. temp_data_p += padding_bytes;
  167. sha512_encode(sha512_ctx->len[1], len_be, 0);
  168. sha512_encode(sha512_ctx->len[0], len_be, 8);
  169. sha512_memcpy(len_be, temp_data_p, 16);
  170. sha512_hash_factory(sha512_ctx, temp_data);
  171. for (i = 0; i < 6; i++) { sha512_encode(sha512_ctx->val[i], output, i * 8); }
  172. for ( ;(i < 8) && (sha512_ctx->is_sha384 == 0); i++) { sha512_encode(sha512_ctx->val[i], output, i * 8); }
  173. }
  174. } // anonymous namespace
  175. void SHA512(void *digest,const void *data,unsigned int len)
  176. {
  177. sha512_ctx_t h;
  178. sha512_init_512(&h,(uint8_t *)data,len);
  179. sha512_stage1(&h);
  180. sha512_stage2(&h,(uint8_t *)digest);
  181. }
  182. void SHA384(void *digest,const void *data,unsigned int len)
  183. {
  184. sha512_ctx_t h;
  185. sha512_init_384(&h,(uint8_t *)data,len);
  186. sha512_stage1(&h);
  187. sha512_stage2(&h,(uint8_t *)digest);
  188. }
  189. } // namespace ZeroTier
  190. #endif // !ZT_HAVE_NATIVE_SHA512