You can not select more than 25 topics Topics must start with a letter or number, can include dashes ('-') and can be up to 35 characters long.

sha256.cpp 11KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190191192193194195196197198199200201202203204205206207208209210211212213214215216217218219220221222223224225226227228229230231232233234235236237238239240241242243244245246247248249
  1. // Copyright (c) 2014 The Starwels developers
  2. // Distributed under the MIT software license, see the accompanying
  3. // file COPYING or http://www.opensource.org/licenses/mit-license.php.
  4. #include "crypto/sha256.h"
  5. #include "crypto/common.h"
  6. #include <assert.h>
  7. #include <string.h>
  8. #include <atomic>
  9. #if defined(__x86_64__) || defined(__amd64__)
  10. #if defined(EXPERIMENTAL_ASM)
  11. #include <cpuid.h>
  12. namespace sha256_sse4
  13. {
  14. void Transform(uint32_t* s, const unsigned char* chunk, size_t blocks);
  15. }
  16. #endif
  17. #endif
  18. // Internal implementation code.
  19. namespace
  20. {
  21. /// Internal SHA-256 implementation.
  22. namespace sha256
  23. {
  24. uint32_t inline Ch(uint32_t x, uint32_t y, uint32_t z) { return z ^ (x & (y ^ z)); }
  25. uint32_t inline Maj(uint32_t x, uint32_t y, uint32_t z) { return (x & y) | (z & (x | y)); }
  26. uint32_t inline Sigma0(uint32_t x) { return (x >> 2 | x << 30) ^ (x >> 13 | x << 19) ^ (x >> 22 | x << 10); }
  27. uint32_t inline Sigma1(uint32_t x) { return (x >> 6 | x << 26) ^ (x >> 11 | x << 21) ^ (x >> 25 | x << 7); }
  28. uint32_t inline sigma0(uint32_t x) { return (x >> 7 | x << 25) ^ (x >> 18 | x << 14) ^ (x >> 3); }
  29. uint32_t inline sigma1(uint32_t x) { return (x >> 17 | x << 15) ^ (x >> 19 | x << 13) ^ (x >> 10); }
  30. /** One round of SHA-256. */
  31. void inline Round(uint32_t a, uint32_t b, uint32_t c, uint32_t& d, uint32_t e, uint32_t f, uint32_t g, uint32_t& h, uint32_t k, uint32_t w)
  32. {
  33. uint32_t t1 = h + Sigma1(e) + Ch(e, f, g) + k + w;
  34. uint32_t t2 = Sigma0(a) + Maj(a, b, c);
  35. d += t1;
  36. h = t1 + t2;
  37. }
  38. /** Initialize SHA-256 state. */
  39. void inline Initialize(uint32_t* s)
  40. {
  41. s[0] = 0x6a09e667ul;
  42. s[1] = 0xbb67ae85ul;
  43. s[2] = 0x3c6ef372ul;
  44. s[3] = 0xa54ff53aul;
  45. s[4] = 0x510e527ful;
  46. s[5] = 0x9b05688cul;
  47. s[6] = 0x1f83d9abul;
  48. s[7] = 0x5be0cd19ul;
  49. }
  50. /** Perform a number of SHA-256 transformations, processing 64-byte chunks. */
  51. void Transform(uint32_t* s, const unsigned char* chunk, size_t blocks)
  52. {
  53. while (blocks--) {
  54. uint32_t a = s[0], b = s[1], c = s[2], d = s[3], e = s[4], f = s[5], g = s[6], h = s[7];
  55. uint32_t w0, w1, w2, w3, w4, w5, w6, w7, w8, w9, w10, w11, w12, w13, w14, w15;
  56. Round(a, b, c, d, e, f, g, h, 0x428a2f98, w0 = ReadBE32(chunk + 0));
  57. Round(h, a, b, c, d, e, f, g, 0x71374491, w1 = ReadBE32(chunk + 4));
  58. Round(g, h, a, b, c, d, e, f, 0xb5c0fbcf, w2 = ReadBE32(chunk + 8));
  59. Round(f, g, h, a, b, c, d, e, 0xe9b5dba5, w3 = ReadBE32(chunk + 12));
  60. Round(e, f, g, h, a, b, c, d, 0x3956c25b, w4 = ReadBE32(chunk + 16));
  61. Round(d, e, f, g, h, a, b, c, 0x59f111f1, w5 = ReadBE32(chunk + 20));
  62. Round(c, d, e, f, g, h, a, b, 0x923f82a4, w6 = ReadBE32(chunk + 24));
  63. Round(b, c, d, e, f, g, h, a, 0xab1c5ed5, w7 = ReadBE32(chunk + 28));
  64. Round(a, b, c, d, e, f, g, h, 0xd807aa98, w8 = ReadBE32(chunk + 32));
  65. Round(h, a, b, c, d, e, f, g, 0x12835b01, w9 = ReadBE32(chunk + 36));
  66. Round(g, h, a, b, c, d, e, f, 0x243185be, w10 = ReadBE32(chunk + 40));
  67. Round(f, g, h, a, b, c, d, e, 0x550c7dc3, w11 = ReadBE32(chunk + 44));
  68. Round(e, f, g, h, a, b, c, d, 0x72be5d74, w12 = ReadBE32(chunk + 48));
  69. Round(d, e, f, g, h, a, b, c, 0x80deb1fe, w13 = ReadBE32(chunk + 52));
  70. Round(c, d, e, f, g, h, a, b, 0x9bdc06a7, w14 = ReadBE32(chunk + 56));
  71. Round(b, c, d, e, f, g, h, a, 0xc19bf174, w15 = ReadBE32(chunk + 60));
  72. Round(a, b, c, d, e, f, g, h, 0xe49b69c1, w0 += sigma1(w14) + w9 + sigma0(w1));
  73. Round(h, a, b, c, d, e, f, g, 0xefbe4786, w1 += sigma1(w15) + w10 + sigma0(w2));
  74. Round(g, h, a, b, c, d, e, f, 0x0fc19dc6, w2 += sigma1(w0) + w11 + sigma0(w3));
  75. Round(f, g, h, a, b, c, d, e, 0x240ca1cc, w3 += sigma1(w1) + w12 + sigma0(w4));
  76. Round(e, f, g, h, a, b, c, d, 0x2de92c6f, w4 += sigma1(w2) + w13 + sigma0(w5));
  77. Round(d, e, f, g, h, a, b, c, 0x4a7484aa, w5 += sigma1(w3) + w14 + sigma0(w6));
  78. Round(c, d, e, f, g, h, a, b, 0x5cb0a9dc, w6 += sigma1(w4) + w15 + sigma0(w7));
  79. Round(b, c, d, e, f, g, h, a, 0x76f988da, w7 += sigma1(w5) + w0 + sigma0(w8));
  80. Round(a, b, c, d, e, f, g, h, 0x983e5152, w8 += sigma1(w6) + w1 + sigma0(w9));
  81. Round(h, a, b, c, d, e, f, g, 0xa831c66d, w9 += sigma1(w7) + w2 + sigma0(w10));
  82. Round(g, h, a, b, c, d, e, f, 0xb00327c8, w10 += sigma1(w8) + w3 + sigma0(w11));
  83. Round(f, g, h, a, b, c, d, e, 0xbf597fc7, w11 += sigma1(w9) + w4 + sigma0(w12));
  84. Round(e, f, g, h, a, b, c, d, 0xc6e00bf3, w12 += sigma1(w10) + w5 + sigma0(w13));
  85. Round(d, e, f, g, h, a, b, c, 0xd5a79147, w13 += sigma1(w11) + w6 + sigma0(w14));
  86. Round(c, d, e, f, g, h, a, b, 0x06ca6351, w14 += sigma1(w12) + w7 + sigma0(w15));
  87. Round(b, c, d, e, f, g, h, a, 0x14292967, w15 += sigma1(w13) + w8 + sigma0(w0));
  88. Round(a, b, c, d, e, f, g, h, 0x27b70a85, w0 += sigma1(w14) + w9 + sigma0(w1));
  89. Round(h, a, b, c, d, e, f, g, 0x2e1b2138, w1 += sigma1(w15) + w10 + sigma0(w2));
  90. Round(g, h, a, b, c, d, e, f, 0x4d2c6dfc, w2 += sigma1(w0) + w11 + sigma0(w3));
  91. Round(f, g, h, a, b, c, d, e, 0x53380d13, w3 += sigma1(w1) + w12 + sigma0(w4));
  92. Round(e, f, g, h, a, b, c, d, 0x650a7354, w4 += sigma1(w2) + w13 + sigma0(w5));
  93. Round(d, e, f, g, h, a, b, c, 0x766a0abb, w5 += sigma1(w3) + w14 + sigma0(w6));
  94. Round(c, d, e, f, g, h, a, b, 0x81c2c92e, w6 += sigma1(w4) + w15 + sigma0(w7));
  95. Round(b, c, d, e, f, g, h, a, 0x92722c85, w7 += sigma1(w5) + w0 + sigma0(w8));
  96. Round(a, b, c, d, e, f, g, h, 0xa2bfe8a1, w8 += sigma1(w6) + w1 + sigma0(w9));
  97. Round(h, a, b, c, d, e, f, g, 0xa81a664b, w9 += sigma1(w7) + w2 + sigma0(w10));
  98. Round(g, h, a, b, c, d, e, f, 0xc24b8b70, w10 += sigma1(w8) + w3 + sigma0(w11));
  99. Round(f, g, h, a, b, c, d, e, 0xc76c51a3, w11 += sigma1(w9) + w4 + sigma0(w12));
  100. Round(e, f, g, h, a, b, c, d, 0xd192e819, w12 += sigma1(w10) + w5 + sigma0(w13));
  101. Round(d, e, f, g, h, a, b, c, 0xd6990624, w13 += sigma1(w11) + w6 + sigma0(w14));
  102. Round(c, d, e, f, g, h, a, b, 0xf40e3585, w14 += sigma1(w12) + w7 + sigma0(w15));
  103. Round(b, c, d, e, f, g, h, a, 0x106aa070, w15 += sigma1(w13) + w8 + sigma0(w0));
  104. Round(a, b, c, d, e, f, g, h, 0x19a4c116, w0 += sigma1(w14) + w9 + sigma0(w1));
  105. Round(h, a, b, c, d, e, f, g, 0x1e376c08, w1 += sigma1(w15) + w10 + sigma0(w2));
  106. Round(g, h, a, b, c, d, e, f, 0x2748774c, w2 += sigma1(w0) + w11 + sigma0(w3));
  107. Round(f, g, h, a, b, c, d, e, 0x34b0bcb5, w3 += sigma1(w1) + w12 + sigma0(w4));
  108. Round(e, f, g, h, a, b, c, d, 0x391c0cb3, w4 += sigma1(w2) + w13 + sigma0(w5));
  109. Round(d, e, f, g, h, a, b, c, 0x4ed8aa4a, w5 += sigma1(w3) + w14 + sigma0(w6));
  110. Round(c, d, e, f, g, h, a, b, 0x5b9cca4f, w6 += sigma1(w4) + w15 + sigma0(w7));
  111. Round(b, c, d, e, f, g, h, a, 0x682e6ff3, w7 += sigma1(w5) + w0 + sigma0(w8));
  112. Round(a, b, c, d, e, f, g, h, 0x748f82ee, w8 += sigma1(w6) + w1 + sigma0(w9));
  113. Round(h, a, b, c, d, e, f, g, 0x78a5636f, w9 += sigma1(w7) + w2 + sigma0(w10));
  114. Round(g, h, a, b, c, d, e, f, 0x84c87814, w10 += sigma1(w8) + w3 + sigma0(w11));
  115. Round(f, g, h, a, b, c, d, e, 0x8cc70208, w11 += sigma1(w9) + w4 + sigma0(w12));
  116. Round(e, f, g, h, a, b, c, d, 0x90befffa, w12 += sigma1(w10) + w5 + sigma0(w13));
  117. Round(d, e, f, g, h, a, b, c, 0xa4506ceb, w13 += sigma1(w11) + w6 + sigma0(w14));
  118. Round(c, d, e, f, g, h, a, b, 0xbef9a3f7, w14 + sigma1(w12) + w7 + sigma0(w15));
  119. Round(b, c, d, e, f, g, h, a, 0xc67178f2, w15 + sigma1(w13) + w8 + sigma0(w0));
  120. s[0] += a;
  121. s[1] += b;
  122. s[2] += c;
  123. s[3] += d;
  124. s[4] += e;
  125. s[5] += f;
  126. s[6] += g;
  127. s[7] += h;
  128. chunk += 64;
  129. }
  130. }
  131. } // namespace sha256
  132. typedef void (*TransformType)(uint32_t*, const unsigned char*, size_t);
  133. bool SelfTest(TransformType tr) {
  134. static const unsigned char in1[65] = {0, 0x80};
  135. static const unsigned char in2[129] = {
  136. 0,
  137. 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32,
  138. 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32,
  139. 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,
  140. 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, 2, 0
  141. };
  142. static const uint32_t init[8] = {0x6a09e667ul, 0xbb67ae85ul, 0x3c6ef372ul, 0xa54ff53aul, 0x510e527ful, 0x9b05688cul, 0x1f83d9abul, 0x5be0cd19ul};
  143. static const uint32_t out1[8] = {0xe3b0c442ul, 0x98fc1c14ul, 0x9afbf4c8ul, 0x996fb924ul, 0x27ae41e4ul, 0x649b934cul, 0xa495991bul, 0x7852b855ul};
  144. static const uint32_t out2[8] = {0xce4153b0ul, 0x147c2a86ul, 0x3ed4298eul, 0xe0676bc8ul, 0x79fc77a1ul, 0x2abe1f49ul, 0xb2b055dful, 0x1069523eul};
  145. uint32_t buf[8];
  146. memcpy(buf, init, sizeof(buf));
  147. // Process nothing, and check we remain in the initial state.
  148. tr(buf, nullptr, 0);
  149. if (memcmp(buf, init, sizeof(buf))) return false;
  150. // Process the padded empty string (unaligned)
  151. tr(buf, in1 + 1, 1);
  152. if (memcmp(buf, out1, sizeof(buf))) return false;
  153. // Process 64 spaces (unaligned)
  154. memcpy(buf, init, sizeof(buf));
  155. tr(buf, in2 + 1, 2);
  156. if (memcmp(buf, out2, sizeof(buf))) return false;
  157. return true;
  158. }
  159. TransformType Transform = sha256::Transform;
  160. } // namespace
  161. std::string SHA256AutoDetect()
  162. {
  163. #if defined(EXPERIMENTAL_ASM) && (defined(__x86_64__) || defined(__amd64__))
  164. uint32_t eax, ebx, ecx, edx;
  165. if (__get_cpuid(1, &eax, &ebx, &ecx, &edx) && (ecx >> 19) & 1) {
  166. Transform = sha256_sse4::Transform;
  167. assert(SelfTest(Transform));
  168. return "sse4";
  169. }
  170. #endif
  171. assert(SelfTest(Transform));
  172. return "standard";
  173. }
  174. ////// SHA-256
  175. CSHA256::CSHA256() : bytes(0)
  176. {
  177. sha256::Initialize(s);
  178. }
  179. CSHA256& CSHA256::Write(const unsigned char* data, size_t len)
  180. {
  181. const unsigned char* end = data + len;
  182. size_t bufsize = bytes % 64;
  183. if (bufsize && bufsize + len >= 64) {
  184. // Fill the buffer, and process it.
  185. memcpy(buf + bufsize, data, 64 - bufsize);
  186. bytes += 64 - bufsize;
  187. data += 64 - bufsize;
  188. Transform(s, buf, 1);
  189. bufsize = 0;
  190. }
  191. if (end - data >= 64) {
  192. size_t blocks = (end - data) / 64;
  193. Transform(s, data, blocks);
  194. data += 64 * blocks;
  195. bytes += 64 * blocks;
  196. }
  197. if (end > data) {
  198. // Fill the buffer with what remains.
  199. memcpy(buf + bufsize, data, end - data);
  200. bytes += end - data;
  201. }
  202. return *this;
  203. }
  204. void CSHA256::Finalize(unsigned char hash[OUTPUT_SIZE])
  205. {
  206. static const unsigned char pad[64] = {0x80};
  207. unsigned char sizedesc[8];
  208. WriteBE64(sizedesc, bytes << 3);
  209. Write(pad, 1 + ((119 - (bytes % 64)) % 64));
  210. Write(sizedesc, 8);
  211. WriteBE32(hash, s[0]);
  212. WriteBE32(hash + 4, s[1]);
  213. WriteBE32(hash + 8, s[2]);
  214. WriteBE32(hash + 12, s[3]);
  215. WriteBE32(hash + 16, s[4]);
  216. WriteBE32(hash + 20, s[5]);
  217. WriteBE32(hash + 24, s[6]);
  218. WriteBE32(hash + 28, s[7]);
  219. }
  220. CSHA256& CSHA256::Reset()
  221. {
  222. bytes = 0;
  223. sha256::Initialize(s);
  224. return *this;
  225. }