Hello,
I was computing some correlations while I noticed that I can easily make an implementation significantly faster than coffcoef by replacing sum() with count_nonzero(). Therefore, I think whenever sum() is called on a boolean array, it should use count_nonzero or such instead of addition.