選択できるのは25トピックまでです。 トピックは、先頭が英数字で、英数字とダッシュ('-')を使用した35文字以内のものにしてください。

keccakf.go 10KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190191192193194195196197198199200201202203204205206207208209210211212213214215216217218219220221222223224225226227228229230231232233234235236237238239240241242243244245246247248249250251252253254255256257258259260261262263264265266267268269270271272273274275276277278279280281282283284285286287288289290291292293294295296297298299300301302303304305306307308309310311312313314315316317318319320321322323324325326327328329330331332333334335336337338339340341342343344345346347348349350351352353354355356357358359360361362363364365366367368369370371372373374375376377378379380381382383384385386387388389390391392393394395396397398399400401402403404405406407408409410411412413414415
  1. // Copyright 2014 The Go Authors. All rights reserved.
  2. // Use of this source code is governed by a BSD-style
  3. // license that can be found in the LICENSE file.
  4. //go:build !amd64 || purego || !gc
  5. // +build !amd64 purego !gc
  6. package sha3
  7. import "math/bits"
  8. // rc stores the round constants for use in the ι step.
  9. var rc = [24]uint64{
  10. 0x0000000000000001,
  11. 0x0000000000008082,
  12. 0x800000000000808A,
  13. 0x8000000080008000,
  14. 0x000000000000808B,
  15. 0x0000000080000001,
  16. 0x8000000080008081,
  17. 0x8000000000008009,
  18. 0x000000000000008A,
  19. 0x0000000000000088,
  20. 0x0000000080008009,
  21. 0x000000008000000A,
  22. 0x000000008000808B,
  23. 0x800000000000008B,
  24. 0x8000000000008089,
  25. 0x8000000000008003,
  26. 0x8000000000008002,
  27. 0x8000000000000080,
  28. 0x000000000000800A,
  29. 0x800000008000000A,
  30. 0x8000000080008081,
  31. 0x8000000000008080,
  32. 0x0000000080000001,
  33. 0x8000000080008008,
  34. }
  35. // keccakF1600 applies the Keccak permutation to a 1600b-wide
  36. // state represented as a slice of 25 uint64s.
  37. func keccakF1600(a *[25]uint64) {
  38. // Implementation translated from Keccak-inplace.c
  39. // in the keccak reference code.
  40. var t, bc0, bc1, bc2, bc3, bc4, d0, d1, d2, d3, d4 uint64
  41. for i := 0; i < 24; i += 4 {
  42. // Combines the 5 steps in each round into 2 steps.
  43. // Unrolls 4 rounds per loop and spreads some steps across rounds.
  44. // Round 1
  45. bc0 = a[0] ^ a[5] ^ a[10] ^ a[15] ^ a[20]
  46. bc1 = a[1] ^ a[6] ^ a[11] ^ a[16] ^ a[21]
  47. bc2 = a[2] ^ a[7] ^ a[12] ^ a[17] ^ a[22]
  48. bc3 = a[3] ^ a[8] ^ a[13] ^ a[18] ^ a[23]
  49. bc4 = a[4] ^ a[9] ^ a[14] ^ a[19] ^ a[24]
  50. d0 = bc4 ^ (bc1<<1 | bc1>>63)
  51. d1 = bc0 ^ (bc2<<1 | bc2>>63)
  52. d2 = bc1 ^ (bc3<<1 | bc3>>63)
  53. d3 = bc2 ^ (bc4<<1 | bc4>>63)
  54. d4 = bc3 ^ (bc0<<1 | bc0>>63)
  55. bc0 = a[0] ^ d0
  56. t = a[6] ^ d1
  57. bc1 = bits.RotateLeft64(t, 44)
  58. t = a[12] ^ d2
  59. bc2 = bits.RotateLeft64(t, 43)
  60. t = a[18] ^ d3
  61. bc3 = bits.RotateLeft64(t, 21)
  62. t = a[24] ^ d4
  63. bc4 = bits.RotateLeft64(t, 14)
  64. a[0] = bc0 ^ (bc2 &^ bc1) ^ rc[i]
  65. a[6] = bc1 ^ (bc3 &^ bc2)
  66. a[12] = bc2 ^ (bc4 &^ bc3)
  67. a[18] = bc3 ^ (bc0 &^ bc4)
  68. a[24] = bc4 ^ (bc1 &^ bc0)
  69. t = a[10] ^ d0
  70. bc2 = bits.RotateLeft64(t, 3)
  71. t = a[16] ^ d1
  72. bc3 = bits.RotateLeft64(t, 45)
  73. t = a[22] ^ d2
  74. bc4 = bits.RotateLeft64(t, 61)
  75. t = a[3] ^ d3
  76. bc0 = bits.RotateLeft64(t, 28)
  77. t = a[9] ^ d4
  78. bc1 = bits.RotateLeft64(t, 20)
  79. a[10] = bc0 ^ (bc2 &^ bc1)
  80. a[16] = bc1 ^ (bc3 &^ bc2)
  81. a[22] = bc2 ^ (bc4 &^ bc3)
  82. a[3] = bc3 ^ (bc0 &^ bc4)
  83. a[9] = bc4 ^ (bc1 &^ bc0)
  84. t = a[20] ^ d0
  85. bc4 = bits.RotateLeft64(t, 18)
  86. t = a[1] ^ d1
  87. bc0 = bits.RotateLeft64(t, 1)
  88. t = a[7] ^ d2
  89. bc1 = bits.RotateLeft64(t, 6)
  90. t = a[13] ^ d3
  91. bc2 = bits.RotateLeft64(t, 25)
  92. t = a[19] ^ d4
  93. bc3 = bits.RotateLeft64(t, 8)
  94. a[20] = bc0 ^ (bc2 &^ bc1)
  95. a[1] = bc1 ^ (bc3 &^ bc2)
  96. a[7] = bc2 ^ (bc4 &^ bc3)
  97. a[13] = bc3 ^ (bc0 &^ bc4)
  98. a[19] = bc4 ^ (bc1 &^ bc0)
  99. t = a[5] ^ d0
  100. bc1 = bits.RotateLeft64(t, 36)
  101. t = a[11] ^ d1
  102. bc2 = bits.RotateLeft64(t, 10)
  103. t = a[17] ^ d2
  104. bc3 = bits.RotateLeft64(t, 15)
  105. t = a[23] ^ d3
  106. bc4 = bits.RotateLeft64(t, 56)
  107. t = a[4] ^ d4
  108. bc0 = bits.RotateLeft64(t, 27)
  109. a[5] = bc0 ^ (bc2 &^ bc1)
  110. a[11] = bc1 ^ (bc3 &^ bc2)
  111. a[17] = bc2 ^ (bc4 &^ bc3)
  112. a[23] = bc3 ^ (bc0 &^ bc4)
  113. a[4] = bc4 ^ (bc1 &^ bc0)
  114. t = a[15] ^ d0
  115. bc3 = bits.RotateLeft64(t, 41)
  116. t = a[21] ^ d1
  117. bc4 = bits.RotateLeft64(t, 2)
  118. t = a[2] ^ d2
  119. bc0 = bits.RotateLeft64(t, 62)
  120. t = a[8] ^ d3
  121. bc1 = bits.RotateLeft64(t, 55)
  122. t = a[14] ^ d4
  123. bc2 = bits.RotateLeft64(t, 39)
  124. a[15] = bc0 ^ (bc2 &^ bc1)
  125. a[21] = bc1 ^ (bc3 &^ bc2)
  126. a[2] = bc2 ^ (bc4 &^ bc3)
  127. a[8] = bc3 ^ (bc0 &^ bc4)
  128. a[14] = bc4 ^ (bc1 &^ bc0)
  129. // Round 2
  130. bc0 = a[0] ^ a[5] ^ a[10] ^ a[15] ^ a[20]
  131. bc1 = a[1] ^ a[6] ^ a[11] ^ a[16] ^ a[21]
  132. bc2 = a[2] ^ a[7] ^ a[12] ^ a[17] ^ a[22]
  133. bc3 = a[3] ^ a[8] ^ a[13] ^ a[18] ^ a[23]
  134. bc4 = a[4] ^ a[9] ^ a[14] ^ a[19] ^ a[24]
  135. d0 = bc4 ^ (bc1<<1 | bc1>>63)
  136. d1 = bc0 ^ (bc2<<1 | bc2>>63)
  137. d2 = bc1 ^ (bc3<<1 | bc3>>63)
  138. d3 = bc2 ^ (bc4<<1 | bc4>>63)
  139. d4 = bc3 ^ (bc0<<1 | bc0>>63)
  140. bc0 = a[0] ^ d0
  141. t = a[16] ^ d1
  142. bc1 = bits.RotateLeft64(t, 44)
  143. t = a[7] ^ d2
  144. bc2 = bits.RotateLeft64(t, 43)
  145. t = a[23] ^ d3
  146. bc3 = bits.RotateLeft64(t, 21)
  147. t = a[14] ^ d4
  148. bc4 = bits.RotateLeft64(t, 14)
  149. a[0] = bc0 ^ (bc2 &^ bc1) ^ rc[i+1]
  150. a[16] = bc1 ^ (bc3 &^ bc2)
  151. a[7] = bc2 ^ (bc4 &^ bc3)
  152. a[23] = bc3 ^ (bc0 &^ bc4)
  153. a[14] = bc4 ^ (bc1 &^ bc0)
  154. t = a[20] ^ d0
  155. bc2 = bits.RotateLeft64(t, 3)
  156. t = a[11] ^ d1
  157. bc3 = bits.RotateLeft64(t, 45)
  158. t = a[2] ^ d2
  159. bc4 = bits.RotateLeft64(t, 61)
  160. t = a[18] ^ d3
  161. bc0 = bits.RotateLeft64(t, 28)
  162. t = a[9] ^ d4
  163. bc1 = bits.RotateLeft64(t, 20)
  164. a[20] = bc0 ^ (bc2 &^ bc1)
  165. a[11] = bc1 ^ (bc3 &^ bc2)
  166. a[2] = bc2 ^ (bc4 &^ bc3)
  167. a[18] = bc3 ^ (bc0 &^ bc4)
  168. a[9] = bc4 ^ (bc1 &^ bc0)
  169. t = a[15] ^ d0
  170. bc4 = bits.RotateLeft64(t, 18)
  171. t = a[6] ^ d1
  172. bc0 = bits.RotateLeft64(t, 1)
  173. t = a[22] ^ d2
  174. bc1 = bits.RotateLeft64(t, 6)
  175. t = a[13] ^ d3
  176. bc2 = bits.RotateLeft64(t, 25)
  177. t = a[4] ^ d4
  178. bc3 = bits.RotateLeft64(t, 8)
  179. a[15] = bc0 ^ (bc2 &^ bc1)
  180. a[6] = bc1 ^ (bc3 &^ bc2)
  181. a[22] = bc2 ^ (bc4 &^ bc3)
  182. a[13] = bc3 ^ (bc0 &^ bc4)
  183. a[4] = bc4 ^ (bc1 &^ bc0)
  184. t = a[10] ^ d0
  185. bc1 = bits.RotateLeft64(t, 36)
  186. t = a[1] ^ d1
  187. bc2 = bits.RotateLeft64(t, 10)
  188. t = a[17] ^ d2
  189. bc3 = bits.RotateLeft64(t, 15)
  190. t = a[8] ^ d3
  191. bc4 = bits.RotateLeft64(t, 56)
  192. t = a[24] ^ d4
  193. bc0 = bits.RotateLeft64(t, 27)
  194. a[10] = bc0 ^ (bc2 &^ bc1)
  195. a[1] = bc1 ^ (bc3 &^ bc2)
  196. a[17] = bc2 ^ (bc4 &^ bc3)
  197. a[8] = bc3 ^ (bc0 &^ bc4)
  198. a[24] = bc4 ^ (bc1 &^ bc0)
  199. t = a[5] ^ d0
  200. bc3 = bits.RotateLeft64(t, 41)
  201. t = a[21] ^ d1
  202. bc4 = bits.RotateLeft64(t, 2)
  203. t = a[12] ^ d2
  204. bc0 = bits.RotateLeft64(t, 62)
  205. t = a[3] ^ d3
  206. bc1 = bits.RotateLeft64(t, 55)
  207. t = a[19] ^ d4
  208. bc2 = bits.RotateLeft64(t, 39)
  209. a[5] = bc0 ^ (bc2 &^ bc1)
  210. a[21] = bc1 ^ (bc3 &^ bc2)
  211. a[12] = bc2 ^ (bc4 &^ bc3)
  212. a[3] = bc3 ^ (bc0 &^ bc4)
  213. a[19] = bc4 ^ (bc1 &^ bc0)
  214. // Round 3
  215. bc0 = a[0] ^ a[5] ^ a[10] ^ a[15] ^ a[20]
  216. bc1 = a[1] ^ a[6] ^ a[11] ^ a[16] ^ a[21]
  217. bc2 = a[2] ^ a[7] ^ a[12] ^ a[17] ^ a[22]
  218. bc3 = a[3] ^ a[8] ^ a[13] ^ a[18] ^ a[23]
  219. bc4 = a[4] ^ a[9] ^ a[14] ^ a[19] ^ a[24]
  220. d0 = bc4 ^ (bc1<<1 | bc1>>63)
  221. d1 = bc0 ^ (bc2<<1 | bc2>>63)
  222. d2 = bc1 ^ (bc3<<1 | bc3>>63)
  223. d3 = bc2 ^ (bc4<<1 | bc4>>63)
  224. d4 = bc3 ^ (bc0<<1 | bc0>>63)
  225. bc0 = a[0] ^ d0
  226. t = a[11] ^ d1
  227. bc1 = bits.RotateLeft64(t, 44)
  228. t = a[22] ^ d2
  229. bc2 = bits.RotateLeft64(t, 43)
  230. t = a[8] ^ d3
  231. bc3 = bits.RotateLeft64(t, 21)
  232. t = a[19] ^ d4
  233. bc4 = bits.RotateLeft64(t, 14)
  234. a[0] = bc0 ^ (bc2 &^ bc1) ^ rc[i+2]
  235. a[11] = bc1 ^ (bc3 &^ bc2)
  236. a[22] = bc2 ^ (bc4 &^ bc3)
  237. a[8] = bc3 ^ (bc0 &^ bc4)
  238. a[19] = bc4 ^ (bc1 &^ bc0)
  239. t = a[15] ^ d0
  240. bc2 = bits.RotateLeft64(t, 3)
  241. t = a[1] ^ d1
  242. bc3 = bits.RotateLeft64(t, 45)
  243. t = a[12] ^ d2
  244. bc4 = bits.RotateLeft64(t, 61)
  245. t = a[23] ^ d3
  246. bc0 = bits.RotateLeft64(t, 28)
  247. t = a[9] ^ d4
  248. bc1 = bits.RotateLeft64(t, 20)
  249. a[15] = bc0 ^ (bc2 &^ bc1)
  250. a[1] = bc1 ^ (bc3 &^ bc2)
  251. a[12] = bc2 ^ (bc4 &^ bc3)
  252. a[23] = bc3 ^ (bc0 &^ bc4)
  253. a[9] = bc4 ^ (bc1 &^ bc0)
  254. t = a[5] ^ d0
  255. bc4 = bits.RotateLeft64(t, 18)
  256. t = a[16] ^ d1
  257. bc0 = bits.RotateLeft64(t, 1)
  258. t = a[2] ^ d2
  259. bc1 = bits.RotateLeft64(t, 6)
  260. t = a[13] ^ d3
  261. bc2 = bits.RotateLeft64(t, 25)
  262. t = a[24] ^ d4
  263. bc3 = bits.RotateLeft64(t, 8)
  264. a[5] = bc0 ^ (bc2 &^ bc1)
  265. a[16] = bc1 ^ (bc3 &^ bc2)
  266. a[2] = bc2 ^ (bc4 &^ bc3)
  267. a[13] = bc3 ^ (bc0 &^ bc4)
  268. a[24] = bc4 ^ (bc1 &^ bc0)
  269. t = a[20] ^ d0
  270. bc1 = bits.RotateLeft64(t, 36)
  271. t = a[6] ^ d1
  272. bc2 = bits.RotateLeft64(t, 10)
  273. t = a[17] ^ d2
  274. bc3 = bits.RotateLeft64(t, 15)
  275. t = a[3] ^ d3
  276. bc4 = bits.RotateLeft64(t, 56)
  277. t = a[14] ^ d4
  278. bc0 = bits.RotateLeft64(t, 27)
  279. a[20] = bc0 ^ (bc2 &^ bc1)
  280. a[6] = bc1 ^ (bc3 &^ bc2)
  281. a[17] = bc2 ^ (bc4 &^ bc3)
  282. a[3] = bc3 ^ (bc0 &^ bc4)
  283. a[14] = bc4 ^ (bc1 &^ bc0)
  284. t = a[10] ^ d0
  285. bc3 = bits.RotateLeft64(t, 41)
  286. t = a[21] ^ d1
  287. bc4 = bits.RotateLeft64(t, 2)
  288. t = a[7] ^ d2
  289. bc0 = bits.RotateLeft64(t, 62)
  290. t = a[18] ^ d3
  291. bc1 = bits.RotateLeft64(t, 55)
  292. t = a[4] ^ d4
  293. bc2 = bits.RotateLeft64(t, 39)
  294. a[10] = bc0 ^ (bc2 &^ bc1)
  295. a[21] = bc1 ^ (bc3 &^ bc2)
  296. a[7] = bc2 ^ (bc4 &^ bc3)
  297. a[18] = bc3 ^ (bc0 &^ bc4)
  298. a[4] = bc4 ^ (bc1 &^ bc0)
  299. // Round 4
  300. bc0 = a[0] ^ a[5] ^ a[10] ^ a[15] ^ a[20]
  301. bc1 = a[1] ^ a[6] ^ a[11] ^ a[16] ^ a[21]
  302. bc2 = a[2] ^ a[7] ^ a[12] ^ a[17] ^ a[22]
  303. bc3 = a[3] ^ a[8] ^ a[13] ^ a[18] ^ a[23]
  304. bc4 = a[4] ^ a[9] ^ a[14] ^ a[19] ^ a[24]
  305. d0 = bc4 ^ (bc1<<1 | bc1>>63)
  306. d1 = bc0 ^ (bc2<<1 | bc2>>63)
  307. d2 = bc1 ^ (bc3<<1 | bc3>>63)
  308. d3 = bc2 ^ (bc4<<1 | bc4>>63)
  309. d4 = bc3 ^ (bc0<<1 | bc0>>63)
  310. bc0 = a[0] ^ d0
  311. t = a[1] ^ d1
  312. bc1 = bits.RotateLeft64(t, 44)
  313. t = a[2] ^ d2
  314. bc2 = bits.RotateLeft64(t, 43)
  315. t = a[3] ^ d3
  316. bc3 = bits.RotateLeft64(t, 21)
  317. t = a[4] ^ d4
  318. bc4 = bits.RotateLeft64(t, 14)
  319. a[0] = bc0 ^ (bc2 &^ bc1) ^ rc[i+3]
  320. a[1] = bc1 ^ (bc3 &^ bc2)
  321. a[2] = bc2 ^ (bc4 &^ bc3)
  322. a[3] = bc3 ^ (bc0 &^ bc4)
  323. a[4] = bc4 ^ (bc1 &^ bc0)
  324. t = a[5] ^ d0
  325. bc2 = bits.RotateLeft64(t, 3)
  326. t = a[6] ^ d1
  327. bc3 = bits.RotateLeft64(t, 45)
  328. t = a[7] ^ d2
  329. bc4 = bits.RotateLeft64(t, 61)
  330. t = a[8] ^ d3
  331. bc0 = bits.RotateLeft64(t, 28)
  332. t = a[9] ^ d4
  333. bc1 = bits.RotateLeft64(t, 20)
  334. a[5] = bc0 ^ (bc2 &^ bc1)
  335. a[6] = bc1 ^ (bc3 &^ bc2)
  336. a[7] = bc2 ^ (bc4 &^ bc3)
  337. a[8] = bc3 ^ (bc0 &^ bc4)
  338. a[9] = bc4 ^ (bc1 &^ bc0)
  339. t = a[10] ^ d0
  340. bc4 = bits.RotateLeft64(t, 18)
  341. t = a[11] ^ d1
  342. bc0 = bits.RotateLeft64(t, 1)
  343. t = a[12] ^ d2
  344. bc1 = bits.RotateLeft64(t, 6)
  345. t = a[13] ^ d3
  346. bc2 = bits.RotateLeft64(t, 25)
  347. t = a[14] ^ d4
  348. bc3 = bits.RotateLeft64(t, 8)
  349. a[10] = bc0 ^ (bc2 &^ bc1)
  350. a[11] = bc1 ^ (bc3 &^ bc2)
  351. a[12] = bc2 ^ (bc4 &^ bc3)
  352. a[13] = bc3 ^ (bc0 &^ bc4)
  353. a[14] = bc4 ^ (bc1 &^ bc0)
  354. t = a[15] ^ d0
  355. bc1 = bits.RotateLeft64(t, 36)
  356. t = a[16] ^ d1
  357. bc2 = bits.RotateLeft64(t, 10)
  358. t = a[17] ^ d2
  359. bc3 = bits.RotateLeft64(t, 15)
  360. t = a[18] ^ d3
  361. bc4 = bits.RotateLeft64(t, 56)
  362. t = a[19] ^ d4
  363. bc0 = bits.RotateLeft64(t, 27)
  364. a[15] = bc0 ^ (bc2 &^ bc1)
  365. a[16] = bc1 ^ (bc3 &^ bc2)
  366. a[17] = bc2 ^ (bc4 &^ bc3)
  367. a[18] = bc3 ^ (bc0 &^ bc4)
  368. a[19] = bc4 ^ (bc1 &^ bc0)
  369. t = a[20] ^ d0
  370. bc3 = bits.RotateLeft64(t, 41)
  371. t = a[21] ^ d1
  372. bc4 = bits.RotateLeft64(t, 2)
  373. t = a[22] ^ d2
  374. bc0 = bits.RotateLeft64(t, 62)
  375. t = a[23] ^ d3
  376. bc1 = bits.RotateLeft64(t, 55)
  377. t = a[24] ^ d4
  378. bc2 = bits.RotateLeft64(t, 39)
  379. a[20] = bc0 ^ (bc2 &^ bc1)
  380. a[21] = bc1 ^ (bc3 &^ bc2)
  381. a[22] = bc2 ^ (bc4 &^ bc3)
  382. a[23] = bc3 ^ (bc0 &^ bc4)
  383. a[24] = bc4 ^ (bc1 &^ bc0)
  384. }
  385. }