Holy crimminey, I'm totally lost in your Fortran syntax. :-)
What you describe might be a bug in our MPI_IN_PLACE handling for
MPI_ALLREDUCE.
Could you possible make a small test case that a) we can run, and b) uses
straightforward Fortran? (avoid using terms like "assumed shape" and "assumed
Note also that coding the mpi_allreduce as:
call
mpi_allreduce(MPI_IN_PLACE,phim(0,1,1,1,grp),phim_size*im*jm*kmloc(coords(2)+1),mpi_real,mpi_sum,ang_com,ierr)
results in the same freezing behavior in the 60th iteration. (I don't
recall why the arrays were being passed, possibly just a mistak
I am seeing mpi_allreduce operations freeze execution of my code on some
moderately-sized problems. The freeze does not manifest itself in every
problem. In addition, it is in a portion of the code that is repeated many
times. In the problem discussed below, the problem appears in the 60th
itera
