check that constraint tree has at least 4 taxa and print a clear erro…

…r message otherwise (#200, #202)
amkozlov · Nov 21, 2024 · ecdf971 · ecdf971
1 parent e12b46b
commit ecdf971
Show file tree

Hide file tree

Showing 3 changed files with 5 additions and 2 deletions.
diff --git a/libs/coraxlib b/libs/coraxlib
diff --git a/src/main.cpp b/src/main.cpp
@@ -1471,6 +1471,9 @@ void load_constraint(RaxmlInstance& instance)
         (cons_tree.num_tips() == parted_msa.taxon_count() ? "" : "non-") <<
         "comprehensive constraint tree with " << cons_tree.num_tips() << " taxa" << endl << endl;
 
+    if (cons_tree.num_tips() < 4)
+      throw runtime_error("Constraint tree must have at least 4 taxa!");
+
     // check if taxa names are consistent between constraint tree and MSA
     {
       NameList missing_taxa;

diff --git a/src/version.h b/src/version.h
@@ -1,3 +1,3 @@
 #define RAXML_VERSION "2.0.0-dev"
-#define RAXML_DATE    "30.10.2024"
+#define RAXML_DATE    "21.11.2024"
 #define RAXML_INTVER  200