[SciPy-dev] Single precision FFT

Tom Grydeland tom.grydeland@gmail....
Thu Jan 8 02:00:28 CST 2009


On Wed, Jan 7, 2009 at 4:45 PM, David Cournapeau
<david@ar.media.kyoto-u.ac.jp> wrote:

> I got bored during a meeting: single precision is now implemented in
> scipy.fftpack for complex fft (1d and nd) and 'real' fft (that is real
> input -> hermitian output for rfft). It should be transparent, e.g:
>
> from scipy.fftpack import fft, rfft
> import numpy as np
>
> fft(np.random.randn(10)).dtype == np.cdouble
> fft(np.random.randn(10).astype(np.complex64)).dtype == np.complex64

Excellent!

I'd like to bring the the docstrings up to date on this as well.

"If the input is single-precision (real or complex), the output will
be single-precision also" or do you have other suggestions?

> David


-- 
Tom Grydeland
  <Tom.Grydeland@(gmail.com)>


More information about the Scipy-dev mailing list