8 #ifndef BOTAN_MP_MADD_H__ 9 #define BOTAN_MP_MADD_H__ 11 #include <botan/mp_types.h> 15 #if (BOTAN_MP_WORD_BITS != 64) 16 #error The mp_asm64 module requires that BOTAN_MP_WORD_BITS == 64 19 #if defined(BOTAN_TARGET_ARCH_IS_ALPHA) 21 #define BOTAN_WORD_MUL(a,b,z1,z0) do { \ 22 asm("umulh %1,%2,%0" : "=r" (z0) : "r" (a), "r" (b)); \ 26 #elif defined(BOTAN_TARGET_ARCH_IS_IA64) 28 #define BOTAN_WORD_MUL(a,b,z1,z0) do { \ 29 asm("xmpy.hu %0=%1,%2" : "=f" (z0) : "f" (a), "f" (b)); \ 33 #elif defined(BOTAN_TARGET_ARCH_IS_PPC64) 35 #define BOTAN_WORD_MUL(a,b,z1,z0) do { \ 36 asm("mulhdu %0,%1,%2" : "=r" (z0) : "r" (a), "r" (b) : "cc"); \ 40 #elif defined(BOTAN_TARGET_ARCH_IS_MIPS64) 42 #define BOTAN_WORD_MUL(a,b,z1,z0) do { \ 43 typedef unsigned int uint128_t __attribute__((mode(TI))); \ 44 uint128_t r = (uint128_t)a * b; \ 45 z0 = (r >> 64) & 0xFFFFFFFFFFFFFFFF; \ 46 z1 = (r ) & 0xFFFFFFFFFFFFFFFF; \ 56 const size_t MP_HWORD_BITS = BOTAN_MP_WORD_BITS / 2;
57 const word MP_HWORD_MASK = ((word)1 << MP_HWORD_BITS) - 1;
59 const word a_hi = (a >> MP_HWORD_BITS);
60 const word a_lo = (a & MP_HWORD_MASK);
61 const word b_hi = (b >> MP_HWORD_BITS);
62 const word b_lo = (b & MP_HWORD_MASK);
64 word x0 = a_hi * b_hi;
65 word x1 = a_lo * b_hi;
66 word x2 = a_hi * b_lo;
67 word x3 = a_lo * b_lo;
69 x2 += x3 >> (MP_HWORD_BITS);
73 x0 += ((word)1 << MP_HWORD_BITS);
75 *z0 = x0 + (x2 >> MP_HWORD_BITS);
76 *z1 = ((x2 & MP_HWORD_MASK) << MP_HWORD_BITS) + (x3 & MP_HWORD_MASK);
79 #define BOTAN_WORD_MUL(a,b,z1,z0) bigint_2word_mul(a, b, &z1, &z0) void bigint_2word_mul(word a, word b, word *z1, word *z0)
word word_madd3(word a, word b, word c, word *d)
#define BOTAN_WORD_MUL(a, b, z1, z0)
word word_madd2(word a, word b, word *c)