SecP224R1Field.cs 9.1 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190191192193194195196197198199200201202203204205206207208209210211212213214215216217218219220221222223224225226227228229230231232233234235236237238239240241242243244245246247248249250251252253254255256257258259260261262263264265266267268269270271272273274275276277278279280281282283284285286287288289290291292293294295296297298299300301302303304305306307308309310311312313314315316317318319320321322323324325326327328329330331332333334335336337338339340341
  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.Crypto.Utilities;
  6. using BestHTTP.SecureProtocol.Org.BouncyCastle.Math.Raw;
  7. using BestHTTP.SecureProtocol.Org.BouncyCastle.Security;
  8. namespace BestHTTP.SecureProtocol.Org.BouncyCastle.Math.EC.Custom.Sec
  9. {
  10. internal class SecP224R1Field
  11. {
  12. // 2^224 - 2^96 + 1
  13. internal static readonly uint[] P = new uint[]{ 0x00000001, 0x00000000, 0x00000000, 0xFFFFFFFF, 0xFFFFFFFF,
  14. 0xFFFFFFFF, 0xFFFFFFFF };
  15. private static readonly uint[] PExt = new uint[]{ 0x00000001, 0x00000000, 0x00000000, 0xFFFFFFFE, 0xFFFFFFFF,
  16. 0xFFFFFFFF, 0x00000000, 0x00000002, 0x00000000, 0x00000000, 0xFFFFFFFE, 0xFFFFFFFF, 0xFFFFFFFF, 0xFFFFFFFF
  17. };
  18. private static readonly uint[] PExtInv = new uint[]{ 0xFFFFFFFF, 0xFFFFFFFF, 0xFFFFFFFF, 0x00000001,
  19. 0x00000000, 0x00000000, 0xFFFFFFFF, 0xFFFFFFFD, 0xFFFFFFFF, 0xFFFFFFFF, 0x00000001 };
  20. private const uint P6 = 0xFFFFFFFF;
  21. private const uint PExt13 = 0xFFFFFFFF;
  22. public static void Add(uint[] x, uint[] y, uint[] z)
  23. {
  24. uint c = Nat224.Add(x, y, z);
  25. if (c != 0 || (z[6] == P6 && Nat224.Gte(z, P)))
  26. {
  27. AddPInvTo(z);
  28. }
  29. }
  30. public static void AddExt(uint[] xx, uint[] yy, uint[] zz)
  31. {
  32. uint c = Nat.Add(14, xx, yy, zz);
  33. if (c != 0 || (zz[13] == PExt13 && Nat.Gte(14, zz, PExt)))
  34. {
  35. if (Nat.AddTo(PExtInv.Length, PExtInv, zz) != 0)
  36. {
  37. Nat.IncAt(14, zz, PExtInv.Length);
  38. }
  39. }
  40. }
  41. public static void AddOne(uint[] x, uint[] z)
  42. {
  43. uint c = Nat.Inc(7, x, z);
  44. if (c != 0 || (z[6] == P6 && Nat224.Gte(z, P)))
  45. {
  46. AddPInvTo(z);
  47. }
  48. }
  49. public static uint[] FromBigInteger(BigInteger x)
  50. {
  51. uint[] z = Nat.FromBigInteger(224, x);
  52. if (z[6] == P6 && Nat224.Gte(z, P))
  53. {
  54. Nat224.SubFrom(P, z);
  55. }
  56. return z;
  57. }
  58. public static void Half(uint[] x, uint[] z)
  59. {
  60. if ((x[0] & 1) == 0)
  61. {
  62. Nat.ShiftDownBit(7, x, 0, z);
  63. }
  64. else
  65. {
  66. uint c = Nat224.Add(x, P, z);
  67. Nat.ShiftDownBit(7, z, c);
  68. }
  69. }
  70. public static void Inv(uint[] x, uint[] z)
  71. {
  72. Mod.CheckedModOddInverse(P, x, z);
  73. }
  74. public static int IsZero(uint[] x)
  75. {
  76. uint d = 0;
  77. for (int i = 0; i < 7; ++i)
  78. {
  79. d |= x[i];
  80. }
  81. d = (d >> 1) | (d & 1);
  82. return ((int)d - 1) >> 31;
  83. }
  84. public static void Multiply(uint[] x, uint[] y, uint[] z)
  85. {
  86. uint[] tt = Nat224.CreateExt();
  87. Nat224.Mul(x, y, tt);
  88. Reduce(tt, z);
  89. }
  90. public static void MultiplyAddToExt(uint[] x, uint[] y, uint[] zz)
  91. {
  92. uint c = Nat224.MulAddTo(x, y, zz);
  93. if (c != 0 || (zz[13] == PExt13 && Nat.Gte(14, zz, PExt)))
  94. {
  95. if (Nat.AddTo(PExtInv.Length, PExtInv, zz) != 0)
  96. {
  97. Nat.IncAt(14, zz, PExtInv.Length);
  98. }
  99. }
  100. }
  101. public static void Negate(uint[] x, uint[] z)
  102. {
  103. if (0 != IsZero(x))
  104. {
  105. Nat224.Sub(P, P, z);
  106. }
  107. else
  108. {
  109. Nat224.Sub(P, x, z);
  110. }
  111. }
  112. public static void Random(SecureRandom r, uint[] z)
  113. {
  114. byte[] bb = new byte[7 * 4];
  115. do
  116. {
  117. r.NextBytes(bb);
  118. Pack.LE_To_UInt32(bb, 0, z, 0, 7);
  119. }
  120. while (0 == Nat.LessThan(7, z, P));
  121. }
  122. public static void RandomMult(SecureRandom r, uint[] z)
  123. {
  124. do
  125. {
  126. Random(r, z);
  127. }
  128. while (0 != IsZero(z));
  129. }
  130. public static void Reduce(uint[] xx, uint[] z)
  131. {
  132. long xx10 = xx[10], xx11 = xx[11], xx12 = xx[12], xx13 = xx[13];
  133. const long n = 1;
  134. long t0 = (long)xx[7] + xx11 - n;
  135. long t1 = (long)xx[8] + xx12;
  136. long t2 = (long)xx[9] + xx13;
  137. long cc = 0;
  138. cc += (long)xx[0] - t0;
  139. long z0 = (uint)cc;
  140. cc >>= 32;
  141. cc += (long)xx[1] - t1;
  142. z[1] = (uint)cc;
  143. cc >>= 32;
  144. cc += (long)xx[2] - t2;
  145. z[2] = (uint)cc;
  146. cc >>= 32;
  147. cc += (long)xx[3] + t0 - xx10;
  148. long z3 = (uint)cc;
  149. cc >>= 32;
  150. cc += (long)xx[4] + t1 - xx11;
  151. z[4] = (uint)cc;
  152. cc >>= 32;
  153. cc += (long)xx[5] + t2 - xx12;
  154. z[5] = (uint)cc;
  155. cc >>= 32;
  156. cc += (long)xx[6] + xx10 - xx13;
  157. z[6] = (uint)cc;
  158. cc >>= 32;
  159. cc += n;
  160. Debug.Assert(cc >= 0);
  161. z3 += cc;
  162. z0 -= cc;
  163. z[0] = (uint)z0;
  164. cc = z0 >> 32;
  165. if (cc != 0)
  166. {
  167. cc += (long)z[1];
  168. z[1] = (uint)cc;
  169. cc >>= 32;
  170. cc += (long)z[2];
  171. z[2] = (uint)cc;
  172. z3 += cc >> 32;
  173. }
  174. z[3] = (uint)z3;
  175. cc = z3 >> 32;
  176. Debug.Assert(cc == 0 || cc == 1);
  177. if ((cc != 0 && Nat.IncAt(7, z, 4) != 0)
  178. || (z[6] == P6 && Nat224.Gte(z, P)))
  179. {
  180. AddPInvTo(z);
  181. }
  182. }
  183. public static void Reduce32(uint x, uint[] z)
  184. {
  185. long cc = 0;
  186. if (x != 0)
  187. {
  188. long xx07 = x;
  189. cc += (long)z[0] - xx07;
  190. z[0] = (uint)cc;
  191. cc >>= 32;
  192. if (cc != 0)
  193. {
  194. cc += (long)z[1];
  195. z[1] = (uint)cc;
  196. cc >>= 32;
  197. cc += (long)z[2];
  198. z[2] = (uint)cc;
  199. cc >>= 32;
  200. }
  201. cc += (long)z[3] + xx07;
  202. z[3] = (uint)cc;
  203. cc >>= 32;
  204. Debug.Assert(cc == 0 || cc == 1);
  205. }
  206. if ((cc != 0 && Nat.IncAt(7, z, 4) != 0)
  207. || (z[6] == P6 && Nat224.Gte(z, P)))
  208. {
  209. AddPInvTo(z);
  210. }
  211. }
  212. public static void Square(uint[] x, uint[] z)
  213. {
  214. uint[] tt = Nat224.CreateExt();
  215. Nat224.Square(x, tt);
  216. Reduce(tt, z);
  217. }
  218. public static void SquareN(uint[] x, int n, uint[] z)
  219. {
  220. Debug.Assert(n > 0);
  221. uint[] tt = Nat224.CreateExt();
  222. Nat224.Square(x, tt);
  223. Reduce(tt, z);
  224. while (--n > 0)
  225. {
  226. Nat224.Square(z, tt);
  227. Reduce(tt, z);
  228. }
  229. }
  230. public static void Subtract(uint[] x, uint[] y, uint[] z)
  231. {
  232. int c = Nat224.Sub(x, y, z);
  233. if (c != 0)
  234. {
  235. SubPInvFrom(z);
  236. }
  237. }
  238. public static void SubtractExt(uint[] xx, uint[] yy, uint[] zz)
  239. {
  240. int c = Nat.Sub(14, xx, yy, zz);
  241. if (c != 0)
  242. {
  243. if (Nat.SubFrom(PExtInv.Length, PExtInv, zz) != 0)
  244. {
  245. Nat.DecAt(14, zz, PExtInv.Length);
  246. }
  247. }
  248. }
  249. public static void Twice(uint[] x, uint[] z)
  250. {
  251. uint c = Nat.ShiftUpBit(7, x, 0, z);
  252. if (c != 0 || (z[6] == P6 && Nat224.Gte(z, P)))
  253. {
  254. AddPInvTo(z);
  255. }
  256. }
  257. private static void AddPInvTo(uint[] z)
  258. {
  259. long c = (long)z[0] - 1;
  260. z[0] = (uint)c;
  261. c >>= 32;
  262. if (c != 0)
  263. {
  264. c += (long)z[1];
  265. z[1] = (uint)c;
  266. c >>= 32;
  267. c += (long)z[2];
  268. z[2] = (uint)c;
  269. c >>= 32;
  270. }
  271. c += (long)z[3] + 1;
  272. z[3] = (uint)c;
  273. c >>= 32;
  274. if (c != 0)
  275. {
  276. Nat.IncAt(7, z, 4);
  277. }
  278. }
  279. private static void SubPInvFrom(uint[] z)
  280. {
  281. long c = (long)z[0] + 1;
  282. z[0] = (uint)c;
  283. c >>= 32;
  284. if (c != 0)
  285. {
  286. c += (long)z[1];
  287. z[1] = (uint)c;
  288. c >>= 32;
  289. c += (long)z[2];
  290. z[2] = (uint)c;
  291. c >>= 32;
  292. }
  293. c += (long)z[3] - 1;
  294. z[3] = (uint)c;
  295. c >>= 32;
  296. if (c != 0)
  297. {
  298. Nat.DecAt(7, z, 4);
  299. }
  300. }
  301. }
  302. }
  303. #pragma warning restore
  304. #endif