blob: 1283c97e0a442a9475f39ffc5bcc66ab90a574d3 (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
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) + ((0x2 * ((uint128_t)x7 * x19)) + ((0x2 * ((uint128_t)x9 * x17)) + ((0x2 * ((uint128_t)x11 * x15)) + ((uint128_t)x10 * x13)))));
{ uint128_t x21 = ((((uint128_t)x5 * x19) + ((0x2 * ((uint128_t)x7 * x17)) + ((0x2 * ((uint128_t)x9 * x15)) + ((uint128_t)x11 * x13)))) + (0x13 * ((uint128_t)x10 * x18)));
{ uint128_t x22 = ((((uint128_t)x5 * x17) + ((0x2 * ((uint128_t)x7 * x15)) + ((uint128_t)x9 * x13))) + (0x13 * (((uint128_t)x11 * x18) + ((uint128_t)x10 * x19))));
{ uint128_t x23 = ((((uint128_t)x5 * x15) + ((uint128_t)x7 * x13)) + (0x13 * (((uint128_t)x9 * x18) + (((uint128_t)x11 * x19) + ((uint128_t)x10 * x17)))));
{ uint128_t x24 = (((uint128_t)x5 * x13) + (0x13 * ((0x2 * ((uint128_t)x7 * x18)) + ((0x2 * ((uint128_t)x9 * x19)) + ((0x2 * ((uint128_t)x11 * x17)) + (0x2 * ((uint128_t)x10 * x15)))))));
{ uint64_t x25 = (uint64_t) (x24 >> 0x27);
{ uint64_t x26 = ((uint64_t)x24 & 0x7fffffffff);
{ uint128_t x27 = (x25 + x23);
{ uint64_t x28 = (uint64_t) (x27 >> 0x26);
{ uint64_t x29 = ((uint64_t)x27 & 0x3fffffffff);
{ uint128_t x30 = (x28 + x22);
{ uint64_t x31 = (uint64_t) (x30 >> 0x26);
{ uint64_t x32 = ((uint64_t)x30 & 0x3fffffffff);
{ uint128_t x33 = (x31 + x21);
{ uint64_t x34 = (uint64_t) (x33 >> 0x26);
{ uint64_t x35 = ((uint64_t)x33 & 0x3fffffffff);
{ uint128_t x36 = (x34 + x20);
{ uint64_t x37 = (uint64_t) (x36 >> 0x26);
{ uint64_t x38 = ((uint64_t)x36 & 0x3fffffffff);
{ uint64_t x39 = (x26 + (0x13 * x37));
{ uint64_t x40 = (x39 >> 0x27);
{ uint64_t x41 = (x39 & 0x7fffffffff);
{ uint64_t x42 = (x40 + x29);
{ uint64_t x43 = (x42 >> 0x26);
{ uint64_t x44 = (x42 & 0x3fffffffff);
out[0] = x38;
out[1] = x35;
out[2] = x43 + x32;
out[3] = x44;
out[4] = x41;
}}}}}}}}}}}}}}}}}}}}}}}}}
// caller: uint64_t out[5];
|
