Skip to content

Commit e8686a5

Browse files
committed
Add additional path scaffolding iterations
1 parent c8a1356 commit e8686a5

File tree

1 file changed

+10
-1
lines changed
  • assembler/src/common/modules/path_extend/pipeline

1 file changed

+10
-1
lines changed

assembler/src/common/modules/path_extend/pipeline/launcher.cpp

Lines changed: 10 additions & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -850,7 +850,16 @@ void PathExtendLauncher::Launch() {
850850
const double cluster_length_percentile = params_.pe_cfg.read_cloud.scaff_con.cluster_length_percentile;
851851
size_t length_upper_bound = length_bound_estimator.EstimateUpperBound(cluster_statistics_extractor,
852852
cluster_length_percentile);
853-
ScaffoldPaths(polished_paths);
853+
size_t max_scaffolding_iterations = 5;
854+
size_t prev_path_number = polished_paths.size();
855+
size_t curr_path_number = 0;
856+
size_t curr_scaffolding_iter = 1;
857+
while (prev_path_number != curr_path_number and curr_scaffolding_iter < max_scaffolding_iterations) {
858+
INFO("Starting scaffolding iteration " << curr_scaffolding_iter << " for " << prev_path_number << " paths");
859+
ScaffoldPaths(polished_paths);
860+
curr_path_number = polished_paths.size();
861+
++curr_scaffolding_iter;
862+
}
854863

855864
PolishPaths(polished_paths, contig_paths, polished_map);
856865
DebugOutputPaths(contig_paths, "merged_polished_paths");

0 commit comments

Comments
 (0)