use build_global() to fix issue with using inappropriate number of th…

…reads

Cargo.toml

@@ -1,6 +1,6 @@
 [package]
 name = "STRdust"
-version = "0.8.4"
+version = "0.8.5"
 edition = "2021"
 
 

src/call.rs

@@ -29,7 +29,7 @@ pub fn genotype_repeats(args: Cli) {
         // output is sorted by chrom, start and end
         rayon::ThreadPoolBuilder::new()
             .num_threads(args.threads)
-            .build()
+            .build_global()
             .expect("Failed to create threadpool");
         // genotypes contains the output of the genotyping, a struct instance
         let genotypes = Mutex::new(Vec::new());
@@ -59,7 +59,7 @@ pub fn genotype_repeats(args: Cli) {
         // this does not use the BufWriter (as that doesn't work across threads), but writes directly to stdout
         rayon::ThreadPoolBuilder::new()
             .num_threads(args.threads)
-            .build()
+            .build_global()
             .expect("Failed to create threadpool");
         let num_intervals = repeats.num_intervals;
         repeats

src/genotype.rs

@@ -76,6 +76,7 @@ fn genotype_repeat(
     // Create an index for minimap2 alignment to the artificial reference
     let aligner = minimap2::Aligner::builder()
         .map_ont()
+        .with_index_threads(1)
         .with_cigar()
         .with_seq(&repeat_compressed_reference)
         .unwrap_or_else(|err| panic!("Unable to build index:\n{err}"));