git.sur5r.net Git - freertos/blob - FreeRTOS-Plus/Source/CyaSSL/ctaocrypt/src/fp_mul_comba_12.i

   1 /* fp_mul_comba_12.i
   2  *
   3  * Copyright (C) 2006-2011 Sawtooth Consulting Ltd.
   4  *
   5  * This file is part of CyaSSL.
   6  *
   7  * CyaSSL 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  * CyaSSL 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., 59 Temple Place - Suite 330, Boston, MA 02111-1307, USA
  20  */
  21
  22
  23 #ifdef TFM_MUL12
  24 void fp_mul_comba12(fp_int *A, fp_int *B, fp_int *C)
  25 {
  26    fp_digit c0, c1, c2, at[24];
  27
  28    memcpy(at, A->dp, 12 * sizeof(fp_digit));
  29    memcpy(at+12, B->dp, 12 * sizeof(fp_digit));
  30    COMBA_START;
  31
  32    COMBA_CLEAR;
  33    /* 0 */
  34    MULADD(at[0], at[12]);
  35    COMBA_STORE(C->dp[0]);
  36    /* 1 */
  37    COMBA_FORWARD;
  38    MULADD(at[0], at[13]);    MULADD(at[1], at[12]);
  39    COMBA_STORE(C->dp[1]);
  40    /* 2 */
  41    COMBA_FORWARD;
  42    MULADD(at[0], at[14]);    MULADD(at[1], at[13]);    MULADD(at[2], at[12]);
  43    COMBA_STORE(C->dp[2]);
  44    /* 3 */
  45    COMBA_FORWARD;
  46    MULADD(at[0], at[15]);    MULADD(at[1], at[14]);    MULADD(at[2], at[13]);    MULADD(at[3], at[12]);
  47    COMBA_STORE(C->dp[3]);
  48    /* 4 */
  49    COMBA_FORWARD;
  50    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]);
  51    COMBA_STORE(C->dp[4]);
  52    /* 5 */
  53    COMBA_FORWARD;
  54    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]);
  55    COMBA_STORE(C->dp[5]);
  56    /* 6 */
  57    COMBA_FORWARD;
  58    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]);
  59    COMBA_STORE(C->dp[6]);
  60    /* 7 */
  61    COMBA_FORWARD;
  62    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]);
  63    COMBA_STORE(C->dp[7]);
  64    /* 8 */
  65    COMBA_FORWARD;
  66    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]);
  67    COMBA_STORE(C->dp[8]);
  68    /* 9 */
  69    COMBA_FORWARD;
  70    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]);
  71    COMBA_STORE(C->dp[9]);
  72    /* 10 */
  73    COMBA_FORWARD;
  74    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]);
  75    COMBA_STORE(C->dp[10]);
  76    /* 11 */
  77    COMBA_FORWARD;
  78    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]);
  79    COMBA_STORE(C->dp[11]);
  80    /* 12 */
  81    COMBA_FORWARD;
  82    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]);
  83    COMBA_STORE(C->dp[12]);
  84    /* 13 */
  85    COMBA_FORWARD;
  86    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]);
  87    COMBA_STORE(C->dp[13]);
  88    /* 14 */
  89    COMBA_FORWARD;
  90    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]);
  91    COMBA_STORE(C->dp[14]);
  92    /* 15 */
  93    COMBA_FORWARD;
  94    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]);
  95    COMBA_STORE(C->dp[15]);
  96    /* 16 */
  97    COMBA_FORWARD;
  98    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]);
  99    COMBA_STORE(C->dp[16]);
 100    /* 17 */
 101    COMBA_FORWARD;
 102    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]);
 103    COMBA_STORE(C->dp[17]);
 104    /* 18 */
 105    COMBA_FORWARD;
 106    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]);
 107    COMBA_STORE(C->dp[18]);
 108    /* 19 */
 109    COMBA_FORWARD;
 110    MULADD(at[8], at[23]);    MULADD(at[9], at[22]);    MULADD(at[10], at[21]);    MULADD(at[11], at[20]);
 111    COMBA_STORE(C->dp[19]);
 112    /* 20 */
 113    COMBA_FORWARD;
 114    MULADD(at[9], at[23]);    MULADD(at[10], at[22]);    MULADD(at[11], at[21]);
 115    COMBA_STORE(C->dp[20]);
 116    /* 21 */
 117    COMBA_FORWARD;
 118    MULADD(at[10], at[23]);    MULADD(at[11], at[22]);
 119    COMBA_STORE(C->dp[21]);
 120    /* 22 */
 121    COMBA_FORWARD;
 122    MULADD(at[11], at[23]);
 123    COMBA_STORE(C->dp[22]);
 124    COMBA_STORE2(C->dp[23]);
 125    C->used = 24;
 126    C->sign = A->sign ^ B->sign;
 127    fp_clamp(C);
 128    COMBA_FINI;
 129 }
 130 #endif