********> bugfix.30 Author: Mike Crowley Date: Sept. 1998 Programs: sander Severity: PME runs using MPI will fail for some implementations of MPI: seems to affect Solaris and recent SGI MPI's Description: A communication array used in the FFT transpose was dimensioned incorrectly; there is no work-around Fix: *** ew_fft.f Fri Sep 18 12:15:46 1998 --- ew_fft.f.50bug Fri Sep 18 12:11:21 1998 *************** *** 988,994 **** _REAL_ targ(*), src(*), tmp(*), tmp1(*) integer n3, ldx, ktask, num_recv, numval, numtrans, numtries integer i,j,k,ks,k0,k00,jjtask,jtask,ibuff,ifoo,rtask ! integer maxval,tmpval,ii,ishift #include "parallel.h" #include "mpif.h" --- 988,994 ---- _REAL_ targ(*), src(*), tmp(*), tmp1(*) integer n3, ldx, ktask, num_recv, numval, numtrans, numtries integer i,j,k,ks,k0,k00,jjtask,jtask,ibuff,ifoo,rtask ! integer maxval,tmpval,ii,ishift,ireq,isnd_stat #include "parallel.h" #include "mpif.h" *************** *** 996,1002 **** #ifdef TERRA # include "tlibs/Tcommlib_f.h" #endif - integer ireq,isnd_stat(MPI_STATUS_SIZE) rtask = mytaskid ishift = 1 --- 996,1001 ---- *************** *** 1116,1122 **** c _REAL_ targ(*), src(*), tmp(*), tmp1(*) integer n3, ldx, ktask, num_recv, numval, numtrans, numtries integer i,j,k,ks,k0,k00,jjtask,jtask, ibuff,ifoo,rtask ! integer maxval, tmpval, ii, ishift #include "parallel.h" #include "mpif.h" --- 1115,1121 ---- c _REAL_ targ(*), src(*), tmp(*), tmp1(*) integer n3, ldx, ktask, num_recv, numval, numtrans, numtries integer i,j,k,ks,k0,k00,jjtask,jtask, ibuff,ifoo,rtask ! integer maxval, tmpval, ii, ishift, ireq, isnd_stat #include "parallel.h" #include "mpif.h" *************** *** 1124,1131 **** #ifdef TERRA # include "tlibs/Tcommlib_f.h" #endif - integer ireq,isnd_stat(MPI_STATUS_SIZE) - rtask = mytaskid ishift = 1 do jjtask = mytaskid 1, mytaskid numtasks-1 --- 1123,1128 ----