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

   1 /* fp_sqr_comba_17.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_SQR17
  24 void fp_sqr_comba17(fp_int *A, fp_int *B)
  25 {
  26    fp_digit *a, b[34], c0, c1, c2, sc0, sc1, sc2;
  27
  28    a = A->dp;
  29    COMBA_START;
  30
  31    /* clear carries */
  32    CLEAR_CARRY;
  33
  34    /* output 0 */
  35    SQRADD(a[0],a[0]);
  36    COMBA_STORE(b[0]);
  37
  38    /* output 1 */
  39    CARRY_FORWARD;
  40    SQRADD2(a[0], a[1]);
  41    COMBA_STORE(b[1]);
  42
  43    /* output 2 */
  44    CARRY_FORWARD;
  45    SQRADD2(a[0], a[2]); SQRADD(a[1], a[1]);
  46    COMBA_STORE(b[2]);
  47
  48    /* output 3 */
  49    CARRY_FORWARD;
  50    SQRADD2(a[0], a[3]); SQRADD2(a[1], a[2]);
  51    COMBA_STORE(b[3]);
  52
  53    /* output 4 */
  54    CARRY_FORWARD;
  55    SQRADD2(a[0], a[4]); SQRADD2(a[1], a[3]); SQRADD(a[2], a[2]);
  56    COMBA_STORE(b[4]);
  57
  58    /* output 5 */
  59    CARRY_FORWARD;
  60    SQRADDSC(a[0], a[5]); SQRADDAC(a[1], a[4]); SQRADDAC(a[2], a[3]); SQRADDDB;
  61    COMBA_STORE(b[5]);
  62
  63    /* output 6 */
  64    CARRY_FORWARD;
  65    SQRADDSC(a[0], a[6]); SQRADDAC(a[1], a[5]); SQRADDAC(a[2], a[4]); SQRADDDB; SQRADD(a[3], a[3]);
  66    COMBA_STORE(b[6]);
  67
  68    /* output 7 */
  69    CARRY_FORWARD;
  70    SQRADDSC(a[0], a[7]); SQRADDAC(a[1], a[6]); SQRADDAC(a[2], a[5]); SQRADDAC(a[3], a[4]); SQRADDDB;
  71    COMBA_STORE(b[7]);
  72
  73    /* output 8 */
  74    CARRY_FORWARD;
  75    SQRADDSC(a[0], a[8]); SQRADDAC(a[1], a[7]); SQRADDAC(a[2], a[6]); SQRADDAC(a[3], a[5]); SQRADDDB; SQRADD(a[4], a[4]);
  76    COMBA_STORE(b[8]);
  77
  78    /* output 9 */
  79    CARRY_FORWARD;
  80    SQRADDSC(a[0], a[9]); SQRADDAC(a[1], a[8]); SQRADDAC(a[2], a[7]); SQRADDAC(a[3], a[6]); SQRADDAC(a[4], a[5]); SQRADDDB;
  81    COMBA_STORE(b[9]);
  82
  83    /* output 10 */
  84    CARRY_FORWARD;
  85    SQRADDSC(a[0], a[10]); SQRADDAC(a[1], a[9]); SQRADDAC(a[2], a[8]); SQRADDAC(a[3], a[7]); SQRADDAC(a[4], a[6]); SQRADDDB; SQRADD(a[5], a[5]);
  86    COMBA_STORE(b[10]);
  87
  88    /* output 11 */
  89    CARRY_FORWARD;
  90    SQRADDSC(a[0], a[11]); SQRADDAC(a[1], a[10]); SQRADDAC(a[2], a[9]); SQRADDAC(a[3], a[8]); SQRADDAC(a[4], a[7]); SQRADDAC(a[5], a[6]); SQRADDDB;
  91    COMBA_STORE(b[11]);
  92
  93    /* output 12 */
  94    CARRY_FORWARD;
  95    SQRADDSC(a[0], a[12]); SQRADDAC(a[1], a[11]); SQRADDAC(a[2], a[10]); SQRADDAC(a[3], a[9]); SQRADDAC(a[4], a[8]); SQRADDAC(a[5], a[7]); SQRADDDB; SQRADD(a[6], a[6]);
  96    COMBA_STORE(b[12]);
  97
  98    /* output 13 */
  99    CARRY_FORWARD;
 100    SQRADDSC(a[0], a[13]); SQRADDAC(a[1], a[12]); SQRADDAC(a[2], a[11]); SQRADDAC(a[3], a[10]); SQRADDAC(a[4], a[9]); SQRADDAC(a[5], a[8]); SQRADDAC(a[6], a[7]); SQRADDDB;
 101    COMBA_STORE(b[13]);
 102
 103    /* output 14 */
 104    CARRY_FORWARD;
 105    SQRADDSC(a[0], a[14]); SQRADDAC(a[1], a[13]); SQRADDAC(a[2], a[12]); SQRADDAC(a[3], a[11]); SQRADDAC(a[4], a[10]); SQRADDAC(a[5], a[9]); SQRADDAC(a[6], a[8]); SQRADDDB; SQRADD(a[7], a[7]);
 106    COMBA_STORE(b[14]);
 107
 108    /* output 15 */
 109    CARRY_FORWARD;
 110    SQRADDSC(a[0], a[15]); SQRADDAC(a[1], a[14]); SQRADDAC(a[2], a[13]); SQRADDAC(a[3], a[12]); SQRADDAC(a[4], a[11]); SQRADDAC(a[5], a[10]); SQRADDAC(a[6], a[9]); SQRADDAC(a[7], a[8]); SQRADDDB;
 111    COMBA_STORE(b[15]);
 112
 113    /* output 16 */
 114    CARRY_FORWARD;
 115    SQRADDSC(a[0], a[16]); SQRADDAC(a[1], a[15]); SQRADDAC(a[2], a[14]); SQRADDAC(a[3], a[13]); SQRADDAC(a[4], a[12]); SQRADDAC(a[5], a[11]); SQRADDAC(a[6], a[10]); SQRADDAC(a[7], a[9]); SQRADDDB; SQRADD(a[8], a[8]);
 116    COMBA_STORE(b[16]);
 117
 118    /* output 17 */
 119    CARRY_FORWARD;
 120    SQRADDSC(a[1], a[16]); SQRADDAC(a[2], a[15]); SQRADDAC(a[3], a[14]); SQRADDAC(a[4], a[13]); SQRADDAC(a[5], a[12]); SQRADDAC(a[6], a[11]); SQRADDAC(a[7], a[10]); SQRADDAC(a[8], a[9]); SQRADDDB;
 121    COMBA_STORE(b[17]);
 122
 123    /* output 18 */
 124    CARRY_FORWARD;
 125    SQRADDSC(a[2], a[16]); SQRADDAC(a[3], a[15]); SQRADDAC(a[4], a[14]); SQRADDAC(a[5], a[13]); SQRADDAC(a[6], a[12]); SQRADDAC(a[7], a[11]); SQRADDAC(a[8], a[10]); SQRADDDB; SQRADD(a[9], a[9]);
 126    COMBA_STORE(b[18]);
 127
 128    /* output 19 */
 129    CARRY_FORWARD;
 130    SQRADDSC(a[3], a[16]); SQRADDAC(a[4], a[15]); SQRADDAC(a[5], a[14]); SQRADDAC(a[6], a[13]); SQRADDAC(a[7], a[12]); SQRADDAC(a[8], a[11]); SQRADDAC(a[9], a[10]); SQRADDDB;
 131    COMBA_STORE(b[19]);
 132
 133    /* output 20 */
 134    CARRY_FORWARD;
 135    SQRADDSC(a[4], a[16]); SQRADDAC(a[5], a[15]); SQRADDAC(a[6], a[14]); SQRADDAC(a[7], a[13]); SQRADDAC(a[8], a[12]); SQRADDAC(a[9], a[11]); SQRADDDB; SQRADD(a[10], a[10]);
 136    COMBA_STORE(b[20]);
 137
 138    /* output 21 */
 139    CARRY_FORWARD;
 140    SQRADDSC(a[5], a[16]); SQRADDAC(a[6], a[15]); SQRADDAC(a[7], a[14]); SQRADDAC(a[8], a[13]); SQRADDAC(a[9], a[12]); SQRADDAC(a[10], a[11]); SQRADDDB;
 141    COMBA_STORE(b[21]);
 142
 143    /* output 22 */
 144    CARRY_FORWARD;
 145    SQRADDSC(a[6], a[16]); SQRADDAC(a[7], a[15]); SQRADDAC(a[8], a[14]); SQRADDAC(a[9], a[13]); SQRADDAC(a[10], a[12]); SQRADDDB; SQRADD(a[11], a[11]);
 146    COMBA_STORE(b[22]);
 147
 148    /* output 23 */
 149    CARRY_FORWARD;
 150    SQRADDSC(a[7], a[16]); SQRADDAC(a[8], a[15]); SQRADDAC(a[9], a[14]); SQRADDAC(a[10], a[13]); SQRADDAC(a[11], a[12]); SQRADDDB;
 151    COMBA_STORE(b[23]);
 152
 153    /* output 24 */
 154    CARRY_FORWARD;
 155    SQRADDSC(a[8], a[16]); SQRADDAC(a[9], a[15]); SQRADDAC(a[10], a[14]); SQRADDAC(a[11], a[13]); SQRADDDB; SQRADD(a[12], a[12]);
 156    COMBA_STORE(b[24]);
 157
 158    /* output 25 */
 159    CARRY_FORWARD;
 160    SQRADDSC(a[9], a[16]); SQRADDAC(a[10], a[15]); SQRADDAC(a[11], a[14]); SQRADDAC(a[12], a[13]); SQRADDDB;
 161    COMBA_STORE(b[25]);
 162
 163    /* output 26 */
 164    CARRY_FORWARD;
 165    SQRADDSC(a[10], a[16]); SQRADDAC(a[11], a[15]); SQRADDAC(a[12], a[14]); SQRADDDB; SQRADD(a[13], a[13]);
 166    COMBA_STORE(b[26]);
 167
 168    /* output 27 */
 169    CARRY_FORWARD;
 170    SQRADDSC(a[11], a[16]); SQRADDAC(a[12], a[15]); SQRADDAC(a[13], a[14]); SQRADDDB;
 171    COMBA_STORE(b[27]);
 172
 173    /* output 28 */
 174    CARRY_FORWARD;
 175    SQRADD2(a[12], a[16]); SQRADD2(a[13], a[15]); SQRADD(a[14], a[14]);
 176    COMBA_STORE(b[28]);
 177
 178    /* output 29 */
 179    CARRY_FORWARD;
 180    SQRADD2(a[13], a[16]); SQRADD2(a[14], a[15]);
 181    COMBA_STORE(b[29]);
 182
 183    /* output 30 */
 184    CARRY_FORWARD;
 185    SQRADD2(a[14], a[16]); SQRADD(a[15], a[15]);
 186    COMBA_STORE(b[30]);
 187
 188    /* output 31 */
 189    CARRY_FORWARD;
 190    SQRADD2(a[15], a[16]);
 191    COMBA_STORE(b[31]);
 192
 193    /* output 32 */
 194    CARRY_FORWARD;
 195    SQRADD(a[16], a[16]);
 196    COMBA_STORE(b[32]);
 197    COMBA_STORE2(b[33]);
 198    COMBA_FINI;
 199
 200    B->used = 34;
 201    B->sign = FP_ZPOS;
 202    memcpy(B->dp, b, 34 * sizeof(fp_digit));
 203    fp_clamp(B);
 204 }
 205 #endif
 206
 207