124 lines
3.3 KiB
C
124 lines
3.3 KiB
C
// keccak.c
|
|
// 19-Nov-11 Markku-Juhani O. Saarinen <mjos@iki.fi>
|
|
// A baseline Keccak (3rd round) implementation.
|
|
|
|
#include <stdint.h>
|
|
#include <string.h>
|
|
|
|
#define HASH_DATA_AREA 136
|
|
#define KECCAK_ROUNDS 24
|
|
|
|
#ifndef ROTL64
|
|
#define ROTL64(x, y) (((x) << (y)) | ((x) >> (64 - (y))))
|
|
#endif
|
|
|
|
const uint64_t keccakf_rndc[24] =
|
|
{
|
|
0x0000000000000001, 0x0000000000008082, 0x800000000000808a,
|
|
0x8000000080008000, 0x000000000000808b, 0x0000000080000001,
|
|
0x8000000080008081, 0x8000000000008009, 0x000000000000008a,
|
|
0x0000000000000088, 0x0000000080008009, 0x000000008000000a,
|
|
0x000000008000808b, 0x800000000000008b, 0x8000000000008089,
|
|
0x8000000000008003, 0x8000000000008002, 0x8000000000000080,
|
|
0x000000000000800a, 0x800000008000000a, 0x8000000080008081,
|
|
0x8000000000008080, 0x0000000080000001, 0x8000000080008008
|
|
};
|
|
|
|
const int keccakf_rotc[24] =
|
|
{
|
|
1, 3, 6, 10, 15, 21, 28, 36, 45, 55, 2, 14,
|
|
27, 41, 56, 8, 25, 43, 62, 18, 39, 61, 20, 44
|
|
};
|
|
|
|
const int keccakf_piln[24] =
|
|
{
|
|
10, 7, 11, 17, 18, 3, 5, 16, 8, 21, 24, 4,
|
|
15, 23, 19, 13, 12, 2, 20, 14, 22, 9, 6, 1
|
|
};
|
|
|
|
// update the state with given number of rounds
|
|
|
|
void keccakf(uint64_t st[25], int rounds)
|
|
{
|
|
int i, j, round;
|
|
uint64_t t, bc[5];
|
|
|
|
for (round = 0; round < rounds; ++round) {
|
|
|
|
// Theta
|
|
bc[0] = st[0] ^ st[5] ^ st[10] ^ st[15] ^ st[20];
|
|
bc[1] = st[1] ^ st[6] ^ st[11] ^ st[16] ^ st[21];
|
|
bc[2] = st[2] ^ st[7] ^ st[12] ^ st[17] ^ st[22];
|
|
bc[3] = st[3] ^ st[8] ^ st[13] ^ st[18] ^ st[23];
|
|
bc[4] = st[4] ^ st[9] ^ st[14] ^ st[19] ^ st[24];
|
|
|
|
for (i = 0; i < 5; ++i) {
|
|
t = bc[(i + 4) % 5] ^ ROTL64(bc[(i + 1) % 5], 1);
|
|
st[i ] ^= t;
|
|
st[i + 5] ^= t;
|
|
st[i + 10] ^= t;
|
|
st[i + 15] ^= t;
|
|
st[i + 20] ^= t;
|
|
}
|
|
|
|
// Rho Pi
|
|
t = st[1];
|
|
for (i = 0; i < 24; ++i) {
|
|
bc[0] = st[keccakf_piln[i]];
|
|
st[keccakf_piln[i]] = ROTL64(t, keccakf_rotc[i]);
|
|
t = bc[0];
|
|
}
|
|
|
|
// Chi
|
|
for (j = 0; j < 25; j += 5) {
|
|
bc[0] = st[j ];
|
|
bc[1] = st[j + 1];
|
|
bc[2] = st[j + 2];
|
|
bc[3] = st[j + 3];
|
|
bc[4] = st[j + 4];
|
|
st[j ] ^= (~bc[1]) & bc[2];
|
|
st[j + 1] ^= (~bc[2]) & bc[3];
|
|
st[j + 2] ^= (~bc[3]) & bc[4];
|
|
st[j + 3] ^= (~bc[4]) & bc[0];
|
|
st[j + 4] ^= (~bc[0]) & bc[1];
|
|
}
|
|
|
|
// Iota
|
|
st[0] ^= keccakf_rndc[round];
|
|
}
|
|
}
|
|
|
|
// compute a keccak hash (md) of given byte length from "in"
|
|
typedef uint64_t state_t[25];
|
|
|
|
void keccak(const uint8_t *in, int inlen, uint8_t *md, int mdlen)
|
|
{
|
|
state_t st;
|
|
uint8_t temp[144];
|
|
int i, rsiz, rsizw;
|
|
|
|
rsiz = sizeof(state_t) == mdlen ? HASH_DATA_AREA : 200 - 2 * mdlen;
|
|
rsizw = rsiz / 8;
|
|
|
|
memset(st, 0, sizeof(st));
|
|
|
|
for ( ; inlen >= rsiz; inlen -= rsiz, in += rsiz) {
|
|
for (i = 0; i < rsizw; i++)
|
|
st[i] ^= ((uint64_t *) in)[i];
|
|
keccakf(st, KECCAK_ROUNDS);
|
|
}
|
|
|
|
// last block and padding
|
|
memcpy(temp, in, inlen);
|
|
temp[inlen++] = 1;
|
|
memset(temp + inlen, 0, rsiz - inlen);
|
|
temp[rsiz - 1] |= 0x80;
|
|
|
|
for (i = 0; i < rsizw; i++)
|
|
st[i] ^= ((uint64_t *) temp)[i];
|
|
|
|
keccakf(st, KECCAK_ROUNDS);
|
|
|
|
memcpy(md, st, mdlen);
|
|
}
|