<br><br><div class="gmail_quote">On Sat, Mar 14, 2009 at 2:16 PM, Sturla Molden <span dir="ltr">&lt;<a href="mailto:sturla@molden.no">sturla@molden.no</a>&gt;</span> wrote:<br><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
<br>
<br>
1) I have noticed that fftpack_litemodule.c does not release the GIL<br>
around calls to functions in fftpack.c. I cannot se any obvious reason for<br>
this. As far as I can tell, the functions in fftpack.c are re-entrant.<br>
<br>
2) If fftpack_lite did release the GIL, it would allow functions in<br>
numpy.fft to use multithreading for multiple FFTs in parallel<br>
(threading.Thread are ok, not special compilation needed).<br>
<br>
3) Is there any reason numpy.fft does not have dct? If not, I&#39;d suggest<br>
addition of numpy.fft.dct and numpy.fft.idct.<br>
<br>
4) Regarding ticket #400: Cython now makes this easy. NumPy&#39;s FFTs should<br>
be exposed to C extesions without calling back to Python.<br>
<br>
<br>
Can I open a ticket for this and take care of it? At least 1, 2 and 4<br>
should only take me an hour or so to write, so it might even be ready for<br>
1.3.0.<br>
</blockquote><div><br>Give it a shot. Note that the fft transforms also use int instead of intp, which limits the maximum transform size to 32 bits. Fixing that is somewhere on my todo list but I would be happy to leave it to you ;) Although I expect transforms &gt; 2GB aren&#39;t all that common.<br>
<br>Chuck <br></div><br></div><br>