fp_mul_comba_12.i 6.4 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147
  1. /* fp_mul_comba_12.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_MUL12
  22. int fp_mul_comba12(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[24];
  27. #else
  28. fp_digit *at;
  29. #endif
  30. #ifdef WOLFSSL_SMALL_STACK
  31. at = (fp_digit*)XMALLOC(sizeof(fp_digit) * 24, NULL, DYNAMIC_TYPE_TMP_BUFFER);
  32. if (at == NULL)
  33. return FP_MEM;
  34. #endif
  35. XMEMCPY(at, A->dp, 12 * sizeof(fp_digit));
  36. XMEMCPY(at+12, B->dp, 12 * sizeof(fp_digit));
  37. COMBA_START;
  38. COMBA_CLEAR;
  39. /* 0 */
  40. MULADD(at[0], at[12]);
  41. COMBA_STORE(C->dp[0]);
  42. /* 1 */
  43. COMBA_FORWARD;
  44. MULADD(at[0], at[13]); MULADD(at[1], at[12]);
  45. COMBA_STORE(C->dp[1]);
  46. /* 2 */
  47. COMBA_FORWARD;
  48. MULADD(at[0], at[14]); MULADD(at[1], at[13]); MULADD(at[2], at[12]);
  49. COMBA_STORE(C->dp[2]);
  50. /* 3 */
  51. COMBA_FORWARD;
  52. MULADD(at[0], at[15]); MULADD(at[1], at[14]); MULADD(at[2], at[13]); MULADD(at[3], at[12]);
  53. COMBA_STORE(C->dp[3]);
  54. /* 4 */
  55. COMBA_FORWARD;
  56. MULADD(at[0], at[16]); MULADD(at[1], at[15]); MULADD(at[2], at[14]); MULADD(at[3], at[13]); MULADD(at[4], at[12]);
  57. COMBA_STORE(C->dp[4]);
  58. /* 5 */
  59. COMBA_FORWARD;
  60. MULADD(at[0], at[17]); MULADD(at[1], at[16]); MULADD(at[2], at[15]); MULADD(at[3], at[14]); MULADD(at[4], at[13]); MULADD(at[5], at[12]);
  61. COMBA_STORE(C->dp[5]);
  62. /* 6 */
  63. COMBA_FORWARD;
  64. MULADD(at[0], at[18]); MULADD(at[1], at[17]); MULADD(at[2], at[16]); MULADD(at[3], at[15]); MULADD(at[4], at[14]); MULADD(at[5], at[13]); MULADD(at[6], at[12]);
  65. COMBA_STORE(C->dp[6]);
  66. /* 7 */
  67. COMBA_FORWARD;
  68. MULADD(at[0], at[19]); MULADD(at[1], at[18]); MULADD(at[2], at[17]); MULADD(at[3], at[16]); MULADD(at[4], at[15]); MULADD(at[5], at[14]); MULADD(at[6], at[13]); MULADD(at[7], at[12]);
  69. COMBA_STORE(C->dp[7]);
  70. /* 8 */
  71. COMBA_FORWARD;
  72. MULADD(at[0], at[20]); MULADD(at[1], at[19]); MULADD(at[2], at[18]); MULADD(at[3], at[17]); MULADD(at[4], at[16]); MULADD(at[5], at[15]); MULADD(at[6], at[14]); MULADD(at[7], at[13]); MULADD(at[8], at[12]);
  73. COMBA_STORE(C->dp[8]);
  74. /* 9 */
  75. COMBA_FORWARD;
  76. 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]); MULADD(at[5], at[16]); MULADD(at[6], at[15]); MULADD(at[7], at[14]); MULADD(at[8], at[13]); MULADD(at[9], at[12]);
  77. COMBA_STORE(C->dp[9]);
  78. /* 10 */
  79. COMBA_FORWARD;
  80. 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]); MULADD(at[6], at[16]); MULADD(at[7], at[15]); MULADD(at[8], at[14]); MULADD(at[9], at[13]); MULADD(at[10], at[12]);
  81. COMBA_STORE(C->dp[10]);
  82. /* 11 */
  83. COMBA_FORWARD;
  84. 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]); MULADD(at[7], at[16]); MULADD(at[8], at[15]); MULADD(at[9], at[14]); MULADD(at[10], at[13]); MULADD(at[11], at[12]);
  85. COMBA_STORE(C->dp[11]);
  86. /* 12 */
  87. COMBA_FORWARD;
  88. 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]); MULADD(at[8], at[16]); MULADD(at[9], at[15]); MULADD(at[10], at[14]); MULADD(at[11], at[13]);
  89. COMBA_STORE(C->dp[12]);
  90. /* 13 */
  91. COMBA_FORWARD;
  92. 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]); MULADD(at[9], at[16]); MULADD(at[10], at[15]); MULADD(at[11], at[14]);
  93. COMBA_STORE(C->dp[13]);
  94. /* 14 */
  95. COMBA_FORWARD;
  96. 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]); MULADD(at[10], at[16]); MULADD(at[11], at[15]);
  97. COMBA_STORE(C->dp[14]);
  98. /* 15 */
  99. COMBA_FORWARD;
  100. 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]); MULADD(at[11], at[16]);
  101. COMBA_STORE(C->dp[15]);
  102. /* 16 */
  103. COMBA_FORWARD;
  104. 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]);
  105. COMBA_STORE(C->dp[16]);
  106. /* 17 */
  107. COMBA_FORWARD;
  108. 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]);
  109. COMBA_STORE(C->dp[17]);
  110. /* 18 */
  111. COMBA_FORWARD;
  112. 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]);
  113. COMBA_STORE(C->dp[18]);
  114. /* 19 */
  115. COMBA_FORWARD;
  116. MULADD(at[8], at[23]); MULADD(at[9], at[22]); MULADD(at[10], at[21]); MULADD(at[11], at[20]);
  117. COMBA_STORE(C->dp[19]);
  118. /* 20 */
  119. COMBA_FORWARD;
  120. MULADD(at[9], at[23]); MULADD(at[10], at[22]); MULADD(at[11], at[21]);
  121. COMBA_STORE(C->dp[20]);
  122. /* 21 */
  123. COMBA_FORWARD;
  124. MULADD(at[10], at[23]); MULADD(at[11], at[22]);
  125. COMBA_STORE(C->dp[21]);
  126. /* 22 */
  127. COMBA_FORWARD;
  128. MULADD(at[11], at[23]);
  129. COMBA_STORE(C->dp[22]);
  130. COMBA_STORE2(C->dp[23]);
  131. C->used = 24;
  132. C->sign = A->sign ^ B->sign;
  133. fp_clamp(C);
  134. COMBA_FINI;
  135. #ifdef WOLFSSL_SMALL_STACK
  136. XFREE(at, NULL, DYNAMIC_TYPE_TMP_BUFFER);
  137. #endif
  138. return FP_OKAY;
  139. }
  140. #endif