fixed species roots atrandomly resolvd multifurcations, but only fixe…

…d for --ml so far
Pas-Kapli · Jul 4, 2020 · b831c9c · b831c9c
1 parent 3fd12fb
commit b831c9c
Showing 1 changed file with 4 additions and 0 deletions.
diff --git a/src/dp.c b/src/dp.c
@@ -79,6 +79,10 @@ static void dp_recurse(rtree_t * node, long method)
       /* if at least one of the two entries is not valid/filled, skip */
       if (!v_vec[j].filled || !w_vec[k].filled) continue;
 
+      /* Fix species roots at randomlyresolved multifurcations */
+      if (node->left->length == 0 && node->left->left && v_vec[j].species_count == 1) continue;
+      if (node->right->length == 0 && node->right->left && w_vec[k].species_count == 1) continue;
+
       int i = j + k + u_edge_count;
 
       /* set the number of species */