spmd_fvb_apar.F File Reference

#include "implicit_f.inc"
#include "spmd.inc"
#include "task_c.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	spmd_fvb_apar (nelcut, nbf, nesct, nerwl, nesrg, nesmd1, fvpbuf)

Function/Subroutine Documentation

spmd_fvb_apar()

subroutine spmd_fvb_apar	(	integer	nelcut,
		integer	nbf,
		integer	nesct,
		integer	nerwl,
		integer	nesrg,
		integer	nesmd1,
		integer, dimension(*)	fvpbuf )

Definition at line 32 of file spmd_fvb_apar.F.

C-----------------------------------------------
C   M o d u l e s 
C-----------------------------------------------
      USE fvbag_mod
C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
        USE spmd_comm_world_mod, ONLY : spmd_comm_world
#include      "implicit_f.inc"
C-----------------------------------------------------------------
C   M e s s a g e   P a s s i n g
C-----------------------------------------------
#include "spmd.inc"
C-----------------------------------------------
C   C o m m o n   B l o c k s
C-----------------------------------------------
#include      "task_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER NELCUT, NBF, NESCT, NERWL, NESRG, NESMD1, FVPBUF(*)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
#ifdef MPI
      INTEGER FVIAD, I, PMAIN, J, K, KK, N, NN, ITAG, MSGOFF,MSGOFF2,
     .        NPOLH_ANIM, IERR, STAT(MPI_STATUS_SIZE), IAD, IADP
C-----------------------------------------------
      INTEGER, DIMENSION(:), ALLOCATABLE :: ITAGT, IBUF
C-----------------------------------------------
      DATA msgoff/7043/
      DATA msgoff2/7044/
C
      iadp=0
      fviad=nelcut+nbf+nesct+nerwl+nesrg+nesmd1
      DO i=1,nfvbag
         pmain=fvspmd(i)%PMAIN
         IF (ispmd==0) THEN
            IF (ispmd==pmain-1) THEN
               ALLOCATE(itagt(fvdata(i)%NNTR))
               DO j=1,fvdata(i)%NNTR
                  itagt(j)=0
               ENDDO
C
               DO j=1,fvdata(i)%NPOLH_ANIM
                  DO k=fvdata(i)%IFVPADR_ANIM(j),
     .                        fvdata(i)%IFVPADR_ANIM(j+1)-1
                     kk=fvdata(i)%IFVPOLH_ANIM(k)
                     DO n=fvdata(i)%IFVTADR_ANIM(kk),
     .                          fvdata(i)%IFVTADR_ANIM(kk+1)-1
                        nn=fvdata(i)%IFVPOLY_ANIM(n)
                        IF (itagt(nn)==0) THEN
                           fviad=fviad+1
                           itagt(nn)=1
                        ENDIF
                     ENDDO
                  ENDDO
                  fvpbuf(iadp+j)=fviad
               ENDDO
               iadp=iadp+fvdata(i)%NPOLH_ANIM
C
               DEALLOCATE(itagt)
            ELSE
               itag=msgoff
               CALL mpi_recv(npolh_anim, 1, mpi_integer, 
     .                       it_spmd(pmain), itag, spmd_comm_world,
     .                       stat, ierr)
C
               ALLOCATE(ibuf(npolh_anim))
               itag=msgoff2
               CALL mpi_recv(ibuf, npolh_anim, mpi_integer,
     .                       it_spmd(pmain), itag, spmd_comm_world,
     .                       stat, ierr)
C
               DO j=1,npolh_anim
                  fvpbuf(iadp+j)=fviad+ibuf(j)
               ENDDO
               fviad=fviad+ibuf(npolh_anim)
               iadp=iadp+npolh_anim
               DEALLOCATE(ibuf)
            ENDIF
         ELSE
            IF (ispmd==pmain-1) THEN
               itag=msgoff
               CALL mpi_send(fvdata(i)%NPOLH_ANIM, 1, mpi_integer,
     .                       it_spmd(1), itag, spmd_comm_world,
     .                       ierr)
C
               ALLOCATE(ibuf(fvdata(i)%NPOLH_ANIM), 
     .                       itagt(fvdata(i)%NNTR))
               iad=0
C
               DO j=1,fvdata(i)%NNTR
                  itagt(j)=0
               ENDDO
C
               DO j=1,fvdata(i)%NPOLH_ANIM
                  DO k=fvdata(i)%IFVPADR_ANIM(j),
     .                        fvdata(i)%IFVPADR_ANIM(j+1)-1
                     kk=fvdata(i)%IFVPOLH_ANIM(k)
                     DO n=fvdata(i)%IFVTADR_ANIM(kk),
     .                          fvdata(i)%IFVTADR_ANIM(kk+1)-1
                        nn=fvdata(i)%IFVPOLY_ANIM(n)
                        IF (itagt(nn)==0) THEN
                           iad=iad+1
                           itagt(nn)=1
                        ENDIF
                     ENDDO
                  ENDDO
                  ibuf(j)=iad
               ENDDO
               itag=msgoff2
               CALL mpi_send(ibuf, fvdata(i)%NPOLH_ANIM, mpi_integer,
     .                       it_spmd(1), itag, spmd_comm_world,
     .                        ierr)
C
               DEALLOCATE(ibuf, itagt)                              
            ENDIF
         ENDIF
      ENDDO
C
#endif
      RETURN