public static void sqrt(long[] x, long[] z) { long u0, u1; u0 = Interleave.unshuffle(x[0]); u1 = Interleave.unshuffle(x[1]); long e0 = (u0 & 0x00000000FFFFFFFFL) | (u1 << 32); long c0 = (u0 >>> 32) | (u1 & 0xFFFFFFFF00000000L); u0 = Interleave.unshuffle(x[2]); long e1 = (u0 & 0x00000000FFFFFFFFL) ^ (x[3] << 32); long c1 = (u0 >>> 32); z[0] = e0 ^ (c0 << 8); z[1] = e1 ^ (c1 << 8) ^ (c0 >>> 56) ^ (c0 << 33); z[2] = (c1 >>> 56) ^ (c1 << 33) ^ (c0 >>> 31); z[3] = (c1 >>> 31); }
public static void sqrt(long[] x, long[] z) { long u0, u1; u0 = Interleave.unshuffle(x[0]); u1 = Interleave.unshuffle(x[1]); long e0 = (u0 & 0x00000000FFFFFFFFL) | (u1 << 32); long c0 = (u0 >>> 32) | (u1 & 0xFFFFFFFF00000000L); u0 = Interleave.unshuffle(x[2]); u1 = Interleave.unshuffle(x[3]); long e1 = (u0 & 0x00000000FFFFFFFFL) | (u1 << 32); long c1 = (u0 >>> 32) | (u1 & 0xFFFFFFFF00000000L); u0 = Interleave.unshuffle(x[4]); u1 = Interleave.unshuffle(x[5]); long e2 = (u0 & 0x00000000FFFFFFFFL) | (u1 << 32); long c2 = (u0 >>> 32) | (u1 & 0xFFFFFFFF00000000L); u0 = Interleave.unshuffle(x[6]); long e3 = (u0 & 0x00000000FFFFFFFFL); long c3 = (u0 >>> 32); z[0] = e0 ^ (c0 << 44); z[1] = e1 ^ (c1 << 44) ^ (c0 >>> 20); z[2] = e2 ^ (c2 << 44) ^ (c1 >>> 20); z[3] = e3 ^ (c3 << 44) ^ (c2 >>> 20) ^ (c0 << 13); z[4] = (c3 >>> 20) ^ (c1 << 13) ^ (c0 >>> 51); z[5] = (c2 << 13) ^ (c1 >>> 51); z[6] = (c3 << 13) ^ (c2 >>> 51); // assert (c3 >>> 51) == 0; }
public static void sqrt(long[] x, long[] z) { long[] odd = Nat320.create64(); long u0, u1; u0 = Interleave.unshuffle(x[0]); u1 = Interleave.unshuffle(x[1]); long e0 = (u0 & 0x00000000FFFFFFFFL) | (u1 << 32); odd[0] = (u0 >>> 32) | (u1 & 0xFFFFFFFF00000000L); u0 = Interleave.unshuffle(x[2]); u1 = Interleave.unshuffle(x[3]); long e1 = (u0 & 0x00000000FFFFFFFFL) | (u1 << 32); odd[1] = (u0 >>> 32) | (u1 & 0xFFFFFFFF00000000L); u0 = Interleave.unshuffle(x[4]); long e2 = (u0 & 0x00000000FFFFFFFFL); odd[2] = (u0 >>> 32); multiply(odd, ROOT_Z, z); z[0] ^= e0; z[1] ^= e1; z[2] ^= e2; }
public static void sqrt(long[] x, long[] z) { long[] odd = Nat192.create64(); long u0, u1; u0 = Interleave.unshuffle(x[0]); u1 = Interleave.unshuffle(x[1]); long e0 = (u0 & 0x00000000FFFFFFFFL) | (u1 << 32); odd[0] = (u0 >>> 32) | (u1 & 0xFFFFFFFF00000000L); u0 = Interleave.unshuffle(x[2]); long e1 = (u0 & 0x00000000FFFFFFFFL); odd[1] = (u0 >>> 32); multiply(odd, ROOT_Z, z); z[0] ^= e0; z[1] ^= e1; }
public static void sqrt(long[] x, long[] z) { long[] odd = Nat192.create64(); long u0, u1; u0 = Interleave.unshuffle(x[0]); u1 = Interleave.unshuffle(x[1]); long e0 = (u0 & 0x00000000FFFFFFFFL) | (u1 << 32); odd[0] = (u0 >>> 32) | (u1 & 0xFFFFFFFF00000000L); u0 = Interleave.unshuffle(x[2]); long e1 = (u0 & 0x00000000FFFFFFFFL); odd[1] = (u0 >>> 32); multiply(odd, ROOT_Z, z); z[0] ^= e0; z[1] ^= e1; }
public static void sqrt(long[] x, long[] z) u0 = Interleave.unshuffle(x[0]); u1 = Interleave.unshuffle(x[1]); long e0 = (u0 & 0x00000000FFFFFFFFL) | (u1 << 32); long c0 = (u0 >>> 32) | (u1 & 0xFFFFFFFF00000000L); u0 = Interleave.unshuffle(x[2]); u1 = Interleave.unshuffle(x[3]); long e1 = (u0 & 0x00000000FFFFFFFFL) | (u1 << 32); long c1 = (u0 >>> 32) | (u1 & 0xFFFFFFFF00000000L);
public static void sqrt(long[] x, long[] z) u0 = Interleave.unshuffle(x[0]); u1 = Interleave.unshuffle(x[1]); long e0 = (u0 & 0x00000000FFFFFFFFL) | (u1 << 32); long c0 = (u0 >>> 32) | (u1 & 0xFFFFFFFF00000000L); u0 = Interleave.unshuffle(x[2]); u1 = Interleave.unshuffle(x[3]); long e1 = (u0 & 0x00000000FFFFFFFFL) | (u1 << 32); long c1 = (u0 >>> 32) | (u1 & 0xFFFFFFFF00000000L);
public static void sqrt(long[] x, long[] z) { long[] evn = Nat576.create64(), odd = Nat576.create64(); int pos = 0; for (int i = 0; i < 4; ++i) { long u0 = Interleave.unshuffle(x[pos++]); long u1 = Interleave.unshuffle(x[pos++]); evn[i] = (u0 & 0x00000000FFFFFFFFL) | (u1 << 32); odd[i] = (u0 >>> 32) | (u1 & 0xFFFFFFFF00000000L); } { long u0 = Interleave.unshuffle(x[pos]); evn[4] = (u0 & 0x00000000FFFFFFFFL); odd[4] = (u0 >>> 32); } multiply(odd, ROOT_Z, z); add(z, evn, z); }