endpoint fix

jeromekelleher · Sep 11, 2023 · 0afc2b5 · 0afc2b5
1 parent 5d60e6f
commit 0afc2b5
Showing 1 changed file with 6 additions and 4 deletions.
diff --git a/c/examples/multichrom_wright_fisher.c b/c/examples/multichrom_wright_fisher.c
@@ -226,14 +226,16 @@ simulate(
                         right_parent, child, NULL, 0);
                     check_tsk_error(ret);
                 } else {
-                    assert(chunk_left < breakpoint && breakpoint < chunk_right);
+                    assert(chunk_left < breakpoint && breakpoint <= chunk_right);
                     /* Breakpoint is within this chunk so need two edges */
                     ret = tsk_edge_table_add_row(&tcs[k].edges, chunk_left, breakpoint,
                         left_parent, child, NULL, 0);
                     check_tsk_error(ret);
-                    ret = tsk_edge_table_add_row(&tcs[k].edges, breakpoint, chunk_right,
-                        right_parent, child, NULL, 0);
-                    check_tsk_error(ret);
+                    if (breakpoint < chunk_right) {
+                        ret = tsk_edge_table_add_row(&tcs[k].edges, breakpoint, chunk_right,
+                            right_parent, child, NULL, 0);
+                        check_tsk_error(ret);
+                    }
                 }
                 chunk_left = chunk_right;
             }