private static void RP(int[] nc, int[] d1, int[] e1, int[] f1, int[] t) { Nat224.copy(nc, f1); int[] d0 = Nat224.create(); int[] e0 = Nat224.create(); for (int i = 0; i < 7; ++i) { Nat224.copy(d1, d0); Nat224.copy(e1, e0); int j = 1 << i; while (--j >= 0) { RS(d1, e1, f1, t); } RM(nc, d0, e0, d1, e1, f1, t); } }
private static void RP(int[] nc, int[] d1, int[] e1, int[] f1, int[] t) { Nat224.copy(nc, f1); int[] d0 = Nat224.create(); int[] e0 = Nat224.create(); for (int i = 0; i < 7; ++i) { Nat224.copy(d1, d0); Nat224.copy(e1, e0); int j = 1 << i; while (--j >= 0) { RS(d1, e1, f1, t); } RM(nc, d0, e0, d1, e1, f1, t); } }
private static boolean isSquare(int[] x) { int[] t1 = Nat224.create(); int[] t2 = Nat224.create(); Nat224.copy(x, t1); for (int i = 0; i < 7; ++i) { Nat224.copy(t1, t2); SecP224R1Field.squareN(t1, 1 << i, t1); SecP224R1Field.multiply(t1, t2, t1); } SecP224R1Field.squareN(t1, 95, t1); return Nat224.isOne(t1); }
private static boolean isSquare(int[] x) { int[] t1 = Nat224.create(); int[] t2 = Nat224.create(); Nat224.copy(x, t1); for (int i = 0; i < 7; ++i) { Nat224.copy(t1, t2); SecP224R1Field.squareN(t1, 1 << i, t1); SecP224R1Field.multiply(t1, t2, t1); } SecP224R1Field.squareN(t1, 95, t1); return Nat224.isOne(t1); }
private static boolean trySqrt(int[] nc, int[] r, int[] t) { int[] d1 = Nat224.create(); Nat224.copy(r, d1); int[] e1 = Nat224.create(); e1[0] = 1; int[] f1 = Nat224.create(); RP(nc, d1, e1, f1, t); int[] d0 = Nat224.create(); int[] e0 = Nat224.create(); for (int k = 1; k < 96; ++k) { Nat224.copy(d1, d0); Nat224.copy(e1, e0); RS(d1, e1, f1, t); if (Nat224.isZero(d1)) { Mod.invert(SecP224R1Field.P, e0, t); SecP224R1Field.multiply(t, d0, t); return true; } } return false; } }
private static boolean trySqrt(int[] nc, int[] r, int[] t) { int[] d1 = Nat224.create(); Nat224.copy(r, d1); int[] e1 = Nat224.create(); e1[0] = 1; int[] f1 = Nat224.create(); RP(nc, d1, e1, f1, t); int[] d0 = Nat224.create(); int[] e0 = Nat224.create(); for (int k = 1; k < 96; ++k) { Nat224.copy(d1, d0); Nat224.copy(e1, e0); RS(d1, e1, f1, t); if (Nat224.isZero(d1)) { Mod.invert(SecP224R1Field.P, e0, t); SecP224R1Field.multiply(t, d0, t); return true; } } return false; } }
private static void RM(int[] nc, int[] d0, int[] e0, int[] d1, int[] e1, int[] f1, int[] t) { SecP224R1Field.multiply(e1, e0, t); SecP224R1Field.multiply(t, nc, t); SecP224R1Field.multiply(d1, d0, f1); SecP224R1Field.add(f1, t, f1); SecP224R1Field.multiply(d1, e0, t); Nat224.copy(f1, d1); SecP224R1Field.multiply(e1, d0, e1); SecP224R1Field.add(e1, t, e1); SecP224R1Field.square(e1, f1); SecP224R1Field.multiply(f1, nc, f1); }
private static void RM(int[] nc, int[] d0, int[] e0, int[] d1, int[] e1, int[] f1, int[] t) { SecP224R1Field.multiply(e1, e0, t); SecP224R1Field.multiply(t, nc, t); SecP224R1Field.multiply(d1, d0, f1); SecP224R1Field.add(f1, t, f1); SecP224R1Field.multiply(d1, e0, t); Nat224.copy(f1, d1); SecP224R1Field.multiply(e1, d0, e1); SecP224R1Field.add(e1, t, e1); SecP224R1Field.square(e1, f1); SecP224R1Field.multiply(f1, nc, f1); }