internal static void Square(ulong[] x, ulong[] z)
{
ulong[] t = new ulong[4];
Interleave.Expand64To128Rev(x[0], t, 0);
Interleave.Expand64To128Rev(x[1], t, 2);
ulong z0 = t[0], z1 = t[1], z2 = t[2], z3 = t[3];
z1 ^= z3 ^ (z3 >> 1) ^ (z3 >> 2) ^ (z3 >> 7);
// z2 ^= (z3 << 63) ^ (z3 << 62) ^ (z3 << 57);
z2 ^= (z3 << 62) ^ (z3 << 57);
z0 ^= z2 ^ (z2 >> 1) ^ (z2 >> 2) ^ (z2 >> 7);
z1 ^= (z2 << 63) ^ (z2 << 62) ^ (z2 << 57);
z[0] = z0;
z[1] = z1;
}