[pooma-dev] FFT

[Gasper Tkacik]

Hi Guenther!

Thank you for the advice. I took a look at the code. However, I somehow 
cannot get the correct results when I use Real-to-Complex transform 
(rfftw) on 2D arrays. Is it possible that this is so because the 
algorithm does in-place transform and lays out the values in the 
half-complex packing? Has this been tested on rfftw or has it only been 
tested with complex-to-complex stuff? I am afraid to use complex because 
of the overhead.

Best regards, Gasper.

Richard Guenther wrote:

>On Thu, 6 Feb 2003, Gasper Tkacik wrote:
>
>  
>
>>Hello everybody!
>>
>>Is there any support for FFT operations in Pooma? For example, given a
>>Array<double, 3>, is there any preferred way of doing FFT?
>>If not, is there any way of doing it in parallel with the constructs
>>that Pooma provides?
>>    
>>
>
>Try looking at src/Transform/WrapFFTW.h which seems to make use of libfftw
>to do FFT.
>