fp_mul_comba_17.i 11 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187
  1. /* fp_mul_comba_17.i
  2. *
  3. * Copyright (C) 2006-2023 wolfSSL Inc.
  4. *
  5. * This file is part of wolfSSL.
  6. *
  7. * wolfSSL is free software; you can redistribute it and/or modify
  8. * it under the terms of the GNU General Public License as published by
  9. * the Free Software Foundation; either version 2 of the License, or
  10. * (at your option) any later version.
  11. *
  12. * wolfSSL is distributed in the hope that it will be useful,
  13. * but WITHOUT ANY WARRANTY; without even the implied warranty of
  14. * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
  15. * GNU General Public License for more details.
  16. *
  17. * You should have received a copy of the GNU General Public License
  18. * along with this program; if not, write to the Free Software
  19. * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1335, USA
  20. */
  21. #ifdef TFM_MUL17
  22. int fp_mul_comba17(fp_int *A, fp_int *B, fp_int *C)
  23. {
  24. fp_digit c0, c1, c2;
  25. #ifndef WOLFSSL_SMALL_STACK
  26. fp_digit at[34];
  27. #else
  28. fp_digit *at;
  29. #endif
  30. #ifdef WOLFSSL_SMALL_STACK
  31. at = (fp_digit*)XMALLOC(sizeof(fp_digit) * 34, NULL, DYNAMIC_TYPE_TMP_BUFFER);
  32. if (at == NULL)
  33. return FP_MEM;
  34. #endif
  35. XMEMCPY(at, A->dp, 17 * sizeof(fp_digit));
  36. XMEMCPY(at+17, B->dp, 17 * sizeof(fp_digit));
  37. COMBA_START;
  38. COMBA_CLEAR;
  39. /* 0 */
  40. MULADD(at[0], at[17]);
  41. COMBA_STORE(C->dp[0]);
  42. /* 1 */
  43. COMBA_FORWARD;
  44. MULADD(at[0], at[18]); MULADD(at[1], at[17]);
  45. COMBA_STORE(C->dp[1]);
  46. /* 2 */
  47. COMBA_FORWARD;
  48. MULADD(at[0], at[19]); MULADD(at[1], at[18]); MULADD(at[2], at[17]);
  49. COMBA_STORE(C->dp[2]);
  50. /* 3 */
  51. COMBA_FORWARD;
  52. MULADD(at[0], at[20]); MULADD(at[1], at[19]); MULADD(at[2], at[18]); MULADD(at[3], at[17]);
  53. COMBA_STORE(C->dp[3]);
  54. /* 4 */
  55. COMBA_FORWARD;
  56. MULADD(at[0], at[21]); MULADD(at[1], at[20]); MULADD(at[2], at[19]); MULADD(at[3], at[18]); MULADD(at[4], at[17]);
  57. COMBA_STORE(C->dp[4]);
  58. /* 5 */
  59. COMBA_FORWARD;
  60. MULADD(at[0], at[22]); MULADD(at[1], at[21]); MULADD(at[2], at[20]); MULADD(at[3], at[19]); MULADD(at[4], at[18]); MULADD(at[5], at[17]);
  61. COMBA_STORE(C->dp[5]);
  62. /* 6 */
  63. COMBA_FORWARD;
  64. MULADD(at[0], at[23]); MULADD(at[1], at[22]); MULADD(at[2], at[21]); MULADD(at[3], at[20]); MULADD(at[4], at[19]); MULADD(at[5], at[18]); MULADD(at[6], at[17]);
  65. COMBA_STORE(C->dp[6]);
  66. /* 7 */
  67. COMBA_FORWARD;
  68. MULADD(at[0], at[24]); MULADD(at[1], at[23]); MULADD(at[2], at[22]); MULADD(at[3], at[21]); MULADD(at[4], at[20]); MULADD(at[5], at[19]); MULADD(at[6], at[18]); MULADD(at[7], at[17]);
  69. COMBA_STORE(C->dp[7]);
  70. /* 8 */
  71. COMBA_FORWARD;
  72. MULADD(at[0], at[25]); MULADD(at[1], at[24]); MULADD(at[2], at[23]); MULADD(at[3], at[22]); MULADD(at[4], at[21]); MULADD(at[5], at[20]); MULADD(at[6], at[19]); MULADD(at[7], at[18]); MULADD(at[8], at[17]);
  73. COMBA_STORE(C->dp[8]);
  74. /* 9 */
  75. COMBA_FORWARD;
  76. MULADD(at[0], at[26]); MULADD(at[1], at[25]); MULADD(at[2], at[24]); MULADD(at[3], at[23]); MULADD(at[4], at[22]); MULADD(at[5], at[21]); MULADD(at[6], at[20]); MULADD(at[7], at[19]); MULADD(at[8], at[18]); MULADD(at[9], at[17]);
  77. COMBA_STORE(C->dp[9]);
  78. /* 10 */
  79. COMBA_FORWARD;
  80. MULADD(at[0], at[27]); MULADD(at[1], at[26]); MULADD(at[2], at[25]); MULADD(at[3], at[24]); MULADD(at[4], at[23]); MULADD(at[5], at[22]); MULADD(at[6], at[21]); MULADD(at[7], at[20]); MULADD(at[8], at[19]); MULADD(at[9], at[18]); MULADD(at[10], at[17]);
  81. COMBA_STORE(C->dp[10]);
  82. /* 11 */
  83. COMBA_FORWARD;
  84. MULADD(at[0], at[28]); MULADD(at[1], at[27]); MULADD(at[2], at[26]); MULADD(at[3], at[25]); MULADD(at[4], at[24]); MULADD(at[5], at[23]); MULADD(at[6], at[22]); MULADD(at[7], at[21]); MULADD(at[8], at[20]); MULADD(at[9], at[19]); MULADD(at[10], at[18]); MULADD(at[11], at[17]);
  85. COMBA_STORE(C->dp[11]);
  86. /* 12 */
  87. COMBA_FORWARD;
  88. MULADD(at[0], at[29]); MULADD(at[1], at[28]); MULADD(at[2], at[27]); MULADD(at[3], at[26]); MULADD(at[4], at[25]); MULADD(at[5], at[24]); MULADD(at[6], at[23]); MULADD(at[7], at[22]); MULADD(at[8], at[21]); MULADD(at[9], at[20]); MULADD(at[10], at[19]); MULADD(at[11], at[18]); MULADD(at[12], at[17]);
  89. COMBA_STORE(C->dp[12]);
  90. /* 13 */
  91. COMBA_FORWARD;
  92. MULADD(at[0], at[30]); MULADD(at[1], at[29]); MULADD(at[2], at[28]); MULADD(at[3], at[27]); MULADD(at[4], at[26]); MULADD(at[5], at[25]); MULADD(at[6], at[24]); MULADD(at[7], at[23]); MULADD(at[8], at[22]); MULADD(at[9], at[21]); MULADD(at[10], at[20]); MULADD(at[11], at[19]); MULADD(at[12], at[18]); MULADD(at[13], at[17]);
  93. COMBA_STORE(C->dp[13]);
  94. /* 14 */
  95. COMBA_FORWARD;
  96. MULADD(at[0], at[31]); MULADD(at[1], at[30]); MULADD(at[2], at[29]); MULADD(at[3], at[28]); MULADD(at[4], at[27]); MULADD(at[5], at[26]); MULADD(at[6], at[25]); MULADD(at[7], at[24]); MULADD(at[8], at[23]); MULADD(at[9], at[22]); MULADD(at[10], at[21]); MULADD(at[11], at[20]); MULADD(at[12], at[19]); MULADD(at[13], at[18]); MULADD(at[14], at[17]);
  97. COMBA_STORE(C->dp[14]);
  98. /* 15 */
  99. COMBA_FORWARD;
  100. MULADD(at[0], at[32]); MULADD(at[1], at[31]); MULADD(at[2], at[30]); MULADD(at[3], at[29]); MULADD(at[4], at[28]); MULADD(at[5], at[27]); MULADD(at[6], at[26]); MULADD(at[7], at[25]); MULADD(at[8], at[24]); MULADD(at[9], at[23]); MULADD(at[10], at[22]); MULADD(at[11], at[21]); MULADD(at[12], at[20]); MULADD(at[13], at[19]); MULADD(at[14], at[18]); MULADD(at[15], at[17]);
  101. COMBA_STORE(C->dp[15]);
  102. /* 16 */
  103. COMBA_FORWARD;
  104. MULADD(at[0], at[33]); MULADD(at[1], at[32]); MULADD(at[2], at[31]); MULADD(at[3], at[30]); MULADD(at[4], at[29]); MULADD(at[5], at[28]); MULADD(at[6], at[27]); MULADD(at[7], at[26]); MULADD(at[8], at[25]); MULADD(at[9], at[24]); MULADD(at[10], at[23]); MULADD(at[11], at[22]); MULADD(at[12], at[21]); MULADD(at[13], at[20]); MULADD(at[14], at[19]); MULADD(at[15], at[18]); MULADD(at[16], at[17]);
  105. COMBA_STORE(C->dp[16]);
  106. /* 17 */
  107. COMBA_FORWARD;
  108. MULADD(at[1], at[33]); MULADD(at[2], at[32]); MULADD(at[3], at[31]); MULADD(at[4], at[30]); MULADD(at[5], at[29]); MULADD(at[6], at[28]); MULADD(at[7], at[27]); MULADD(at[8], at[26]); MULADD(at[9], at[25]); MULADD(at[10], at[24]); MULADD(at[11], at[23]); MULADD(at[12], at[22]); MULADD(at[13], at[21]); MULADD(at[14], at[20]); MULADD(at[15], at[19]); MULADD(at[16], at[18]);
  109. COMBA_STORE(C->dp[17]);
  110. /* 18 */
  111. COMBA_FORWARD;
  112. MULADD(at[2], at[33]); MULADD(at[3], at[32]); MULADD(at[4], at[31]); MULADD(at[5], at[30]); MULADD(at[6], at[29]); MULADD(at[7], at[28]); MULADD(at[8], at[27]); MULADD(at[9], at[26]); MULADD(at[10], at[25]); MULADD(at[11], at[24]); MULADD(at[12], at[23]); MULADD(at[13], at[22]); MULADD(at[14], at[21]); MULADD(at[15], at[20]); MULADD(at[16], at[19]);
  113. COMBA_STORE(C->dp[18]);
  114. /* 19 */
  115. COMBA_FORWARD;
  116. MULADD(at[3], at[33]); MULADD(at[4], at[32]); MULADD(at[5], at[31]); MULADD(at[6], at[30]); MULADD(at[7], at[29]); MULADD(at[8], at[28]); MULADD(at[9], at[27]); MULADD(at[10], at[26]); MULADD(at[11], at[25]); MULADD(at[12], at[24]); MULADD(at[13], at[23]); MULADD(at[14], at[22]); MULADD(at[15], at[21]); MULADD(at[16], at[20]);
  117. COMBA_STORE(C->dp[19]);
  118. /* 20 */
  119. COMBA_FORWARD;
  120. MULADD(at[4], at[33]); MULADD(at[5], at[32]); MULADD(at[6], at[31]); MULADD(at[7], at[30]); MULADD(at[8], at[29]); MULADD(at[9], at[28]); MULADD(at[10], at[27]); MULADD(at[11], at[26]); MULADD(at[12], at[25]); MULADD(at[13], at[24]); MULADD(at[14], at[23]); MULADD(at[15], at[22]); MULADD(at[16], at[21]);
  121. COMBA_STORE(C->dp[20]);
  122. /* 21 */
  123. COMBA_FORWARD;
  124. MULADD(at[5], at[33]); MULADD(at[6], at[32]); MULADD(at[7], at[31]); MULADD(at[8], at[30]); MULADD(at[9], at[29]); MULADD(at[10], at[28]); MULADD(at[11], at[27]); MULADD(at[12], at[26]); MULADD(at[13], at[25]); MULADD(at[14], at[24]); MULADD(at[15], at[23]); MULADD(at[16], at[22]);
  125. COMBA_STORE(C->dp[21]);
  126. /* 22 */
  127. COMBA_FORWARD;
  128. MULADD(at[6], at[33]); MULADD(at[7], at[32]); MULADD(at[8], at[31]); MULADD(at[9], at[30]); MULADD(at[10], at[29]); MULADD(at[11], at[28]); MULADD(at[12], at[27]); MULADD(at[13], at[26]); MULADD(at[14], at[25]); MULADD(at[15], at[24]); MULADD(at[16], at[23]);
  129. COMBA_STORE(C->dp[22]);
  130. /* 23 */
  131. COMBA_FORWARD;
  132. MULADD(at[7], at[33]); MULADD(at[8], at[32]); MULADD(at[9], at[31]); MULADD(at[10], at[30]); MULADD(at[11], at[29]); MULADD(at[12], at[28]); MULADD(at[13], at[27]); MULADD(at[14], at[26]); MULADD(at[15], at[25]); MULADD(at[16], at[24]);
  133. COMBA_STORE(C->dp[23]);
  134. /* 24 */
  135. COMBA_FORWARD;
  136. MULADD(at[8], at[33]); MULADD(at[9], at[32]); MULADD(at[10], at[31]); MULADD(at[11], at[30]); MULADD(at[12], at[29]); MULADD(at[13], at[28]); MULADD(at[14], at[27]); MULADD(at[15], at[26]); MULADD(at[16], at[25]);
  137. COMBA_STORE(C->dp[24]);
  138. /* 25 */
  139. COMBA_FORWARD;
  140. MULADD(at[9], at[33]); MULADD(at[10], at[32]); MULADD(at[11], at[31]); MULADD(at[12], at[30]); MULADD(at[13], at[29]); MULADD(at[14], at[28]); MULADD(at[15], at[27]); MULADD(at[16], at[26]);
  141. COMBA_STORE(C->dp[25]);
  142. /* 26 */
  143. COMBA_FORWARD;
  144. MULADD(at[10], at[33]); MULADD(at[11], at[32]); MULADD(at[12], at[31]); MULADD(at[13], at[30]); MULADD(at[14], at[29]); MULADD(at[15], at[28]); MULADD(at[16], at[27]);
  145. COMBA_STORE(C->dp[26]);
  146. /* 27 */
  147. COMBA_FORWARD;
  148. MULADD(at[11], at[33]); MULADD(at[12], at[32]); MULADD(at[13], at[31]); MULADD(at[14], at[30]); MULADD(at[15], at[29]); MULADD(at[16], at[28]);
  149. COMBA_STORE(C->dp[27]);
  150. /* 28 */
  151. COMBA_FORWARD;
  152. MULADD(at[12], at[33]); MULADD(at[13], at[32]); MULADD(at[14], at[31]); MULADD(at[15], at[30]); MULADD(at[16], at[29]);
  153. COMBA_STORE(C->dp[28]);
  154. /* 29 */
  155. COMBA_FORWARD;
  156. MULADD(at[13], at[33]); MULADD(at[14], at[32]); MULADD(at[15], at[31]); MULADD(at[16], at[30]);
  157. COMBA_STORE(C->dp[29]);
  158. /* 30 */
  159. COMBA_FORWARD;
  160. MULADD(at[14], at[33]); MULADD(at[15], at[32]); MULADD(at[16], at[31]);
  161. COMBA_STORE(C->dp[30]);
  162. /* 31 */
  163. COMBA_FORWARD;
  164. MULADD(at[15], at[33]); MULADD(at[16], at[32]);
  165. COMBA_STORE(C->dp[31]);
  166. /* 32 */
  167. COMBA_FORWARD;
  168. MULADD(at[16], at[33]);
  169. COMBA_STORE(C->dp[32]);
  170. COMBA_STORE2(C->dp[33]);
  171. C->used = 34;
  172. C->sign = A->sign ^ B->sign;
  173. fp_clamp(C);
  174. COMBA_FINI;
  175. #ifdef WOLFSSL_SMALL_STACK
  176. XFREE(at, NULL, DYNAMIC_TYPE_TMP_BUFFER);
  177. #endif
  178. return FP_OKAY;
  179. }
  180. #endif