1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
|
/*
* Lowest Level MPI Algorithms
* (C) 1999-2008,2013 Jack Lloyd
* 2006 Luca Piccarreta
*
* Botan is released under the Simplified BSD License (see license.txt)
*/
#ifndef BOTAN_MP_WORD_MULADD_H__
#define BOTAN_MP_WORD_MULADD_H__
#include <botan/mp_types.h>
#include <botan/mul128.h>
namespace Botan {
#if (BOTAN_MP_WORD_BITS == 8)
typedef u16bit dword;
#define BOTAN_HAS_MP_DWORD
#elif (BOTAN_MP_WORD_BITS == 16)
typedef u32bit dword;
#define BOTAN_HAS_MP_DWORD
#elif (BOTAN_MP_WORD_BITS == 32)
typedef u64bit dword;
#define BOTAN_HAS_MP_DWORD
#elif (BOTAN_MP_WORD_BITS == 64)
#if defined(BOTAN_TARGET_HAS_NATIVE_UINT128)
typedef uint128_t dword;
#define BOTAN_HAS_MP_DWORD
#else
// No native 128 bit integer type; use mul64x64_128 instead
#endif
#else
#error BOTAN_MP_WORD_BITS must be 8, 16, 32, or 64
#endif
#if defined(BOTAN_TARGET_ARCH_IS_X86_32) && (BOTAN_MP_WORD_BITS == 32)
#if defined(BOTAN_USE_GCC_INLINE_ASM)
#define BOTAN_MP_USE_X86_32_ASM
#define ASM(x) x "\n\t"
#elif defined(BOTAN_TARGET_COMPILER_IS_MSVC)
#define BOTAN_MP_USE_X86_32_MSVC_ASM
#endif
#elif defined(BOTAN_TARGET_ARCH_IS_X86_64) && (BOTAN_MP_WORD_BITS == 64) && (BOTAN_USE_GCC_INLINE_ASM)
#define BOTAN_MP_USE_X86_64_ASM
#define ASM(x) x "\n\t"
#endif
#if defined(BOTAN_MP_USE_X86_32_ASM) || defined(BOTAN_MP_USE_X86_64_ASM)
#define ASM(x) x "\n\t"
#endif
/*
* Word Multiply/Add
*/
inline word word_madd2(word a, word b, word* c)
{
#if defined(BOTAN_MP_USE_X86_32_ASM)
asm(
ASM("mull %[b]")
ASM("addl %[c],%[a]")
ASM("adcl $0,%[carry]")
: [a]"=a"(a), [b]"=rm"(b), [carry]"=&d"(*c)
: "0"(a), "1"(b), [c]"g"(*c) : "cc");
return a;
#elif defined(BOTAN_MP_USE_X86_64_ASM)
asm(
ASM("mulq %[b]")
ASM("addq %[c],%[a]")
ASM("adcq $0,%[carry]")
: [a]"=a"(a), [b]"=rm"(b), [carry]"=&d"(*c)
: "0"(a), "1"(b), [c]"g"(*c) : "cc");
return a;
#elif defined(BOTAN_HAS_MP_DWORD)
const dword s = static_cast<dword>(a) * b + *c;
*c = static_cast<word>(s >> BOTAN_MP_WORD_BITS);
return static_cast<word>(s);
#else
static_assert(BOTAN_MP_WORD_BITS == 64, "Unexpected word size");
word hi = 0, lo = 0;
mul64x64_128(a, b, &lo, &hi);
lo += *c;
hi += (lo < *c); // carry?
*c = hi;
return lo;
#endif
}
/*
* Word Multiply/Add
*/
inline word word_madd3(word a, word b, word c, word* d)
{
#if defined(BOTAN_MP_USE_X86_32_ASM)
asm(
ASM("mull %[b]")
ASM("addl %[c],%[a]")
ASM("adcl $0,%[carry]")
ASM("addl %[d],%[a]")
ASM("adcl $0,%[carry]")
: [a]"=a"(a), [b]"=rm"(b), [carry]"=&d"(*d)
: "0"(a), "1"(b), [c]"g"(c), [d]"g"(*d) : "cc");
return a;
#elif defined(BOTAN_MP_USE_X86_64_ASM)
asm(
ASM("mulq %[b]")
ASM("addq %[c],%[a]")
ASM("adcq $0,%[carry]")
ASM("addq %[d],%[a]")
ASM("adcq $0,%[carry]")
: [a]"=a"(a), [b]"=rm"(b), [carry]"=&d"(*d)
: "0"(a), "1"(b), [c]"g"(c), [d]"g"(*d) : "cc");
return a;
#elif defined(BOTAN_HAS_MP_DWORD)
const dword s = static_cast<dword>(a) * b + c + *d;
*d = static_cast<word>(s >> BOTAN_MP_WORD_BITS);
return static_cast<word>(s);
#else
static_assert(BOTAN_MP_WORD_BITS == 64, "Unexpected word size");
word hi = 0, lo = 0;
mul64x64_128(a, b, &lo, &hi);
lo += c;
hi += (lo < c); // carry?
lo += *d;
hi += (lo < *d); // carry?
*d = hi;
return lo;
#endif
}
#if defined(ASM)
#undef ASM
#endif
}
#endif
|
