jipolanco / PencilFFTs.jl

Fast Fourier transforms of MPI-distributed Julia arrays
https://jipolanco.github.io/PencilFFTs.jl/dev/
MIT License
77 stars 7 forks source link

In-place real-to-complex FFTs? #21

Open jipolanco opened 3 years ago

jipolanco commented 3 years ago

Not that simple due to the change from a real to a complex element type.

This will require changes to ManyPencilArray in PencilArrays, to allow for variable element types of the aliased arrays. I'm not sure if this can be done efficiently with reinterpret, and unsafe workarounds may be needed.

Also, in-place r2c FFTs are not currently wrapped by FFTW.jl (due to this and other complications), and custom wrappers will likely be needed.