The OPAL Discussion Forum

["Dragos Constantin" <dragos.constantin AT gmail.com>]

Hi Philippe,
This reminds me of the FFT function in Matlab. One has to use a function
called 'fftshift' to obtain a correct matrix representation of the outcome.
From the definition of 'fftshift':

"Y = fftshift(X) rearranges a Fourier transform X by shifting the zero-
frequency component to the center of the array."

The plot of the FFT without applying the fftshift looks like your plots.

I am also getting similar results when I reshape a vector to a matrix using a
wrong sequence for the indices. If the reshape order is not done right the
surface plot will look like your plots.

It might be that your post-processing python script is doing something funny
with the indices. Try to swap x and y when you build your matrices and see if
that helps. I am not sure if a simple x-y swap would fix the np dependence
(np=number of mpi processes).

Please report back if you find the issue.

Thanks,
Dragos