diff --git a/src/rsrs/rsrs_cycle.rs b/src/rsrs/rsrs_cycle.rs
index e0be718..37fee5e 100644
--- a/src/rsrs/rsrs_cycle.rs
+++ b/src/rsrs/rsrs_cycle.rs
@@ -43,7 +43,6 @@ use rustc_hash::FxHashSet;
 use std::{
     collections::HashMap,
     path::{Path, PathBuf},
-    sync::Arc,
     time::{Duration, Instant},
 };
 
@@ -65,7 +64,6 @@ pub struct Rsrs<Item: RlstScalar> {
     pub stats: Stats,
     options: RsrsOptions<Item>,
     anticipated_fixed_rank_samples: Option<FixedRankSampleBudget>,
-    thread_pool: Arc<rayon::ThreadPool>,
 }
 
 #[derive(Debug, Clone, PartialEq, Eq)]
@@ -839,7 +837,6 @@ where
         };
 
         println!("Configured number of threads = {}", options.num_threads);
-        let thread_pool = Arc::new(build_rsrs_thread_pool(options.num_threads));
         let stats = Stats {
             sampling_time: Vec::new(),
             sample_loading_time: 0_u128,
@@ -886,7 +883,6 @@ where
             active_samples: 0,
             options,
             anticipated_fixed_rank_samples,
-            thread_pool,
         }
     }
 
@@ -1415,14 +1411,14 @@ where
         level: usize,
         update_type: &UpdateType<Item>,
     ) -> (u128, u128, Option<u128>) {
-        let thread_pool = &self.thread_pool;
+        let thread_pool = build_rsrs_thread_pool(self.options.num_threads);
         let (mut tot_id_update, mut tot_lu_update) = self.y_data.update_samples(
             update_start,
             samples_to_update,
             level,
             update_type,
             &self.options.fact_type,
-            thread_pool,
+            &thread_pool,
             self.options.num_threads,
         );
 
@@ -1439,7 +1435,7 @@ where
                 level,
                 update_type,
                 &self.options.fact_type,
-                thread_pool,
+                &thread_pool,
                 self.options.num_threads,
             );
             tot_id_update += tot_z_id_update;
@@ -1528,7 +1524,7 @@ where
         let current_level_keys: Vec<MortonKey> =
             self.level_indexing.level_keys.iter().copied().collect();
         // Process all batches *sequentially*, each batch using Rayon internally
-        let thread_pool = Arc::clone(&self.thread_pool);
+        let thread_pool = build_rsrs_thread_pool(self.options.num_threads);
         let batches_res: Vec<_> = thread_pool.install(|| {
             independent_near_fields
                 .into_iter()
@@ -1693,6 +1689,7 @@ where
                                     self.options.sketching.fixed_rank_sampling_mode,
                                     &self.options,
                                 );
+
                                 <Item as Skel<Item, Space>>::lu_step(
                                     &skel_box,
                                     scratch,
@@ -1814,7 +1811,7 @@ where
             )
         });
         let start = Instant::now();
-        let thread_pool = Arc::clone(&self.thread_pool);
+        let thread_pool = build_rsrs_thread_pool(self.options.num_threads);
         let id_level_iteration_res: Vec<_> = thread_pool.install(|| {
             current_box_indices
                 .par_iter()
@@ -1964,7 +1961,7 @@ where
         let lu_step_start: Instant = Instant::now();
 
         let mut inactive_inds = Vec::new();
-        let thread_pool = Arc::clone(&self.thread_pool);
+        let thread_pool = build_rsrs_thread_pool(self.options.num_threads);
         let batches_res: Vec<_> = independent_near_fields
             .into_iter()
             .map(|batch| {