SecP224R1Field.cs 8.1 KB

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