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

Removing parameter q and bitcount? #352

Open
ksahlin opened this issue Oct 2, 2023 · 1 comment
Open

Removing parameter q and bitcount? #352

ksahlin opened this issue Oct 2, 2023 · 1 comment

Comments

@ksahlin
Copy link
Owner

ksahlin commented Oct 2, 2023

Posting from discussion in #350 (comment)

Should we remove q and Bitcount() in the randstrobe linking function? Bitcount() uses popcount which seems computationally expensive compared to other link methods (I can show you some results from another study offline). q and Bitcount() were initially introduced to skew the sampling towards shorter seeds, which positively affected the shortest reads (=< 100ish), but with above settings, read lengths 50-100 does not even need this computation.

Note: For read lengths were a proper window is used (150 and up) this may affect the recent parameter optimisations.

@marcelm
Copy link
Collaborator

marcelm commented Oct 3, 2023

The std::bitset<64>::count() method compiles down to a single POPCNT machine instruction, so it should be quite fast (with the proper compiler options). I can measure it, though.

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

No branches or pull requests

2 participants