API: Move estimate_rank function out of raw
ENH+FIX: Making CUDA resampling more efficient, fixing DC/Nyquist error