Fix copyTree for a taxa_set to work with rooted tree properly (report…

…ed by Suha Naser)
Cibiv · Mar 27, 2020 · db6696e · db6696e
1 parent 714d6cb
commit db6696e
Show file tree

Hide file tree

Showing 2 changed files with 7 additions and 0 deletions.
diff --git a/tree/constrainttree.cpp b/tree/constrainttree.cpp
@@ -103,6 +103,12 @@ void ConstraintTree::initFromTree() {
 
 void ConstraintTree::readConstraint(MTree &src_tree) {
     copyTree(&src_tree);
+    // convert all rooted constraint tree to unrooted
+    if (rooted) {
+        if (verbose_mode >= VB_MED)
+            cout << "Converting rooted constraint tree to unrooted" << endl;
+        convertToUnrooted();
+    }
     initFromTree();
 }
 

diff --git a/tree/mtree.cpp b/tree/mtree.cpp
@@ -142,6 +142,7 @@ void MTree::copyTree(MTree *tree) {
 }
 
 void MTree::copyTree(MTree *tree, string &taxa_set) {
+    rooted = tree->rooted;
     if (rooted) {
         ASSERT(tree->rooted);
         taxa_set.push_back(1);