ab_merkle_tree/

balanced.rs

use crate::{hash_pair, hash_pair_block, hash_pairs};
use ab_blake3::{BLOCK_LEN, OUT_LEN};
#[cfg(feature = "alloc")]
use alloc::boxed::Box;
use array_reshape::Flatten;
use core::iter::TrustedLen;
use core::mem::MaybeUninit;
use core::num::NonZero;

/// Optimal number of blocks for hashing at once to saturate BLAKE3 SIMD on any hardware
const BATCH_HASH_NUM_BLOCKS: usize = 16;
/// Number of leaves that corresponds to [`BATCH_HASH_NUM_BLOCKS`]
const BATCH_HASH_NUM_LEAVES: usize = BATCH_HASH_NUM_BLOCKS * BLOCK_LEN / OUT_LEN;

const TREE_SIZE_WITHOUT_LEAVES<const N: usize>: usize = {
    // Ensuring only supported `N` can be specified for [`BalancedMerkleTree`]
    assert!(
        N.is_power_of_two(),
        "Balanced Merkle Tree must have a number of leaves that is a power of 2"
    );

    assert!(
        N > 1,
        "This Balanced Merkle Tree must have more than one leaf"
    );

    N - 1
};
/// Number of elements in a proof for a tree with `N` leaves
pub const PROOF_ELEMENTS<const N: usize>: usize = N.ilog2() as usize;
const ROOT_ONLY_LARGE_STACK_SIZE<const N: usize>: usize = {
    // For small trees the large stack is not used, so the returned value does not matter as long as
    // it compiles
    if N < BATCH_HASH_NUM_LEAVES {
        1
    } else {
        (N / BATCH_HASH_NUM_LEAVES).ilog2() as usize + 1
    }
};

/// Merkle Tree variant that has hash-sized leaves and is fully balanced according to configured
/// generic parameter.
///
/// This can be considered a general case of [`UnbalancedMerkleTree`]. The root and proofs are
/// identical for both in case the number of leaves is a power of two. For the number of leaves that
/// is a power of two [`UnbalancedMerkleTree`] is useful when a single proof needs to be generated
/// and the number of leaves is very large (it can generate proofs with very little RAM usage
/// compared to this version).
///
/// [`UnbalancedMerkleTree`]: crate::unbalanced::UnbalancedMerkleTree
///
/// This Merkle Tree implementation is best suited for use cases when proofs for all (or most) of
/// the elements need to be generated and the whole tree easily fits into memory. It can also be
/// constructed and proofs can be generated efficiently without heap allocations.
///
/// With all parameters of the tree known statically, it results in the most efficient version of
/// the code being generated for a given set of parameters.
#[derive(Debug)]
pub struct BalancedMerkleTree<'a, const N: usize> {
    leaves: &'a [[u8; OUT_LEN]; N],
    // This tree doesn't include leaves because we have them in `leaves` field
    tree: [[u8; OUT_LEN]; TREE_SIZE_WITHOUT_LEAVES::<N>],
}

// TODO: Optimize by implementing SIMD-accelerated hashing of multiple values:
//  https://github.com/BLAKE3-team/BLAKE3/issues/478
impl<'a, const N: usize> BalancedMerkleTree<'a, N>
where
    [(); TREE_SIZE_WITHOUT_LEAVES::<N>]:,
{
    /// Create a new tree from a fixed set of elements.
    ///
    /// The data structure is statically allocated and might be too large to fit on the stack!
    /// If that is the case, use `new_boxed()` method.
    // TODO: Unlock on RISC-V, it started failing since https://github.com/nazar-pc/abundance/pull/551
    //  for unknown reason
    #[cfg_attr(
        all(feature = "no-panic", not(target_arch = "riscv64")),
        no_panic::no_panic
    )]
    pub fn new(leaves: &'a [[u8; OUT_LEN]; N]) -> Self {
        // TODO: Technically redundant with the where-clause, which is not enforced in generic code
        //  due to https://github.com/rust-lang/rust/issues/164019, remove once fixed
        const { TREE_SIZE_WITHOUT_LEAVES::<N> };

        let mut tree = [MaybeUninit::<[u8; OUT_LEN]>::uninit(); _];

        Self::init_internal(leaves, &mut tree);

        Self {
            leaves,
            // SAFETY: Statically guaranteed for all elements to be initialized
            tree: unsafe { tree.transpose().assume_init() },
        }
    }

    /// Like [`Self::new()`], but used pre-allocated memory for instantiation
    // TODO: Unlock on RISC-V, it started failing since https://github.com/nazar-pc/abundance/pull/551
    //  for unknown reason
    #[cfg_attr(
        all(feature = "no-panic", not(target_arch = "riscv64")),
        no_panic::no_panic
    )]
    pub fn new_in<'b>(
        instance: &'b mut MaybeUninit<Self>,
        leaves: &'a [[u8; OUT_LEN]; N],
    ) -> &'b mut Self {
        // TODO: Technically redundant with the where-clause, which is not enforced in generic code
        //  due to https://github.com/rust-lang/rust/issues/164019, remove once fixed
        const { TREE_SIZE_WITHOUT_LEAVES::<N> };

        let instance_ptr = instance.as_mut_ptr();
        // SAFETY: Valid and correctly aligned non-null pointer
        unsafe {
            (&raw mut (*instance_ptr).leaves).write(leaves);
        }
        let tree = {
            // SAFETY: Valid and correctly aligned non-null pointer
            let tree_ptr = unsafe { &raw mut (*instance_ptr).tree };
            // SAFETY: Allocated and correctly aligned uninitialized data
            unsafe {
                tree_ptr
                    .cast::<[MaybeUninit<[u8; OUT_LEN]>; TREE_SIZE_WITHOUT_LEAVES::<N>]>()
                    .as_mut_unchecked()
            }
        };

        Self::init_internal(leaves, tree);

        // SAFETY: Initialized field by field above
        unsafe { instance.assume_init_mut() }
    }

    /// Like [`Self::new()`], but creates heap-allocated instance, avoiding excessive stack usage
    /// for large trees
    #[cfg(feature = "alloc")]
    pub fn new_boxed(leaves: &'a [[u8; OUT_LEN]; N]) -> Box<Self> {
        let mut instance = Box::<Self>::new_uninit();

        Self::new_in(&mut instance, leaves);

        // SAFETY: Initialized by constructor above
        unsafe { instance.assume_init() }
    }

    // TODO: Unlock on RISC-V, it started failing since https://github.com/nazar-pc/abundance/pull/551
    //  for unknown reason
    #[cfg_attr(
        all(feature = "no-panic", not(target_arch = "riscv64")),
        no_panic::no_panic
    )]
    fn init_internal(
        leaves: &[[u8; OUT_LEN]; N],
        tree: &mut [MaybeUninit<[u8; OUT_LEN]>; TREE_SIZE_WITHOUT_LEAVES::<N>],
    ) {
        let mut tree_hashes = tree.as_mut_slice();
        let mut level_hashes = leaves.as_slice();

        while level_hashes.len() > 1 {
            let num_pairs = level_hashes.len() / 2;
            let parent_hashes;
            // SAFETY: The size of the tree is statically known to match the number of leaves and
            // levels of hashes
            (parent_hashes, tree_hashes) = unsafe { tree_hashes.split_at_mut_unchecked(num_pairs) };

            let (parent_hashes_chunks, parent_hashes_remainder) =
                parent_hashes.as_chunks_mut::<BATCH_HASH_NUM_BLOCKS>();
            if parent_hashes_remainder.is_empty() {
                for (pairs, hashes) in level_hashes
                    .as_chunks::<BATCH_HASH_NUM_LEAVES>()
                    .0
                    .iter()
                    .zip(parent_hashes_chunks)
                {
                    // TODO: This memory copy is unfortunate, make hashing write into this memory
                    //  directly once blake3 API improves
                    *hashes = MaybeUninit::new(hash_pairs(pairs)).into();
                }
            } else {
                for (pair, parent_hash) in level_hashes
                    .as_chunks::<{ BLOCK_LEN / OUT_LEN }>()
                    .0
                    .iter()
                    .zip(parent_hashes.iter_mut())
                {
                    parent_hash.write(hash_pair_block(pair.flatten_ref()));
                }
            }

            // SAFETY: Just initialized
            level_hashes = unsafe { parent_hashes.assume_init_ref() };
        }
    }

    // TODO: Method that generates not only root, but also proof, like Unbalanced Merkle Tree
    /// Compute Merkle Tree root.
    ///
    /// This is functionally equivalent to creating an instance first and calling [`Self::root()`]
    /// method, but is faster and avoids heap allocation when root is the only thing that is needed.
    #[inline]
    #[cfg_attr(feature = "no-panic", no_panic::no_panic)]
    pub fn compute_root_only(leaves: &[[u8; OUT_LEN]; N]) -> [u8; OUT_LEN] {
        // TODO: Technically redundant with the where-clause, which is not enforced in generic code
        //  due to https://github.com/rust-lang/rust/issues/164019, remove once fixed
        const { TREE_SIZE_WITHOUT_LEAVES::<N> };

        // Special case for small trees below optimal SIMD width
        match N {
            2 => {
                let [root] = hash_pairs(leaves);

                return root;
            }
            4 => {
                let hashes = hash_pairs::<2, _>(leaves);
                let [root] = hash_pairs(&hashes);

                return root;
            }
            8 => {
                let hashes = hash_pairs::<4, _>(leaves);
                let hashes = hash_pairs::<2, _>(&hashes);
                let [root] = hash_pairs(&hashes);

                return root;
            }
            16 => {
                let hashes = hash_pairs::<8, _>(leaves);
                let hashes = hash_pairs::<4, _>(&hashes);
                let hashes = hash_pairs::<2, _>(&hashes);
                let [root] = hash_pairs(&hashes);

                return root;
            }
            _ => {
                // We know this is the case
                assert!(N >= BATCH_HASH_NUM_LEAVES);
            }
        }

        // Stack of intermediate nodes per tree level. The logic here is the same as with a small
        // tree above, except we store `BATCH_HASH_NUM_BLOCKS` hashes per level and do a
        // post-processing step at the very end to collapse them into a single root hash.
        let mut stack = [[[0u8; OUT_LEN]; BATCH_HASH_NUM_BLOCKS]; ROOT_ONLY_LARGE_STACK_SIZE::<N>];

        // This variable allows reusing and reducing stack usage instead of having a separate
        // `current` variable
        let mut parent_current = [[0u8; OUT_LEN]; BATCH_HASH_NUM_LEAVES];
        for (num_chunks, chunk_leaves) in leaves
            .as_chunks::<BATCH_HASH_NUM_LEAVES>()
            .0
            .iter()
            .enumerate()
        {
            let current_half = &mut parent_current[BATCH_HASH_NUM_BLOCKS..];

            let current = hash_pairs::<BATCH_HASH_NUM_BLOCKS, _>(chunk_leaves);
            current_half.copy_from_slice(&current);

            // Every bit set to `1` corresponds to an active Merkle Tree level
            let lowest_active_levels = num_chunks.trailing_ones() as usize;
            for parent in &mut stack[..lowest_active_levels] {
                let parent_half = &mut parent_current[..BATCH_HASH_NUM_BLOCKS];
                parent_half.copy_from_slice(parent);

                let current = hash_pairs::<BATCH_HASH_NUM_BLOCKS, _>(&parent_current);

                let current_half = &mut parent_current[BATCH_HASH_NUM_BLOCKS..];
                current_half.copy_from_slice(&current);
            }

            let current_half = &mut parent_current[BATCH_HASH_NUM_BLOCKS..];

            // Place freshly computed 8 hashes into the first inactive level
            stack[lowest_active_levels].copy_from_slice(current_half);
        }

        let hashes = &stack[ROOT_ONLY_LARGE_STACK_SIZE::<N> - 1];
        let hashes = hash_pairs::<{ BATCH_HASH_NUM_BLOCKS / 2 }, _>(hashes);
        let hashes = hash_pairs::<{ BATCH_HASH_NUM_BLOCKS / 4 }, _>(&hashes);
        let hashes = hash_pairs::<{ BATCH_HASH_NUM_BLOCKS / 8 }, _>(&hashes);
        let [root] = hash_pairs::<{ BATCH_HASH_NUM_BLOCKS / 16 }, _>(&hashes);

        root
    }

    /// Get the root of Merkle Tree
    #[inline]
    #[cfg_attr(feature = "no-panic", no_panic::no_panic)]
    pub fn root(&self) -> [u8; OUT_LEN] {
        *self
            .tree
            .last()
            .or(self.leaves.last())
            .expect("There is always at least one leaf hash; qed")
    }

    /// Iterator over proofs in the same order as provided leaf hashes
    #[cfg_attr(feature = "no-panic", no_panic::no_panic)]
    pub fn all_proofs(&self) -> ProofsIterator<'_, N> {
        ProofsIterator {
            leaves: self.leaves,
            tree: &self.tree,
            leaf_index: 0,
        }
    }

    /// Verify previously generated proof
    #[inline]
    #[cfg_attr(feature = "no-panic", no_panic::no_panic)]
    pub fn verify(
        root: &[u8; OUT_LEN],
        proof: &[[u8; OUT_LEN]; PROOF_ELEMENTS::<N>],
        leaf_index: usize,
        leaf: [u8; OUT_LEN],
    ) -> bool {
        // TODO: Technically redundant with the where-clause, which is not enforced in generic code
        //  due to https://github.com/rust-lang/rust/issues/164019, remove once fixed
        const { TREE_SIZE_WITHOUT_LEAVES::<N> };

        if leaf_index >= N {
            return false;
        }

        let mut computed_root = leaf;

        let mut position = leaf_index;
        for hash in proof {
            computed_root = if position.is_multiple_of(2) {
                hash_pair(&computed_root, hash)
            } else {
                hash_pair(hash, &computed_root)
            };

            position /= 2;
        }

        root == &computed_root
    }
}

/// Iterator over proofs for a balanced Merkle tree
#[derive(Debug)]
pub struct ProofsIterator<'a, const N: usize> {
    leaves: &'a [[u8; OUT_LEN]; N],
    tree: &'a [[u8; OUT_LEN]; TREE_SIZE_WITHOUT_LEAVES::<N>],
    /// Index of the next leaf, the iterator is exhausted once it reaches `N`
    leaf_index: usize,
}

impl<const N: usize> ProofsIterator<'_, N> {
    #[inline(always)]
    fn remaining(&self) -> usize {
        N - self.leaf_index
    }
}

impl<const N: usize> Iterator for ProofsIterator<'_, N> {
    type Item = [[u8; OUT_LEN]; PROOF_ELEMENTS::<N>];

    #[cfg_attr(feature = "no-panic", no_panic::no_panic)]
    fn next(&mut self) -> Option<Self::Item> {
        let index = self.leaf_index;

        // The line below is a more efficient branchless version of this:
        // let sibling_index = if index % 2 == 0 {
        //     index + 1
        // } else {
        //     index - 1
        // };
        let sibling_index = index ^ 1;
        // `N` is an even number, so the sibling is missing exactly when `index >= N`, meaning the
        // iterator is exhausted
        let sibling_hash = *self.leaves.get(sibling_index)?;
        self.leaf_index += 1;

        let mut proof = [MaybeUninit::<[u8; OUT_LEN]>::uninit(); _];
        proof[0].write(sibling_hash);

        // Part that is shared between left and right leaf proofs
        let shared_proof = &mut proof[1..];

        let mut tree_hashes = self.tree.as_slice();
        let mut parent_position = index / 2;
        let mut parent_level_size = N / 2;

        for hash in shared_proof {
            let parent_other_position = parent_position ^ 1;

            let other_hash = tree_hashes.get(parent_other_position).expect(
                "`parent_position` is within the current level of the tree, whose size is an \
                    even number; qed",
            );
            hash.write(*other_hash);
            tree_hashes = &tree_hashes[parent_level_size..];

            parent_position /= 2;
            parent_level_size /= 2;
        }

        // SAFETY: Just initialized
        Some(unsafe { proof.transpose().assume_init() })
    }

    #[inline(always)]
    fn size_hint(&self) -> (usize, Option<usize>) {
        let remaining = self.remaining();
        (remaining, Some(remaining))
    }

    #[inline(always)]
    fn count(self) -> usize {
        self.remaining()
    }

    #[cfg_attr(feature = "no-panic", no_panic::no_panic)]
    fn last(mut self) -> Option<Self::Item> {
        if self.leaf_index >= N {
            return None;
        }
        self.leaf_index = N - 1;
        self.next()
    }

    #[inline(always)]
    fn advance_by(&mut self, n: usize) -> Result<(), NonZero<usize>> {
        let advance = n.min(self.remaining());
        self.leaf_index += advance;
        NonZero::new(n - advance).map_or(Ok(()), Err)
    }

    #[inline(always)]
    fn nth(&mut self, n: usize) -> Option<Self::Item> {
        match self.advance_by(n) {
            Ok(()) => self.next(),
            Err(_) => None,
        }
    }
}

impl<const N: usize> ExactSizeIterator for ProofsIterator<'_, N> {
    #[inline(always)]
    fn len(&self) -> usize {
        self.remaining()
    }
}

// SAFETY: size_hint is always exact
unsafe impl<const N: usize> TrustedLen for ProofsIterator<'_, N> {}