ecc_lookup_relation.hpp

Go to the documentation of this file.

// === AUDIT STATUS ===
// internal:    { status: not started, auditors: [], date: YYYY-MM-DD }
// external_1:  { status: not started, auditors: [], date: YYYY-MM-DD }
// external_2:  { status: not started, auditors: [], date: YYYY-MM-DD }
// =====================
 
#pragma once
#include <array>
#include <tuple>
 
#include "barretenberg/common/constexpr_utils.hpp"
#include "barretenberg/polynomials/univariate.hpp"
#include "barretenberg/relations/relation_types.hpp"
 
namespace bb {
 
template <typename FF_> class ECCVMLookupRelationImpl {
  public:
    using FF = FF_;
    static constexpr size_t READ_TERMS = 4;
    static constexpr size_t WRITE_TERMS = 2;
    // 1 + polynomial degree of this relation
    static constexpr size_t LENGTH = READ_TERMS + WRITE_TERMS + 3; // 9
 
    static constexpr std::array<size_t, 2> SUBRELATION_PARTIAL_LENGTHS{
        LENGTH, // grand product construction sub-relation
        LENGTH  // left-shiftable polynomial sub-relation
    };
 
    static constexpr std::array<bool, 2> SUBRELATION_LINEARLY_INDEPENDENT = { true, false };
 
    template <typename AllValues> static bool operation_exists_at_row(const AllValues& row)
 
    {
        return (row.msm_add == 1) || (row.msm_skew == 1) || (row.precompute_select == 1);
    }
 
    template <typename AllEntities> static auto& get_inverse_polynomial(AllEntities& in) { return in.lookup_inverses; }
 
    template <typename Accumulator, typename AllEntities>
    static Accumulator compute_inverse_exists(const AllEntities& in)
    {
        using View = typename Accumulator::View;
 
        const auto row_has_write = View(in.precompute_select);
        const auto row_has_read = View(in.msm_add) + View(in.msm_skew);
        return row_has_write + row_has_read - (row_has_write * row_has_read);
    }
 
    template <typename Accumulator, size_t index, typename AllEntities>
    static Accumulator lookup_read_counts(const AllEntities& in)
    {
        using View = typename Accumulator::View;
 
        if constexpr (index == 0) {
            return Accumulator(View(in.lookup_read_counts_0));
        }
        if constexpr (index == 1) {
            return Accumulator(View(in.lookup_read_counts_1));
        }
        return Accumulator(1);
    }
 
    template <typename Accumulator, size_t read_index, typename AllEntities>
    static Accumulator compute_read_term_predicate(const AllEntities& in)
 
    {
        using View = typename Accumulator::View;
 
        if constexpr (read_index == 0) {
            return Accumulator(View(in.msm_add1));
        }
        if constexpr (read_index == 1) {
            return Accumulator(View(in.msm_add2));
        }
        if constexpr (read_index == 2) {
            return Accumulator(View(in.msm_add3));
        }
        if constexpr (read_index == 3) {
            return Accumulator(View(in.msm_add4));
        }
        return Accumulator(1);
    }
 
    template <typename Accumulator, size_t write_index, typename AllEntities>
    static Accumulator compute_write_term_predicate(const AllEntities& in)
    {
        using View = typename Accumulator::View;
 
        if constexpr (write_index == 0) {
            return Accumulator(View(in.precompute_select));
        }
        if constexpr (write_index == 1) {
            // TODO(https://github.com/AztecProtocol/barretenberg/issues/750) Is this a bug?
            return Accumulator(View(in.precompute_select));
        }
        return Accumulator(1);
    }
    template <typename Accumulator, size_t write_index, typename AllEntities, typename Parameters>
    static Accumulator compute_write_term(const AllEntities& in, const Parameters& params)
    {
        using View = typename Accumulator::View;
 
        static_assert(write_index < WRITE_TERMS);
        // write_index == 0 means our wNAF digit is positive (i.e., ∈{1, 3..., 15}).
        // write_index == 1 means our wNAF digit is negative (i.e., ∈{-15, -13..., -1})
 
        // round starts at 0 and increments to 7
        // point starts at 15[P] and decrements to [P]
        // a slice value of 0 maps to -15[P]
 
        // we have computed `(15 - 2 * round)[P] =: (precompute_tx, precompute_ty)`.
        // `round`∈{0, 1..., 7}
        // if write_index == 0, we want to write (pc, 15 - 2 * round, precompute_tx, precompute_ty)
        // if write_index == 1, we want to write (pc, round, precompute_tx, -precompute_ty)
        // to sum up, both:
        //      (pc, round, precompute_tx, -precompute_ty) _and_
        //      (pc, 15 - 2 * round, precompute_tx, precompute_ty)
        // will be written to the lookup table.
        //
        // therefore, if `pc` corresponds to the elliptic curve point [P], we will write:
        // | pc | 0  | -15[P].x | -15[P].y |
        // | pc | 1  | -13[P].x | -13[P].y |
        // | pc | 2  | -11[P].x | -11[P].y |
        // | pc | 3  | -9[P].x  | -9[P].y  |
        // | pc | 4  | -7[P].x  | -7[P].y  |
        // | pc | 5  | -5[P].x  | -5[P].y  |
        // | pc | 6  | -3[P].x  | -3[P].y  |
        // | pc | 7  | -1[P].x  | -1[P].y  |
        // | pc | 8  |   [P].x  |   [P].y  |
        // | pc | 9  |  3[P].x  |  3[P].y  |
        // | pc | 10 |  5[P].x  |  5[P].y  |
        // | pc | 11 |  7[P].x  |  7[P].y  |
        // | pc | 12 |  9[P].x  |  9[P].y  |
        // | pc | 13 | 11[P].x  | 11[P].y  |
        // | pc | 14 | 13[P].x  | 13[P].y  |
        // | pc | 15 | 15[P].x  | 15[P].y  |
 
        const auto& precompute_pc = View(in.precompute_pc);
        const auto& tx = View(in.precompute_tx);
        const auto& ty = View(in.precompute_ty);
        const auto& precompute_round = View(in.precompute_round);
        const auto& gamma = params.gamma;
        const auto& beta = params.beta;
        const auto& beta_sqr = params.beta_sqr;
        const auto& beta_cube = params.beta_cube;
 
        if constexpr (write_index == 0) {
            const auto positive_slice_value = -(precompute_round) + 15;
            const auto positive_term =
                precompute_pc + gamma + positive_slice_value * beta + tx * beta_sqr + ty * beta_cube;
            return positive_term; // degree 1
        }
        if constexpr (write_index == 1) {
            const auto negative_term = precompute_pc + gamma + precompute_round * beta + tx * beta_sqr - ty * beta_cube;
            return negative_term; // degree 1
        }
        return Accumulator(1);
    }
 
    template <typename Accumulator, size_t read_index, typename AllEntities, typename Parameters>
    static Accumulator compute_read_term(const AllEntities& in, const Parameters& params)
    {
        using View = typename Accumulator::View;
 
        // read term: (pc, compressed_slice, (2 * compressed_slice - 15)[P])
        // (the latter term is of course represented via an x and y coordinate.)
        static_assert(read_index < READ_TERMS);
        const auto& gamma = params.gamma;
        const auto& beta = params.beta;
        const auto& beta_sqr = params.beta_sqr;
        const auto& beta_cube = params.beta_cube;
        const auto& msm_pc = View(in.msm_pc);
        const auto& msm_count = View(in.msm_count);
        const auto& msm_slice1 = View(in.msm_slice1);
        const auto& msm_slice2 = View(in.msm_slice2);
        const auto& msm_slice3 = View(in.msm_slice3);
        const auto& msm_slice4 = View(in.msm_slice4);
        const auto& msm_x1 = View(in.msm_x1);
        const auto& msm_x2 = View(in.msm_x2);
        const auto& msm_x3 = View(in.msm_x3);
        const auto& msm_x4 = View(in.msm_x4);
        const auto& msm_y1 = View(in.msm_y1);
        const auto& msm_y2 = View(in.msm_y2);
        const auto& msm_y3 = View(in.msm_y3);
        const auto& msm_y4 = View(in.msm_y4);
 
        // Recall that `pc` stands for point-counter. We recall how to compute the current pc.
        //
        // row pc = value of pc after msm
        // msm_count = number of (128-bit) multiplications processed so far in current MSM round (NOT INCLUDING current
        // row) current_pc = msm_pc - msm_count next_pc = current_pc - {0, 1, 2, 3}, depending on how many adds are
        // performed in the current row.
        const auto current_pc = msm_pc - msm_count;
 
        if constexpr (read_index == 0) {
            const auto read_term1 = (current_pc) + gamma + msm_slice1 * beta + msm_x1 * beta_sqr + msm_y1 * beta_cube;
            return read_term1; // degree 1
        }
        if constexpr (read_index == 1) {
            const auto read_term2 =
                (current_pc - 1) + gamma + msm_slice2 * beta + msm_x2 * beta_sqr + msm_y2 * beta_cube;
            return read_term2; // degree 1
        }
        if constexpr (read_index == 2) {
            const auto read_term3 =
                (current_pc - 2) + gamma + msm_slice3 * beta + msm_x3 * beta_sqr + msm_y3 * beta_cube;
            return read_term3; // degree 1
        }
        if constexpr (read_index == 3) {
            const auto read_term4 =
                (current_pc - 3) + gamma + msm_slice4 * beta + msm_x4 * beta_sqr + msm_y4 * beta_cube;
            return read_term4; // degree 1
        }
        return Accumulator(1);
    }
 
    template <typename ContainerOverSubrelations, typename AllEntities, typename Parameters>
    static void accumulate(ContainerOverSubrelations& accumulator,
                           const AllEntities& in,
                           const Parameters& params,
                           const FF& scaling_factor);
};
 
template <typename FF> using ECCVMLookupRelation = Relation<ECCVMLookupRelationImpl<FF>>;
 
} // namespace bb