Speed up with C extensions

[jimmymathews]

In #338 we saw that there is an opportunity to speed up very specific computations -- in that case bitmasking in a long loop -- by delegation from our primary Python package to a C extension.

For this issue, implement this C extension (guide) and to replace a certain numpy call used in #338. This can serve as a template for further speedup of other specific computations in the future.