File: simple.cpp

package info (click to toggle)
iqtree 2.0.7%2Bdfsg-2
  • links: PTS, VCS
  • area: main
  • in suites: forky, sid
  • size: 14,700 kB
  • sloc: cpp: 142,571; ansic: 57,789; sh: 275; python: 242; makefile: 95
file content (92 lines) | stat: -rw-r--r-- 4,229 bytes parent folder | download | duplicates (2)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
#include <catch.hpp>
#include <sstream>
#include <terraces/errors.hpp>
#include <terraces/simple.hpp>

using namespace terraces::simple;

namespace terraces {
namespace tests {

TEST_CASE("simple_errors") {
	// malformed nwk
	CHECK_THROWS_AS(is_on_terrace("(s1,", "1 1\n1 s1"), bad_input_error);
	// mismatching sizes
	CHECK_THROWS_AS(is_on_terrace("(s1,s2)", "1 1\n1 s1"), bad_input_error);
	// malformed matrix
	CHECK_THROWS_AS(is_on_terrace("(s1,s2)", "1 1\na s1"), bad_input_error);
	// multifurcating newick tree
	CHECK_THROWS_AS(is_on_terrace("(s1,s2,s3,s4)", "4 1\n1 s1\n1 s2\n1 s3\n1 s4"),
	                bad_input_error);
	// unknown/empty name
	CHECK_THROWS_AS(is_on_terrace("((s1,s2),(,s5))", "4 2\n1 0 s1\n1 0 s2\n0 1 s3\n0 1 s4"),
	                bad_input_error);
	// no root
	CHECK_THROWS_AS(is_on_terrace("((s1,s2),(s3,s4))", "4 2\n1 0 s1\n1 0 s2\n0 1 s3\n0 1 s4"),
	                no_usable_root_error);
}

TEST_CASE("simple_results") {
	CHECK(!is_on_terrace("(s1, (s2, (s3, s4)))", "4 2\n1 0 s1\n1 1 s2\n1 1 s3\n1 1 s4"));
	CHECK(is_on_terrace("(s3, ((s1, s2), (s4, s5)))",
	                    "5 2\n1 0 s1\n1 0 s2\n1 1 s3\n0 1 s4\n0 1 s5"));
	CHECK(get_terrace_size("((s4, (s3, (s2, (s1, s6)))), s5)",
	                       "6 3\n1 0 0 s1\n1 0 0 s2\n0 0 1 s3\n0 1 1 s4\n1 1 1 s5\n0 1 1 s6") ==
	      35);
	std::stringstream ss;
	print_terrace_compressed("((s4, (s3, (s2, (s1, s6)))), s5)",
	                         "6 3\n1 0 0 s1\n1 0 0 s2\n0 0 1 s3\n0 1 1 s4\n1 1 1 s5\n0 1 1 s6",
	                         ss);
	CHECK(ss.str() == "(s5,(s2,((s3,s6),(s1,s4))|(s4,{s1,s3,s6})|((s4,(s3,s6)),s1))|((s3,s6),{"
	                  "s1,s2,s4})|({s2,s3,s6},(s1,s4))|(s4,{s1,s2,s3,s6})|((s2,s4),{s1,s3,s6})|"
	                  "((s4,(s3,s6)),(s1,s2))|(((s3,s6),(s2,s4))|(s4,{s2,s3,s6})|((s4,(s3,s6)),"
	                  "s2),s1))");
	ss.str("");
	print_terrace("((s4, (s3, (s2, (s1, s6)))), s5)",
	              "6 3\n1 0 0 s1\n1 0 0 s2\n0 0 1 s3\n0 1 1 s4\n1 1 1 s5\n0 1 1 s6", ss);
	CHECK(ss.str() == "(s5,(s2,((s3,s6),(s1,s4))));\n"
	                  "(s5,(s2,(s4,(s1,(s3,s6)))));\n"
	                  "(s5,(s2,(s4,(s3,(s1,s6)))));\n"
	                  "(s5,(s2,(s4,((s1,s3),s6))));\n"
	                  "(s5,(s2,((s4,(s3,s6)),s1)));\n"
	                  "(s5,((s3,s6),(s1,(s2,s4))));\n"
	                  "(s5,((s3,s6),(s2,(s1,s4))));\n"
	                  "(s5,((s3,s6),((s1,s2),s4)));\n"
	                  "(s5,((s2,(s3,s6)),(s1,s4)));\n"
	                  "(s5,((s3,(s2,s6)),(s1,s4)));\n"
	                  "(s5,(((s2,s3),s6),(s1,s4)));\n"
	                  "(s5,(s4,(s1,(s2,(s3,s6)))));\n"
	                  "(s5,(s4,(s1,(s3,(s2,s6)))));\n"
	                  "(s5,(s4,(s1,((s2,s3),s6))));\n"
	                  "(s5,(s4,(s2,(s1,(s3,s6)))));\n"
	                  "(s5,(s4,(s2,(s3,(s1,s6)))));\n"
	                  "(s5,(s4,(s2,((s1,s3),s6))));\n"
	                  "(s5,(s4,((s1,s2),(s3,s6))));\n"
	                  "(s5,(s4,(s3,(s1,(s2,s6)))));\n"
	                  "(s5,(s4,(s3,(s2,(s1,s6)))));\n"
	                  "(s5,(s4,(s3,((s1,s2),s6))));\n"
	                  "(s5,(s4,((s1,s3),(s2,s6))));\n"
	                  "(s5,(s4,((s2,s3),(s1,s6))));\n"
	                  "(s5,(s4,((s1,(s2,s3)),s6)));\n"
	                  "(s5,(s4,((s2,(s1,s3)),s6)));\n"
	                  "(s5,(s4,(((s1,s2),s3),s6)));\n"
	                  "(s5,((s2,s4),(s1,(s3,s6))));\n"
	                  "(s5,((s2,s4),(s3,(s1,s6))));\n"
	                  "(s5,((s2,s4),((s1,s3),s6)));\n"
	                  "(s5,((s4,(s3,s6)),(s1,s2)));\n"
	                  "(s5,(((s3,s6),(s2,s4)),s1));\n"
	                  "(s5,((s4,(s2,(s3,s6))),s1));\n"
	                  "(s5,((s4,(s3,(s2,s6))),s1));\n"
	                  "(s5,((s4,((s2,s3),s6)),s1));\n"
	                  "(s5,(((s4,(s3,s6)),s2),s1));\n");
}

TEST_CASE("simple_results_force") {
	CHECK(get_terrace_size("((s4, (s3, (s2, (s1, s6)))), s5)", "6 5\n0 1 0 0 0 s1\n0 1 0 0 0 "
	                                                           "s2\n0 0 0 0 1 s3\n0 0 1 0 1 "
	                                                           "s4\n0 1 1 0 1 s5\n0 0 1 0 1 s6",
	                       true) == 35);
}

} // namespace tests
} // namespace terraces