1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
|
package org.spongycastle.math.raw;
public abstract class Nat512
{
public static void mul(int[] x, int[] y, int[] zz)
{
Nat256.mul(x, y, zz);
Nat256.mul(x, 8, y, 8, zz, 16);
int c24 = Nat256.addToEachOther(zz, 8, zz, 16);
int c16 = c24 + Nat256.addTo(zz, 0, zz, 8, 0);
c24 += Nat256.addTo(zz, 24, zz, 16, c16);
int[] dx = Nat256.create(), dy = Nat256.create();
boolean neg = Nat256.diff(x, 8, x, 0, dx, 0) != Nat256.diff(y, 8, y, 0, dy, 0);
int[] tt = Nat256.createExt();
Nat256.mul(dx, dy, tt);
c24 += neg ? Nat.addTo(16, tt, 0, zz, 8) : Nat.subFrom(16, tt, 0, zz, 8);
Nat.addWordAt(32, c24, zz, 24);
}
public static void square(int[] x, int[] zz)
{
Nat256.square(x, zz);
Nat256.square(x, 8, zz, 16);
int c24 = Nat256.addToEachOther(zz, 8, zz, 16);
int c16 = c24 + Nat256.addTo(zz, 0, zz, 8, 0);
c24 += Nat256.addTo(zz, 24, zz, 16, c16);
int[] dx = Nat256.create();
Nat256.diff(x, 8, x, 0, dx, 0);
int[] tt = Nat256.createExt();
Nat256.square(dx, tt);
c24 += Nat.subFrom(16, tt, 0, zz, 8);
Nat.addWordAt(32, c24, zz, 24);
}
}
|