Skip to content
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

faster min_node_support when support is 0 #4095

Merged
merged 1 commit into from
Sep 22, 2023
Merged

faster min_node_support when support is 0 #4095

merged 1 commit into from
Sep 22, 2023

Conversation

glennhickey
Copy link
Contributor

@glennhickey glennhickey commented Sep 21, 2023

Changelog Entry

To be copied to the draft changelog by merger:

  • vg call patched so that certain problem cases no longer take forever.

Description

I'm not sure if there's a github issue for this, but I think it comes from here:
https://www.biostars.org/p/9572558/

and @jeizenga was kind enough to pass over the data, with which it is easy to reproduce.

But I still don't really understand what's going on. Usually it boils down to one big snarl, but in this case it just seems to be going a bit slow on every single snarl. And whenever I pause it, it's always inside PackedTraversalSupportFinder::get_min_node_support()

I don't know why this function would cause a slowdown on this particular graph, and I suspect there's something going on that I should be figuring out but am not.

The good news is that there's a fairly obvious speedup: quit computing the min() if 0 is ever found. And... somehow that seems to be enough to get it to run through on this input.

@jeizenga
Copy link
Contributor

This Biostars question, actually: https://www.biostars.org/p/9572558/

@jeizenga
Copy link
Contributor

A thought I just had: someone pointed in this comment that this is a very repetitive genome stemming from LTR transposon activity. Maybe it could have something to do with that?

@glennhickey glennhickey merged commit a97a233 into master Sep 22, 2023
1 check passed
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

Successfully merging this pull request may close these issues.

3 participants