1 files changed, 56 insertions, 0 deletions
diff --git a/src/Specific/solinas64_2e243m9/femul.c b/src/Specific/solinas64_2e243m9/femul.c
new file mode 100644
index 000000000..2a5ad5f77
--- /dev/null
+++ b/src/Specific/solinas64_2e243m9/femul.c
@@ -0,0 +1,56 @@
+#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 x12, uint64_t x13, uint64_t x11, uint64_t x9, uint64_t x7, uint64_t x5, uint64_t x22, uint64_t x23, uint64_t x21, uint64_t x19, uint64_t x17, uint64_t x15)
+{  uint128_t x24 = (((uint128_t)x5 * x22) + (((uint128_t)x7 * x23) + (((uint128_t)x9 * x21) + (((uint128_t)x11 * x19) + (((uint128_t)x13 * x17) + ((uint128_t)x12 * x15))))));
+{  uint128_t x25 = ((((uint128_t)x5 * x23) + ((0x2 * ((uint128_t)x7 * x21)) + (((uint128_t)x9 * x19) + ((0x2 * ((uint128_t)x11 * x17)) + ((uint128_t)x13 * x15))))) + (0x9 * (0x2 * ((uint128_t)x12 * x22))));
+{  uint128_t x26 = ((((uint128_t)x5 * x21) + (((uint128_t)x7 * x19) + (((uint128_t)x9 * x17) + ((uint128_t)x11 * x15)))) + (0x9 * (((uint128_t)x13 * x22) + ((uint128_t)x12 * x23))));
+{  uint128_t x27 = ((((uint128_t)x5 * x19) + ((0x2 * ((uint128_t)x7 * x17)) + ((uint128_t)x9 * x15))) + (0x9 * ((0x2 * ((uint128_t)x11 * x22)) + (((uint128_t)x13 * x23) + (0x2 * ((uint128_t)x12 * x21))))));
+{  uint128_t x28 = ((((uint128_t)x5 * x17) + ((uint128_t)x7 * x15)) + (0x9 * (((uint128_t)x9 * x22) + (((uint128_t)x11 * x23) + (((uint128_t)x13 * x21) + ((uint128_t)x12 * x19))))));
+{  uint128_t x29 = (((uint128_t)x5 * x15) + (0x9 * ((0x2 * ((uint128_t)x7 * x22)) + (((uint128_t)x9 * x23) + ((0x2 * ((uint128_t)x11 * x21)) + (((uint128_t)x13 * x19) + (0x2 * ((uint128_t)x12 * x17))))))));
+{  uint64_t x30 = (uint64_t) (x29 >> 0x29);
+{  uint64_t x31 = ((uint64_t)x29 & 0x1ffffffffff);
+{  uint128_t x32 = (x30 + x28);
+{  uint64_t x33 = (uint64_t) (x32 >> 0x28);
+{  uint64_t x34 = ((uint64_t)x32 & 0xffffffffff);
+{  uint128_t x35 = (x33 + x27);
+{  uint64_t x36 = (uint64_t) (x35 >> 0x29);
+{  uint64_t x37 = ((uint64_t)x35 & 0x1ffffffffff);
+{  uint128_t x38 = (x36 + x26);
+{  uint64_t x39 = (uint64_t) (x38 >> 0x28);
+{  uint64_t x40 = ((uint64_t)x38 & 0xffffffffff);
+{  uint128_t x41 = (x39 + x25);
+{  uint64_t x42 = (uint64_t) (x41 >> 0x29);
+{  uint64_t x43 = ((uint64_t)x41 & 0x1ffffffffff);
+{  uint128_t x44 = (x42 + x24);
+{  uint64_t x45 = (uint64_t) (x44 >> 0x28);
+{  uint64_t x46 = ((uint64_t)x44 & 0xffffffffff);
+{  uint64_t x47 = (x31 + (0x9 * x45));
+{  uint64_t x48 = (x47 >> 0x29);
+{  uint64_t x49 = (x47 & 0x1ffffffffff);
+{  uint64_t x50 = (x48 + x34);
+{  uint64_t x51 = (x50 >> 0x28);
+{  uint64_t x52 = (x50 & 0xffffffffff);
+out[0] = x46;
+out[1] = x43;
+out[2] = x40;
+out[3] = x51 + x37;
+out[4] = x52;
+out[5] = x49;
+}}}}}}}}}}}}}}}}}}}}}}}}}}}}}
+// caller: uint64_t out[6];