fp_mul_comba_12.c

00001 #define TFM_DEFINES
00002 #include "fp_mul_comba.c"
00003 
00004 #ifdef TFM_MUL12
00005 void fp_mul_comba12(fp_int *A, fp_int *B, fp_int *C)
00006 {
00007    fp_digit c0, c1, c2, at[24];
00008 
00009    memcpy(at, A->dp, 12 * sizeof(fp_digit));
00010    memcpy(at+12, B->dp, 12 * sizeof(fp_digit));
00011    COMBA_START;
00012 
00013    COMBA_CLEAR;
00014    /* 0 */
00015    MULADD(at[0], at[12]); 
00016    COMBA_STORE(C->dp[0]);
00017    /* 1 */
00018    COMBA_FORWARD;
00019    MULADD(at[0], at[13]);    MULADD(at[1], at[12]); 
00020    COMBA_STORE(C->dp[1]);
00021    /* 2 */
00022    COMBA_FORWARD;
00023    MULADD(at[0], at[14]);    MULADD(at[1], at[13]);    MULADD(at[2], at[12]); 
00024    COMBA_STORE(C->dp[2]);
00025    /* 3 */
00026    COMBA_FORWARD;
00027    MULADD(at[0], at[15]);    MULADD(at[1], at[14]);    MULADD(at[2], at[13]);    MULADD(at[3], at[12]); 
00028    COMBA_STORE(C->dp[3]);
00029    /* 4 */
00030    COMBA_FORWARD;
00031    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]); 
00032    COMBA_STORE(C->dp[4]);
00033    /* 5 */
00034    COMBA_FORWARD;
00035    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]); 
00036    COMBA_STORE(C->dp[5]);
00037    /* 6 */
00038    COMBA_FORWARD;
00039    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]); 
00040    COMBA_STORE(C->dp[6]);
00041    /* 7 */
00042    COMBA_FORWARD;
00043    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]); 
00044    COMBA_STORE(C->dp[7]);
00045    /* 8 */
00046    COMBA_FORWARD;
00047    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]); 
00048    COMBA_STORE(C->dp[8]);
00049    /* 9 */
00050    COMBA_FORWARD;
00051    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]); 
00052    COMBA_STORE(C->dp[9]);
00053    /* 10 */
00054    COMBA_FORWARD;
00055    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]); 
00056    COMBA_STORE(C->dp[10]);
00057    /* 11 */
00058    COMBA_FORWARD;
00059    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]); 
00060    COMBA_STORE(C->dp[11]);
00061    /* 12 */
00062    COMBA_FORWARD;
00063    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]); 
00064    COMBA_STORE(C->dp[12]);
00065    /* 13 */
00066    COMBA_FORWARD;
00067    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]); 
00068    COMBA_STORE(C->dp[13]);
00069    /* 14 */
00070    COMBA_FORWARD;
00071    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]); 
00072    COMBA_STORE(C->dp[14]);
00073    /* 15 */
00074    COMBA_FORWARD;
00075    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]); 
00076    COMBA_STORE(C->dp[15]);
00077    /* 16 */
00078    COMBA_FORWARD;
00079    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]); 
00080    COMBA_STORE(C->dp[16]);
00081    /* 17 */
00082    COMBA_FORWARD;
00083    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]); 
00084    COMBA_STORE(C->dp[17]);
00085    /* 18 */
00086    COMBA_FORWARD;
00087    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]); 
00088    COMBA_STORE(C->dp[18]);
00089    /* 19 */
00090    COMBA_FORWARD;
00091    MULADD(at[8], at[23]);    MULADD(at[9], at[22]);    MULADD(at[10], at[21]);    MULADD(at[11], at[20]); 
00092    COMBA_STORE(C->dp[19]);
00093    /* 20 */
00094    COMBA_FORWARD;
00095    MULADD(at[9], at[23]);    MULADD(at[10], at[22]);    MULADD(at[11], at[21]); 
00096    COMBA_STORE(C->dp[20]);
00097    /* 21 */
00098    COMBA_FORWARD;
00099    MULADD(at[10], at[23]);    MULADD(at[11], at[22]); 
00100    COMBA_STORE(C->dp[21]);
00101    /* 22 */
00102    COMBA_FORWARD;
00103    MULADD(at[11], at[23]); 
00104    COMBA_STORE(C->dp[22]);
00105    COMBA_STORE2(C->dp[23]);
00106    C->used = 24;
00107    C->sign = A->sign ^ B->sign;
00108    fp_clamp(C);
00109    COMBA_FINI;
00110 }
00111 #endif