[SciPy-Dev] Accuracy of single-precision FFT

Charles R Harris charlesr.harris@gmail....
Fri Jun 25 14:39:49 CDT 2010


On Fri, Jun 25, 2010 at 12:50 PM, Pauli Virtanen <pav@iki.fi> wrote:

> Fri, 25 Jun 2010 12:38:44 -0400, Anne Archibald wrote:
> [clip]
> > It seems to me that the right short-term solution is to switch to double
> > precision for difficult sizes (any non-power-of-two is probably fine)
> > and add a warning to the docstring (in particular mentioning the
> > horrible performance as well as the roundoff error, which can be an
> > issue for doubles too). This we can do for 0.8.
>
> Agreed. Done here:
>
>
> http://github.com/pv/scipy-work/compare/0.8.x...disable-single-precision-fft
>
> Seems to work, I think it should go in to the 0.8.x branch ASAP.
>
>
I'd also keep at least power of 5 also so that round decimal numbers work
well.

Chuck
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mail.scipy.org/pipermail/scipy-dev/attachments/20100625/09207761/attachment.html 


More information about the SciPy-Dev mailing list