mgm/mul64.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 "math/big"
  19
  20 const Mul64MaxBit = 64 - 1
  21
  22 var R64 = big.NewInt(0).SetBytes([]byte{
  23         0x00, 0x00, 0x00, 0x00, 0x00, 0x00, 0x00, 0x1b,
  24 })
  25
  26 type mul64 struct {
  27         x   *big.Int
  28         y   *big.Int
  29         z   *big.Int
  30         buf [8]byte
  31 }
  32
  33 func newMul64() *mul64 {
  34         return &mul64{
  35                 x: big.NewInt(0),
  36                 y: big.NewInt(0),
  37                 z: big.NewInt(0),
  38         }
  39 }
  40
  41 func (mul *mul64) Mul(x, y []byte) []byte {
  42         mul.x.SetBytes(x)
  43         mul.y.SetBytes(y)
  44         mul.z.SetInt64(0)
  45         for mul.y.BitLen() != 0 {
  46                 if mul.y.Bit(0) == 1 {
  47                         mul.z.Xor(mul.z, mul.x)
  48                 }
  49                 if mul.x.Bit(Mul64MaxBit) == 1 {
  50                         mul.x.SetBit(mul.x, Mul64MaxBit, 0)
  51                         mul.x.Lsh(mul.x, 1)
  52                         mul.x.Xor(mul.x, R64)
  53                 } else {
  54                         mul.x.Lsh(mul.x, 1)
  55                 }
  56                 mul.y.Rsh(mul.y, 1)
  57         }
  58         zBytes := mul.z.Bytes()
  59         rem := len(x) - len(zBytes)
  60         for i := 0; i < rem; i++ {
  61                 mul.buf[i] = 0
  62         }
  63         copy(mul.buf[rem:], zBytes)
  64         return mul.buf[:]
  65 }