Update display logs and c files

10 files changed, 582 insertions, 0 deletions

diff --git a/src/Specific/montgomery32_2e510m290x2e496m1_16limbs/feadd.c b/src/Specific/montgomery32_2e510m290x2e496m1_16limbs/feadd.c
new file mode 100644
index 000000000..9ad623033
--- /dev/null
+++ b/src/Specific/montgomery32_2e510m290x2e496m1_16limbs/feadd.c
@@ -0,0 +1,100 @@
+static void feadd(uint32_t out[16], const uint32_t in1[16], const uint32_t in2[16]) {
+  { const uint32_t x32 = in1[15];
+  { const uint32_t x33 = in1[14];
+  { const uint32_t x31 = in1[13];
+  { const uint32_t x29 = in1[12];
+  { const uint32_t x27 = in1[11];
+  { const uint32_t x25 = in1[10];
+  { const uint32_t x23 = in1[9];
+  { const uint32_t x21 = in1[8];
+  { const uint32_t x19 = in1[7];
+  { const uint32_t x17 = in1[6];
+  { const uint32_t x15 = in1[5];
+  { const uint32_t x13 = in1[4];
+  { const uint32_t x11 = in1[3];
+  { const uint32_t x9 = in1[2];
+  { const uint32_t x7 = in1[1];
+  { const uint32_t x5 = in1[0];
+  { const uint32_t x62 = in2[15];
+  { const uint32_t x63 = in2[14];
+  { const uint32_t x61 = in2[13];
+  { const uint32_t x59 = in2[12];
+  { const uint32_t x57 = in2[11];
+  { const uint32_t x55 = in2[10];
+  { const uint32_t x53 = in2[9];
+  { const uint32_t x51 = in2[8];
+  { const uint32_t x49 = in2[7];
+  { const uint32_t x47 = in2[6];
+  { const uint32_t x45 = in2[5];
+  { const uint32_t x43 = in2[4];
+  { const uint32_t x41 = in2[3];
+  { const uint32_t x39 = in2[2];
+  { const uint32_t x37 = in2[1];
+  { const uint32_t x35 = in2[0];
+  { uint32_t x65; uint8_t x66 = _addcarryx_u32(0x0, x5, x35, &x65);
+  { uint32_t x68; uint8_t x69 = _addcarryx_u32(x66, x7, x37, &x68);
+  { uint32_t x71; uint8_t x72 = _addcarryx_u32(x69, x9, x39, &x71);
+  { uint32_t x74; uint8_t x75 = _addcarryx_u32(x72, x11, x41, &x74);
+  { uint32_t x77; uint8_t x78 = _addcarryx_u32(x75, x13, x43, &x77);
+  { uint32_t x80; uint8_t x81 = _addcarryx_u32(x78, x15, x45, &x80);
+  { uint32_t x83; uint8_t x84 = _addcarryx_u32(x81, x17, x47, &x83);
+  { uint32_t x86; uint8_t x87 = _addcarryx_u32(x84, x19, x49, &x86);
+  { uint32_t x89; uint8_t x90 = _addcarryx_u32(x87, x21, x51, &x89);
+  { uint32_t x92; uint8_t x93 = _addcarryx_u32(x90, x23, x53, &x92);
+  { uint32_t x95; uint8_t x96 = _addcarryx_u32(x93, x25, x55, &x95);
+  { uint32_t x98; uint8_t x99 = _addcarryx_u32(x96, x27, x57, &x98);
+  { uint32_t x101; uint8_t x102 = _addcarryx_u32(x99, x29, x59, &x101);
+  { uint32_t x104; uint8_t x105 = _addcarryx_u32(x102, x31, x61, &x104);
+  { uint32_t x107; uint8_t x108 = _addcarryx_u32(x105, x33, x63, &x107);
+  { uint32_t x110; uint8_t x111 = _addcarryx_u32(x108, x32, x62, &x110);
+  { uint32_t x113; uint8_t x114 = _subborrow_u32(0x0, x65, 0xffffffff, &x113);
+  { uint32_t x116; uint8_t x117 = _subborrow_u32(x114, x68, 0xffffffff, &x116);
+  { uint32_t x119; uint8_t x120 = _subborrow_u32(x117, x71, 0xffffffff, &x119);
+  { uint32_t x122; uint8_t x123 = _subborrow_u32(x120, x74, 0xffffffff, &x122);
+  { uint32_t x125; uint8_t x126 = _subborrow_u32(x123, x77, 0xffffffff, &x125);
+  { uint32_t x128; uint8_t x129 = _subborrow_u32(x126, x80, 0xffffffff, &x128);
+  { uint32_t x131; uint8_t x132 = _subborrow_u32(x129, x83, 0xffffffff, &x131);
+  { uint32_t x134; uint8_t x135 = _subborrow_u32(x132, x86, 0xffffffff, &x134);
+  { uint32_t x137; uint8_t x138 = _subborrow_u32(x135, x89, 0xffffffff, &x137);
+  { uint32_t x140; uint8_t x141 = _subborrow_u32(x138, x92, 0xffffffff, &x140);
+  { uint32_t x143; uint8_t x144 = _subborrow_u32(x141, x95, 0xffffffff, &x143);
+  { uint32_t x146; uint8_t x147 = _subborrow_u32(x144, x98, 0xffffffff, &x146);
+  { uint32_t x149; uint8_t x150 = _subborrow_u32(x147, x101, 0xffffffff, &x149);
+  { uint32_t x152; uint8_t x153 = _subborrow_u32(x150, x104, 0xffffffff, &x152);
+  { uint32_t x155; uint8_t x156 = _subborrow_u32(x153, x107, 0xffffffff, &x155);
+  { uint32_t x158; uint8_t x159 = _subborrow_u32(x156, x110, 0x3eddffff, &x158);
+  { uint32_t _; uint8_t x162 = _subborrow_u32(x159, x111, 0x0, &_);
+  { uint32_t x163 = cmovznz32(x162, x158, x110);
+  { uint32_t x164 = cmovznz32(x162, x155, x107);
+  { uint32_t x165 = cmovznz32(x162, x152, x104);
+  { uint32_t x166 = cmovznz32(x162, x149, x101);
+  { uint32_t x167 = cmovznz32(x162, x146, x98);
+  { uint32_t x168 = cmovznz32(x162, x143, x95);
+  { uint32_t x169 = cmovznz32(x162, x140, x92);
+  { uint32_t x170 = cmovznz32(x162, x137, x89);
+  { uint32_t x171 = cmovznz32(x162, x134, x86);
+  { uint32_t x172 = cmovznz32(x162, x131, x83);
+  { uint32_t x173 = cmovznz32(x162, x128, x80);
+  { uint32_t x174 = cmovznz32(x162, x125, x77);
+  { uint32_t x175 = cmovznz32(x162, x122, x74);
+  { uint32_t x176 = cmovznz32(x162, x119, x71);
+  { uint32_t x177 = cmovznz32(x162, x116, x68);
+  { uint32_t x178 = cmovznz32(x162, x113, x65);
+  out[0] = x178;
+  out[1] = x177;
+  out[2] = x176;
+  out[3] = x175;
+  out[4] = x174;
+  out[5] = x173;
+  out[6] = x172;
+  out[7] = x171;
+  out[8] = x170;
+  out[9] = x169;
+  out[10] = x168;
+  out[11] = x167;
+  out[12] = x166;
+  out[13] = x165;
+  out[14] = x164;
+  out[15] = x163;
+  }}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}
+}
diff --git a/src/Specific/montgomery32_2e510m290x2e496m1_16limbs/feaddDisplay.log b/src/Specific/montgomery32_2e510m290x2e496m1_16limbs/feaddDisplay.log
new file mode 100644
index 000000000..36c81d199
--- /dev/null
+++ b/src/Specific/montgomery32_2e510m290x2e496m1_16limbs/feaddDisplay.log
@@ -0,0 +1,56 @@
+λ x x0 : word32 * word32 * word32 * word32 * word32 * word32 * word32 * word32 * word32 * word32 * word32 * word32 * word32 * word32 * word32 * word32,
+Interp-η
+(λ var : Syntax.base_type → Type,
+ λ '(x32, x33, x31, x29, x27, x25, x23, x21, x19, x17, x15, x13, x11, x9, x7, x5, (x62, x63, x61, x59, x57, x55, x53, x51, x49, x47, x45, x43, x41, x39, x37, x35))%core,
+ uint32_t x65, uint8_t x66 = addcarryx_u32(0x0, x5, x35);
+ uint32_t x68, uint8_t x69 = addcarryx_u32(x66, x7, x37);
+ uint32_t x71, uint8_t x72 = addcarryx_u32(x69, x9, x39);
+ uint32_t x74, uint8_t x75 = addcarryx_u32(x72, x11, x41);
+ uint32_t x77, uint8_t x78 = addcarryx_u32(x75, x13, x43);
+ uint32_t x80, uint8_t x81 = addcarryx_u32(x78, x15, x45);
+ uint32_t x83, uint8_t x84 = addcarryx_u32(x81, x17, x47);
+ uint32_t x86, uint8_t x87 = addcarryx_u32(x84, x19, x49);
+ uint32_t x89, uint8_t x90 = addcarryx_u32(x87, x21, x51);
+ uint32_t x92, uint8_t x93 = addcarryx_u32(x90, x23, x53);
+ uint32_t x95, uint8_t x96 = addcarryx_u32(x93, x25, x55);
+ uint32_t x98, uint8_t x99 = addcarryx_u32(x96, x27, x57);
+ uint32_t x101, uint8_t x102 = addcarryx_u32(x99, x29, x59);
+ uint32_t x104, uint8_t x105 = addcarryx_u32(x102, x31, x61);
+ uint32_t x107, uint8_t x108 = addcarryx_u32(x105, x33, x63);
+ uint32_t x110, uint8_t x111 = addcarryx_u32(x108, x32, x62);
+ uint32_t x113, uint8_t x114 = subborrow_u32(0x0, x65, 0xffffffff);
+ uint32_t x116, uint8_t x117 = subborrow_u32(x114, x68, 0xffffffff);
+ uint32_t x119, uint8_t x120 = subborrow_u32(x117, x71, 0xffffffff);
+ uint32_t x122, uint8_t x123 = subborrow_u32(x120, x74, 0xffffffff);
+ uint32_t x125, uint8_t x126 = subborrow_u32(x123, x77, 0xffffffff);
+ uint32_t x128, uint8_t x129 = subborrow_u32(x126, x80, 0xffffffff);
+ uint32_t x131, uint8_t x132 = subborrow_u32(x129, x83, 0xffffffff);
+ uint32_t x134, uint8_t x135 = subborrow_u32(x132, x86, 0xffffffff);
+ uint32_t x137, uint8_t x138 = subborrow_u32(x135, x89, 0xffffffff);
+ uint32_t x140, uint8_t x141 = subborrow_u32(x138, x92, 0xffffffff);
+ uint32_t x143, uint8_t x144 = subborrow_u32(x141, x95, 0xffffffff);
+ uint32_t x146, uint8_t x147 = subborrow_u32(x144, x98, 0xffffffff);
+ uint32_t x149, uint8_t x150 = subborrow_u32(x147, x101, 0xffffffff);
+ uint32_t x152, uint8_t x153 = subborrow_u32(x150, x104, 0xffffffff);
+ uint32_t x155, uint8_t x156 = subborrow_u32(x153, x107, 0xffffffff);
+ uint32_t x158, uint8_t x159 = subborrow_u32(x156, x110, 0x3eddffff);
+ uint32_t _, uint8_t x162 = subborrow_u32(x159, x111, 0x0);
+ uint32_t x163 = cmovznz32(x162, x158, x110);
+ uint32_t x164 = cmovznz32(x162, x155, x107);
+ uint32_t x165 = cmovznz32(x162, x152, x104);
+ uint32_t x166 = cmovznz32(x162, x149, x101);
+ uint32_t x167 = cmovznz32(x162, x146, x98);
+ uint32_t x168 = cmovznz32(x162, x143, x95);
+ uint32_t x169 = cmovznz32(x162, x140, x92);
+ uint32_t x170 = cmovznz32(x162, x137, x89);
+ uint32_t x171 = cmovznz32(x162, x134, x86);
+ uint32_t x172 = cmovznz32(x162, x131, x83);
+ uint32_t x173 = cmovznz32(x162, x128, x80);
+ uint32_t x174 = cmovznz32(x162, x125, x77);
+ uint32_t x175 = cmovznz32(x162, x122, x74);
+ uint32_t x176 = cmovznz32(x162, x119, x71);
+ uint32_t x177 = cmovznz32(x162, x116, x68);
+ uint32_t x178 = cmovznz32(x162, x113, x65);
+ return (x163, x164, x165, x166, x167, x168, x169, x170, x171, x172, x173, x174, x175, x176, x177, x178))
+(x, x0)%core
+     : word32 * word32 * word32 * word32 * word32 * word32 * word32 * word32 * word32 * word32 * word32 * word32 * word32 * word32 * word32 * word32 → word32 * word32 * word32 * word32 * word32 * word32 * word32 * word32 * word32 * word32 * word32 * word32 * word32 * word32 * word32 * word32 → ReturnType (uint32_t * uint32_t * uint32_t * uint32_t * uint32_t * uint32_t * uint32_t * uint32_t * uint32_t * uint32_t * uint32_t * uint32_t * uint32_t * uint32_t * uint32_t * uint32_t)
diff --git a/src/Specific/montgomery32_2e510m290x2e496m1_16limbs/fenz.c b/src/Specific/montgomery32_2e510m290x2e496m1_16limbs/fenz.c
new file mode 100644
index 000000000..5c5d21b0e
--- /dev/null
+++ b/src/Specific/montgomery32_2e510m290x2e496m1_16limbs/fenz.c
@@ -0,0 +1,35 @@
+static void fenz(ReturnType uint32_t out[1], const uint32_t in1[16]) {
+  { const uint32_t x29 = in1[15];
+  { const uint32_t x30 = in1[14];
+  { const uint32_t x28 = in1[13];
+  { const uint32_t x26 = in1[12];
+  { const uint32_t x24 = in1[11];
+  { const uint32_t x22 = in1[10];
+  { const uint32_t x20 = in1[9];
+  { const uint32_t x18 = in1[8];
+  { const uint32_t x16 = in1[7];
+  { const uint32_t x14 = in1[6];
+  { const uint32_t x12 = in1[5];
+  { const uint32_t x10 = in1[4];
+  { const uint32_t x8 = in1[3];
+  { const uint32_t x6 = in1[2];
+  { const uint32_t x4 = in1[1];
+  { const uint32_t x2 = in1[0];
+  { uint32_t x31 = (x30 | x29);
+  { uint32_t x32 = (x28 | x31);
+  { uint32_t x33 = (x26 | x32);
+  { uint32_t x34 = (x24 | x33);
+  { uint32_t x35 = (x22 | x34);
+  { uint32_t x36 = (x20 | x35);
+  { uint32_t x37 = (x18 | x36);
+  { uint32_t x38 = (x16 | x37);
+  { uint32_t x39 = (x14 | x38);
+  { uint32_t x40 = (x12 | x39);
+  { uint32_t x41 = (x10 | x40);
+  { uint32_t x42 = (x8 | x41);
+  { uint32_t x43 = (x6 | x42);
+  { uint32_t x44 = (x4 | x43);
+  { uint32_t x45 = (x2 | x44);
+  out[0] = x45;
+  }}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}
+}
diff --git a/src/Specific/montgomery32_2e510m290x2e496m1_16limbs/fenzDisplay.log b/src/Specific/montgomery32_2e510m290x2e496m1_16limbs/fenzDisplay.log
new file mode 100644
index 000000000..eb0b6ea33
--- /dev/null
+++ b/src/Specific/montgomery32_2e510m290x2e496m1_16limbs/fenzDisplay.log
@@ -0,0 +1,22 @@
+λ x : word32 * word32 * word32 * word32 * word32 * word32 * word32 * word32 * word32 * word32 * word32 * word32 * word32 * word32 * word32 * word32,
+Interp-η
+(λ var : Syntax.base_type → Type,
+ λ '(x29, x30, x28, x26, x24, x22, x20, x18, x16, x14, x12, x10, x8, x6, x4, x2)%core,
+ uint32_t x31 = (x30 | x29);
+ uint32_t x32 = (x28 | x31);
+ uint32_t x33 = (x26 | x32);
+ uint32_t x34 = (x24 | x33);
+ uint32_t x35 = (x22 | x34);
+ uint32_t x36 = (x20 | x35);
+ uint32_t x37 = (x18 | x36);
+ uint32_t x38 = (x16 | x37);
+ uint32_t x39 = (x14 | x38);
+ uint32_t x40 = (x12 | x39);
+ uint32_t x41 = (x10 | x40);
+ uint32_t x42 = (x8 | x41);
+ uint32_t x43 = (x6 | x42);
+ uint32_t x44 = (x4 | x43);
+ uint32_t x45 = (x2 | x44);
+ return x45)
+x
+     : word32 * word32 * word32 * word32 * word32 * word32 * word32 * word32 * word32 * word32 * word32 * word32 * word32 * word32 * word32 * word32 → ReturnType uint32_t
diff --git a/src/Specific/montgomery32_2e512m491x2e496m1_16limbs/feadd.c b/src/Specific/montgomery32_2e512m491x2e496m1_16limbs/feadd.c
new file mode 100644
index 000000000..b107d1c5e
--- /dev/null
+++ b/src/Specific/montgomery32_2e512m491x2e496m1_16limbs/feadd.c
@@ -0,0 +1,100 @@
+static void feadd(uint32_t out[16], const uint32_t in1[16], const uint32_t in2[16]) {
+  { const uint32_t x32 = in1[15];
+  { const uint32_t x33 = in1[14];
+  { const uint32_t x31 = in1[13];
+  { const uint32_t x29 = in1[12];
+  { const uint32_t x27 = in1[11];
+  { const uint32_t x25 = in1[10];
+  { const uint32_t x23 = in1[9];
+  { const uint32_t x21 = in1[8];
+  { const uint32_t x19 = in1[7];
+  { const uint32_t x17 = in1[6];
+  { const uint32_t x15 = in1[5];
+  { const uint32_t x13 = in1[4];
+  { const uint32_t x11 = in1[3];
+  { const uint32_t x9 = in1[2];
+  { const uint32_t x7 = in1[1];
+  { const uint32_t x5 = in1[0];
+  { const uint32_t x62 = in2[15];
+  { const uint32_t x63 = in2[14];
+  { const uint32_t x61 = in2[13];
+  { const uint32_t x59 = in2[12];
+  { const uint32_t x57 = in2[11];
+  { const uint32_t x55 = in2[10];
+  { const uint32_t x53 = in2[9];
+  { const uint32_t x51 = in2[8];
+  { const uint32_t x49 = in2[7];
+  { const uint32_t x47 = in2[6];
+  { const uint32_t x45 = in2[5];
+  { const uint32_t x43 = in2[4];
+  { const uint32_t x41 = in2[3];
+  { const uint32_t x39 = in2[2];
+  { const uint32_t x37 = in2[1];
+  { const uint32_t x35 = in2[0];
+  { uint32_t x65; uint8_t x66 = _addcarryx_u32(0x0, x5, x35, &x65);
+  { uint32_t x68; uint8_t x69 = _addcarryx_u32(x66, x7, x37, &x68);
+  { uint32_t x71; uint8_t x72 = _addcarryx_u32(x69, x9, x39, &x71);
+  { uint32_t x74; uint8_t x75 = _addcarryx_u32(x72, x11, x41, &x74);
+  { uint32_t x77; uint8_t x78 = _addcarryx_u32(x75, x13, x43, &x77);
+  { uint32_t x80; uint8_t x81 = _addcarryx_u32(x78, x15, x45, &x80);
+  { uint32_t x83; uint8_t x84 = _addcarryx_u32(x81, x17, x47, &x83);
+  { uint32_t x86; uint8_t x87 = _addcarryx_u32(x84, x19, x49, &x86);
+  { uint32_t x89; uint8_t x90 = _addcarryx_u32(x87, x21, x51, &x89);
+  { uint32_t x92; uint8_t x93 = _addcarryx_u32(x90, x23, x53, &x92);
+  { uint32_t x95; uint8_t x96 = _addcarryx_u32(x93, x25, x55, &x95);
+  { uint32_t x98; uint8_t x99 = _addcarryx_u32(x96, x27, x57, &x98);
+  { uint32_t x101; uint8_t x102 = _addcarryx_u32(x99, x29, x59, &x101);
+  { uint32_t x104; uint8_t x105 = _addcarryx_u32(x102, x31, x61, &x104);
+  { uint32_t x107; uint8_t x108 = _addcarryx_u32(x105, x33, x63, &x107);
+  { uint32_t x110; uint8_t x111 = _addcarryx_u32(x108, x32, x62, &x110);
+  { uint32_t x113; uint8_t x114 = _subborrow_u32(0x0, x65, 0xffffffff, &x113);
+  { uint32_t x116; uint8_t x117 = _subborrow_u32(x114, x68, 0xffffffff, &x116);
+  { uint32_t x119; uint8_t x120 = _subborrow_u32(x117, x71, 0xffffffff, &x119);
+  { uint32_t x122; uint8_t x123 = _subborrow_u32(x120, x74, 0xffffffff, &x122);
+  { uint32_t x125; uint8_t x126 = _subborrow_u32(x123, x77, 0xffffffff, &x125);
+  { uint32_t x128; uint8_t x129 = _subborrow_u32(x126, x80, 0xffffffff, &x128);
+  { uint32_t x131; uint8_t x132 = _subborrow_u32(x129, x83, 0xffffffff, &x131);
+  { uint32_t x134; uint8_t x135 = _subborrow_u32(x132, x86, 0xffffffff, &x134);
+  { uint32_t x137; uint8_t x138 = _subborrow_u32(x135, x89, 0xffffffff, &x137);
+  { uint32_t x140; uint8_t x141 = _subborrow_u32(x138, x92, 0xffffffff, &x140);
+  { uint32_t x143; uint8_t x144 = _subborrow_u32(x141, x95, 0xffffffff, &x143);
+  { uint32_t x146; uint8_t x147 = _subborrow_u32(x144, x98, 0xffffffff, &x146);
+  { uint32_t x149; uint8_t x150 = _subborrow_u32(x147, x101, 0xffffffff, &x149);
+  { uint32_t x152; uint8_t x153 = _subborrow_u32(x150, x104, 0xffffffff, &x152);
+  { uint32_t x155; uint8_t x156 = _subborrow_u32(x153, x107, 0xffffffff, &x155);
+  { uint32_t x158; uint8_t x159 = _subborrow_u32(x156, x110, 0xfe14ffff, &x158);
+  { uint32_t _; uint8_t x162 = _subborrow_u32(x159, x111, 0x0, &_);
+  { uint32_t x163 = cmovznz32(x162, x158, x110);
+  { uint32_t x164 = cmovznz32(x162, x155, x107);
+  { uint32_t x165 = cmovznz32(x162, x152, x104);
+  { uint32_t x166 = cmovznz32(x162, x149, x101);
+  { uint32_t x167 = cmovznz32(x162, x146, x98);
+  { uint32_t x168 = cmovznz32(x162, x143, x95);
+  { uint32_t x169 = cmovznz32(x162, x140, x92);
+  { uint32_t x170 = cmovznz32(x162, x137, x89);
+  { uint32_t x171 = cmovznz32(x162, x134, x86);
+  { uint32_t x172 = cmovznz32(x162, x131, x83);
+  { uint32_t x173 = cmovznz32(x162, x128, x80);
+  { uint32_t x174 = cmovznz32(x162, x125, x77);
+  { uint32_t x175 = cmovznz32(x162, x122, x74);
+  { uint32_t x176 = cmovznz32(x162, x119, x71);
+  { uint32_t x177 = cmovznz32(x162, x116, x68);
+  { uint32_t x178 = cmovznz32(x162, x113, x65);
+  out[0] = x178;
+  out[1] = x177;
+  out[2] = x176;
+  out[3] = x175;
+  out[4] = x174;
+  out[5] = x173;
+  out[6] = x172;
+  out[7] = x171;
+  out[8] = x170;
+  out[9] = x169;
+  out[10] = x168;
+  out[11] = x167;
+  out[12] = x166;
+  out[13] = x165;
+  out[14] = x164;
+  out[15] = x163;
+  }}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}
+}
diff --git a/src/Specific/montgomery32_2e512m491x2e496m1_16limbs/feaddDisplay.log b/src/Specific/montgomery32_2e512m491x2e496m1_16limbs/feaddDisplay.log
new file mode 100644
index 000000000..e0fa0330d
--- /dev/null
+++ b/src/Specific/montgomery32_2e512m491x2e496m1_16limbs/feaddDisplay.log
@@ -0,0 +1,56 @@
+λ x x0 : word32 * word32 * word32 * word32 * word32 * word32 * word32 * word32 * word32 * word32 * word32 * word32 * word32 * word32 * word32 * word32,
+Interp-η
+(λ var : Syntax.base_type → Type,
+ λ '(x32, x33, x31, x29, x27, x25, x23, x21, x19, x17, x15, x13, x11, x9, x7, x5, (x62, x63, x61, x59, x57, x55, x53, x51, x49, x47, x45, x43, x41, x39, x37, x35))%core,
+ uint32_t x65, uint8_t x66 = addcarryx_u32(0x0, x5, x35);
+ uint32_t x68, uint8_t x69 = addcarryx_u32(x66, x7, x37);
+ uint32_t x71, uint8_t x72 = addcarryx_u32(x69, x9, x39);
+ uint32_t x74, uint8_t x75 = addcarryx_u32(x72, x11, x41);
+ uint32_t x77, uint8_t x78 = addcarryx_u32(x75, x13, x43);
+ uint32_t x80, uint8_t x81 = addcarryx_u32(x78, x15, x45);
+ uint32_t x83, uint8_t x84 = addcarryx_u32(x81, x17, x47);
+ uint32_t x86, uint8_t x87 = addcarryx_u32(x84, x19, x49);
+ uint32_t x89, uint8_t x90 = addcarryx_u32(x87, x21, x51);
+ uint32_t x92, uint8_t x93 = addcarryx_u32(x90, x23, x53);
+ uint32_t x95, uint8_t x96 = addcarryx_u32(x93, x25, x55);
+ uint32_t x98, uint8_t x99 = addcarryx_u32(x96, x27, x57);
+ uint32_t x101, uint8_t x102 = addcarryx_u32(x99, x29, x59);
+ uint32_t x104, uint8_t x105 = addcarryx_u32(x102, x31, x61);
+ uint32_t x107, uint8_t x108 = addcarryx_u32(x105, x33, x63);
+ uint32_t x110, uint8_t x111 = addcarryx_u32(x108, x32, x62);
+ uint32_t x113, uint8_t x114 = subborrow_u32(0x0, x65, 0xffffffff);
+ uint32_t x116, uint8_t x117 = subborrow_u32(x114, x68, 0xffffffff);
+ uint32_t x119, uint8_t x120 = subborrow_u32(x117, x71, 0xffffffff);
+ uint32_t x122, uint8_t x123 = subborrow_u32(x120, x74, 0xffffffff);
+ uint32_t x125, uint8_t x126 = subborrow_u32(x123, x77, 0xffffffff);
+ uint32_t x128, uint8_t x129 = subborrow_u32(x126, x80, 0xffffffff);
+ uint32_t x131, uint8_t x132 = subborrow_u32(x129, x83, 0xffffffff);
+ uint32_t x134, uint8_t x135 = subborrow_u32(x132, x86, 0xffffffff);
+ uint32_t x137, uint8_t x138 = subborrow_u32(x135, x89, 0xffffffff);
+ uint32_t x140, uint8_t x141 = subborrow_u32(x138, x92, 0xffffffff);
+ uint32_t x143, uint8_t x144 = subborrow_u32(x141, x95, 0xffffffff);
+ uint32_t x146, uint8_t x147 = subborrow_u32(x144, x98, 0xffffffff);
+ uint32_t x149, uint8_t x150 = subborrow_u32(x147, x101, 0xffffffff);
+ uint32_t x152, uint8_t x153 = subborrow_u32(x150, x104, 0xffffffff);
+ uint32_t x155, uint8_t x156 = subborrow_u32(x153, x107, 0xffffffff);
+ uint32_t x158, uint8_t x159 = subborrow_u32(x156, x110, 0xfe14ffff);
+ uint32_t _, uint8_t x162 = subborrow_u32(x159, x111, 0x0);
+ uint32_t x163 = cmovznz32(x162, x158, x110);
+ uint32_t x164 = cmovznz32(x162, x155, x107);
+ uint32_t x165 = cmovznz32(x162, x152, x104);
+ uint32_t x166 = cmovznz32(x162, x149, x101);
+ uint32_t x167 = cmovznz32(x162, x146, x98);
+ uint32_t x168 = cmovznz32(x162, x143, x95);
+ uint32_t x169 = cmovznz32(x162, x140, x92);
+ uint32_t x170 = cmovznz32(x162, x137, x89);
+ uint32_t x171 = cmovznz32(x162, x134, x86);
+ uint32_t x172 = cmovznz32(x162, x131, x83);
+ uint32_t x173 = cmovznz32(x162, x128, x80);
+ uint32_t x174 = cmovznz32(x162, x125, x77);
+ uint32_t x175 = cmovznz32(x162, x122, x74);
+ uint32_t x176 = cmovznz32(x162, x119, x71);
+ uint32_t x177 = cmovznz32(x162, x116, x68);
+ uint32_t x178 = cmovznz32(x162, x113, x65);
+ return (x163, x164, x165, x166, x167, x168, x169, x170, x171, x172, x173, x174, x175, x176, x177, x178))
+(x, x0)%core
+     : word32 * word32 * word32 * word32 * word32 * word32 * word32 * word32 * word32 * word32 * word32 * word32 * word32 * word32 * word32 * word32 → word32 * word32 * word32 * word32 * word32 * word32 * word32 * word32 * word32 * word32 * word32 * word32 * word32 * word32 * word32 * word32 → ReturnType (uint32_t * uint32_t * uint32_t * uint32_t * uint32_t * uint32_t * uint32_t * uint32_t * uint32_t * uint32_t * uint32_t * uint32_t * uint32_t * uint32_t * uint32_t * uint32_t)
diff --git a/src/Specific/montgomery32_2e512m491x2e496m1_16limbs/fenz.c b/src/Specific/montgomery32_2e512m491x2e496m1_16limbs/fenz.c
new file mode 100644
index 000000000..5c5d21b0e
--- /dev/null
+++ b/src/Specific/montgomery32_2e512m491x2e496m1_16limbs/fenz.c
@@ -0,0 +1,35 @@
+static void fenz(ReturnType uint32_t out[1], const uint32_t in1[16]) {
+  { const uint32_t x29 = in1[15];
+  { const uint32_t x30 = in1[14];
+  { const uint32_t x28 = in1[13];
+  { const uint32_t x26 = in1[12];
+  { const uint32_t x24 = in1[11];
+  { const uint32_t x22 = in1[10];
+  { const uint32_t x20 = in1[9];
+  { const uint32_t x18 = in1[8];
+  { const uint32_t x16 = in1[7];
+  { const uint32_t x14 = in1[6];
+  { const uint32_t x12 = in1[5];
+  { const uint32_t x10 = in1[4];
+  { const uint32_t x8 = in1[3];
+  { const uint32_t x6 = in1[2];
+  { const uint32_t x4 = in1[1];
+  { const uint32_t x2 = in1[0];
+  { uint32_t x31 = (x30 | x29);
+  { uint32_t x32 = (x28 | x31);
+  { uint32_t x33 = (x26 | x32);
+  { uint32_t x34 = (x24 | x33);
+  { uint32_t x35 = (x22 | x34);
+  { uint32_t x36 = (x20 | x35);
+  { uint32_t x37 = (x18 | x36);
+  { uint32_t x38 = (x16 | x37);
+  { uint32_t x39 = (x14 | x38);
+  { uint32_t x40 = (x12 | x39);
+  { uint32_t x41 = (x10 | x40);
+  { uint32_t x42 = (x8 | x41);
+  { uint32_t x43 = (x6 | x42);
+  { uint32_t x44 = (x4 | x43);
+  { uint32_t x45 = (x2 | x44);
+  out[0] = x45;
+  }}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}
+}
diff --git a/src/Specific/montgomery32_2e512m491x2e496m1_16limbs/fenzDisplay.log b/src/Specific/montgomery32_2e512m491x2e496m1_16limbs/fenzDisplay.log
new file mode 100644
index 000000000..eb0b6ea33
--- /dev/null
+++ b/src/Specific/montgomery32_2e512m491x2e496m1_16limbs/fenzDisplay.log
@@ -0,0 +1,22 @@
+λ x : word32 * word32 * word32 * word32 * word32 * word32 * word32 * word32 * word32 * word32 * word32 * word32 * word32 * word32 * word32 * word32,
+Interp-η
+(λ var : Syntax.base_type → Type,
+ λ '(x29, x30, x28, x26, x24, x22, x20, x18, x16, x14, x12, x10, x8, x6, x4, x2)%core,
+ uint32_t x31 = (x30 | x29);
+ uint32_t x32 = (x28 | x31);
+ uint32_t x33 = (x26 | x32);
+ uint32_t x34 = (x24 | x33);
+ uint32_t x35 = (x22 | x34);
+ uint32_t x36 = (x20 | x35);
+ uint32_t x37 = (x18 | x36);
+ uint32_t x38 = (x16 | x37);
+ uint32_t x39 = (x14 | x38);
+ uint32_t x40 = (x12 | x39);
+ uint32_t x41 = (x10 | x40);
+ uint32_t x42 = (x8 | x41);
+ uint32_t x43 = (x6 | x42);
+ uint32_t x44 = (x4 | x43);
+ uint32_t x45 = (x2 | x44);
+ return x45)
+x
+     : word32 * word32 * word32 * word32 * word32 * word32 * word32 * word32 * word32 * word32 * word32 * word32 * word32 * word32 * word32 * word32 → ReturnType uint32_t
diff --git a/src/Specific/montgomery32_2e512m569_16limbs/feadd.c b/src/Specific/montgomery32_2e512m569_16limbs/feadd.c
new file mode 100644
index 000000000..60bca6599
--- /dev/null
+++ b/src/Specific/montgomery32_2e512m569_16limbs/feadd.c
@@ -0,0 +1,100 @@
+static void feadd(uint32_t out[16], const uint32_t in1[16], const uint32_t in2[16]) {
+  { const uint32_t x32 = in1[15];
+  { const uint32_t x33 = in1[14];
+  { const uint32_t x31 = in1[13];
+  { const uint32_t x29 = in1[12];
+  { const uint32_t x27 = in1[11];
+  { const uint32_t x25 = in1[10];
+  { const uint32_t x23 = in1[9];
+  { const uint32_t x21 = in1[8];
+  { const uint32_t x19 = in1[7];
+  { const uint32_t x17 = in1[6];
+  { const uint32_t x15 = in1[5];
+  { const uint32_t x13 = in1[4];
+  { const uint32_t x11 = in1[3];
+  { const uint32_t x9 = in1[2];
+  { const uint32_t x7 = in1[1];
+  { const uint32_t x5 = in1[0];
+  { const uint32_t x62 = in2[15];
+  { const uint32_t x63 = in2[14];
+  { const uint32_t x61 = in2[13];
+  { const uint32_t x59 = in2[12];
+  { const uint32_t x57 = in2[11];
+  { const uint32_t x55 = in2[10];
+  { const uint32_t x53 = in2[9];
+  { const uint32_t x51 = in2[8];
+  { const uint32_t x49 = in2[7];
+  { const uint32_t x47 = in2[6];
+  { const uint32_t x45 = in2[5];
+  { const uint32_t x43 = in2[4];
+  { const uint32_t x41 = in2[3];
+  { const uint32_t x39 = in2[2];
+  { const uint32_t x37 = in2[1];
+  { const uint32_t x35 = in2[0];
+  { uint32_t x65; uint8_t x66 = _addcarryx_u32(0x0, x5, x35, &x65);
+  { uint32_t x68; uint8_t x69 = _addcarryx_u32(x66, x7, x37, &x68);
+  { uint32_t x71; uint8_t x72 = _addcarryx_u32(x69, x9, x39, &x71);
+  { uint32_t x74; uint8_t x75 = _addcarryx_u32(x72, x11, x41, &x74);
+  { uint32_t x77; uint8_t x78 = _addcarryx_u32(x75, x13, x43, &x77);
+  { uint32_t x80; uint8_t x81 = _addcarryx_u32(x78, x15, x45, &x80);
+  { uint32_t x83; uint8_t x84 = _addcarryx_u32(x81, x17, x47, &x83);
+  { uint32_t x86; uint8_t x87 = _addcarryx_u32(x84, x19, x49, &x86);
+  { uint32_t x89; uint8_t x90 = _addcarryx_u32(x87, x21, x51, &x89);
+  { uint32_t x92; uint8_t x93 = _addcarryx_u32(x90, x23, x53, &x92);
+  { uint32_t x95; uint8_t x96 = _addcarryx_u32(x93, x25, x55, &x95);
+  { uint32_t x98; uint8_t x99 = _addcarryx_u32(x96, x27, x57, &x98);
+  { uint32_t x101; uint8_t x102 = _addcarryx_u32(x99, x29, x59, &x101);
+  { uint32_t x104; uint8_t x105 = _addcarryx_u32(x102, x31, x61, &x104);
+  { uint32_t x107; uint8_t x108 = _addcarryx_u32(x105, x33, x63, &x107);
+  { uint32_t x110; uint8_t x111 = _addcarryx_u32(x108, x32, x62, &x110);
+  { uint32_t x113; uint8_t x114 = _subborrow_u32(0x0, x65, 0xfffffdc7, &x113);
+  { uint32_t x116; uint8_t x117 = _subborrow_u32(x114, x68, 0xffffffff, &x116);
+  { uint32_t x119; uint8_t x120 = _subborrow_u32(x117, x71, 0xffffffff, &x119);
+  { uint32_t x122; uint8_t x123 = _subborrow_u32(x120, x74, 0xffffffff, &x122);
+  { uint32_t x125; uint8_t x126 = _subborrow_u32(x123, x77, 0xffffffff, &x125);
+  { uint32_t x128; uint8_t x129 = _subborrow_u32(x126, x80, 0xffffffff, &x128);
+  { uint32_t x131; uint8_t x132 = _subborrow_u32(x129, x83, 0xffffffff, &x131);
+  { uint32_t x134; uint8_t x135 = _subborrow_u32(x132, x86, 0xffffffff, &x134);
+  { uint32_t x137; uint8_t x138 = _subborrow_u32(x135, x89, 0xffffffff, &x137);
+  { uint32_t x140; uint8_t x141 = _subborrow_u32(x138, x92, 0xffffffff, &x140);
+  { uint32_t x143; uint8_t x144 = _subborrow_u32(x141, x95, 0xffffffff, &x143);
+  { uint32_t x146; uint8_t x147 = _subborrow_u32(x144, x98, 0xffffffff, &x146);
+  { uint32_t x149; uint8_t x150 = _subborrow_u32(x147, x101, 0xffffffff, &x149);
+  { uint32_t x152; uint8_t x153 = _subborrow_u32(x150, x104, 0xffffffff, &x152);
+  { uint32_t x155; uint8_t x156 = _subborrow_u32(x153, x107, 0xffffffff, &x155);
+  { uint32_t x158; uint8_t x159 = _subborrow_u32(x156, x110, 0xffffffff, &x158);
+  { uint32_t _; uint8_t x162 = _subborrow_u32(x159, x111, 0x0, &_);
+  { uint32_t x163 = cmovznz32(x162, x158, x110);
+  { uint32_t x164 = cmovznz32(x162, x155, x107);
+  { uint32_t x165 = cmovznz32(x162, x152, x104);
+  { uint32_t x166 = cmovznz32(x162, x149, x101);
+  { uint32_t x167 = cmovznz32(x162, x146, x98);
+  { uint32_t x168 = cmovznz32(x162, x143, x95);
+  { uint32_t x169 = cmovznz32(x162, x140, x92);
+  { uint32_t x170 = cmovznz32(x162, x137, x89);
+  { uint32_t x171 = cmovznz32(x162, x134, x86);
+  { uint32_t x172 = cmovznz32(x162, x131, x83);
+  { uint32_t x173 = cmovznz32(x162, x128, x80);
+  { uint32_t x174 = cmovznz32(x162, x125, x77);
+  { uint32_t x175 = cmovznz32(x162, x122, x74);
+  { uint32_t x176 = cmovznz32(x162, x119, x71);
+  { uint32_t x177 = cmovznz32(x162, x116, x68);
+  { uint32_t x178 = cmovznz32(x162, x113, x65);
+  out[0] = x178;
+  out[1] = x177;
+  out[2] = x176;
+  out[3] = x175;
+  out[4] = x174;
+  out[5] = x173;
+  out[6] = x172;
+  out[7] = x171;
+  out[8] = x170;
+  out[9] = x169;
+  out[10] = x168;
+  out[11] = x167;
+  out[12] = x166;
+  out[13] = x165;
+  out[14] = x164;
+  out[15] = x163;
+  }}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}
+}
diff --git a/src/Specific/montgomery32_2e512m569_16limbs/feaddDisplay.log b/src/Specific/montgomery32_2e512m569_16limbs/feaddDisplay.log
new file mode 100644
index 000000000..7b83b4a4f
--- /dev/null
+++ b/src/Specific/montgomery32_2e512m569_16limbs/feaddDisplay.log
@@ -0,0 +1,56 @@
+λ x x0 : word32 * word32 * word32 * word32 * word32 * word32 * word32 * word32 * word32 * word32 * word32 * word32 * word32 * word32 * word32 * word32,
+Interp-η
+(λ var : Syntax.base_type → Type,
+ λ '(x32, x33, x31, x29, x27, x25, x23, x21, x19, x17, x15, x13, x11, x9, x7, x5, (x62, x63, x61, x59, x57, x55, x53, x51, x49, x47, x45, x43, x41, x39, x37, x35))%core,
+ uint32_t x65, uint8_t x66 = addcarryx_u32(0x0, x5, x35);
+ uint32_t x68, uint8_t x69 = addcarryx_u32(x66, x7, x37);
+ uint32_t x71, uint8_t x72 = addcarryx_u32(x69, x9, x39);
+ uint32_t x74, uint8_t x75 = addcarryx_u32(x72, x11, x41);
+ uint32_t x77, uint8_t x78 = addcarryx_u32(x75, x13, x43);
+ uint32_t x80, uint8_t x81 = addcarryx_u32(x78, x15, x45);
+ uint32_t x83, uint8_t x84 = addcarryx_u32(x81, x17, x47);
+ uint32_t x86, uint8_t x87 = addcarryx_u32(x84, x19, x49);
+ uint32_t x89, uint8_t x90 = addcarryx_u32(x87, x21, x51);
+ uint32_t x92, uint8_t x93 = addcarryx_u32(x90, x23, x53);
+ uint32_t x95, uint8_t x96 = addcarryx_u32(x93, x25, x55);
+ uint32_t x98, uint8_t x99 = addcarryx_u32(x96, x27, x57);
+ uint32_t x101, uint8_t x102 = addcarryx_u32(x99, x29, x59);
+ uint32_t x104, uint8_t x105 = addcarryx_u32(x102, x31, x61);
+ uint32_t x107, uint8_t x108 = addcarryx_u32(x105, x33, x63);
+ uint32_t x110, uint8_t x111 = addcarryx_u32(x108, x32, x62);
+ uint32_t x113, uint8_t x114 = subborrow_u32(0x0, x65, 0xfffffdc7);
+ uint32_t x116, uint8_t x117 = subborrow_u32(x114, x68, 0xffffffff);
+ uint32_t x119, uint8_t x120 = subborrow_u32(x117, x71, 0xffffffff);
+ uint32_t x122, uint8_t x123 = subborrow_u32(x120, x74, 0xffffffff);
+ uint32_t x125, uint8_t x126 = subborrow_u32(x123, x77, 0xffffffff);
+ uint32_t x128, uint8_t x129 = subborrow_u32(x126, x80, 0xffffffff);
+ uint32_t x131, uint8_t x132 = subborrow_u32(x129, x83, 0xffffffff);
+ uint32_t x134, uint8_t x135 = subborrow_u32(x132, x86, 0xffffffff);
+ uint32_t x137, uint8_t x138 = subborrow_u32(x135, x89, 0xffffffff);
+ uint32_t x140, uint8_t x141 = subborrow_u32(x138, x92, 0xffffffff);
+ uint32_t x143, uint8_t x144 = subborrow_u32(x141, x95, 0xffffffff);
+ uint32_t x146, uint8_t x147 = subborrow_u32(x144, x98, 0xffffffff);
+ uint32_t x149, uint8_t x150 = subborrow_u32(x147, x101, 0xffffffff);
+ uint32_t x152, uint8_t x153 = subborrow_u32(x150, x104, 0xffffffff);
+ uint32_t x155, uint8_t x156 = subborrow_u32(x153, x107, 0xffffffff);
+ uint32_t x158, uint8_t x159 = subborrow_u32(x156, x110, 0xffffffff);
+ uint32_t _, uint8_t x162 = subborrow_u32(x159, x111, 0x0);
+ uint32_t x163 = cmovznz32(x162, x158, x110);
+ uint32_t x164 = cmovznz32(x162, x155, x107);
+ uint32_t x165 = cmovznz32(x162, x152, x104);
+ uint32_t x166 = cmovznz32(x162, x149, x101);
+ uint32_t x167 = cmovznz32(x162, x146, x98);
+ uint32_t x168 = cmovznz32(x162, x143, x95);
+ uint32_t x169 = cmovznz32(x162, x140, x92);
+ uint32_t x170 = cmovznz32(x162, x137, x89);
+ uint32_t x171 = cmovznz32(x162, x134, x86);
+ uint32_t x172 = cmovznz32(x162, x131, x83);
+ uint32_t x173 = cmovznz32(x162, x128, x80);
+ uint32_t x174 = cmovznz32(x162, x125, x77);
+ uint32_t x175 = cmovznz32(x162, x122, x74);
+ uint32_t x176 = cmovznz32(x162, x119, x71);
+ uint32_t x177 = cmovznz32(x162, x116, x68);
+ uint32_t x178 = cmovznz32(x162, x113, x65);
+ return (x163, x164, x165, x166, x167, x168, x169, x170, x171, x172, x173, x174, x175, x176, x177, x178))
+(x, x0)%core
+     : word32 * word32 * word32 * word32 * word32 * word32 * word32 * word32 * word32 * word32 * word32 * word32 * word32 * word32 * word32 * word32 → word32 * word32 * word32 * word32 * word32 * word32 * word32 * word32 * word32 * word32 * word32 * word32 * word32 * word32 * word32 * word32 → ReturnType (uint32_t * uint32_t * uint32_t * uint32_t * uint32_t * uint32_t * uint32_t * uint32_t * uint32_t * uint32_t * uint32_t * uint32_t * uint32_t * uint32_t * uint32_t * uint32_t)