Add missing C test

Author: benjeffery

c/tests/test_trees.c

@@ -8776,6 +8776,51 @@ test_init_compute_mutation_parents(void)
     tsk_treeseq_free(&ts);
 }
 
+static void
+test_init_compute_mutation_parents_errors(void)
+{
+    int ret;
+    tsk_id_t row_ret;
+    tsk_table_collection_t tables;
+    tsk_treeseq_t ts;
+    const char *sites = "0.5       0\n"
+                        "0         0\n";
+
+    ret = tsk_table_collection_init(&tables, 0);
+    CU_ASSERT_EQUAL_FATAL(ret, 0);
+
+    tables.sequence_length = 1;
+    parse_nodes(single_tree_ex_nodes, &tables.nodes);
+    CU_ASSERT_EQUAL_FATAL(tables.nodes.num_rows, 7);
+    parse_edges(single_tree_ex_edges, &tables.edges);
+    CU_ASSERT_EQUAL_FATAL(tables.edges.num_rows, 6);
+    parse_sites(sites, &tables.sites);
+    CU_ASSERT_EQUAL_FATAL(tables.sites.num_rows, 2);
+    tables.sequence_length = 1.0;
+
+    ret = tsk_treeseq_init(
+        &ts, &tables, TSK_TS_INIT_BUILD_INDEXES | TSK_TS_INIT_COMPUTE_MUTATION_PARENTS);
+    CU_ASSERT_EQUAL_FATAL(ret, TSK_ERR_UNSORTED_SITES);
+    tsk_treeseq_free(&ts);
+
+    tsk_site_table_clear(&tables.sites);
+    row_ret = tsk_site_table_add_row(&tables.sites, 0.5, "A", 1, NULL, 0);
+    CU_ASSERT_EQUAL_FATAL(row_ret, 0);
+    row_ret = tsk_mutation_table_add_row(
+        &tables.mutations, 0, 0, TSK_NULL, TSK_UNKNOWN_TIME, "A", 1, NULL, 0);
+    CU_ASSERT_EQUAL_FATAL(row_ret, 0);
+    row_ret = tsk_mutation_table_add_row(
+        &tables.mutations, 0, 4, TSK_NULL, TSK_UNKNOWN_TIME, "A", 1, NULL, 0);
+    CU_ASSERT_EQUAL_FATAL(row_ret, 1);
+
+    ret = tsk_treeseq_init(
+        &ts, &tables, TSK_TS_INIT_BUILD_INDEXES | TSK_TS_INIT_COMPUTE_MUTATION_PARENTS);
+    CU_ASSERT_EQUAL_FATAL(ret, TSK_ERR_MUTATION_PARENT_AFTER_CHILD);
+    tsk_treeseq_free(&ts);
+
+    tsk_table_collection_free(&tables);
+}
+
 int
 main(int argc, char **argv)
 {
@@ -8985,6 +9030,8 @@ main(int argc, char **argv)
         { "test_init_take_ownership_no_edge_metadata",
             test_init_take_ownership_no_edge_metadata },
         { "test_init_compute_mutation_parents", test_init_compute_mutation_parents },
+        { "test_init_compute_mutation_parents_errors",
+            test_init_compute_mutation_parents_errors },
         { NULL, NULL },
     };