@@ -301,10 +301,7 @@ pub mod simulation {
301301
302302 // state variables for site/mutation tables
303303 let num_sites = f64:: from ( seqlen. into ( ) ) as usize ;
304- let mut site_last_mutation_order = vec ! [ 0usize ; num_sites] ;
305-
306- let mut site_last_mutation_tables = vec ! [ MutationId :: NULL ; num_sites] ;
307- let mut site_last_mutation_tr_tbls = vec ! [ MutationId :: NULL ; num_sites] ;
304+ let site_last_mutation_order = vec ! [ 0usize ; num_sites] ;
308305
309306 let mut site_id_map_tables = vec ! [ SiteId :: NULL ; num_sites] ;
310307 let mut site_id_map_tr_tbls = vec ! [ SiteId :: NULL ; num_sites] ;
@@ -387,14 +384,10 @@ pub mod simulation {
387384 tables. add_site ( mut_pos as f64 , Some ( b"a" ) ) . unwrap ( ) ;
388385 }
389386 // add mutation
390- let parent_mut = site_last_mutation_tables[ mut_pos] ;
391387 let site = site_id_map_tables[ mut_pos] ;
392- let new_mutation = tables
388+ let _ = tables
393389 . add_mutation ( site, c, MutationId :: NULL , t, Some ( derived_state) )
394390 . unwrap ( ) ;
395-
396- site_last_mutation_tables[ mut_pos] = new_mutation;
397- site_last_mutation_order[ mut_pos] += 1 ;
398391 }
399392
400393 find_overlaps ( s, e, & intervals, & mut buffer) ;
@@ -410,12 +403,10 @@ pub mod simulation {
410403 tr_tbls. add_site ( mut_pos as f64 , Some ( b"a" ) ) . unwrap ( ) ;
411404 }
412405 // add mutation
413- let parent_mut = site_last_mutation_tr_tbls[ mut_pos] ;
414406 let site = site_id_map_tr_tbls[ mut_pos] ;
415- let new_mutation = tr_tbls
407+ let _ = tr_tbls
416408 . add_mutation ( site, c, MutationId :: NULL , t, Some ( derived_state) )
417409 . unwrap ( ) ;
418- site_last_mutation_tr_tbls[ mut_pos] = new_mutation;
419410 }
420411 }
421412 }
@@ -447,7 +438,7 @@ pub mod simulation {
447438 // build indices
448439 tables. build_index ( ) . unwrap ( ) ;
449440 tr_tbls. build_index ( ) . unwrap ( ) ;
450-
441+
451442 // Calculate mutation parents
452443 // TODO: need safe API for this.
453444 let code = unsafe {
0 commit comments