MPIExchangeAsync_MappedScalar3D_t Subroutine

public subroutine MPIExchangeAsync_MappedScalar3D_t(this, mesh)

Arguments

TypeIntentOptionalAttributesName
class(MappedScalar3D_t), intent(inout) :: this
type(Mesh3D), intent(inout) :: mesh

Calls

proc~~mpiexchangeasync_mappedscalar3d_t~~CallsGraph proc~mpiexchangeasync_mappedscalar3d_t MPIExchangeAsync_MappedScalar3D_t mpi_irecv mpi_irecv proc~mpiexchangeasync_mappedscalar3d_t->mpi_irecv mpi_isend mpi_isend proc~mpiexchangeasync_mappedscalar3d_t->mpi_isend

Contents


Source Code

  subroutine MPIExchangeAsync_MappedScalar3D_t(this,mesh)
    implicit none
    class(MappedScalar3D_t),intent(inout) :: this
    type(Mesh3D),intent(inout) :: mesh
    ! Local
    integer :: e1,s1,e2,s2,ivar
    integer :: globalSideId,r2,tag
    integer :: iError
    integer :: msgCount

    msgCount = 0

    do ivar = 1,this%nvar
      do e1 = 1,this%nElem
        do s1 = 1,6

          e2 = mesh%sideInfo(3,s1,e1) ! Neighbor Element
          if(e2 > 0) then
            r2 = mesh%decomp%elemToRank(e2) ! Neighbor Rank

            if(r2 /= mesh%decomp%rankId) then

              s2 = mesh%sideInfo(4,s1,e1)/10
              globalSideId = abs(mesh%sideInfo(2,s1,e1))
              tag = globalsideid+mesh%nUniqueSides*(ivar-1)

              msgCount = msgCount+1
              call MPI_IRECV(this%extBoundary(:,:,s1,e1,ivar), &
                             (this%interp%N+1)*(this%interp%N+1), &
                             mesh%decomp%mpiPrec, &
                             r2,globalSideId, &
                             mesh%decomp%mpiComm, &
                             mesh%decomp%requests(msgCount),iError)

              msgCount = msgCount+1
              call MPI_ISEND(this%boundary(:,:,s1,e1,ivar), &
                             (this%interp%N+1)*(this%interp%N+1), &
                             mesh%decomp%mpiPrec, &
                             r2,globalSideId, &
                             mesh%decomp%mpiComm, &
                             mesh%decomp%requests(msgCount),iError)
            endif
          endif

        enddo
      enddo
    enddo

    mesh%decomp%msgCount = msgCount

  endsubroutine MPIExchangeAsync_MappedScalar3D_t