-
Notifications
You must be signed in to change notification settings - Fork 89
add keep_unary for pedigrees and dtwf #2145
New issue
Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.
By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.
Already on GitHub? Sign in to your account
add keep_unary for pedigrees and dtwf #2145
Conversation
|
Nice! @sgravel I think this should more-or-less do what we want with the pedigree sims, so we could start thinking about trying it out and seeing what we get? |
|
Perfect! Presumably this does not record the ancestors above the MRCA yet, right? |
Yes, that's correct @sgravel. This would require an additional flag indicating the stopping point of the simulation. We have been discussing this exact issue though while thinking about potential performance improvements (see #2121). |
|
Ok, we'll do some testing! |
|
Let's update this and get it merged I think @GertjanBisschop - I guess we should try to register some issues for any remaining TODOs and corner cases so we don't forget them. We can update to the new API later. |
8c464f3 to
3490343
Compare
Codecov Report
@@ Coverage Diff @@
## main #2145 +/- ##
==========================================
- Coverage 91.45% 91.38% -0.08%
==========================================
Files 20 20
Lines 11170 11199 +29
Branches 2268 2280 +12
==========================================
+ Hits 10216 10234 +18
- Misses 521 529 +8
- Partials 433 436 +3
Flags with carried forward coverage won't be shown. Click here to find out more.
Continue to review full report at Codecov.
|
|
Not sure how I might have reduced code coverage that significantly. What do I do to fix this? |
|
Right I guess this means I'll have to write some additional tests in C. |
b2d17ad to
a000f35
Compare
|
Let's have a look through this tomorrow @GertjanBisschop - it might be simpler to push through the final API rather than getting the C test coverage up. |
|
Implemented in #2176. |
This pull request stores all nodes through which some ancestral material passes for pedigrees and dtwf simulations when the flag
). Is this as simple as simply using the next available number in line (
record_unary=Trueis set. This is a first pass to allow for testing (see #2132).This pull request requires defining a
MSP_NODE_IS_PASS_THROUGH_EVENT(see#define MSP_NODE_IS_PT_EVENT (1u << 22))?As with #2130 unit tests are still basic and require further work.