pub fn calculate_scan_occurrences_gaussian_par( times: &[f64], means: Vec<f64>, sigmas: Vec<f64>, target_p: f64, step_size: f64, n_lower_start: f64, n_upper_start: f64, num_threads: usize, ) -> Vec<Vec<i32>>