Changeset 999 for codes/icosagcm/trunk/src/parallel/transfert_requests.F90

Timestamp:

01/10/20 12:20:51 (4 years ago)

Author:

adurocher

Message:

transfert_mpi : aggregate messages

Only one message is sent per MPI process pair.
One message was sent for each pair of tiles, and each MPI process can have multiple tiles.

File:

• codes/icosagcm/trunk/src/parallel/transfert_requests.F90 (modified) (2 diffs)

codes/icosagcm/trunk/src/parallel/transfert_requests.F90

@@ -218 +218 @@
       integer :: k, i, last, i_min, tmp
 
-      displs = (points%i-1) + (points%j-1)*iim
-      if(allocated(points%elt)) displs = displs + (points%elt-1)*iim*jjm
+      displs = (points%i(1:points%npoints)-1) + (points%j(1:points%npoints)-1)*iim
+      if(allocated(points%elt)) displs = displs + (points%elt(1:points%npoints)-1)*iim*jjm
 
       last=0
@@ -226 +226 @@
         tmp = displs(i_min)
         displs(i_min) = displs(i)
-        if(last==0 .or. displs(i_min) /= displs(last)) last=last+1
+        if( last==0 ) then
+          last=last+1
+        else if( displs(i_min) /= displs(last) ) then 
+          last=last+1
+        endif
         displs(last) = tmp
       end do