git.sur5r.net Git - freertos/blob - FreeRTOS-Plus/Source/WolfSSL/wolfcrypt/src/fp_sqr_comba_12.i

   1 /* fp_sqr_comba_12.i
   2  *
   3  * Copyright (C) 2006-2015 wolfSSL Inc.
   4  *
   5  * This file is part of wolfSSL. (formerly known as CyaSSL)
   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-1301, USA
  20  */
  21
  22
  23 #ifdef TFM_SQR12
  24 void fp_sqr_comba12(fp_int *A, fp_int *B)
  25 {
  26    fp_digit *a, b[24], c0, c1, c2, sc0, sc1, sc2;
  27 #ifdef TFM_ISO
  28    fp_word tt;
  29 #endif
  30    a = A->dp;
  31    COMBA_START;
  32
  33    /* clear carries */
  34    CLEAR_CARRY;
  35
  36    /* output 0 */
  37    SQRADD(a[0],a[0]);
  38    COMBA_STORE(b[0]);
  39
  40    /* output 1 */
  41    CARRY_FORWARD;
  42    SQRADD2(a[0], a[1]);
  43    COMBA_STORE(b[1]);
  44
  45    /* output 2 */
  46    CARRY_FORWARD;
  47    SQRADD2(a[0], a[2]); SQRADD(a[1], a[1]);
  48    COMBA_STORE(b[2]);
  49
  50    /* output 3 */
  51    CARRY_FORWARD;
  52    SQRADD2(a[0], a[3]); SQRADD2(a[1], a[2]);
  53    COMBA_STORE(b[3]);
  54
  55    /* output 4 */
  56    CARRY_FORWARD;
  57    SQRADD2(a[0], a[4]); SQRADD2(a[1], a[3]); SQRADD(a[2], a[2]);
  58    COMBA_STORE(b[4]);
  59
  60    /* output 5 */
  61    CARRY_FORWARD;
  62    SQRADDSC(a[0], a[5]); SQRADDAC(a[1], a[4]); SQRADDAC(a[2], a[3]); SQRADDDB;
  63    COMBA_STORE(b[5]);
  64
  65    /* output 6 */
  66    CARRY_FORWARD;
  67    SQRADDSC(a[0], a[6]); SQRADDAC(a[1], a[5]); SQRADDAC(a[2], a[4]); SQRADDDB; SQRADD(a[3], a[3]);
  68    COMBA_STORE(b[6]);
  69
  70    /* output 7 */
  71    CARRY_FORWARD;
  72    SQRADDSC(a[0], a[7]); SQRADDAC(a[1], a[6]); SQRADDAC(a[2], a[5]); SQRADDAC(a[3], a[4]); SQRADDDB;
  73    COMBA_STORE(b[7]);
  74
  75    /* output 8 */
  76    CARRY_FORWARD;
  77    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]);
  78    COMBA_STORE(b[8]);
  79
  80    /* output 9 */
  81    CARRY_FORWARD;
  82    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;
  83    COMBA_STORE(b[9]);
  84
  85    /* output 10 */
  86    CARRY_FORWARD;
  87    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]);
  88    COMBA_STORE(b[10]);
  89
  90    /* output 11 */
  91    CARRY_FORWARD;
  92    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;
  93    COMBA_STORE(b[11]);
  94
  95    /* output 12 */
  96    CARRY_FORWARD;
  97    SQRADDSC(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]);
  98    COMBA_STORE(b[12]);
  99
 100    /* output 13 */
 101    CARRY_FORWARD;
 102    SQRADDSC(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;
 103    COMBA_STORE(b[13]);
 104
 105    /* output 14 */
 106    CARRY_FORWARD;
 107    SQRADDSC(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]);
 108    COMBA_STORE(b[14]);
 109
 110    /* output 15 */
 111    CARRY_FORWARD;
 112    SQRADDSC(a[4], a[11]); SQRADDAC(a[5], a[10]); SQRADDAC(a[6], a[9]); SQRADDAC(a[7], a[8]); SQRADDDB;
 113    COMBA_STORE(b[15]);
 114
 115    /* output 16 */
 116    CARRY_FORWARD;
 117    SQRADDSC(a[5], a[11]); SQRADDAC(a[6], a[10]); SQRADDAC(a[7], a[9]); SQRADDDB; SQRADD(a[8], a[8]);
 118    COMBA_STORE(b[16]);
 119
 120    /* output 17 */
 121    CARRY_FORWARD;
 122    SQRADDSC(a[6], a[11]); SQRADDAC(a[7], a[10]); SQRADDAC(a[8], a[9]); SQRADDDB;
 123    COMBA_STORE(b[17]);
 124
 125    /* output 18 */
 126    CARRY_FORWARD;
 127    SQRADD2(a[7], a[11]); SQRADD2(a[8], a[10]); SQRADD(a[9], a[9]);
 128    COMBA_STORE(b[18]);
 129
 130    /* output 19 */
 131    CARRY_FORWARD;
 132    SQRADD2(a[8], a[11]); SQRADD2(a[9], a[10]);
 133    COMBA_STORE(b[19]);
 134
 135    /* output 20 */
 136    CARRY_FORWARD;
 137    SQRADD2(a[9], a[11]); SQRADD(a[10], a[10]);
 138    COMBA_STORE(b[20]);
 139
 140    /* output 21 */
 141    CARRY_FORWARD;
 142    SQRADD2(a[10], a[11]);
 143    COMBA_STORE(b[21]);
 144
 145    /* output 22 */
 146    CARRY_FORWARD;
 147    SQRADD(a[11], a[11]);
 148    COMBA_STORE(b[22]);
 149    COMBA_STORE2(b[23]);
 150    COMBA_FINI;
 151
 152    B->used = 24;
 153    B->sign = FP_ZPOS;
 154    memcpy(B->dp, b, 24 * sizeof(fp_digit));
 155    fp_clamp(B);
 156 }
 157 #endif
 158
 159