SecP256R1Field.cs 8.6 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190191192193194195196197198199200201202203204205206207208209210211212213214215216217218219220221222223224225226227228229230231232233234235236237238239240241242243244245246247248249250251252253254255256257258259260261262263264265266267268269270271272273274275276277278279280281282283284285286287288289290291292293294295296297298299300301302303304305306307308309310311312313314315316
  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.Sec
  7. {
  8. internal class SecP256R1Field
  9. {
  10. // 2^256 - 2^224 + 2^192 + 2^96 - 1
  11. internal static readonly uint[] P = new uint[]{ 0xFFFFFFFF, 0xFFFFFFFF, 0xFFFFFFFF, 0x00000000, 0x00000000, 0x00000000,
  12. 0x00000001, 0xFFFFFFFF };
  13. internal static readonly uint[] PExt = new uint[]{ 0x00000001, 0x00000000, 0x00000000, 0xFFFFFFFE, 0xFFFFFFFF,
  14. 0xFFFFFFFF, 0xFFFFFFFE, 0x00000001, 0xFFFFFFFE, 0x00000001, 0xFFFFFFFE, 0x00000001, 0x00000001, 0xFFFFFFFE,
  15. 0x00000002, 0xFFFFFFFE };
  16. internal const uint P7 = 0xFFFFFFFF;
  17. internal const uint PExt15 = 0xFFFFFFFE;
  18. public static void Add(uint[] x, uint[] y, uint[] z)
  19. {
  20. uint c = Nat256.Add(x, y, z);
  21. if (c != 0 || (z[7] == P7 && Nat256.Gte(z, P)))
  22. {
  23. AddPInvTo(z);
  24. }
  25. }
  26. public static void AddExt(uint[] xx, uint[] yy, uint[] zz)
  27. {
  28. uint c = Nat.Add(16, xx, yy, zz);
  29. if (c != 0 || (zz[15] >= PExt15 && Nat.Gte(16, zz, PExt)))
  30. {
  31. Nat.SubFrom(16, PExt, zz);
  32. }
  33. }
  34. public static void AddOne(uint[] x, uint[] z)
  35. {
  36. uint c = Nat.Inc(8, x, z);
  37. if (c != 0 || (z[7] == P7 && Nat256.Gte(z, P)))
  38. {
  39. AddPInvTo(z);
  40. }
  41. }
  42. public static uint[] FromBigInteger(BigInteger x)
  43. {
  44. uint[] z = Nat256.FromBigInteger(x);
  45. if (z[7] == P7 && 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. uint c = Nat256.Add(x, P, z);
  60. Nat.ShiftDownBit(8, z, c);
  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. uint c = Nat256.MulAddTo(x, y, zz);
  72. if (c != 0 || (zz[15] >= PExt15 && Nat.Gte(16, zz, PExt)))
  73. {
  74. Nat.SubFrom(16, PExt, 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. long xx08 = xx[8], xx09 = xx[9], xx10 = xx[10], xx11 = xx[11];
  91. long xx12 = xx[12], xx13 = xx[13], xx14 = xx[14], xx15 = xx[15];
  92. const long n = 6;
  93. xx08 -= n;
  94. long t0 = xx08 + xx09;
  95. long t1 = xx09 + xx10;
  96. long t2 = xx10 + xx11 - xx15;
  97. long t3 = xx11 + xx12;
  98. long t4 = xx12 + xx13;
  99. long t5 = xx13 + xx14;
  100. long t6 = xx14 + xx15;
  101. long t7 = t5 - t0;
  102. long cc = 0;
  103. cc += (long)xx[0] - t3 - t7;
  104. z[0] = (uint)cc;
  105. cc >>= 32;
  106. cc += (long)xx[1] + t1 - t4 - t6;
  107. z[1] = (uint)cc;
  108. cc >>= 32;
  109. cc += (long)xx[2] + t2 - t5;
  110. z[2] = (uint)cc;
  111. cc >>= 32;
  112. cc += (long)xx[3] + (t3 << 1) + t7 - t6;
  113. z[3] = (uint)cc;
  114. cc >>= 32;
  115. cc += (long)xx[4] + (t4 << 1) + xx14 - t1;
  116. z[4] = (uint)cc;
  117. cc >>= 32;
  118. cc += (long)xx[5] + (t5 << 1) - t2;
  119. z[5] = (uint)cc;
  120. cc >>= 32;
  121. cc += (long)xx[6] + (t6 << 1) + t7;
  122. z[6] = (uint)cc;
  123. cc >>= 32;
  124. cc += (long)xx[7] + (xx15 << 1) + xx08 - t2 - t4;
  125. z[7] = (uint)cc;
  126. cc >>= 32;
  127. cc += n;
  128. Debug.Assert(cc >= 0);
  129. Reduce32((uint)cc, z);
  130. }
  131. public static void Reduce32(uint x, uint[] z)
  132. {
  133. long cc = 0;
  134. if (x != 0)
  135. {
  136. long xx08 = x;
  137. cc += (long)z[0] + xx08;
  138. z[0] = (uint)cc;
  139. cc >>= 32;
  140. if (cc != 0)
  141. {
  142. cc += (long)z[1];
  143. z[1] = (uint)cc;
  144. cc >>= 32;
  145. cc += (long)z[2];
  146. z[2] = (uint)cc;
  147. cc >>= 32;
  148. }
  149. cc += (long)z[3] - xx08;
  150. z[3] = (uint)cc;
  151. cc >>= 32;
  152. if (cc != 0)
  153. {
  154. cc += (long)z[4];
  155. z[4] = (uint)cc;
  156. cc >>= 32;
  157. cc += (long)z[5];
  158. z[5] = (uint)cc;
  159. cc >>= 32;
  160. }
  161. cc += (long)z[6] - xx08;
  162. z[6] = (uint)cc;
  163. cc >>= 32;
  164. cc += (long)z[7] + xx08;
  165. z[7] = (uint)cc;
  166. cc >>= 32;
  167. Debug.Assert(cc == 0 || cc == 1);
  168. }
  169. if (cc != 0 || (z[7] == P7 && Nat256.Gte(z, P)))
  170. {
  171. AddPInvTo(z);
  172. }
  173. }
  174. public static void Square(uint[] x, uint[] z)
  175. {
  176. uint[] tt = Nat256.CreateExt();
  177. Nat256.Square(x, tt);
  178. Reduce(tt, z);
  179. }
  180. public static void SquareN(uint[] x, int n, uint[] z)
  181. {
  182. Debug.Assert(n > 0);
  183. uint[] tt = Nat256.CreateExt();
  184. Nat256.Square(x, tt);
  185. Reduce(tt, z);
  186. while (--n > 0)
  187. {
  188. Nat256.Square(z, tt);
  189. Reduce(tt, z);
  190. }
  191. }
  192. public static void Subtract(uint[] x, uint[] y, uint[] z)
  193. {
  194. int c = Nat256.Sub(x, y, z);
  195. if (c != 0)
  196. {
  197. SubPInvFrom(z);
  198. }
  199. }
  200. public static void SubtractExt(uint[] xx, uint[] yy, uint[] zz)
  201. {
  202. int c = Nat.Sub(16, xx, yy, zz);
  203. if (c != 0)
  204. {
  205. Nat.AddTo(16, PExt, zz);
  206. }
  207. }
  208. public static void Twice(uint[] x, uint[] z)
  209. {
  210. uint c = Nat.ShiftUpBit(8, x, 0, z);
  211. if (c != 0 || (z[7] == P7 && Nat256.Gte(z, P)))
  212. {
  213. AddPInvTo(z);
  214. }
  215. }
  216. private static void AddPInvTo(uint[] z)
  217. {
  218. long c = (long)z[0] + 1;
  219. z[0] = (uint)c;
  220. c >>= 32;
  221. if (c != 0)
  222. {
  223. c += (long)z[1];
  224. z[1] = (uint)c;
  225. c >>= 32;
  226. c += (long)z[2];
  227. z[2] = (uint)c;
  228. c >>= 32;
  229. }
  230. c += (long)z[3] - 1;
  231. z[3] = (uint)c;
  232. c >>= 32;
  233. if (c != 0)
  234. {
  235. c += (long)z[4];
  236. z[4] = (uint)c;
  237. c >>= 32;
  238. c += (long)z[5];
  239. z[5] = (uint)c;
  240. c >>= 32;
  241. }
  242. c += (long)z[6] - 1;
  243. z[6] = (uint)c;
  244. c >>= 32;
  245. c += (long)z[7] + 1;
  246. z[7] = (uint)c;
  247. //c >>= 32;
  248. }
  249. private static void SubPInvFrom(uint[] z)
  250. {
  251. long c = (long)z[0] - 1;
  252. z[0] = (uint)c;
  253. c >>= 32;
  254. if (c != 0)
  255. {
  256. c += (long)z[1];
  257. z[1] = (uint)c;
  258. c >>= 32;
  259. c += (long)z[2];
  260. z[2] = (uint)c;
  261. c >>= 32;
  262. }
  263. c += (long)z[3] + 1;
  264. z[3] = (uint)c;
  265. c >>= 32;
  266. if (c != 0)
  267. {
  268. c += (long)z[4];
  269. z[4] = (uint)c;
  270. c >>= 32;
  271. c += (long)z[5];
  272. z[5] = (uint)c;
  273. c >>= 32;
  274. }
  275. c += (long)z[6] + 1;
  276. z[6] = (uint)c;
  277. c >>= 32;
  278. c += (long)z[7] - 1;
  279. z[7] = (uint)c;
  280. //c >>= 32;
  281. }
  282. }
  283. }
  284. #pragma warning restore
  285. #endif