diff options
Diffstat (limited to 'src/Specific/solinas64_2e255m765/femul.c')
-rw-r--r-- | src/Specific/solinas64_2e255m765/femul.c | 51 |
1 files changed, 51 insertions, 0 deletions
diff --git a/src/Specific/solinas64_2e255m765/femul.c b/src/Specific/solinas64_2e255m765/femul.c new file mode 100644 index 000000000..010cb8eba --- /dev/null +++ b/src/Specific/solinas64_2e255m765/femul.c @@ -0,0 +1,51 @@ +#include <stdint.h> +#include <stdbool.h> +#include <x86intrin.h> +#include "liblow.h" + +#include "femul.h" + +typedef unsigned int uint128_t __attribute__((mode(TI))); + +#if (defined(__GNUC__) || defined(__GNUG__)) && !(defined(__clang__)||defined(__INTEL_COMPILER)) +// https://gcc.gnu.org/bugzilla/show_bug.cgi?id=81294 +#define _subborrow_u32 __builtin_ia32_sbb_u32 +#define _subborrow_u64 __builtin_ia32_sbb_u64 +#endif + +#undef force_inline +#define force_inline __attribute__((always_inline)) + +void force_inline femul(uint64_t* out, uint64_t x10, uint64_t x11, uint64_t x9, uint64_t x7, uint64_t x5, uint64_t x18, uint64_t x19, uint64_t x17, uint64_t x15, uint64_t x13) +{ uint128_t x20 = (((uint128_t)x5 * x18) + (((uint128_t)x7 * x19) + (((uint128_t)x9 * x17) + (((uint128_t)x11 * x15) + ((uint128_t)x10 * x13))))); +{ uint128_t x21 = ((((uint128_t)x5 * x19) + (((uint128_t)x7 * x17) + (((uint128_t)x9 * x15) + ((uint128_t)x11 * x13)))) + (0x2fd * ((uint128_t)x10 * x18))); +{ uint128_t x22 = ((((uint128_t)x5 * x17) + (((uint128_t)x7 * x15) + ((uint128_t)x9 * x13))) + (0x2fd * (((uint128_t)x11 * x18) + ((uint128_t)x10 * x19)))); +{ uint128_t x23 = ((((uint128_t)x5 * x15) + ((uint128_t)x7 * x13)) + (0x2fd * (((uint128_t)x9 * x18) + (((uint128_t)x11 * x19) + ((uint128_t)x10 * x17))))); +{ uint128_t x24 = (((uint128_t)x5 * x13) + (0x2fd * (((uint128_t)x7 * x18) + (((uint128_t)x9 * x19) + (((uint128_t)x11 * x17) + ((uint128_t)x10 * x15)))))); +{ uint128_t x25 = (x24 >> 0x33); +{ uint64_t x26 = ((uint64_t)x24 & 0x7ffffffffffff); +{ uint128_t x27 = (x25 + x23); +{ uint128_t x28 = (x27 >> 0x33); +{ uint64_t x29 = ((uint64_t)x27 & 0x7ffffffffffff); +{ uint128_t x30 = (x28 + x22); +{ uint128_t x31 = (x30 >> 0x33); +{ uint64_t x32 = ((uint64_t)x30 & 0x7ffffffffffff); +{ uint128_t x33 = (x31 + x21); +{ uint128_t x34 = (x33 >> 0x33); +{ uint64_t x35 = ((uint64_t)x33 & 0x7ffffffffffff); +{ uint128_t x36 = (x34 + x20); +{ uint64_t x37 = (uint64_t) (x36 >> 0x33); +{ uint64_t x38 = ((uint64_t)x36 & 0x7ffffffffffff); +{ uint128_t x39 = (x26 + ((uint128_t)0x2fd * x37)); +{ uint64_t x40 = (uint64_t) (x39 >> 0x33); +{ uint64_t x41 = ((uint64_t)x39 & 0x7ffffffffffff); +{ uint64_t x42 = (x40 + x29); +{ uint64_t x43 = (x42 >> 0x33); +{ uint64_t x44 = (x42 & 0x7ffffffffffff); +out[0] = x38; +out[1] = x35; +out[2] = x43 + x32; +out[3] = x44; +out[4] = x41; +}}}}}}}}}}}}}}}}}}}}}}}}} +// caller: uint64_t out[5]; |