poly1305-donna32.c 4.7 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190191192193194195196197198199200201202203204205
  1. // SPDX-License-Identifier: GPL-2.0 OR MIT
  2. /*
  3. * Copyright (C) 2015-2019 Jason A. Donenfeld <[email protected]>. All Rights Reserved.
  4. *
  5. * This is based in part on Andrew Moon's poly1305-donna, which is in the
  6. * public domain.
  7. */
  8. #include <linux/kernel.h>
  9. #include <asm/unaligned.h>
  10. #include <crypto/internal/poly1305.h>
  11. void poly1305_core_setkey(struct poly1305_core_key *key,
  12. const u8 raw_key[POLY1305_BLOCK_SIZE])
  13. {
  14. /* r &= 0xffffffc0ffffffc0ffffffc0fffffff */
  15. key->key.r[0] = (get_unaligned_le32(&raw_key[0])) & 0x3ffffff;
  16. key->key.r[1] = (get_unaligned_le32(&raw_key[3]) >> 2) & 0x3ffff03;
  17. key->key.r[2] = (get_unaligned_le32(&raw_key[6]) >> 4) & 0x3ffc0ff;
  18. key->key.r[3] = (get_unaligned_le32(&raw_key[9]) >> 6) & 0x3f03fff;
  19. key->key.r[4] = (get_unaligned_le32(&raw_key[12]) >> 8) & 0x00fffff;
  20. /* s = 5*r */
  21. key->precomputed_s.r[0] = key->key.r[1] * 5;
  22. key->precomputed_s.r[1] = key->key.r[2] * 5;
  23. key->precomputed_s.r[2] = key->key.r[3] * 5;
  24. key->precomputed_s.r[3] = key->key.r[4] * 5;
  25. }
  26. EXPORT_SYMBOL(poly1305_core_setkey);
  27. void poly1305_core_blocks(struct poly1305_state *state,
  28. const struct poly1305_core_key *key, const void *src,
  29. unsigned int nblocks, u32 hibit)
  30. {
  31. const u8 *input = src;
  32. u32 r0, r1, r2, r3, r4;
  33. u32 s1, s2, s3, s4;
  34. u32 h0, h1, h2, h3, h4;
  35. u64 d0, d1, d2, d3, d4;
  36. u32 c;
  37. if (!nblocks)
  38. return;
  39. hibit <<= 24;
  40. r0 = key->key.r[0];
  41. r1 = key->key.r[1];
  42. r2 = key->key.r[2];
  43. r3 = key->key.r[3];
  44. r4 = key->key.r[4];
  45. s1 = key->precomputed_s.r[0];
  46. s2 = key->precomputed_s.r[1];
  47. s3 = key->precomputed_s.r[2];
  48. s4 = key->precomputed_s.r[3];
  49. h0 = state->h[0];
  50. h1 = state->h[1];
  51. h2 = state->h[2];
  52. h3 = state->h[3];
  53. h4 = state->h[4];
  54. do {
  55. /* h += m[i] */
  56. h0 += (get_unaligned_le32(&input[0])) & 0x3ffffff;
  57. h1 += (get_unaligned_le32(&input[3]) >> 2) & 0x3ffffff;
  58. h2 += (get_unaligned_le32(&input[6]) >> 4) & 0x3ffffff;
  59. h3 += (get_unaligned_le32(&input[9]) >> 6) & 0x3ffffff;
  60. h4 += (get_unaligned_le32(&input[12]) >> 8) | hibit;
  61. /* h *= r */
  62. d0 = ((u64)h0 * r0) + ((u64)h1 * s4) +
  63. ((u64)h2 * s3) + ((u64)h3 * s2) +
  64. ((u64)h4 * s1);
  65. d1 = ((u64)h0 * r1) + ((u64)h1 * r0) +
  66. ((u64)h2 * s4) + ((u64)h3 * s3) +
  67. ((u64)h4 * s2);
  68. d2 = ((u64)h0 * r2) + ((u64)h1 * r1) +
  69. ((u64)h2 * r0) + ((u64)h3 * s4) +
  70. ((u64)h4 * s3);
  71. d3 = ((u64)h0 * r3) + ((u64)h1 * r2) +
  72. ((u64)h2 * r1) + ((u64)h3 * r0) +
  73. ((u64)h4 * s4);
  74. d4 = ((u64)h0 * r4) + ((u64)h1 * r3) +
  75. ((u64)h2 * r2) + ((u64)h3 * r1) +
  76. ((u64)h4 * r0);
  77. /* (partial) h %= p */
  78. c = (u32)(d0 >> 26);
  79. h0 = (u32)d0 & 0x3ffffff;
  80. d1 += c;
  81. c = (u32)(d1 >> 26);
  82. h1 = (u32)d1 & 0x3ffffff;
  83. d2 += c;
  84. c = (u32)(d2 >> 26);
  85. h2 = (u32)d2 & 0x3ffffff;
  86. d3 += c;
  87. c = (u32)(d3 >> 26);
  88. h3 = (u32)d3 & 0x3ffffff;
  89. d4 += c;
  90. c = (u32)(d4 >> 26);
  91. h4 = (u32)d4 & 0x3ffffff;
  92. h0 += c * 5;
  93. c = (h0 >> 26);
  94. h0 = h0 & 0x3ffffff;
  95. h1 += c;
  96. input += POLY1305_BLOCK_SIZE;
  97. } while (--nblocks);
  98. state->h[0] = h0;
  99. state->h[1] = h1;
  100. state->h[2] = h2;
  101. state->h[3] = h3;
  102. state->h[4] = h4;
  103. }
  104. EXPORT_SYMBOL(poly1305_core_blocks);
  105. void poly1305_core_emit(const struct poly1305_state *state, const u32 nonce[4],
  106. void *dst)
  107. {
  108. u8 *mac = dst;
  109. u32 h0, h1, h2, h3, h4, c;
  110. u32 g0, g1, g2, g3, g4;
  111. u64 f;
  112. u32 mask;
  113. /* fully carry h */
  114. h0 = state->h[0];
  115. h1 = state->h[1];
  116. h2 = state->h[2];
  117. h3 = state->h[3];
  118. h4 = state->h[4];
  119. c = h1 >> 26;
  120. h1 = h1 & 0x3ffffff;
  121. h2 += c;
  122. c = h2 >> 26;
  123. h2 = h2 & 0x3ffffff;
  124. h3 += c;
  125. c = h3 >> 26;
  126. h3 = h3 & 0x3ffffff;
  127. h4 += c;
  128. c = h4 >> 26;
  129. h4 = h4 & 0x3ffffff;
  130. h0 += c * 5;
  131. c = h0 >> 26;
  132. h0 = h0 & 0x3ffffff;
  133. h1 += c;
  134. /* compute h + -p */
  135. g0 = h0 + 5;
  136. c = g0 >> 26;
  137. g0 &= 0x3ffffff;
  138. g1 = h1 + c;
  139. c = g1 >> 26;
  140. g1 &= 0x3ffffff;
  141. g2 = h2 + c;
  142. c = g2 >> 26;
  143. g2 &= 0x3ffffff;
  144. g3 = h3 + c;
  145. c = g3 >> 26;
  146. g3 &= 0x3ffffff;
  147. g4 = h4 + c - (1UL << 26);
  148. /* select h if h < p, or h + -p if h >= p */
  149. mask = (g4 >> ((sizeof(u32) * 8) - 1)) - 1;
  150. g0 &= mask;
  151. g1 &= mask;
  152. g2 &= mask;
  153. g3 &= mask;
  154. g4 &= mask;
  155. mask = ~mask;
  156. h0 = (h0 & mask) | g0;
  157. h1 = (h1 & mask) | g1;
  158. h2 = (h2 & mask) | g2;
  159. h3 = (h3 & mask) | g3;
  160. h4 = (h4 & mask) | g4;
  161. /* h = h % (2^128) */
  162. h0 = ((h0) | (h1 << 26)) & 0xffffffff;
  163. h1 = ((h1 >> 6) | (h2 << 20)) & 0xffffffff;
  164. h2 = ((h2 >> 12) | (h3 << 14)) & 0xffffffff;
  165. h3 = ((h3 >> 18) | (h4 << 8)) & 0xffffffff;
  166. if (likely(nonce)) {
  167. /* mac = (h + nonce) % (2^128) */
  168. f = (u64)h0 + nonce[0];
  169. h0 = (u32)f;
  170. f = (u64)h1 + nonce[1] + (f >> 32);
  171. h1 = (u32)f;
  172. f = (u64)h2 + nonce[2] + (f >> 32);
  173. h2 = (u32)f;
  174. f = (u64)h3 + nonce[3] + (f >> 32);
  175. h3 = (u32)f;
  176. }
  177. put_unaligned_le32(h0, &mac[0]);
  178. put_unaligned_le32(h1, &mac[4]);
  179. put_unaligned_le32(h2, &mac[8]);
  180. put_unaligned_le32(h3, &mac[12]);
  181. }
  182. EXPORT_SYMBOL(poly1305_core_emit);