Hello,
i'm studying the parallelized version of a solving 2D heat equation code in
order to understand cartesian topology and the famous "MPI_CART_SHIFT".
Here's my problem at this part of the code :
-------------------------------------------------------------------------------------------------------------------------------------------
call MPI_INIT(infompi)
comm = MPI_COMM_WORLD
call MPI_COMM_SIZE(comm,nproc,infompi)
call MPI_COMM_RANK(comm,me,infompi)
!
......
! Create 2D cartesian grid
periods(:) = .false.
ndims = 2
dims(1)=x_domains
dims(2)=y_domains
CALL MPI_CART_CREATE(MPI_COMM_WORLD, ndims, dims, periods, &
reorganisation,comm2d,infompi)
!
! Identify neighbors
!
NeighBor(:) = MPI_PROC_NULL
! Left/West and right/Est neigbors
CALL MPI_CART_SHIFT(comm2d,0,1,NeighBor(W),NeighBor(E),infompi)
print *,'rank=', me
print *, 'here first mpi_cart_shift : neighbor(w)=',NeighBor(W)
print *, 'here first mpi_cart_shift : neighbor(e)=',NeighBor(E)
...
-------------------------------------------------------------------------------------------------------------------------------------------
with x_domains=y_domains=2
and i get at the execution :" mpirun -np 4 ./explicitPar"
rank= 0
here first mpi_cart_shift : neighbor(w)= -1
here first mpi_cart_shift : neighbor(e)= 2
rank= 3
here first mpi_cart_shift : neighbor(w)= 1
here first mpi_cart_shift : neighbor(e)= -1
rank= 2
here first mpi_cart_shift : neighbor(w)= 0
here first mpi_cart_shift : neighbor(e)= -1
rank= 1
here first mpi_cart_shift : neighbor(w)= -1
here first mpi_cart_shift : neighbor(e)= 3
I saw that if the rank is out of the topology and wihtout periodicity, the
rank should be equal to MPI_UNDEFINED whis is assigned to -32766 in "mpif.h"
. So, why have i got the value "-1" ?
On my Macbook pro, i get the value "-2".
Thanks in advance.
|
![[Flag]](../../../../../community/mirrors/flag-us.gif){kind=small}