This function should be refactored to behave like fwdpy11's. We need a lookup based on position for selected mutations, and only insert them into the fixations if the position doesn't already exist.
This issue does not affect simulation correctness. Rather, what happens is that fixation times of selected variants are lost when using this function.
(This does not affect fwdpy11, which uses custom functions for this.)
This function should be refactored to behave like fwdpy11's. We need a lookup based on position for selected mutations, and only insert them into the fixations if the position doesn't already exist.
This issue does not affect simulation correctness. Rather, what happens is that fixation times of selected variants are lost when using this function.
(This does not affect fwdpy11, which uses custom functions for this.)