search

molpopgen / fwdpp

fwdpp is a C++ template library for implementing efficient forward-time population genetic simulations
http://fwdpp.readthedocs.io
GNU General Public License v3.0
27 stars 11 forks source link

Simplification return value #167

Closed molpopgen closed 5 years ago

molpopgen commented 5 years ago

This PR is a WIP to see if the ideas in #163 are helpful.

There is a logic error in the issue as described. The number of indexes of mutations that are "simplified out" is not the same as the length of a recycling queue. The correct procedure is something like:

  1. generate a list of mutation indexes that are preserved after simplification
  2. generate a list of all possible mutation indexes
  3. The recycling queue is the set difference of these two lists.