1 files changed, 51 insertions, 0 deletions
diff --git a/src/Specific/solinas64_2e285m9/fesquare.c b/src/Specific/solinas64_2e285m9/fesquare.c
new file mode 100644
index 000000000..a15f688d5
--- /dev/null
+++ b/src/Specific/solinas64_2e285m9/fesquare.c
@@ -0,0 +1,51 @@
+#include <stdint.h>
+#include <stdbool.h>
+#include <x86intrin.h>
+#include "liblow.h"
+
+#include "fesquare.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 fesquare(uint64_t* out, uint64_t x7, uint64_t x8, uint64_t x6, uint64_t x4, uint64_t x2)
+{  uint128_t x9 = (((uint128_t)x2 * x7) + (((uint128_t)x4 * x8) + (((uint128_t)x6 * x6) + (((uint128_t)x8 * x4) + ((uint128_t)x7 * x2)))));
+{  uint128_t x10 = ((((uint128_t)x2 * x8) + (((uint128_t)x4 * x6) + (((uint128_t)x6 * x4) + ((uint128_t)x8 * x2)))) + (0x9 * ((uint128_t)x7 * x7)));
+{  uint128_t x11 = ((((uint128_t)x2 * x6) + (((uint128_t)x4 * x4) + ((uint128_t)x6 * x2))) + (0x9 * (((uint128_t)x8 * x7) + ((uint128_t)x7 * x8))));
+{  uint128_t x12 = ((((uint128_t)x2 * x4) + ((uint128_t)x4 * x2)) + (0x9 * (((uint128_t)x6 * x7) + (((uint128_t)x8 * x8) + ((uint128_t)x7 * x6)))));
+{  uint128_t x13 = (((uint128_t)x2 * x2) + (0x9 * (((uint128_t)x4 * x7) + (((uint128_t)x6 * x8) + (((uint128_t)x8 * x6) + ((uint128_t)x7 * x4))))));
+{  uint128_t x14 = (x13 >> 0x39);
+{  uint64_t x15 = ((uint64_t)x13 & 0x1ffffffffffffff);
+{  uint128_t x16 = (x14 + x12);
+{  uint128_t x17 = (x16 >> 0x39);
+{  uint64_t x18 = ((uint64_t)x16 & 0x1ffffffffffffff);
+{  uint128_t x19 = (x17 + x11);
+{  uint128_t x20 = (x19 >> 0x39);
+{  uint64_t x21 = ((uint64_t)x19 & 0x1ffffffffffffff);
+{  uint128_t x22 = (x20 + x10);
+{  uint128_t x23 = (x22 >> 0x39);
+{  uint64_t x24 = ((uint64_t)x22 & 0x1ffffffffffffff);
+{  uint128_t x25 = (x23 + x9);
+{  uint64_t x26 = (uint64_t) (x25 >> 0x39);
+{  uint64_t x27 = ((uint64_t)x25 & 0x1ffffffffffffff);
+{  uint128_t x28 = (x15 + ((uint128_t)0x9 * x26));
+{  uint64_t x29 = (uint64_t) (x28 >> 0x39);
+{  uint64_t x30 = ((uint64_t)x28 & 0x1ffffffffffffff);
+{  uint64_t x31 = (x29 + x18);
+{  uint64_t x32 = (x31 >> 0x39);
+{  uint64_t x33 = (x31 & 0x1ffffffffffffff);
+out[0] = x27;
+out[1] = x24;
+out[2] = x32 + x21;
+out[3] = x33;
+out[4] = x30;
+}}}}}}}}}}}}}}}}}}}}}}}}}
+// caller: uint64_t out[5];