keccak.hpp

Go to the documentation of this file.

// === AUDIT STATUS ===
// internal:    { status: Complete, auditors: [Nishat], commit: }
// external_1:  { status: not started, auditors: [], commit: }
// external_2:  { status: not started, auditors: [], commit: }
// =====================
 
#pragma once
#include "barretenberg/stdlib/primitives/byte_array/byte_array.hpp"
#include <array>
 
namespace bb::stdlib {
 
template <typename Builder> class keccak {
  public:
    using witness_ct = stdlib::witness_t<Builder>;
    using field_ct = stdlib::field_t<Builder>;
    using bool_ct = stdlib::bool_t<Builder>;
    using byte_array_ct = stdlib::byte_array<Builder>;
 
    // base of extended representation we use for efficient logic operations
    static constexpr uint256_t BASE = 11;
 
    static constexpr size_t NUM_KECCAK_ROUNDS = 24;
 
    // 1 "lane" = 64 bits. Instead of interpreting the keccak sponge as 1,600 bits, it's easier to work over 64-bit
    // "lanes". 1,600 / 64 = 25.
    static constexpr size_t NUM_KECCAK_LANES = 25;
    static constexpr size_t KECCAK_LANE_SIZE = 64;
 
    // round constants. Used in IOTA round
    static constexpr std::array<uint64_t, NUM_KECCAK_ROUNDS> RC = {
        0x0000000000000001, 0x0000000000008082, 0x800000000000808a, 0x8000000080008000, 0x000000000000808b,
        0x0000000080000001, 0x8000000080008081, 0x8000000000008009, 0x000000000000008a, 0x0000000000000088,
        0x0000000080008009, 0x000000008000000a, 0x000000008000808b, 0x800000000000008b, 0x8000000000008089,
        0x8000000000008003, 0x8000000000008002, 0x8000000000000080, 0x000000000000800a, 0x800000008000000a,
        0x8000000080008081, 0x8000000000008080, 0x0000000080000001, 0x8000000080008008
    };
 
    // Rotation offsets, y vertically, x horizontally: r[y * 5 + x]
    static constexpr std::array<size_t, NUM_KECCAK_LANES> ROTATIONS = {
        0, 1, 62, 28, 27, 36, 44, 6, 55, 20, 3, 10, 43, 25, 39, 41, 45, 15, 21, 8, 18, 2, 61, 56, 14,
    };
 
    static constexpr uint256_t convert_to_sparse(uint256_t input)
    {
        std::array<uint64_t, 64> out_bits;
        size_t count = 0;
        while (input > 0) {
            uint64_t bit = static_cast<uint64_t>(input & 1);
            out_bits[count++] = bit;
            input = input >> 1;
        }
        uint256_t output = 0;
        for (size_t i = 0; i < count; ++i) {
            output *= BASE;
            output += out_bits[count - 1 - i];
        }
        return output;
    };
 
    static constexpr uint256_t normalize_sparse(uint256_t input)
    {
        std::array<uint64_t, 64> out_bits;
        size_t count = 0;
        while (input > 0) {
            const auto [quotient, slice] = input.divmod(BASE);
            uint64_t bit = static_cast<uint64_t>(slice) & 1;
            out_bits[count++] = bit;
            input = quotient;
        }
        uint256_t out;
        for (size_t i = 0; i < count; ++i) {
            out *= BASE;
            out += out_bits[count - 1 - i];
        }
        return out;
    }
 
    static constexpr std::array<uint256_t, NUM_KECCAK_ROUNDS> get_sparse_round_constants()
    {
        std::array<uint256_t, NUM_KECCAK_ROUNDS> output;
        for (size_t i = 0; i < NUM_KECCAK_ROUNDS; ++i) {
            output[i] = convert_to_sparse(RC[i]);
        }
        return output;
    }
    static constexpr std::array<uint256_t, NUM_KECCAK_ROUNDS> SPARSE_RC = get_sparse_round_constants();
 
    static constexpr uint256_t get_chi_offset()
    {
        uint256_t result = 0;
        for (size_t i = 0; i < 64; ++i) {
            result *= 11;
            result += 1;
        }
        return result;
    }
    static constexpr uint256_t CHI_OFFSET = get_chi_offset();
 
    struct keccak_state {
        std::array<field_ct, NUM_KECCAK_LANES> state;
        std::array<field_ct, NUM_KECCAK_LANES> state_msb;
        std::array<field_ct, NUM_KECCAK_LANES> twisted_state;
        Builder* context;
    };
 
    template <size_t lane_index> static field_t<Builder> normalize_and_rotate(const field_ct& limb, field_ct& msb);
    static void compute_twisted_state(keccak_state& internal);
    static void theta(keccak_state& state);
    static void rho(keccak_state& state);
    static void pi(keccak_state& state);
    static void chi(keccak_state& state);
    static void iota(keccak_state& state, size_t round);
 
    static void keccakf1600(keccak_state& state);
 
    static std::vector<uint8_t> hash_native(const std::vector<uint8_t>& data)
    {
        auto hash_result = ethash_keccak256(&data[0], data.size());
 
        std::vector<uint8_t> output;
        output.resize(32);
 
        memcpy((void*)&output[0], (void*)&hash_result.word64s[0], 32);
        return output;
    }
 
    // exposing keccak f1600 permutation
 
    static std::array<field_ct, NUM_KECCAK_LANES> permutation_opcode(std::array<field_ct, NUM_KECCAK_LANES> state,
                                                                     Builder* context);
    static std::array<field_ct, NUM_KECCAK_LANES> extended_2_normal(keccak_state& internal);
};
 
} // namespace bb::stdlib