aboutsummaryrefslogtreecommitdiff
path: root/src/Specific/solinas64_2e174m3/femul.c
diff options
context:
space:
mode:
authorGravatar Jason Gross <jgross@mit.edu>2017-11-02 01:36:04 -0400
committerGravatar Jason Gross <jgross@mit.edu>2017-11-02 01:36:04 -0400
commit5b45aa93efa19d54e50299ff0df8ee46c60f891b (patch)
tree9cea98ab783f557b5cbca0456a4dc3b4498e43c5 /src/Specific/solinas64_2e174m3/femul.c
parente8bda9b779d5762c5868cd09c85142151655d5ca (diff)
Update display logs and c files
Diffstat (limited to 'src/Specific/solinas64_2e174m3/femul.c')
-rw-r--r--src/Specific/solinas64_2e174m3/femul.c41
1 files changed, 41 insertions, 0 deletions
diff --git a/src/Specific/solinas64_2e174m3/femul.c b/src/Specific/solinas64_2e174m3/femul.c
new file mode 100644
index 000000000..4b25d07dd
--- /dev/null
+++ b/src/Specific/solinas64_2e174m3/femul.c
@@ -0,0 +1,41 @@
+#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 x6, uint64_t x7, uint64_t x5, uint64_t x10, uint64_t x11, uint64_t x9)
+{ uint128_t x12 = (((uint128_t)x5 * x10) + (((uint128_t)x7 * x11) + ((uint128_t)x6 * x9)));
+{ uint128_t x13 = ((((uint128_t)x5 * x11) + ((uint128_t)x7 * x9)) + (0x3 * ((uint128_t)x6 * x10)));
+{ uint128_t x14 = (((uint128_t)x5 * x9) + (0x3 * (((uint128_t)x7 * x10) + ((uint128_t)x6 * x11))));
+{ uint128_t x15 = (x14 >> 0x3a);
+{ uint64_t x16 = ((uint64_t)x14 & 0x3ffffffffffffff);
+{ uint128_t x17 = (x15 + x13);
+{ uint64_t x18 = (uint64_t) (x17 >> 0x3a);
+{ uint64_t x19 = ((uint64_t)x17 & 0x3ffffffffffffff);
+{ uint128_t x20 = (x18 + x12);
+{ uint64_t x21 = (uint64_t) (x20 >> 0x3a);
+{ uint64_t x22 = ((uint64_t)x20 & 0x3ffffffffffffff);
+{ uint128_t x23 = (x16 + ((uint128_t)0x3 * x21));
+{ uint64_t x24 = (uint64_t) (x23 >> 0x3a);
+{ uint64_t x25 = ((uint64_t)x23 & 0x3ffffffffffffff);
+{ uint64_t x26 = (x24 + x19);
+{ uint64_t x27 = (x26 >> 0x3a);
+{ uint64_t x28 = (x26 & 0x3ffffffffffffff);
+out[0] = x27 + x22;
+out[1] = x28;
+out[2] = x25;
+}}}}}}}}}}}}}}}}}
+// caller: uint64_t out[3];