summaryrefslogtreecommitdiff
path: root/src/add_python
diff options
context:
space:
mode:
Diffstat (limited to 'src/add_python')
-rw-r--r--src/add_python/lilliput/__init__.py31
-rw-r--r--src/add_python/lilliput/constants.py40
-rw-r--r--src/add_python/lilliput/helpers.py92
-rw-r--r--src/add_python/lilliput/lilliput_ae_1.py155
-rw-r--r--src/add_python/lilliput/lilliput_ae_2.py127
-rw-r--r--src/add_python/lilliput/multiplications.py186
-rw-r--r--src/add_python/lilliput/tbc.py178
7 files changed, 809 insertions, 0 deletions
diff --git a/src/add_python/lilliput/__init__.py b/src/add_python/lilliput/__init__.py
new file mode 100644
index 0000000..5fbc0de
--- /dev/null
+++ b/src/add_python/lilliput/__init__.py
@@ -0,0 +1,31 @@
+from . import lilliput_ae_1
+from . import lilliput_ae_2
+from .constants import NONCE_BYTES
+
+
+_AE_MODES = {
+ 1: lilliput_ae_1,
+ 2: lilliput_ae_2
+}
+
+
+def _check_inputs(key, mode, nonce):
+ valid_key_lengths = (128, 192, 256)
+ if len(key)*8 not in valid_key_lengths:
+ raise ValueError('invalid key size: {} not in {}'.format(len(key)*8, valid_key_lengths))
+
+ if mode not in _AE_MODES:
+ raise ValueError('invalid mode: {} not in {}'.format(mode, tuple(_AE_MODES)))
+
+ if len(nonce) != NONCE_BYTES:
+ raise ValueError('invalid nonce size: expecting {}, have {}'.format(NONCE_BYTES, len(nonce)))
+
+
+def encrypt(plaintext, adata, key, nonce, mode):
+ _check_inputs(key, mode, nonce)
+ return _AE_MODES[mode].encrypt(adata, plaintext, nonce, key)
+
+
+def decrypt(ciphertext, tag, adata, key, nonce, mode):
+ _check_inputs(key, mode, nonce)
+ return _AE_MODES[mode].decrypt(adata, ciphertext, nonce, tag, key)
diff --git a/src/add_python/lilliput/constants.py b/src/add_python/lilliput/constants.py
new file mode 100644
index 0000000..0c9b89f
--- /dev/null
+++ b/src/add_python/lilliput/constants.py
@@ -0,0 +1,40 @@
+BLOCK_BITS = 128
+BLOCK_BYTES = BLOCK_BITS//8
+NONCE_BYTES = 15
+TAG_BYTES = 16
+
+
+Sbox = [
+ 0x20, 0x00, 0xb2, 0x85, 0x3b, 0x35, 0xa6, 0xa4,
+ 0x30, 0xe4, 0x6a, 0x2c, 0xff, 0x59, 0xe2, 0x0e,
+ 0xf8, 0x1e, 0x7a, 0x80, 0x15, 0xbd, 0x3e, 0xb1,
+ 0xe8, 0xf3, 0xa2, 0xc2, 0xda, 0x51, 0x2a, 0x10,
+ 0x21, 0x01, 0x23, 0x78, 0x5c, 0x24, 0x27, 0xb5,
+ 0x37, 0xc7, 0x2b, 0x1f, 0xae, 0x0a, 0x77, 0x5f,
+ 0x6f, 0x09, 0x9d, 0x81, 0x04, 0x5a, 0x29, 0xdc,
+ 0x39, 0x9c, 0x05, 0x57, 0x97, 0x74, 0x79, 0x17,
+ 0x44, 0xc6, 0xe6, 0xe9, 0xdd, 0x41, 0xf2, 0x8a,
+ 0x54, 0xca, 0x6e, 0x4a, 0xe1, 0xad, 0xb6, 0x88,
+ 0x1c, 0x98, 0x7e, 0xce, 0x63, 0x49, 0x3a, 0x5d,
+ 0x0c, 0xef, 0xf6, 0x34, 0x56, 0x25, 0x2e, 0xd6,
+ 0x67, 0x75, 0x55, 0x76, 0xb8, 0xd2, 0x61, 0xd9,
+ 0x71, 0x8b, 0xcd, 0x0b, 0x72, 0x6c, 0x31, 0x4b,
+ 0x69, 0xfd, 0x7b, 0x6d, 0x60, 0x3c, 0x2f, 0x62,
+ 0x3f, 0x22, 0x73, 0x13, 0xc9, 0x82, 0x7f, 0x53,
+ 0x32, 0x12, 0xa0, 0x7c, 0x02, 0x87, 0x84, 0x86,
+ 0x93, 0x4e, 0x68, 0x46, 0x8d, 0xc3, 0xdb, 0xec,
+ 0x9b, 0xb7, 0x89, 0x92, 0xa7, 0xbe, 0x3d, 0xd8,
+ 0xea, 0x50, 0x91, 0xf1, 0x33, 0x38, 0xe0, 0xa9,
+ 0xa3, 0x83, 0xa1, 0x1b, 0xcf, 0x06, 0x95, 0x07,
+ 0x9e, 0xed, 0xb9, 0xf5, 0x4c, 0xc0, 0xf4, 0x2d,
+ 0x16, 0xfa, 0xb4, 0x03, 0x26, 0xb3, 0x90, 0x4f,
+ 0xab, 0x65, 0xfc, 0xfe, 0x14, 0xf7, 0xe3, 0x94,
+ 0xee, 0xac, 0x8c, 0x1a, 0xde, 0xcb, 0x28, 0x40,
+ 0x7d, 0xc8, 0xc4, 0x48, 0x6b, 0xdf, 0xa5, 0x52,
+ 0xe5, 0xfb, 0xd7, 0x64, 0xf9, 0xf0, 0xd3, 0x5e,
+ 0x66, 0x96, 0x8f, 0x1d, 0x45, 0x36, 0xcc, 0xc5,
+ 0x4d, 0x9f, 0xbf, 0x0f, 0xd1, 0x08, 0xeb, 0x43,
+ 0x42, 0x19, 0xe7, 0x99, 0xa8, 0x8e, 0x58, 0xc1,
+ 0x9a, 0xd4, 0x18, 0x47, 0xaa, 0xaf, 0xbc, 0x5b,
+ 0xd5, 0x11, 0xd0, 0xb0, 0x70, 0xbb, 0x0d, 0xba
+]
diff --git a/src/add_python/lilliput/helpers.py b/src/add_python/lilliput/helpers.py
new file mode 100644
index 0000000..8677f06
--- /dev/null
+++ b/src/add_python/lilliput/helpers.py
@@ -0,0 +1,92 @@
+from .constants import BLOCK_BITS, BLOCK_BYTES
+from . import tbc
+
+
+def ArrayToBlockbytesMatrix(array):
+ vector = list(array)
+
+ blocks_nb = len(vector)//BLOCK_BYTES
+
+ block_starts = (
+ i*BLOCK_BYTES for i in range(blocks_nb)
+ )
+
+ matrix = [
+ vector[start:start+BLOCK_BYTES] for start in block_starts
+ ]
+
+ padding_len = len(vector)%BLOCK_BYTES
+
+ if padding_len > 0:
+ padding = vector[-padding_len:]
+ matrix.append(padding)
+
+ return matrix
+
+
+def BlockbytesMatrixToBytes(matrix):
+ return bytes(byte for block in matrix for byte in block)
+
+
+def XorState(state1, state2):
+ return [s1^s2 for (s1, s2) in zip(state1, state2)]
+
+
+def Padding10LSB(X):
+ zeroes = [0] * (BLOCK_BYTES-len(X)-1)
+ return zeroes + [0b10000000] + X
+
+
+def _tweakAssociatedData(t, i, padded):
+ t_bytes = t//8
+ tweak = [0]*(t_bytes)
+
+ mask = 0xff
+ for byte in range(t_bytes-1):
+ tweak[byte] = (i & mask) >> (byte * 8)
+ mask = mask << 8
+
+ mask = (0xf << (8 * t_bytes-1))
+ tweak[-1] = (i & mask) >> ((t_bytes-1)*8)
+ if not padded:
+ tweak[-1] |= 0x20
+ else:
+ tweak[-1] |= 0x60
+
+ return tweak
+
+
+def BuildAuth(t, A, key):
+ Auth = [0 for byte in range(0, BLOCK_BYTES)]
+ l_a = len(A)//BLOCK_BYTES
+ need_padding = len(A)%BLOCK_BYTES > 0
+
+ A = ArrayToBlockbytesMatrix(A)
+
+ for i in range(0, l_a):
+ tweak = _tweakAssociatedData(t, i, padded=False)
+ enc = tbc.encrypt(tweak, key, A[i])
+ Auth = XorState(Auth, enc)
+
+ if not need_padding:
+ return Auth
+
+ tweak = _tweakAssociatedData(t, l_a, padded=True)
+ ad_padded = Padding10LSB(A[l_a])
+ enc = tbc.encrypt(tweak, key, ad_padded)
+ Auth = XorState(Auth, enc)
+
+ return Auth
+
+
+class TagValidationError(Exception):
+ def __init__(self, announced, computed):
+ msg = '\n'.join((
+ 'Invalid tag:',
+ announced.hex().upper()+' (announced)',
+ computed.hex().upper()+' (computed)'
+ ))
+
+ super().__init__(msg)
+ self._announced = announced
+ self._computed = computed
diff --git a/src/add_python/lilliput/lilliput_ae_1.py b/src/add_python/lilliput/lilliput_ae_1.py
new file mode 100644
index 0000000..0da2a95
--- /dev/null
+++ b/src/add_python/lilliput/lilliput_ae_1.py
@@ -0,0 +1,155 @@
+"""
+ OCB 3 for lilliput ae i
+"""
+
+from enum import Enum
+
+from .constants import BLOCK_BYTES, NONCE_BYTES
+from .helpers import (
+ ArrayToBlockbytesMatrix,
+ BlockbytesMatrixToBytes,
+ BuildAuth,
+ Padding10LSB,
+ TagValidationError,
+ XorState
+)
+from . import tbc
+
+
+TWEAK_BITS = 192
+TWEAK_BYTES = TWEAK_BITS//8
+
+
+def LowPart(array, number_bits):
+ shifted = 0
+ for byte in range(0, len(array)):
+ shifted |= (array[byte] << (8 * byte))
+
+ mask = 0
+ for bit in range(0, number_bits):
+ mask |= (0x1 << bit)
+
+ lower_part = shifted & mask
+
+ will_pad = 0
+ if number_bits % 8 != 0:
+ will_pad = 1
+
+ lower_part_byte = []
+ nb_bytes = number_bits//8 + will_pad
+ for byte in range(nb_bytes):
+ lower_part_byte.append(lower_part & 0xff)
+ lower_part = lower_part >> 8
+
+ return lower_part_byte
+
+
+class _MessageTweak(Enum):
+ BLOCK = 0b000
+ NO_PADDING = 0b0001
+ PAD = 0b0100
+ FINAL = 0b0101
+
+
+def TweakMessage(N, j, padding):
+ tweak = [0 for byte in range(0, TWEAK_BYTES)]
+ for byte in range(NONCE_BYTES-1, -1, -1):
+ tweak[byte + (TWEAK_BYTES-NONCE_BYTES)] |= (N[byte] & 0xf0) >> 4
+ tweak[byte + (TWEAK_BYTES-NONCE_BYTES-1)] |= (N[byte] & 0x0f) << 4
+
+ tweak[TWEAK_BYTES-NONCE_BYTES-1] |= ((j >> 64) & 0xf)
+ for byte in range(TWEAK_BYTES-NONCE_BYTES-2, -1, -1):
+ tweak[byte] = (j >> (8 * byte)) & 0xff
+
+ tweak[-1] |= padding.value<<4
+
+ return tweak
+
+
+def TreatMessageEnc(M, N, key):
+ checksum = [0 for byte in range(0, BLOCK_BYTES)]
+
+ l = len(M)//BLOCK_BYTES
+ padding_bytes = len(M)%BLOCK_BYTES
+
+ M = ArrayToBlockbytesMatrix(M)
+ C = []
+
+ for j in range(0, l):
+ checksum = XorState(checksum, M[j])
+ tweak = TweakMessage(N, j, _MessageTweak.BLOCK)
+ C.append(tbc.encrypt(tweak, key, M[j]))
+
+ if padding_bytes == 0:
+ tweak = TweakMessage(N, l, _MessageTweak.NO_PADDING)
+ Final = tbc.encrypt(tweak, key, checksum)
+
+ else:
+ m_padded = Padding10LSB(M[l])
+ checksum = XorState(checksum, m_padded)
+ tweak = TweakMessage(N, l, _MessageTweak.PAD)
+ pad = tbc.encrypt(tweak, key, [0 for byte in range(0, BLOCK_BYTES)])
+
+ lower_part = LowPart(pad, padding_bytes*8)
+ C.append(XorState(M[l], lower_part))
+ tweak_final = TweakMessage(N, l+1, _MessageTweak.FINAL)
+ Final = tbc.encrypt(tweak_final, key, checksum)
+
+ return (Final, C)
+
+
+def TreatMessageDec(C, N, key):
+ checksum = [0 for byte in range(0, BLOCK_BYTES)]
+
+ l = len(C)//BLOCK_BYTES
+ padding_bytes = len(C)%BLOCK_BYTES
+
+ C = ArrayToBlockbytesMatrix(C)
+ M = []
+
+ for j in range(0, l):
+ tweak = TweakMessage(N, j, _MessageTweak.BLOCK)
+ M.append(tbc.decrypt(tweak, key, C[j]))
+ checksum = XorState(checksum, M[j])
+
+ if padding_bytes == 0:
+ tweak = TweakMessage(N, l, _MessageTweak.NO_PADDING)
+ Final = tbc.encrypt(tweak, key, checksum)
+
+ else:
+ tweak = TweakMessage(N, l, _MessageTweak.PAD)
+ pad = tbc.encrypt(tweak, key, [0 for byte in range(0, BLOCK_BYTES)])
+ lower_part = LowPart(pad, padding_bytes*8)
+ M.append(XorState(C[l], lower_part))
+
+ m_padded = Padding10LSB(M[l])
+ checksum = XorState(checksum, m_padded)
+ tweak_final = TweakMessage(N, l+1, _MessageTweak.FINAL)
+ Final = tbc.encrypt(tweak_final, key, checksum)
+
+ return (Final, M)
+
+
+################################################################################
+def encrypt(A, M, N, key):
+ K = list(key)
+
+ Auth = BuildAuth(TWEAK_BITS, A, K)
+ (Final, C) = TreatMessageEnc(M, N, K)
+ tag = XorState(Auth, Final)
+
+ return BlockbytesMatrixToBytes(C), bytes(tag)
+
+
+def decrypt(A, C, N, tag, key):
+ K = list(key)
+ tag = list(tag)
+
+ Auth = BuildAuth(TWEAK_BITS, A, K)
+ (Final, M) = TreatMessageDec(C, N, K)
+ tag2 = XorState(Auth, Final)
+
+ if tag != tag2:
+ raise TagValidationError(tag, tag2)
+
+ return BlockbytesMatrixToBytes(M)
diff --git a/src/add_python/lilliput/lilliput_ae_2.py b/src/add_python/lilliput/lilliput_ae_2.py
new file mode 100644
index 0000000..61aa86e
--- /dev/null
+++ b/src/add_python/lilliput/lilliput_ae_2.py
@@ -0,0 +1,127 @@
+"""
+ SCT 2 for lilliput ae 2
+"""
+
+from .constants import BLOCK_BYTES
+from .helpers import (
+ ArrayToBlockbytesMatrix,
+ BlockbytesMatrixToBytes,
+ BuildAuth,
+ Padding10LSB,
+ TagValidationError,
+ XorState
+)
+from . import tbc
+
+
+TWEAK_BITS = 128
+TWEAK_BYTES = TWEAK_BITS//8
+
+
+def TweakTag(j, padded):
+ tweak = [0 for byte in range(0, TWEAK_BYTES)]
+
+ tweak[TWEAK_BYTES - 1] |= ((j >> 120) & 0xf)
+ for byte in range(TWEAK_BYTES - 2, -1, -1):
+ tweak[byte] = (j >> (8 * byte)) & 0xff
+
+ if padded:
+ tweak[TWEAK_BYTES - 1] |= 0x40
+
+ return tweak
+
+
+def TweakTagEnd(N):
+ tweak = [0 for byte in range(0, TWEAK_BYTES)]
+
+ for byte in range(0, TWEAK_BYTES - 1):
+ tweak[byte] = N[byte]
+ tweak[TWEAK_BYTES - 1] = 0x10
+
+ return tweak
+
+
+def AddTagJ(tag, j):
+ array_j = [0 for byte in range(0, TWEAK_BYTES)]
+ for byte in range(0, TWEAK_BYTES):
+ array_j[byte] = (j >> (byte * 8))
+
+ xorr = XorState(tag, array_j)
+
+ xorr[TWEAK_BYTES - 1] |= 0x80
+
+ return xorr
+
+
+def MesssageAuthTag(M, N, Auth, key):
+ l = len(M)//BLOCK_BYTES
+ need_padding = len(M)%BLOCK_BYTES > 0
+
+ tag = list(Auth)
+ M = ArrayToBlockbytesMatrix(M)
+
+ for j in range(0, l):
+ tweak = TweakTag(j, False)
+ encryption = tbc.encrypt(tweak, key, M[j])
+ tag = XorState(tag, encryption)
+
+ if need_padding:
+ tweak = TweakTag(l, True)
+ m_padded = Padding10LSB(M[l])
+ encryption = tbc.encrypt(tweak, key, m_padded)
+ tag = XorState(tag, encryption)
+
+ tweak = TweakTagEnd(N)
+ encryption = tbc.encrypt(tweak, key, tag)
+ tag = encryption
+
+ return tag
+
+
+def MessageEncryption(M, N, tag, key):
+ l = len(M)//BLOCK_BYTES
+ need_padding = len(M)%BLOCK_BYTES > 0
+
+ M = ArrayToBlockbytesMatrix(M)
+ C = []
+
+ for j in range(0, l):
+ tweak = AddTagJ(tag, j)
+ padded_nonce = list(N) + [0x00]
+ encryption = tbc.encrypt(tweak, key, padded_nonce)
+ C.append(XorState(M[j], encryption))
+
+ if need_padding:
+ tweak = AddTagJ(tag, l)
+ padded_nonce = list(N) + [0x00]
+ encryption = tbc.encrypt(tweak, key, padded_nonce)
+ C.append(XorState(M[l], encryption))
+
+ return C
+
+
+################################################################################
+def encrypt(A, M, N, key):
+ K = list(key)
+
+ Auth = BuildAuth(TWEAK_BITS, A, K)
+ tag = MesssageAuthTag(M, N, Auth, K)
+ C = MessageEncryption(M, N, tag, K)
+
+ return BlockbytesMatrixToBytes(C), bytes(tag)
+
+
+def decrypt(A, C, N, tag, key):
+ K = list(key)
+ tag = list(tag)
+
+ M = BlockbytesMatrixToBytes(
+ MessageEncryption(C, N, tag, K)
+ )
+ Auth = BuildAuth(TWEAK_BITS, A, K)
+ tag2 = MesssageAuthTag(M, N, Auth, K)
+
+ if tag != tag2:
+ raise TagValidationError(tag, tag2)
+
+ return M
diff --git a/src/add_python/lilliput/multiplications.py b/src/add_python/lilliput/multiplications.py
new file mode 100644
index 0000000..c5f1e44
--- /dev/null
+++ b/src/add_python/lilliput/multiplications.py
@@ -0,0 +1,186 @@
+
+# Multiply by matrix M
+def _multiplyM(lane):
+ multiplied_lane = [lane[(byte-1) % 8] for byte in range(0, 8)]
+
+ multiplied_lane[2] ^= ((lane[6] << 2) & 0xff)
+ multiplied_lane[4] ^= ((lane[4] >> 3) & 0xff)
+ multiplied_lane[5] ^= ((lane[5] << 3) & 0xff)
+
+ return multiplied_lane
+
+def _multiplyM2(lane):
+ multiplied_lane = [lane[(byte-2) % 8] for byte in range(0, 8)]
+
+ multiplied_lane[2] ^= ((lane[5] << 2) & 0xff)
+ multiplied_lane[3] ^= ((lane[6] << 2) & 0xff)
+ multiplied_lane[4] ^= ((lane[3] >> 3) & 0xff) ^ ((lane[4] >> 6) & 0xff)
+ multiplied_lane[5] ^= ((lane[5] << 6) & 0xff)
+ multiplied_lane[6] ^= ((lane[5] << 3) & 0xff)
+
+ # binary matrix M1
+ multi_mat_l4_m1 = 0
+ l4 = lane[4]
+ multi_mat_l4_m1 ^= ((l4 & 0x8) >> 3)
+ multi_mat_l4_m1 ^= ((l4 & 0x10) >> 3)
+ multi_mat_l4_m1 ^= ((l4 & 0x20) >> 3)
+ multi_mat_l4_m1 ^= ((l4 & 0x40) >> 3) ^ ((l4 & 0x1) << 3)
+ multi_mat_l4_m1 ^= ((l4 & 0x80) >> 3) ^ ((l4 & 0x2) << 3)
+ multi_mat_l4_m1 ^= ((l4 & 0x04) << 3)
+ multi_mat_l4_m1 ^= ((l4 & 0x08) << 3)
+ multi_mat_l4_m1 ^= ((l4 & 0x10) << 3)
+
+ multiplied_lane[5] ^= multi_mat_l4_m1
+
+ return multiplied_lane
+
+
+def _multiplyM3(lane):
+ multiplied_lane = [lane[(byte-3) % 8] for byte in range(0, 8)]
+
+ multiplied_lane[2] ^= ((lane[4] << 2) & 0xff) ^ ((lane[5] << 5) & 0xff)
+ multiplied_lane[3] ^= ((lane[5] << 2) & 0xff)
+ multiplied_lane[4] ^= ((lane[2] >> 3) & 0xff) ^ ((lane[3] >> 6) & 0xff) ^ ((lane[6] << 2) & 0xff)
+ multiplied_lane[6] ^= ((lane[5] << 6) & 0xff)
+ multiplied_lane[7] ^= ((lane[5] << 3) & 0xff)
+
+ # binary matrix M1
+ multi_mat_l3_m1 = 0
+ l3 = lane[3]
+ multi_mat_l3_m1 ^= ((l3 & 0x8) >> 3)
+ multi_mat_l3_m1 ^= ((l3 & 0x10) >> 3)
+ multi_mat_l3_m1 ^= ((l3 & 0x20) >> 3)
+ multi_mat_l3_m1 ^= ((l3 & 0x40) >> 3) ^ ((l3 & 0x1) << 3)
+ multi_mat_l3_m1 ^= ((l3 & 0x80) >> 3) ^ ((l3 & 0x2) << 3)
+ multi_mat_l3_m1 ^= ((l3 & 0x04) << 3)
+ multi_mat_l3_m1 ^= ((l3 & 0x08) << 3)
+ multi_mat_l3_m1 ^= ((l3 & 0x10) << 3)
+
+ # binary matrix M1
+ multi_mat_l4_m1 = 0
+ l4 = lane[4]
+ multi_mat_l4_m1 ^= ((l4 & 0x8) >> 3)
+ multi_mat_l4_m1 ^= ((l4 & 0x10) >> 3)
+ multi_mat_l4_m1 ^= ((l4 & 0x20) >> 3)
+ multi_mat_l4_m1 ^= ((l4 & 0x40) >> 3) ^ ((l4 & 0x1) << 3)
+ multi_mat_l4_m1 ^= ((l4 & 0x80) >> 3) ^ ((l4 & 0x2) << 3)
+ multi_mat_l4_m1 ^= ((l4 & 0x04) << 3)
+ multi_mat_l4_m1 ^= ((l4 & 0x08) << 3)
+ multi_mat_l4_m1 ^= ((l4 & 0x10) << 3)
+
+ # binary matrix M2
+ multi_mat_l4_m2 = 0
+ l4 = lane[4]
+ multi_mat_l4_m2 ^= ((l4 & 0x40) >> 6)
+ multi_mat_l4_m2 ^= ((l4 & 0x80) >> 6)
+ multi_mat_l4_m2 ^= (l4 & 0x08)
+ multi_mat_l4_m2 ^= (l4 & 0x10)
+ multi_mat_l4_m2 ^= (l4 & 0x20)
+ multi_mat_l4_m2 ^= (l4 & 0x40) ^ ((l4 & 0x1) << 6)
+ multi_mat_l4_m2 ^= (l4 & 0x80) ^ ((l4 & 0x2) << 6)
+
+
+ multiplied_lane[5] ^= multi_mat_l3_m1 ^ multi_mat_l4_m2
+ multiplied_lane[6] ^= multi_mat_l4_m1
+
+ return multiplied_lane
+
+
+def _multiplyMR(lane):
+ multiplied_lane = [lane[(byte+1) % 8] for byte in range(0, 8)]
+
+ multiplied_lane[2] ^= ((lane[4] >> 3) & 0xff)
+ multiplied_lane[4] ^= ((lane[6] << 3) & 0xff)
+ multiplied_lane[5] ^= ((lane[3] << 2) & 0xff)
+
+ return multiplied_lane
+
+
+def _multiplyMR2(lane):
+ multiplied_lane = [lane[(byte+2) % 8] for byte in range(0, 8)]
+
+ multiplied_lane[1] ^= ((lane[4] >> 3) & 0xff)
+ multiplied_lane[2] ^= ((lane[5] >> 3) & 0xff)
+ multiplied_lane[3] ^= ((lane[6] << 3) & 0xff)
+ multiplied_lane[4] ^= ((lane[3] << 2) & 0xff) ^ ((lane[7] << 3) & 0xff)
+ multiplied_lane[5] ^= ((lane[4] << 2) & 0xff)
+
+
+ # binary matrix m3
+ multi_mat_l6_m3 = 0
+ l6 = lane[6]
+ multi_mat_l6_m3 ^= (l6 & 0x1)
+ multi_mat_l6_m3 ^= (l6 & 0x2)
+ multi_mat_l6_m3 ^= (l6 & 0x4)
+ multi_mat_l6_m3 ^= (l6 & 0x8)
+ multi_mat_l6_m3 ^= (l6 & 0x10)
+
+
+ multiplied_lane[2] ^= multi_mat_l6_m3
+
+ return multiplied_lane
+
+def _multiplyMR3(lane):
+ multiplied_lane = [lane[(byte+3) % 8] for byte in range(0, 8)]
+
+ multiplied_lane[0] ^= ((lane[4] >> 3) & 0xff)
+ multiplied_lane[1] ^= ((lane[5] >> 3) & 0xff)
+ multiplied_lane[3] ^= ((lane[3] << 2) & 0xff) ^ ((lane[7] << 3) & 0xff)
+ multiplied_lane[4] ^= ((lane[0] << 3) & 0xff) ^ ((lane[4] << 2) & 0xff)
+ multiplied_lane[5] ^= ((lane[5] << 2) & 0xff) ^ ((lane[6] << 5) & 0xff)
+
+ # binary matrix m3
+ multi_mat_l6_m3 = 0
+ l6 = lane[6]
+ multi_mat_l6_m3 ^= (l6 & 0x1)
+ multi_mat_l6_m3 ^= (l6 & 0x2)
+ multi_mat_l6_m3 ^= (l6 & 0x4)
+ multi_mat_l6_m3 ^= (l6 & 0x8)
+ multi_mat_l6_m3 ^= (l6 & 0x10)
+
+ # binary matrix m3
+ multi_mat_l7_m3 = 0
+ l7 = lane[7]
+ multi_mat_l7_m3 ^= (l7 & 0x1)
+ multi_mat_l7_m3 ^= (l7 & 0x2)
+ multi_mat_l7_m3 ^= (l7 & 0x4)
+ multi_mat_l7_m3 ^= (l7 & 0x8)
+ multi_mat_l7_m3 ^= (l7 & 0x10)
+
+ # binary matrix m4
+ multi_mat_l3_m4 = 0
+ l3 = lane[3]
+ multi_mat_l3_m4 ^= ((l3 & 0x2) >> 1)
+ multi_mat_l3_m4 ^= ((l3 & 0x4) >> 1)
+ multi_mat_l3_m4 ^= ((l3 & 0x8) >> 1)
+ multi_mat_l3_m4 ^= ((l3 & 0x10) >> 1)
+ multi_mat_l3_m4 ^= ((l3 & 0x20) >> 1)
+
+ # binary matrix m1 for MR
+ multi_mat_l6_m1 = 0
+ l6 = lane[6]
+ multi_mat_l6_m1 ^= ((l6 & 0x8) >> 3)
+ multi_mat_l6_m1 ^= ((l6 & 0x10) >> 3)
+ multi_mat_l6_m1 ^= ((l6 & 0x20) >> 3)
+ multi_mat_l6_m1 ^= ((l6 & 0x40) >> 3) ^ ((l6 & 0x1) << 3)
+ multi_mat_l6_m1 ^= ((l6 & 0x80) >> 3) ^ ((l6 & 0x2) << 3)
+ multi_mat_l6_m1 ^= ((l6 & 0x4) << 3)
+ multi_mat_l6_m1 ^= ((l6 & 0x8) << 3)
+ multi_mat_l6_m1 ^= ((l6 & 0x10) << 3)
+
+
+ multiplied_lane[1] ^= multi_mat_l6_m3
+ multiplied_lane[2] ^= multi_mat_l3_m4 ^ multi_mat_l6_m1 ^ multi_mat_l7_m3
+
+ return multiplied_lane
+
+
+ALPHAS = (
+ list, # Identity.
+ _multiplyM,
+ _multiplyM2,
+ _multiplyM3,
+ _multiplyMR,
+ _multiplyMR2,
+ _multiplyMR3
+)
diff --git a/src/add_python/lilliput/tbc.py b/src/add_python/lilliput/tbc.py
new file mode 100644
index 0000000..5291994
--- /dev/null
+++ b/src/add_python/lilliput/tbc.py
@@ -0,0 +1,178 @@
+"""
+ Lilliput TBC
+"""
+from .constants import BLOCK_BYTES, Sbox
+from .multiplications import ALPHAS
+
+
+permutation = [14, 11, 12, 10, 8, 9, 13, 15, 3, 1, 4, 5, 6, 0, 2, 7]
+permutationInv = [13, 9, 14, 8, 10, 11, 12, 15, 4, 5, 3, 1, 2, 6 ,0 ,7]
+
+################################################################################
+
+def BuildTweakey(tweak, key):
+ return tweak+key
+
+#############################
+
+def _lane(TK, j):
+ return TK[j*8:(j+1)*8]
+
+def RoundTweakeySchedule(tweakey):
+ p = len(tweakey)//8
+
+ multiplied_lanes = (
+ ALPHAS[j](_lane(tweakey, j)) for j in range(p)
+ )
+
+ return [byte for lane in multiplied_lanes for byte in lane]
+
+
+def SubTweakeyExtract(tweakey, Ci):
+ RTKi = [0]*8
+
+ for j, byte in enumerate(tweakey):
+ RTKi[j%8] ^= byte
+
+ RTKi[0] ^= Ci
+
+ return RTKi
+
+
+def TweakeyScheduleWhole(tweakey, r):
+ # store main tweakey in TKs[0]
+ # and corresponding round tweakey in RTKs[0]
+ TKs = [tweakey]
+ RTKs = [SubTweakeyExtract(TKs[0], 0)]
+
+ for i in range(1, r):
+ TKs.append(RoundTweakeySchedule(TKs[i-1]))
+ RTKs.append(SubTweakeyExtract(TKs[i], i))
+
+ return RTKs
+
+
+################################################################################
+
+def NonLinearLayer(state, subtweakey):
+
+ variables_xored = [0 for byte in range(0, 8)]
+ for byte in range(0,8):
+ variables_xored[byte] = state[byte] ^ subtweakey[byte]
+
+ variables_sboxed = [0 for byte in range(0, 8)]
+ for byte in range(0, 8):
+ variables_sboxed[byte] = Sbox[variables_xored[byte]]
+
+ state_output = [0 for byte in range(0, BLOCK_BYTES)]
+ for byte in range(0,BLOCK_BYTES):
+ state_output[byte] = state[byte]
+ for byte in range(0, 8):
+ state_output[15 - byte] ^= variables_sboxed[byte]
+
+ return state_output
+
+
+def LinearLayer(state):
+ state_output = [0 for byte in range(0, BLOCK_BYTES)]
+ for byte in range(0, BLOCK_BYTES):
+ state_output[byte] = state[byte]
+
+ for byte in range(1, 8):
+ state_output[15] ^= state[byte]
+
+ for byte in range(9, 15):
+ state_output[byte] ^= state[7]
+
+ return state_output
+
+
+def PermutationLayerEnc(state):
+ state_output = [0 for byte in range(0, BLOCK_BYTES)]
+ for byte in range(0, BLOCK_BYTES):
+ state_output[byte] = state[permutation[byte]]
+
+ return state_output
+
+def PermutationLayerDec(state):
+ state_output = [0 for byte in range(0, BLOCK_BYTES)]
+ for byte in range(0, BLOCK_BYTES):
+ state_output[byte] = state[permutationInv[byte]]
+
+ return state_output
+
+
+def OneRoundEGFNEnc(state, subtweakey):
+ state_non_linear = NonLinearLayer(state, subtweakey)
+ state_linear = LinearLayer(state_non_linear)
+ state_permutation = PermutationLayerEnc(state_linear)
+
+ return state_permutation
+
+def LastRoundEGFN(state, subtweakey):
+ state_non_linear = NonLinearLayer(state, subtweakey)
+ state_linear = LinearLayer(state_non_linear)
+
+ return state_linear
+
+
+def OneRoundEGFNDec(state, subtweakey):
+ state_non_linear = NonLinearLayer(state, subtweakey)
+ state_linear = LinearLayer(state_non_linear)
+ state_permutation = PermutationLayerDec(state_linear)
+
+ return state_permutation
+
+
+def _rounds(key_bytes):
+ rounds = {
+ 128: 32,
+ 192: 36,
+ 256: 42
+ }
+ return rounds[key_bytes*8]
+
+
+################################################################################
+# Lilliput TBC
+
+def encrypt(tweak, key, message):
+ r = _rounds(len(key))
+
+ tweakey = BuildTweakey(tweak, key)
+ RTKs = TweakeyScheduleWhole(tweakey, r)
+
+ state = [0 for byte in range(0, BLOCK_BYTES)]
+ for byte in range(0, BLOCK_BYTES):
+ state[byte] = message[byte]
+
+ for i in range(0, r-1):
+ state_output = OneRoundEGFNEnc(state, RTKs[i])
+
+ for byte in range(0, BLOCK_BYTES):
+ state[byte] = state_output[byte]
+
+ state_output = LastRoundEGFN(state, RTKs[r-1])
+
+ return state_output
+
+
+def decrypt(tweak, key, cipher):
+ r = _rounds(len(key))
+
+ tweakey = BuildTweakey(tweak, key)
+ RTKs = TweakeyScheduleWhole(tweakey, r)
+
+ state = [0 for byte in range(0, BLOCK_BYTES)]
+ for byte in range(0, BLOCK_BYTES):
+ state[byte] = cipher[byte]
+
+ for i in range(0, r-1):
+ state_output = OneRoundEGFNDec(state, RTKs[r-i-1])
+
+ for byte in range(0, BLOCK_BYTES):
+ state[byte] = state_output[byte]
+
+ state_output = LastRoundEGFN(state, RTKs[0])
+
+ return state_output
generated by cgit v1.2.3 (git 2.39.1) at 2025-12-06 16:45:42 +0000