add powerspectrum and cic.

These two functions are extracted from gaepsi.

cic is not particularly fast.
powerspectrum tries to use pyfftw, if not found will use numpy's
slower fft.

powerspectrum is in GADGET normalization.