2 #include "fp_mul_comba.c"
4 #if defined(TFM_MUL4) && FP_SIZE >= 8
5 void fp_mul_comba4(fp_int *A, fp_int *B, fp_int *C)
7 fp_digit c0, c1, c2, at[8];
9 memcpy(at, A->dp, 4 * sizeof(fp_digit));
10 memcpy(at+4, B->dp, 4 * sizeof(fp_digit));
16 COMBA_STORE(C->dp[0]);
19 MULADD(at[0], at[5]); MULADD(at[1], at[4]);
20 COMBA_STORE(C->dp[1]);
23 MULADD(at[0], at[6]); MULADD(at[1], at[5]); MULADD(at[2], at[4]);
24 COMBA_STORE(C->dp[2]);
27 MULADD(at[0], at[7]); MULADD(at[1], at[6]); MULADD(at[2], at[5]); MULADD(at[3], at[4]);
28 COMBA_STORE(C->dp[3]);
31 MULADD(at[1], at[7]); MULADD(at[2], at[6]); MULADD(at[3], at[5]);
32 COMBA_STORE(C->dp[4]);
35 MULADD(at[2], at[7]); MULADD(at[3], at[6]);
36 COMBA_STORE(C->dp[5]);
40 COMBA_STORE(C->dp[6]);
41 COMBA_STORE2(C->dp[7]);
43 C->sign = A->sign ^ B->sign;