if a sampling size larger than the number of paths is provided, the w…

…hole panel is used for genotyping. Paths are no longer duplicated to match the sampling size
eblerjana · Jan 31, 2024 · 76478a5 · 76478a5
1 parent 7c171f3
commit 76478a5
Showing 1 changed file with 6 additions and 0 deletions.
diff --git a/src/commands.cpp b/src/commands.cpp
@@ -497,6 +497,9 @@ int run_single_command(string precomputed_prefix, string readfile, string reffil
 			} else {
 				sampling_size = nr_paths;		
 			}
+		} else if (sampling_size > nr_paths) {
+			// make sure that sampling size does not exceed nr_paths in panel
+			sampling_size = nr_paths;
 		}
 
 		PathSampler path_sampler(nr_paths);
@@ -937,6 +940,9 @@ int run_genotype_command(string precomputed_prefix, string readfile, string outn
 				} else {
 					sampling_size = nr_paths;		
 				}
+			} else if (sampling_size > nr_paths) {
+				// make sure that sampling size does not exceed nr_paths in panel
+				sampling_size = nr_paths;
 			}
 
 			PathSampler path_sampler(nr_paths);