Page MenuHomec4science

Non-Trivial Memory Layout for FFT arrays
Closed, ResolvedPublic

Description

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

Event Timeline

junge created this task.Jan 3 2018, 14:33
junge created this object in space S1 c4science.
junge created this object with visibility "µSpectre (Project)".
junge created this object with edit policy "Subscribers".
junge closed this task as Resolved.Jan 3 2018, 18:49

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

junge moved this task from Backlog to Done on the µSpectre board.Jan 16 2018, 10:13
junge changed the visibility from "µSpectre (Project)" to "Public (No Login Required)".Jan 31 2018, 13:03
junge changed the edit policy from "Subscribers" to "µSpectre (Project)".