diff --git a/benches/smt-mutations.rs b/benches/smt-mutations.rs
index d922ca05..77db4cce 100644
--- a/benches/smt-mutations.rs
+++ b/benches/smt-mutations.rs
@@ -7,6 +7,37 @@ use rand::{prelude::SliceRandom, rngs::StdRng, Rng, SeedableRng};
 
 const DEPTH: u8 = 64;
 
+fn benchmark_compute_mutations(c: &mut Criterion) {
+    let mut group = c.benchmark_group("compute_mutations");
+    group.sample_size(10);
+    group.measurement_time(std::time::Duration::from_secs(300));
+
+    // Fixed seed for reproducibility
+    let rng_seed = 42;
+    let mut rng = StdRng::seed_from_u64(rng_seed);
+
+    // Benchmark for various mutation set sizes
+    for &mutation_count in &[10_000] {
+        group.bench_with_input(
+            BenchmarkId::new("SimpleSmt: compute_mutations", mutation_count),
+            &mutation_count,
+            |b, &mutation_count| {
+                // Batch-based benchmarking
+                b.iter_batched(
+                    || generate_tree_and_mutations(&mut rng, mutation_count),
+                    |(smt, mutations)| {
+                        // Compute mutations in the benchmark to measure execution time
+                        let _mutations = smt.compute_mutations(mutations);
+                    },
+                    BatchSize::SmallInput,
+                );
+            },
+        );
+    }
+
+    group.finish();
+}
+
 fn benchmark_apply_mutations(c: &mut Criterion) {
     let mut group = c.benchmark_group("apply_mutations");
     group.sample_size(10);
@@ -17,7 +48,7 @@ fn benchmark_apply_mutations(c: &mut Criterion) {
     let mut rng = StdRng::seed_from_u64(rng_seed);
 
     // Benchmark for various mutation set sizes
-    for &mutation_count in &[1_000, 100_000] {
+    for &mutation_count in &[10_000] {
         group.bench_with_input(
             BenchmarkId::new("SimpleSmt: apply_mutations", mutation_count),
             &mutation_count,
@@ -25,43 +56,11 @@ fn benchmark_apply_mutations(c: &mut Criterion) {
                 // Batch-based benchmarking
                 b.iter_batched(
                     || {
-                        const REMOVAL_PROBABILITY: f64 = 0.2;
-
-                        // Fill tree with 10x more initial elements
-                        let initial_fill_count = mutation_count * 10;
-
-                        // Initialize the tree with initial random values
-                        let initial_kv_pairs = generate_kv_pairs(&mut rng, initial_fill_count);
-                        let smt = SimpleSmt::<DEPTH>::with_leaves(initial_kv_pairs.iter().cloned())
-                            .unwrap();
-
-                        // Select and change a half of pairs from the filled tree (values to be
-                        // updated or removed with given probability)
-                        let mut mutation_kv_pairs: Vec<_> = initial_kv_pairs
-                            .choose_multiple(&mut rng, mutation_count / 2)
-                            .cloned()
-                            .map(|(key, _value)| {
-                                let value = if rng.gen_bool(REMOVAL_PROBABILITY) {
-                                    EMPTY_WORD
-                                } else {
-                                    generate_word(&mut rng)
-                                };
-
-                                (key, value)
-                            })
-                            .collect();
-
-                        // Append another half of new values (values to be added)
-                        for _ in 0..mutation_count / 2 {
-                            mutation_kv_pairs.push((rng.gen(), generate_word(&mut rng)));
-                        }
+                        let (smt, mutation_kv_pairs) =
+                            generate_tree_and_mutations(&mut rng, mutation_count);
 
                         // Compute mutations
-                        let mutations = smt.compute_mutations(
-                            mutation_kv_pairs
-                                .into_iter()
-                                .map(|(key, value)| (LeafIndex::new(key).unwrap(), value)),
-                        );
+                        let mutations = smt.compute_mutations(mutation_kv_pairs);
 
                         (smt, mutations)
                     },
@@ -78,22 +77,60 @@ fn benchmark_apply_mutations(c: &mut Criterion) {
     group.finish();
 }
 
-criterion_group!(benches, benchmark_apply_mutations);
+criterion_group!(benches, benchmark_compute_mutations, benchmark_apply_mutations);
 criterion_main!(benches);
 
 // HELPER FUNCTIONS
 // =================================================================================================
 
-/// Helper function to generate random `(u64, Word)` key-value pairs
-fn generate_kv_pairs(rng: &mut StdRng, count: usize) -> Vec<(u64, Word)> {
-    (0..count)
+/// Helper function to generate initial `SimpleSmt` tree with random keys and values and mutation
+/// key-value pairs
+fn generate_tree_and_mutations(
+    rng: &mut StdRng,
+    mutation_count: usize,
+) -> (SimpleSmt<DEPTH>, Vec<(LeafIndex<DEPTH>, Word)>) {
+    const INITIAL_FILL_COUNT: usize = 1_000_000;
+    const REMOVAL_PROBABILITY: f64 = 0.2;
+
+    // Initialize the tree with initial random values
+    let initial_kv_pairs: Vec<_> = (0..INITIAL_FILL_COUNT)
         .map(|_| {
-            let key = rng.gen();
+            let key: u64 = rng.gen();
             let value = generate_word(rng);
 
             (key, value)
         })
-        .collect()
+        .collect();
+
+    // Select and change a half of pairs from the filled tree (values to be
+    // updated or removed with given probability)
+    let mut mutation_kv_pairs: Vec<_> = initial_kv_pairs
+        .choose_multiple(rng, mutation_count / 2)
+        .cloned()
+        .map(|(key, _value)| {
+            let value = if rng.gen_bool(REMOVAL_PROBABILITY) {
+                EMPTY_WORD
+            } else {
+                generate_word(rng)
+            };
+
+            (key, value)
+        })
+        .collect();
+
+    // Append another half of new values (values to be added)
+    for _ in 0..mutation_count / 2 {
+        mutation_kv_pairs.push((rng.gen(), generate_word(rng)));
+    }
+
+    let smt = SimpleSmt::<DEPTH>::with_leaves(initial_kv_pairs).unwrap();
+
+    let mutations: Vec<_> = mutation_kv_pairs
+        .into_iter()
+        .map(|(key, value)| (LeafIndex::new(key).unwrap(), value))
+        .collect();
+
+    (smt, mutations)
 }
 
 /// Helper function to generate random `Word`