inc_gatherHyb.f

Go to the documentation of this file.

      do i = 1, noofassites
         call mpi_gather(asobssites(i).age, 1, mpi_real8,
     *       ager, 1, mpi_real8, 0, mpi_comm_world, icon)
         call mpi_gather(asobssites(i).esize, 1, mpi_real8,
     *       esizer, 1, mpi_real8, 0, mpi_comm_world, icon)
         if(mpirank .eq. 0) then
            wage = 0.
            wsize = 0.
            do j = 1, mpisize
               wage = wage + ager(j)*esizer(j)
               wsize = wsize + esizer(j)
            enddo
            if(wsize .gt. 0.) then
               asobssites(i).age =wage/wsize
               asobssites(i).esize = wsize
            else
               asobssites(i).age =0.
               asobssites(i).esize =0.
            endif
         endif
      enddo
!       MPI_IN_PLACE dose not work at KEKB
      if( mpirank .eq. 0 ) then
         do i = 1, noofsites
            nga(i) =0.
            nea(i) =0.
            nmua(i) = 0.
            nhada(i) = 0.
            sumelossa(i) = 0.
         enddo
      endif

!      if(mpirank .eq. 0) then
        call mpi_reduce(ng, nga, noofsites,
     *   mpi_real8, mpi_sum, 0, mpi_comm_world, icon)
!         call mpi_reduce(MPI_IN_PLACE, Ng, NoOfSites,
!     *   MPI_REAL8, MPI_SUM, 0, MPI_COMM_WORLD, icon)

        call mpi_reduce(ne, nea, noofsites,
     *   mpi_real8, mpi_sum, 0, mpi_comm_world, icon)
!         call mpi_reduce(MPI_IN_PLACE, Ne, NoOfSites,
!     *   MPI_REAL8, MPI_SUM, 0, MPI_COMM_WORLD, icon)

       call mpi_reduce(nmu, nmua, noofsites,
     *   mpi_real8, mpi_sum, 0, mpi_comm_world, icon)
!         call mpi_reduce(MPI_IN_PLACE, Nmu, NoOfSites,
!     *   MPI_REAL8, MPI_SUM, 0, MPI_COMM_WORLD, icon)

       call mpi_reduce(nhad, nhada, noofsites,
     *   mpi_real8, mpi_sum, 0, mpi_comm_world, icon)
!         call mpi_reduce(MPI_IN_PLACE, Nhad, NoOfSites,
!     *   MPI_REAL8, MPI_SUM, 0, MPI_COMM_WORLD, icon)

       call mpi_reduce(sumeloss, sumelossa, noofsites,
     *   mpi_real8, mpi_sum, 0, mpi_comm_world, icon)
!         call mpi_reduce(MPI_IN_PLACE, SumEloss, NoOfSites,
!     *   MPI_REAL8, MPI_SUM, 0, MPI_COMM_WORLD, icon)

!     else

!         call mpi_reduce(Ng, Ng, NoOfSites,
!     *   MPI_REAL8, MPI_SUM, 0, MPI_COMM_WORLD, icon)

!         call mpi_reduce(Ne, Ne, NoOfSites,
!     *   MPI_REAL8, MPI_SUM, 0, MPI_COMM_WORLD, icon)

!         call mpi_reduce(Nmu, Nmu, NoOfSites,
!     *   MPI_REAL8, MPI_SUM, 0, MPI_COMM_WORLD, icon)

!         call mpi_reduce(Nhad, Nhad, NoOfSites,
!     *   MPI_REAL8, MPI_SUM, 0, MPI_COMM_WORLD, icon)

!         call mpi_reduce(SumEloss, SumEloss, NoOfSites,
!     *   MPI_REAL8, MPI_SUM, 0, MPI_COMM_WORLD, icon)

!      endif