gather of 3d array broken

[mcarmesin]

I'm running IGG with ParallelStencils on four GPUs, such that the grid is split in x and y but not in z.
Initializing the local array

ψ.= Data.Array([ (x>-.1) && (y>-.1) && (z>-.1) for x in X_loc, y in Y_loc, z in Z_loc]) 

and gathering with

ψ_nohalo = zeros(ComplexF64,nx-2, ny-2, nz-2);
ψ_nohalo .= Array(ψ[2:end-1,2:end-1,2:end-1])
Ψ = me==0 ? Array{ComplexF64,3}(undef, (nx-2)*dims[1],(ny-2)*dims[2],(nz-2)*dims[3]) : nothing
gather!(ψ_nohalo,Ψ)

I obtain the results of the plots

In the yz plane the data are corrupted like there is some indexing wrong when putting together the data .

[luraess]

Thanks for reporting - we'll have a look.

[omlins]

Hi @mcarmesin , sorry for the delay in reply. Could you please try with ImplicitGlobalGrid@0.14.0? After that, a new gather version was introduced (the current version).

[mcarmesin]

I found the problem: in the y-z plot the the data were erroneously transposed, such that Plots.jl distorted the data while plotting.

Sorry for bothering with my stupidity.

[luraess]

Glad you figured out what was causing the issue. Thanks

[omlins]

I found the problem: in the y-z plot the the data were erroneously transposed, such that Plots.jl distorted the data while plotting.

Sorry for bothering with my stupidity.

No problem. Happy to hear that it works now.