Mercurial > dropbear
comparison tomsfastmath/src/mul/fp_mul_comba_12.c @ 643:a362b62d38b2 dropbear-tfm
Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a
with Makefile.in renamed
author | Matt Johnston <matt@ucc.asn.au> |
---|---|
date | Wed, 23 Nov 2011 18:10:20 +0700 |
parents | |
children |
comparison
equal
deleted
inserted
replaced
642:33fd2f3499d2 | 643:a362b62d38b2 |
---|---|
1 #define TFM_DEFINES | |
2 #include "fp_mul_comba.c" | |
3 | |
4 #ifdef TFM_MUL12 | |
5 void fp_mul_comba12(fp_int *A, fp_int *B, fp_int *C) | |
6 { | |
7 fp_digit c0, c1, c2, at[24]; | |
8 | |
9 memcpy(at, A->dp, 12 * sizeof(fp_digit)); | |
10 memcpy(at+12, B->dp, 12 * sizeof(fp_digit)); | |
11 COMBA_START; | |
12 | |
13 COMBA_CLEAR; | |
14 /* 0 */ | |
15 MULADD(at[0], at[12]); | |
16 COMBA_STORE(C->dp[0]); | |
17 /* 1 */ | |
18 COMBA_FORWARD; | |
19 MULADD(at[0], at[13]); MULADD(at[1], at[12]); | |
20 COMBA_STORE(C->dp[1]); | |
21 /* 2 */ | |
22 COMBA_FORWARD; | |
23 MULADD(at[0], at[14]); MULADD(at[1], at[13]); MULADD(at[2], at[12]); | |
24 COMBA_STORE(C->dp[2]); | |
25 /* 3 */ | |
26 COMBA_FORWARD; | |
27 MULADD(at[0], at[15]); MULADD(at[1], at[14]); MULADD(at[2], at[13]); MULADD(at[3], at[12]); | |
28 COMBA_STORE(C->dp[3]); | |
29 /* 4 */ | |
30 COMBA_FORWARD; | |
31 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]); | |
32 COMBA_STORE(C->dp[4]); | |
33 /* 5 */ | |
34 COMBA_FORWARD; | |
35 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]); | |
36 COMBA_STORE(C->dp[5]); | |
37 /* 6 */ | |
38 COMBA_FORWARD; | |
39 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]); | |
40 COMBA_STORE(C->dp[6]); | |
41 /* 7 */ | |
42 COMBA_FORWARD; | |
43 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]); | |
44 COMBA_STORE(C->dp[7]); | |
45 /* 8 */ | |
46 COMBA_FORWARD; | |
47 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]); | |
48 COMBA_STORE(C->dp[8]); | |
49 /* 9 */ | |
50 COMBA_FORWARD; | |
51 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]); | |
52 COMBA_STORE(C->dp[9]); | |
53 /* 10 */ | |
54 COMBA_FORWARD; | |
55 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]); | |
56 COMBA_STORE(C->dp[10]); | |
57 /* 11 */ | |
58 COMBA_FORWARD; | |
59 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]); | |
60 COMBA_STORE(C->dp[11]); | |
61 /* 12 */ | |
62 COMBA_FORWARD; | |
63 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]); | |
64 COMBA_STORE(C->dp[12]); | |
65 /* 13 */ | |
66 COMBA_FORWARD; | |
67 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]); | |
68 COMBA_STORE(C->dp[13]); | |
69 /* 14 */ | |
70 COMBA_FORWARD; | |
71 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]); | |
72 COMBA_STORE(C->dp[14]); | |
73 /* 15 */ | |
74 COMBA_FORWARD; | |
75 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]); | |
76 COMBA_STORE(C->dp[15]); | |
77 /* 16 */ | |
78 COMBA_FORWARD; | |
79 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]); | |
80 COMBA_STORE(C->dp[16]); | |
81 /* 17 */ | |
82 COMBA_FORWARD; | |
83 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]); | |
84 COMBA_STORE(C->dp[17]); | |
85 /* 18 */ | |
86 COMBA_FORWARD; | |
87 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]); | |
88 COMBA_STORE(C->dp[18]); | |
89 /* 19 */ | |
90 COMBA_FORWARD; | |
91 MULADD(at[8], at[23]); MULADD(at[9], at[22]); MULADD(at[10], at[21]); MULADD(at[11], at[20]); | |
92 COMBA_STORE(C->dp[19]); | |
93 /* 20 */ | |
94 COMBA_FORWARD; | |
95 MULADD(at[9], at[23]); MULADD(at[10], at[22]); MULADD(at[11], at[21]); | |
96 COMBA_STORE(C->dp[20]); | |
97 /* 21 */ | |
98 COMBA_FORWARD; | |
99 MULADD(at[10], at[23]); MULADD(at[11], at[22]); | |
100 COMBA_STORE(C->dp[21]); | |
101 /* 22 */ | |
102 COMBA_FORWARD; | |
103 MULADD(at[11], at[23]); | |
104 COMBA_STORE(C->dp[22]); | |
105 COMBA_STORE2(C->dp[23]); | |
106 C->used = 24; | |
107 C->sign = A->sign ^ B->sign; | |
108 fp_clamp(C); | |
109 COMBA_FINI; | |
110 } | |
111 #endif |