search

luntergroup / smcsmc

Demographic inference from whole genomes
GNU General Public License v3.0
11 stars 4 forks source link

assertion failed in scrm #28

Closed shajoezhu closed 9 years ago

shajoezhu commented 9 years ago 
./smcsmc_dbg -nsam 3 -Np 1000 -r 130 1000000 -seed 11

===== BUILDING NEXT GENEALOGY =====
        scrm Segment Nr: 2 | Sequence position: 1612.501
╦ ╦ ╦ -       0 - 1:0x7f6db8e10010(0) 3:0x7f6db8e10160(0) 2:0x7f6db8e10080(0) 
║ ╚╦╝ - 1666.48 - 0x7f6db8e101d0(0) 
╚╦═╝  - 113916.4 - 0x7f6db8e100f0(0) 
        scrm * Recombination at height 77736.96 (above 0x7f6db8e101d0)
        scrm * Cutting subtree below recombination 
        scrm * * New leaf of local tree: 0x7f6db8e10240
        scrm * * New root of subtree: 0x7f6db8e102b0
        scrm * * Done
╦  ╦ ╦  -       0 - 1:0x7f6db8e10010(0) 3:0x7f6db8e10160(0) 2:0x7f6db8e10080(0) 
║  ╚╦╝  - 1666.48 - 0x7f6db8e101d0(0) 
║ ┬ ╚═╦ - 77736.96 - 0x7f6db8e10240(0) 0x7f6db8e102b0(0) 
╚╦┘     - 113916.4 - 0x7f6db8e100f0(0) 
        scrm * Starting coalescence
smcsmc_dbg: src/scrm/src/time_interval.cc:180: void TimeIntervalIterator::searchContemporariesBottomUp(Node*, bool): Assertion `start_node != __null' failed.
shajoezhu commented 9 years ago

this seems to work after scrm commit https://github.com/luntergroup/scrm_jz_stable_branch/commit/e6f7eee1ef9286c053d2c1a389f515f6b472055c