-
Notifications
You must be signed in to change notification settings - Fork 72
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
allow genetic_relatedness_vector
to do windows that don't span the whole genome
#3003
Conversation
c7a784e
to
85e7531
Compare
Codecov ReportAttention: Patch coverage is
Additional details and impacted files@@ Coverage Diff @@
## main #3003 +/- ##
=======================================
Coverage 89.80% 89.81%
=======================================
Files 29 29
Lines 31931 31979 +48
Branches 6179 6190 +11
=======================================
+ Hits 28675 28721 +46
- Misses 1860 1861 +1
- Partials 1396 1397 +1
Flags with carried forward coverage won't be shown. Click here to find out more.
|
c26689e
to
e6bbbb6
Compare
@@ -9977,6 +9984,20 @@ tsk_matvec_calculator_init(tsk_matvec_calculator_t *self, const tsk_treeseq_t *t | |||
tsk_memset(result, 0, num_windows * num_samples * num_weights * sizeof(*result)); | |||
tsk_memset(self->parent, TSK_NULL, num_nodes * sizeof(*self->parent)); | |||
|
|||
ret = tsk_tree_position_init(&self->tree_pos, ts, 0); | |||
if (ret != 0) { | |||
goto out; |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
this function always returns 0 currently
} | ||
ret = tsk_tree_position_seek_forward(&self->tree_pos, index); | ||
if (ret != 0) { | ||
goto out; |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
same here
I think this is ready to go - the only uncovered lines are not cover-able. |
I've also put in the code here to let this work with arbitrary nodes. To see those changes separately look at the commits. |
367d075
to
3e558c0
Compare
This looks good to go! (One test will need to be re-run but it failed for other reasons.) |
3e558c0
to
39e8145
Compare
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
LGTM! I haven't gone through the fine details, but generally look good.
Pretty straightforward; and this allows parallelization across windows.
Edit: I've added into this the ability to pass in a different set of nodes (added here since it's on top of these changes).