Curve25519Field.cs 6.6 KB

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