PrivateKey.cpp 8.6 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190191192193194195196197198199200201202203204205206207208209210211212213214215216217218219220221222223224225226227228229230231232233234235236237238239240241242243244245246247248249250251252253254255256257258259260261262263264265266267268269270271272273274275276277278279280281282283284285286287288289290291292293294295296297298299300301302303304305306307308309310311312313314315316317318319320321322323
  1. #include "PrivateKey.hpp"
  2. extern const curvepoint_fp_t bn_curvegen;
  3. PrivateKey::PrivateKey(const Fp& a, const Fp& b, const Fp& c, const Fp& d, const Fp& e, const Fp& f, const Fp& g, const Fp& h)
  4. {
  5. set(a, b, c, d, e, f, g, h);
  6. }
  7. void PrivateKey::set(const Fp& a, const Fp& b, const Fp& c, const Fp& d, const Fp& e, const Fp& f, const Fp& g, const Fp& h)
  8. {
  9. i1 = a;
  10. j1 = b;
  11. k1 = c;
  12. l1 = d;
  13. i2 = e;
  14. j2 = f;
  15. k2 = g;
  16. l2 = h;
  17. }
  18. void PrivateKey::set(string which, const Fp& input);
  19. {
  20. if (which.length() != 2)
  21. return;
  22. bool var_is_1 = false;
  23. switch (which[1])
  24. {
  25. case '1':
  26. var_is_1 = true;
  27. break;
  28. case '2':
  29. break;
  30. default:
  31. return;
  32. }
  33. switch (which[0])
  34. {
  35. case 'i':
  36. if (var_is_1)
  37. i1 = input;
  38. else i2 = input;
  39. break;
  40. case 'j':
  41. if (var_is_1)
  42. j1 = input;
  43. else j2 = input;
  44. break;
  45. case 'k':
  46. if (var_is_1)
  47. k1 = input;
  48. else k2 = input;
  49. break;
  50. case 'l':
  51. if (var_is_1)
  52. l1 = input;
  53. else l2 = input;
  54. break;
  55. default:
  56. break;
  57. }
  58. }
  59. Fp PrivateKey::get(string which) const
  60. {
  61. if (which.length() != 2)
  62. return Fp();
  63. bool var_is_1 = false;
  64. switch (which[1])
  65. {
  66. case '1':
  67. var_is_1 = true;
  68. break;
  69. case '2':
  70. break;
  71. default:
  72. return Fp();
  73. }
  74. switch (which[0])
  75. {
  76. case 'i':
  77. return var_is_1 ? i1 : i2;
  78. case 'j':
  79. return var_is_1 ? j1 : j2;
  80. case 'k':
  81. return var_is_1 ? k1 : k2;
  82. case 'l':
  83. return var_is_1 ? l1 : l2;
  84. default:
  85. return Fp();
  86. }
  87. }
  88. Bipoint<curvepoint_fp_t> PrivateKey::pi_1(const Bipoint<curvepoint_fp_t>& input) const
  89. {
  90. Bipoint<curvepoint_fp_t> retval;
  91. curvepoint_fp_t temp1, temp2;
  92. const scalar_t i1_s = i1.to_scalar();
  93. const scalar_t j1_s = j1.to_scalar();
  94. const scalar_t k1_s = k1.to_scalar();
  95. const scalar_t l1_s = l1.to_scalar();
  96. curvepoint_fp_scalarmult_vartime(temp1, input[0], j1_s);
  97. curvepoint_fp_scalarmult_vartime(temp1, temp1, k1_s);
  98. curvepoint_fp_neg(temp1, temp1);
  99. curvepoint_fp_makeaffine(temp1);
  100. curvepoint_fp_scalarmult_vartime(temp2, input[1], i1_s);
  101. curvepoint_fp_scalarmult_vartime(temp2, temp2, k1_s);
  102. curvepoint_fp_makeaffine(temp2);
  103. curvepoint_fp_add_vartime(temp1, temp1, temp2);
  104. curvepoint_fp_makeaffine(temp1);
  105. curvepoint_fp_set(retval[0], temp1);
  106. curvepoint_fp_scalarmult_vartime(temp1, input[0], j1_s);
  107. curvepoint_fp_scalarmult_vartime(temp1, temp1, l1_s);
  108. curvepoint_fp_neg(temp1, temp1);
  109. curvepoint_fp_makeaffine(temp1);
  110. curvepoint_fp_scalarmult_vartime(temp2, input[1], i1_s);
  111. curvepoint_fp_scalarmult_vartime(temp2, temp2, l1_s);
  112. curvepoint_fp_makeaffine(temp2);
  113. curvepoint_fp_add_vartime(temp1, temp1, temp2);
  114. curvepoint_fp_makeaffine(temp1);
  115. curvepoint_fp_set(retval[1], temp1);
  116. return retval;
  117. }
  118. Bipoint<twistpoint_fp2_t> PrivateKey::pi_2(const Bipoint<twistpoint_fp2_t>& input) const
  119. {
  120. Bipoint<twistpoint_fp2_t> retval;
  121. twistpoint_fp2_t temp1, temp2;
  122. const scalar_t i2_s = i2.to_scalar();
  123. const scalar_t j2_s = j2.to_scalar();
  124. const scalar_t k2_s = k2.to_scalar();
  125. const scalar_t l2_s = l2.to_scalar();
  126. twistpoint_fp2_scalarmult_vartime(temp1, input[0], j2_s);
  127. twistpoint_fp2_scalarmult_vartime(temp1, temp1, k2_s);
  128. twistpoint_fp2_neg(temp1,temp1);
  129. twistpoint_fp2_makeaffine(temp1);
  130. twistpoint_fp2_scalarmult_vartime(temp2, input[1], i2_s);
  131. twistpoint_fp2_scalarmult_vartime(temp2, temp2, k2_s);
  132. twistpoint_fp2_makeaffine(temp2);
  133. twistpoint_fp2_add_vartime(temp1, temp1, temp2);
  134. twistpoint_fp2_makeaffine(temp1);
  135. twistpoint_fp2_set(retval[0], temp1);
  136. twistpoint_fp2_scalarmult_vartime(temp1, input[0], j2_s);
  137. twistpoint_fp2_scalarmult_vartime(temp1, temp1, l2_s);
  138. twistpoint_fp2_neg(temp1, temp1);
  139. twistpoint_fp2_makeaffine(temp1);
  140. twistpoint_fp2_scalarmult_vartime(temp2, input[1], i2_s);
  141. twistpoint_fp2_scalarmult_vartime(temp2, temp2, l2_s);
  142. twistpoint_fp2_makeaffine(temp2);
  143. twistpoint_fp2_add_vartime(temp1, temp1, temp2);
  144. twistpoint_fp2_makeaffine(temp1);
  145. twistpoint_fp2_set(retval[1], temp1);
  146. return retval;
  147. }
  148. Quadripoint PrivateKey::pi_T(const Quadripoint& input) const
  149. {
  150. Quadripoint retval;
  151. fp12e_t temp1, temp2, temp3, temp4;
  152. const scalar_t i1_s = i1.to_scalar();
  153. const scalar_t j1_s = j1.to_scalar();
  154. const scalar_t k1_s = k1.to_scalar();
  155. const scalar_t l1_s = l1.to_scalar();
  156. const scalar_t i2_s = i2.to_scalar();
  157. const scalar_t j2_s = j2.to_scalar();
  158. const scalar_t k2_s = k2.to_scalar();
  159. const scalar_t l2_s = l2.to_scalar();
  160. fp12e_pow_vartime(temp1, input[0], j1_s);
  161. fp12e_pow_vartime(temp1, temp1, k1_s);
  162. fp12e_pow_vartime(temp1, temp1, j2_s);
  163. fp12e_pow_vartime(temp1, temp1, k2_s);
  164. fp12e_invert(temp2, input[1]);
  165. fp12e_pow_vartime(temp2, temp2, j1_s);
  166. fp12e_pow_vartime(temp2, temp2, k1_s);
  167. fp12e_pow_vartime(temp2, temp2, i2_s);
  168. fp12e_pow_vartime(temp2, temp2, k2_s);
  169. fp12e_invert(temp3, input[2]);
  170. fp12e_pow_vartime(temp3, temp3, i1_s);
  171. fp12e_pow_vartime(temp3, temp3, k1_s);
  172. fp12e_pow_vartime(temp3, temp3, j2_s);
  173. fp12e_pow_vartime(temp3, temp3, k2_s);
  174. fp12e_pow_vartime(temp4, input[3], i1_s);
  175. fp12e_pow_vartime(temp4, temp4, k1_s);
  176. fp12e_pow_vartime(temp4, temp4, i2_s);
  177. fp12e_pow_vartime(temp4, temp4, k2_s);
  178. fp12e_mul(temp1, temp1, temp2);
  179. fp12e_mul(temp2, temp3, temp4);
  180. fp12e_mul(temp1, temp1, temp2);
  181. retval.set(temp1, 0);
  182. fp12e_pow_vartime(temp1, input[0], j1_s);
  183. fp12e_pow_vartime(temp1, temp1, k1_s);
  184. fp12e_pow_vartime(temp1, temp1, j2_s);
  185. fp12e_pow_vartime(temp1, temp1, l2_s);
  186. fp12e_pow_vartime(temp2, input[1], j1_s);
  187. fp12e_invert(temp2, temp2);
  188. fp12e_pow_vartime(temp2, temp2, k1_s);
  189. fp12e_pow_vartime(temp2, temp2, i2_s);
  190. fp12e_pow_vartime(temp2, temp2, l2_s);
  191. fp12e_pow_vartime(temp3, input[2], i1_s);
  192. fp12e_pow_vartime(temp3, temp3, k1_s);
  193. fp12e_pow_vartime(temp3, temp3, j2_s);
  194. fp12e_invert(temp3, temp3);
  195. fp12e_pow_vartime(temp3, temp3, l2_s);
  196. fp12e_pow_vartime(temp4, input[3], i1_s);
  197. fp12e_pow_vartime(temp4, temp4, k1_s);
  198. fp12e_pow_vartime(temp4, temp4, i2_s);
  199. fp12e_pow_vartime(temp4, temp4, l2_s);
  200. fp12e_mul(temp1, temp1, temp2);
  201. fp12e_mul(temp2, temp3, temp4);
  202. fp12e_mul(temp1, temp1, temp2);
  203. retval.set(temp1, 1);
  204. fp12e_pow_vartime(temp1, input[0], j1_s);
  205. fp12e_pow_vartime(temp1, temp1, l1_s);
  206. fp12e_pow_vartime(temp1, temp1, j2_s);
  207. fp12e_pow_vartime(temp1, temp1, k2_s);
  208. fp12e_pow_vartime(temp2, input[1], j1_s);
  209. fp12e_invert(temp2, temp2);
  210. fp12e_pow_vartime(temp2, temp2, l1_s);
  211. fp12e_pow_vartime(temp2, temp2, i2_s);
  212. fp12e_pow_vartime(temp2, temp2, k2_s);
  213. fp12e_pow_vartime(temp3, input[2], i1_s);
  214. fp12e_pow_vartime(temp3, temp3, l1_s);
  215. fp12e_pow_vartime(temp3, temp3, j2_s);
  216. fp12e_invert(temp3, temp3);
  217. fp12e_pow_vartime(temp3, temp3, k2_s);
  218. fp12e_pow_vartime(temp4, input[3], i1_s);
  219. fp12e_pow_vartime(temp4, temp4, l1_s);
  220. fp12e_pow_vartime(temp4, temp4, i2_s);
  221. fp12e_pow_vartime(temp4, temp4, k2_s);
  222. fp12e_mul(temp1, temp1, temp2);
  223. fp12e_mul(temp2, temp3, temp4);
  224. fp12e_mul(temp1, temp1, temp2);
  225. retval.set(temp1, 2);
  226. fp12e_pow_vartime(temp1, input[0], j1_s);
  227. fp12e_pow_vartime(temp1, temp1, l1_s);
  228. fp12e_pow_vartime(temp1, temp1, j2_s);
  229. fp12e_pow_vartime(temp1, temp1, l2_s);
  230. fp12e_pow_vartime(temp2, input[1], j1_s);
  231. fp12e_invert(temp2, temp2);
  232. fp12e_pow_vartime(temp2, temp2, l1_s);
  233. fp12e_pow_vartime(temp2, temp2, i2_s);
  234. fp12e_pow_vartime(temp2, temp2, l2_s);
  235. fp12e_pow_vartime(temp3, input[2], i1_s);
  236. fp12e_pow_vartime(temp3, temp3, l1_s);
  237. fp12e_pow_vartime(temp3, temp3, j2_s);
  238. fp12e_invert(temp3, temp3);
  239. fp12e_pow_vartime(temp3, temp3, l2_s);
  240. fp12e_pow_vartime(temp4, input[3], i1_s);
  241. fp12e_pow_vartime(temp4, temp4, l1_s);
  242. fp12e_pow_vartime(temp4, temp4, i2_s);
  243. fp12e_pow_vartime(temp4, temp4, l2_s);
  244. fp12e_mul(temp1, temp1, temp2);
  245. fp12e_mul(temp2, temp3, temp4);
  246. fp12e_mul(temp1, temp1, temp2);
  247. retval.set(temp1, 3);
  248. return retval;
  249. }