[SciPy-dev] Single precision FFT
Tom Grydeland
tom.grydeland at gmail.com
Thu Jan 8 03:00:28 EST 2009
On Wed, Jan 7, 2009 at 4:45 PM, David Cournapeau
<david at 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