Lots of API clean up, bug fixes, etc. Code is now slower as I've
realized that there was some missed dependence between R_star and
p_se which can't be segmented out by spatial location. Moving forward,
this should be the target for optimization.