search

spectralDNS / mpiFFT4py

Parallel FFT in 3D or 2D using MPI for Python. Slab or pencil decomposition possible in 3D. Note this rep is being deprecated in favour of mpi4py-fft (https://bitbucket.org/mpi4py/mpi4py-fft)
GNU Lesser General Public License v3.0
21 stars 10 forks source link

demo/transforms_realdata.py fails out of the box #12

Open apatlpo opened 6 years ago

apatlpo commented 6 years ago

The error is typically:

    fft = PFFT(MPI.COMM_WORLD, N, collapse=False, slab=2)
  File "/Users/aponte/.miniconda3/envs/shenfun/lib/python3.6/site-packages/mpi4py_fft/mpifft.py", line 106, in __init__
    assert axes[-1] != axis

slab=1 seems to fix the issue. Where can I find a description of what options corresponds to? thanks

mikaem commented 6 years ago

Hi Aurélien The "error" you're mentioning is in mpi4py-fft, not mpiFFT4py. The mistake is understandable, considering the two names:-)

The PFFT class has a keyword attribute slab, but it is only boolean, True or False. To use slab decomposition or not. Default is pencil, so then slab=False.

apatlpo commented 6 years ago

ok, thanks !