Curve25519Field.cs 6.5 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190191192193194195196197198199200201202203204205206207208209210211212213214215216217218219220221222223224225226227228229230231232233234235236237238239240241242243244245246247248249250251252253254255256257
  1. #if !BESTHTTP_DISABLE_ALTERNATE_SSL && (!UNITY_WEBGL || UNITY_EDITOR)
  2. using System;
  3. using System.Diagnostics;
  4. using Org.BouncyCastle.Math.Raw;
  5. namespace Org.BouncyCastle.Math.EC.Custom.Djb
  6. {
  7. internal class Curve25519Field
  8. {
  9. // 2^255 - 2^4 - 2^1 - 1
  10. internal static readonly uint[] P = new uint[]{ 0xFFFFFFED, 0xFFFFFFFF, 0xFFFFFFFF, 0xFFFFFFFF, 0xFFFFFFFF, 0xFFFFFFFF,
  11. 0xFFFFFFFF, 0x7FFFFFFF };
  12. private const uint P7 = 0x7FFFFFFF;
  13. private static readonly uint[] PExt = new uint[]{ 0x00000169, 0x00000000, 0x00000000, 0x00000000, 0x00000000,
  14. 0x00000000, 0x00000000, 0x00000000, 0xFFFFFFED, 0xFFFFFFFF, 0xFFFFFFFF, 0xFFFFFFFF, 0xFFFFFFFF, 0xFFFFFFFF,
  15. 0xFFFFFFFF, 0x3FFFFFFF };
  16. private const uint PInv = 0x13;
  17. public static void Add(uint[] x, uint[] y, uint[] z)
  18. {
  19. Nat256.Add(x, y, z);
  20. if (Nat256.Gte(z, P))
  21. {
  22. SubPFrom(z);
  23. }
  24. }
  25. public static void AddExt(uint[] xx, uint[] yy, uint[] zz)
  26. {
  27. Nat.Add(16, xx, yy, zz);
  28. if (Nat.Gte(16, zz, PExt))
  29. {
  30. SubPExtFrom(zz);
  31. }
  32. }
  33. public static void AddOne(uint[] x, uint[] z)
  34. {
  35. Nat.Inc(8, x, z);
  36. if (Nat256.Gte(z, P))
  37. {
  38. SubPFrom(z);
  39. }
  40. }
  41. public static uint[] FromBigInteger(BigInteger x)
  42. {
  43. uint[] z = Nat256.FromBigInteger(x);
  44. while (Nat256.Gte(z, P))
  45. {
  46. Nat256.SubFrom(P, z);
  47. }
  48. return z;
  49. }
  50. public static void Half(uint[] x, uint[] z)
  51. {
  52. if ((x[0] & 1) == 0)
  53. {
  54. Nat.ShiftDownBit(8, x, 0, z);
  55. }
  56. else
  57. {
  58. Nat256.Add(x, P, z);
  59. Nat.ShiftDownBit(8, z, 0);
  60. }
  61. }
  62. public static void Multiply(uint[] x, uint[] y, uint[] z)
  63. {
  64. uint[] tt = Nat256.CreateExt();
  65. Nat256.Mul(x, y, tt);
  66. Reduce(tt, z);
  67. }
  68. public static void MultiplyAddToExt(uint[] x, uint[] y, uint[] zz)
  69. {
  70. Nat256.MulAddTo(x, y, zz);
  71. if (Nat.Gte(16, zz, PExt))
  72. {
  73. SubPExtFrom(zz);
  74. }
  75. }
  76. public static void Negate(uint[] x, uint[] z)
  77. {
  78. if (Nat256.IsZero(x))
  79. {
  80. Nat256.Zero(z);
  81. }
  82. else
  83. {
  84. Nat256.Sub(P, x, z);
  85. }
  86. }
  87. public static void Reduce(uint[] xx, uint[] z)
  88. {
  89. Debug.Assert(xx[15] >> 30 == 0);
  90. uint xx07 = xx[7];
  91. Nat.ShiftUpBit(8, xx, 8, xx07, z, 0);
  92. uint c = Nat256.MulByWordAddTo(PInv, xx, z) << 1;
  93. uint z7 = z[7];
  94. c += (z7 >> 31) - (xx07 >> 31);
  95. z7 &= P7;
  96. z7 += Nat.AddWordTo(7, c * PInv, z);
  97. z[7] = z7;
  98. if (z7 >= P7 && Nat256.Gte(z, P))
  99. {
  100. SubPFrom(z);
  101. }
  102. }
  103. public static void Reduce27(uint x, uint[] z)
  104. {
  105. Debug.Assert(x >> 26 == 0);
  106. uint z7 = z[7];
  107. uint c = (x << 1 | z7 >> 31);
  108. z7 &= P7;
  109. z7 += Nat.AddWordTo(7, c * PInv, z);
  110. z[7] = z7;
  111. if (z7 >= P7 && Nat256.Gte(z, P))
  112. {
  113. SubPFrom(z);
  114. }
  115. }
  116. public static void Square(uint[] x, uint[] z)
  117. {
  118. uint[] tt = Nat256.CreateExt();
  119. Nat256.Square(x, tt);
  120. Reduce(tt, z);
  121. }
  122. public static void SquareN(uint[] x, int n, uint[] z)
  123. {
  124. Debug.Assert(n > 0);
  125. uint[] tt = Nat256.CreateExt();
  126. Nat256.Square(x, tt);
  127. Reduce(tt, z);
  128. while (--n > 0)
  129. {
  130. Nat256.Square(z, tt);
  131. Reduce(tt, z);
  132. }
  133. }
  134. public static void Subtract(uint[] x, uint[] y, uint[] z)
  135. {
  136. int c = Nat256.Sub(x, y, z);
  137. if (c != 0)
  138. {
  139. AddPTo(z);
  140. }
  141. }
  142. public static void SubtractExt(uint[] xx, uint[] yy, uint[] zz)
  143. {
  144. int c = Nat.Sub(16, xx, yy, zz);
  145. if (c != 0)
  146. {
  147. AddPExtTo(zz);
  148. }
  149. }
  150. public static void Twice(uint[] x, uint[] z)
  151. {
  152. Nat.ShiftUpBit(8, x, 0, z);
  153. if (Nat256.Gte(z, P))
  154. {
  155. SubPFrom(z);
  156. }
  157. }
  158. private static uint AddPTo(uint[] z)
  159. {
  160. long c = (long)z[0] - PInv;
  161. z[0] = (uint)c;
  162. c >>= 32;
  163. if (c != 0)
  164. {
  165. c = Nat.DecAt(7, z, 1);
  166. }
  167. c += (long)z[7] + (P7 + 1);
  168. z[7] = (uint)c;
  169. c >>= 32;
  170. return (uint)c;
  171. }
  172. private static uint AddPExtTo(uint[] zz)
  173. {
  174. long c = (long)zz[0] + PExt[0];
  175. zz[0] = (uint)c;
  176. c >>= 32;
  177. if (c != 0)
  178. {
  179. c = Nat.IncAt(8, zz, 1);
  180. }
  181. c += (long)zz[8] - PInv;
  182. zz[8] = (uint)c;
  183. c >>= 32;
  184. if (c != 0)
  185. {
  186. c = Nat.DecAt(15, zz, 9);
  187. }
  188. c += (long)zz[15] + (PExt[15] + 1);
  189. zz[15] = (uint)c;
  190. c >>= 32;
  191. return (uint)c;
  192. }
  193. private static int SubPFrom(uint[] z)
  194. {
  195. long c = (long)z[0] + PInv;
  196. z[0] = (uint)c;
  197. c >>= 32;
  198. if (c != 0)
  199. {
  200. c = Nat.IncAt(7, z, 1);
  201. }
  202. c += (long)z[7] - (P7 + 1);
  203. z[7] = (uint)c;
  204. c >>= 32;
  205. return (int)c;
  206. }
  207. private static int SubPExtFrom(uint[] zz)
  208. {
  209. long c = (long)zz[0] - PExt[0];
  210. zz[0] = (uint)c;
  211. c >>= 32;
  212. if (c != 0)
  213. {
  214. c = Nat.DecAt(8, zz, 1);
  215. }
  216. c += (long)zz[8] + PInv;
  217. zz[8] = (uint)c;
  218. c >>= 32;
  219. if (c != 0)
  220. {
  221. c = Nat.IncAt(15, zz, 9);
  222. }
  223. c += (long)zz[15] - (PExt[15] + 1);
  224. zz[15] = (uint)c;
  225. c >>= 32;
  226. return (int)c;
  227. }
  228. }
  229. }
  230. #endif