NF: add benchmarking for sphere local maxima

So I don't make it slower by refactoring
RF: refactor local_maxima for index dtype

Refactor the recspeed code to allocate memory via numpy, and use
npy_intp type instead of long (because these can differ).
Merge pull request #135 from matthew-brett/refactor-local-maxima

Refactor local maxima