[SciPy-Dev] Accuracy of single-precision FFT

Pauli Virtanen pav@iki...
Fri Jun 25 13:50:40 CDT 2010


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.

-- 
Pauli Virtanen



More information about the SciPy-Dev mailing list