[WIP] FFT based equilibration detection.

So I'm finding that the FFT-based approach only leads to performance improvements for certain input sizes.  This suggests that we should be conservative, leaving the old code path as default.  

With this conservative approach, I have also made `statsmodels` an *optional* dependency with a delayed import.