tests pass

Author: molpopgen

src/test_fixtures.rs

@@ -390,7 +390,7 @@ pub mod simulation {
                         let parent_mut = site_last_mutation_tables[mut_pos];
                         let site = site_id_map_tables[mut_pos];
                         let new_mutation = tables
-                            .add_mutation(site, c, parent_mut, t, Some(derived_state))
+                            .add_mutation(site, c, MutationId::NULL, t, Some(derived_state))
                             .unwrap();
 
                         site_last_mutation_tables[mut_pos] = new_mutation;
@@ -413,7 +413,7 @@ pub mod simulation {
                             let parent_mut = site_last_mutation_tr_tbls[mut_pos];
                             let site = site_id_map_tr_tbls[mut_pos];
                             let new_mutation = tr_tbls
-                                .add_mutation(site, c, parent_mut, t, Some(derived_state))
+                                .add_mutation(site, c, MutationId::NULL, t, Some(derived_state))
                                 .unwrap();
                             site_last_mutation_tr_tbls[mut_pos] = new_mutation;
                         }
@@ -447,6 +447,23 @@ pub mod simulation {
         // build indices
         tables.build_index().unwrap();
         tr_tbls.build_index().unwrap();
+        
+        // Calculate mutation parents
+        // TODO: need safe API for this.
+        let code = unsafe {
+            crate::sys::bindings::tsk_table_collection_compute_mutation_parents(
+                tables.as_mut_ptr(),
+                0,
+            )
+        };
+        assert_eq!(code, 0);
+        let code = unsafe {
+            crate::sys::bindings::tsk_table_collection_compute_mutation_parents(
+                tr_tbls.as_mut_ptr(),
+                0,
+            )
+        };
+        assert_eq!(code, 0);
 
         // to tree sequences
         let treeseq_opts = TreeSequenceFlags::default();