forked from cloudflare/circl
/
fp_generic.go
97 lines (73 loc) · 2.41 KB
/
fp_generic.go
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
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
//go:build go1.12
// +build go1.12
package fourq
import (
"encoding/binary"
"math/bits"
)
func fpModGeneric(c *Fp) { fpSubGeneric(c, c, &modulusP) }
func fpAddGeneric(c, a, b *Fp) {
a0 := binary.LittleEndian.Uint64(a[0*8 : 1*8])
a1 := binary.LittleEndian.Uint64(a[1*8 : 2*8])
b0 := binary.LittleEndian.Uint64(b[0*8 : 1*8])
b1 := binary.LittleEndian.Uint64(b[1*8 : 2*8])
c0, x := bits.Add64(a0, b0, 0)
c1, _ := bits.Add64(a1, b1, x)
c1, x = bits.Add64(c1, c1, 0)
c0, x = bits.Add64(c0, 0, x)
c1, _ = bits.Add64(c1>>1, 0, x)
binary.LittleEndian.PutUint64(c[0*8:1*8], c0)
binary.LittleEndian.PutUint64(c[1*8:2*8], c1)
}
func fpSubGeneric(c, a, b *Fp) {
a0 := binary.LittleEndian.Uint64(a[0*8 : 1*8])
a1 := binary.LittleEndian.Uint64(a[1*8 : 2*8])
b0 := binary.LittleEndian.Uint64(b[0*8 : 1*8])
b1 := binary.LittleEndian.Uint64(b[1*8 : 2*8])
c0, x := bits.Sub64(a0, b0, 0)
c1, _ := bits.Sub64(a1, b1, x)
c1, x = bits.Add64(c1, c1, 0)
c0, x = bits.Sub64(c0, 0, x)
c1, _ = bits.Sub64(c1>>1, 0, x)
binary.LittleEndian.PutUint64(c[0*8:1*8], c0)
binary.LittleEndian.PutUint64(c[1*8:2*8], c1)
}
func fpMulGeneric(c, a, b *Fp) {
a0 := binary.LittleEndian.Uint64(a[0*8 : 1*8])
a1 := binary.LittleEndian.Uint64(a[1*8 : 2*8])
b0 := binary.LittleEndian.Uint64(b[0*8 : 1*8])
b1 := binary.LittleEndian.Uint64(b[1*8 : 2*8])
c1, c0 := bits.Mul64(a0, b0)
hi, lo := bits.Mul64(a0, b1)
c0, x := bits.Add64(c0, hi<<1, 0)
c1, x = bits.Add64(c1, lo, x)
c2, _ := bits.Add64(0, 0, x)
hi, lo = bits.Mul64(a1, b0)
c0, x = bits.Add64(c0, hi<<1, 0)
c1, x = bits.Add64(c1, lo, x)
c2, _ = bits.Add64(c2, 0, x)
hi, lo = bits.Mul64(a1, b1)
lo, x = bits.Add64(lo, lo, 0)
hi, _ = bits.Add64(hi, hi, x)
c0, x = bits.Add64(c0, lo, 0)
c1, x = bits.Add64(c1, hi, x)
c2, _ = bits.Add64(c2, 0, x)
c1, x = bits.Add64(c1, c1, 0)
c0, x = bits.Add64(c0, c2<<1, x)
c1, _ = bits.Add64(c1>>1, 0, x)
c1, x = bits.Add64(c1, c1, 0)
c0, x = bits.Add64(c0, 0, x)
c1, _ = bits.Add64(c1>>1, 0, x)
binary.LittleEndian.PutUint64(c[0*8:1*8], c0)
binary.LittleEndian.PutUint64(c[1*8:2*8], c1)
}
func fpSqrGeneric(c, a *Fp) { fpMulGeneric(c, a, a) }
func fpHlfGeneric(c, a *Fp) {
a0 := binary.LittleEndian.Uint64(a[0*8 : 1*8])
a1 := binary.LittleEndian.Uint64(a[1*8 : 2*8])
hlf := a0 & 0x1
c0 := (a1 << 63) | (a0 >> 1)
c1 := (hlf << 62) | (a1 >> 1)
binary.LittleEndian.PutUint64(c[0*8:1*8], c0)
binary.LittleEndian.PutUint64(c[1*8:2*8], c1)
}