]> Cypherpunks.ru repositories - gogost.git/blob - mgm/mul128.go
Optimized GF^128
[gogost.git] / mgm / mul128.go
1 // GoGOST -- Pure Go GOST cryptographic functions library
2 // Copyright (C) 2015-2021 Sergey Matveev <stargrave@stargrave.org>
3 //
4 // This program is free software: you can redistribute it and/or modify
5 // it under the terms of the GNU General Public License as published by
6 // the Free Software Foundation, version 3 of the License.
7 //
8 // This program is distributed in the hope that it will be useful,
9 // but WITHOUT ANY WARRANTY; without even the implied warranty of
10 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
11 // GNU General Public License for more details.
12 //
13 // You should have received a copy of the GNU General Public License
14 // along with this program.  If not, see <http://www.gnu.org/licenses/>.
15
16 package mgm
17
18 import "encoding/binary"
19
20 type mul128 struct{ buf [16]byte }
21
22 func newMul128() *mul128 {
23         return &mul128{}
24 }
25
26 func gf128half(n int, t, x0, x1, z0, z1 uint64) (uint64, uint64, uint64, uint64, uint64) {
27         var sign bool
28         for i := 0; i < n; i++ {
29                 if t&1 > 0 {
30                         z0, z1 = z0^x0, z1^x1
31                 }
32                 t >>= 1
33                 sign = x1>>63 > 0
34                 x1 = (x1 << 1) ^ (x0 >> 63)
35                 x0 <<= 1
36                 if sign {
37                         x0 ^= 0x87
38                 }
39         }
40         return t, x0, x1, z0, z1
41 }
42
43 func (mul *mul128) Mul(x, y []byte) []byte {
44         x1 := binary.BigEndian.Uint64(x[:8])
45         x0 := binary.BigEndian.Uint64(x[8:])
46         y1 := binary.BigEndian.Uint64(y[:8])
47         y0 := binary.BigEndian.Uint64(y[8:])
48         t, x0, x1, z0, z1 := gf128half(64, y0, x0, x1, 0, 0)
49         t, x0, x1, z0, z1 = gf128half(63, y1, x0, x1, z0, z1)
50         if t&1 > 0 {
51                 z0, z1 = z0^x0, z1^x1
52         }
53         binary.BigEndian.PutUint64(mul.buf[:8], z1)
54         binary.BigEndian.PutUint64(mul.buf[8:], z0)
55         return mul.buf[:]
56 }