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