diff options
Diffstat (limited to 'core/src/main/java/org/bouncycastle/math/ec/custom/sec/Nat192.java')
| -rw-r--r-- | core/src/main/java/org/bouncycastle/math/ec/custom/sec/Nat192.java | 102 |
1 files changed, 72 insertions, 30 deletions
diff --git a/core/src/main/java/org/bouncycastle/math/ec/custom/sec/Nat192.java b/core/src/main/java/org/bouncycastle/math/ec/custom/sec/Nat192.java index f080cc27..c622350a 100644 --- a/core/src/main/java/org/bouncycastle/math/ec/custom/sec/Nat192.java +++ b/core/src/main/java/org/bouncycastle/math/ec/custom/sec/Nat192.java @@ -596,8 +596,9 @@ public abstract class Nat192 long x_0 = x[0] & M; long zz_1; + int c = 0, w; { - int c = 0, i = 5, j = 12; + int i = 5, j = 12; do { long xVal = (x[i--] & M); @@ -611,7 +612,8 @@ public abstract class Nat192 { long p = x_0 * x_0; zz_1 = ((c << 31) & M) | (p >>> 33); - zz[0] = (int)(p >>> 1); + zz[0] = (int)p; + c = (int)(p >>> 32) & 1; } } @@ -620,7 +622,9 @@ public abstract class Nat192 { zz_1 += x_1 * x_0; - zz[1] = (int)zz_1; + w = (int)zz_1; + zz[1] = (w << 1) | c; + c = w >>> 31; zz_2 += zz_1 >>> 32; } @@ -629,7 +633,9 @@ public abstract class Nat192 long zz_4 = zz[4] & M; { zz_2 += x_2 * x_0; - zz[2] = (int)zz_2; + w = (int)zz_2; + zz[2] = (w << 1) | c; + c = w >>> 31; zz_3 += (zz_2 >>> 32) + x_2 * x_1; zz_4 += zz_3 >>> 32; zz_3 &= M; @@ -640,7 +646,9 @@ public abstract class Nat192 long zz_6 = zz[6] & M; { zz_3 += x_3 * x_0; - zz[3] = (int)zz_3; + w = (int)zz_3; + zz[3] = (w << 1) | c; + c = w >>> 31; zz_4 += (zz_3 >>> 32) + x_3 * x_1; zz_5 += (zz_4 >>> 32) + x_3 * x_2; zz_4 &= M; @@ -653,7 +661,9 @@ public abstract class Nat192 long zz_8 = zz[8] & M; { zz_4 += x_4 * x_0; - zz[4] = (int)zz_4; + w = (int)zz_4; + zz[4] = (w << 1) | c; + c = w >>> 31; zz_5 += (zz_4 >>> 32) + x_4 * x_1; zz_6 += (zz_5 >>> 32) + x_4 * x_2; zz_5 &= M; @@ -668,7 +678,9 @@ public abstract class Nat192 long zz_10 = zz[10] & M; { zz_5 += x_5 * x_0; - zz[5] = (int)zz_5; + w = (int)zz_5; + zz[5] = (w << 1) | c; + c = w >>> 31; zz_6 += (zz_5 >>> 32) + x_5 * x_1; zz_7 += (zz_6 >>> 32) + x_5 * x_2; zz_8 += (zz_7 >>> 32) + x_5 * x_3; @@ -676,14 +688,23 @@ public abstract class Nat192 zz_10 += zz_9 >>> 32; } - zz[6] = (int)zz_6; - zz[7] = (int)zz_7; - zz[8] = (int)zz_8; - zz[9] = (int)zz_9; - zz[10] = (int)zz_10; - zz[11] += (int)(zz_10 >>> 32); - - Nat.shiftUpBit(12, zz, (int)x_0 << 31); + w = (int)zz_6; + zz[6] = (w << 1) | c; + c = w >>> 31; + w = (int)zz_7; + zz[7] = (w << 1) | c; + c = w >>> 31; + w = (int)zz_8; + zz[8] = (w << 1) | c; + c = w >>> 31; + w = (int)zz_9; + zz[9] = (w << 1) | c; + c = w >>> 31; + w = (int)zz_10; + zz[10] = (w << 1) | c; + c = w >>> 31; + w = zz[11] + (int)(zz_10 >> 32); + zz[11] = (w << 1) | c; } public static void square(int[] x, int xOff, int[] zz, int zzOff) @@ -691,8 +712,9 @@ public abstract class Nat192 long x_0 = x[xOff + 0] & M; long zz_1; + int c = 0, w; { - int c = 0, i = 5, j = 12; + int i = 5, j = 12; do { long xVal = (x[xOff + i--] & M); @@ -706,7 +728,8 @@ public abstract class Nat192 { long p = x_0 * x_0; zz_1 = ((c << 31) & M) | (p >>> 33); - zz[zzOff + 0] = (int)(p >>> 1); + zz[zzOff + 0] = (int)p; + c = (int)(p >>> 32) & 1; } } @@ -715,7 +738,9 @@ public abstract class Nat192 { zz_1 += x_1 * x_0; - zz[zzOff + 1] = (int)zz_1; + w = (int)zz_1; + zz[zzOff + 1] = (w << 1) | c; + c = w >>> 31; zz_2 += zz_1 >>> 32; } @@ -724,7 +749,9 @@ public abstract class Nat192 long zz_4 = zz[zzOff + 4] & M; { zz_2 += x_2 * x_0; - zz[zzOff + 2] = (int)zz_2; + w = (int)zz_2; + zz[zzOff + 2] = (w << 1) | c; + c = w >>> 31; zz_3 += (zz_2 >>> 32) + x_2 * x_1; zz_4 += zz_3 >>> 32; zz_3 &= M; @@ -735,7 +762,9 @@ public abstract class Nat192 long zz_6 = zz[zzOff + 6] & M; { zz_3 += x_3 * x_0; - zz[zzOff + 3] = (int)zz_3; + w = (int)zz_3; + zz[zzOff + 3] = (w << 1) | c; + c = w >>> 31; zz_4 += (zz_3 >>> 32) + x_3 * x_1; zz_5 += (zz_4 >>> 32) + x_3 * x_2; zz_4 &= M; @@ -748,7 +777,9 @@ public abstract class Nat192 long zz_8 = zz[zzOff + 8] & M; { zz_4 += x_4 * x_0; - zz[zzOff + 4] = (int)zz_4; + w = (int)zz_4; + zz[zzOff + 4] = (w << 1) | c; + c = w >>> 31; zz_5 += (zz_4 >>> 32) + x_4 * x_1; zz_6 += (zz_5 >>> 32) + x_4 * x_2; zz_5 &= M; @@ -763,7 +794,9 @@ public abstract class Nat192 long zz_10 = zz[zzOff + 10] & M; { zz_5 += x_5 * x_0; - zz[zzOff + 5] = (int)zz_5; + w = (int)zz_5; + zz[zzOff + 5] = (w << 1) | c; + c = w >>> 31; zz_6 += (zz_5 >>> 32) + x_5 * x_1; zz_7 += (zz_6 >>> 32) + x_5 * x_2; zz_8 += (zz_7 >>> 32) + x_5 * x_3; @@ -771,14 +804,23 @@ public abstract class Nat192 zz_10 += zz_9 >>> 32; } - zz[zzOff + 6] = (int)zz_6; - zz[zzOff + 7] = (int)zz_7; - zz[zzOff + 8] = (int)zz_8; - zz[zzOff + 9] = (int)zz_9; - zz[zzOff + 10] = (int)zz_10; - zz[zzOff + 11] += (int)(zz_10 >>> 32); - - Nat.shiftUpBit(12, zz, zzOff, (int)x_0 << 31); + w = (int)zz_6; + zz[zzOff + 6] = (w << 1) | c; + c = w >>> 31; + w = (int)zz_7; + zz[zzOff + 7] = (w << 1) | c; + c = w >>> 31; + w = (int)zz_8; + zz[zzOff + 8] = (w << 1) | c; + c = w >>> 31; + w = (int)zz_9; + zz[zzOff + 9] = (w << 1) | c; + c = w >>> 31; + w = (int)zz_10; + zz[zzOff + 10] = (w << 1) | c; + c = w >>> 31; + w = zz[zzOff + 11] + (int)(zz_10 >> 32); + zz[zzOff + 11] = (w << 1) | c; } public static int sub(int[] x, int[] y, int[] z) |