SecP521R1Field.cs 4.4 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159
  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 SecP521R1Field
  9. {
  10. // 2^521 - 1
  11. internal static readonly uint[] P = new uint[]{ 0xFFFFFFFF, 0xFFFFFFFF, 0xFFFFFFFF, 0xFFFFFFFF, 0xFFFFFFFF, 0xFFFFFFFF, 0xFFFFFFFF,
  12. 0xFFFFFFFF, 0xFFFFFFFF, 0xFFFFFFFF, 0xFFFFFFFF, 0xFFFFFFFF, 0xFFFFFFFF, 0xFFFFFFFF, 0xFFFFFFFF, 0xFFFFFFFF, 0x1FF };
  13. private const int P16 = 0x1FF;
  14. public static void Add(uint[] x, uint[] y, uint[] z)
  15. {
  16. uint c = Nat.Add(16, x, y, z) + x[16] + y[16];
  17. if (c > P16 || (c == P16 && Nat.Eq(16, z, P)))
  18. {
  19. c += Nat.Inc(16, z);
  20. c &= P16;
  21. }
  22. z[16] = c;
  23. }
  24. public static void AddOne(uint[] x, uint[] z)
  25. {
  26. uint c = Nat.Inc(16, x, z) + x[16];
  27. if (c > P16 || (c == P16 && Nat.Eq(16, z, P)))
  28. {
  29. c += Nat.Inc(16, z);
  30. c &= P16;
  31. }
  32. z[16] = c;
  33. }
  34. public static uint[] FromBigInteger(BigInteger x)
  35. {
  36. uint[] z = Nat.FromBigInteger(521, x);
  37. if (Nat.Eq(17, z, P))
  38. {
  39. Nat.Zero(17, z);
  40. }
  41. return z;
  42. }
  43. public static void Half(uint[] x, uint[] z)
  44. {
  45. uint x16 = x[16];
  46. uint c = Nat.ShiftDownBit(16, x, x16, z);
  47. z[16] = (x16 >> 1) | (c >> 23);
  48. }
  49. public static void Multiply(uint[] x, uint[] y, uint[] z)
  50. {
  51. uint[] tt = Nat.Create(33);
  52. ImplMultiply(x, y, tt);
  53. Reduce(tt, z);
  54. }
  55. public static void Negate(uint[] x, uint[] z)
  56. {
  57. if (Nat.IsZero(17, x))
  58. {
  59. Nat.Zero(17, z);
  60. }
  61. else
  62. {
  63. Nat.Sub(17, P, x, z);
  64. }
  65. }
  66. public static void Reduce(uint[] xx, uint[] z)
  67. {
  68. Debug.Assert(xx[32] >> 18 == 0);
  69. uint xx32 = xx[32];
  70. uint c = Nat.ShiftDownBits(16, xx, 16, 9, xx32, z, 0) >> 23;
  71. c += xx32 >> 9;
  72. c += Nat.AddTo(16, xx, z);
  73. if (c > P16 || (c == P16 && Nat.Eq(16, z, P)))
  74. {
  75. c += Nat.Inc(16, z);
  76. c &= P16;
  77. }
  78. z[16] = c;
  79. }
  80. public static void Reduce23(uint[] z)
  81. {
  82. uint z16 = z[16];
  83. uint c = Nat.AddWordTo(16, z16 >> 9, z) + (z16 & P16);
  84. if (c > P16 || (c == P16 && Nat.Eq(16, z, P)))
  85. {
  86. c += Nat.Inc(16, z);
  87. c &= P16;
  88. }
  89. z[16] = c;
  90. }
  91. public static void Square(uint[] x, uint[] z)
  92. {
  93. uint[] tt = Nat.Create(33);
  94. ImplSquare(x, tt);
  95. Reduce(tt, z);
  96. }
  97. public static void SquareN(uint[] x, int n, uint[] z)
  98. {
  99. Debug.Assert(n > 0);
  100. uint[] tt = Nat.Create(33);
  101. ImplSquare(x, tt);
  102. Reduce(tt, z);
  103. while (--n > 0)
  104. {
  105. ImplSquare(z, tt);
  106. Reduce(tt, z);
  107. }
  108. }
  109. public static void Subtract(uint[] x, uint[] y, uint[] z)
  110. {
  111. int c = Nat.Sub(16, x, y, z) + (int)(x[16] - y[16]);
  112. if (c < 0)
  113. {
  114. c += Nat.Dec(16, z);
  115. c &= P16;
  116. }
  117. z[16] = (uint)c;
  118. }
  119. public static void Twice(uint[] x, uint[] z)
  120. {
  121. uint x16 = x[16];
  122. uint c = Nat.ShiftUpBit(16, x, x16 << 23, z) | (x16 << 1);
  123. z[16] = c & P16;
  124. }
  125. protected static void ImplMultiply(uint[] x, uint[] y, uint[] zz)
  126. {
  127. Nat512.Mul(x, y, zz);
  128. uint x16 = x[16], y16 = y[16];
  129. zz[32] = Nat.Mul31BothAdd(16, x16, y, y16, x, zz, 16) + (x16 * y16);
  130. }
  131. protected static void ImplSquare(uint[] x, uint[] zz)
  132. {
  133. Nat512.Square(x, zz);
  134. uint x16 = x[16];
  135. zz[32] = Nat.MulWordAddTo(16, x16 << 1, x, 0, zz, 16) + (x16 * x16);
  136. }
  137. }
  138. }
  139. #pragma warning restore
  140. #endif