When using r2c and c2r transformations, the FFTW library pads the contiguous dimension of the complex array in the case of an even number of entries. In the case of MPI-parallel r2c and c2r transformations, even the real-space arrays need padding. The current implementation of field does not consider this, which will inevitable lead to a bug
Description
Description
Status | Assigned | Task | ||
---|---|---|---|---|
Open | junge | T1843 µSpectre Release 0.9 β | ||
Resolved | junge | T1747 Non-Trivial Memory Layout for FFT arrays |
Event Timeline
Comment Actions
It seems that this problem is better worked around by locally copying data and then using in-place transforms. The alternative kills the eigen interface to fields, which is important in the design of the solvers