16 & ( n, neibpetot, neibpe, stack_import, nod_import, &
17 & stack_export, nod_export, &
18 & ws, wr, x, solver_comm,my_rank)
21 implicit real*8 (a-h,o-z)
25 integer(kind=kint ) ,
intent(in) :: N
26 integer(kind=kint ) ,
intent(in) :: NEIBPETOT
27 integer(kind=kint ),
pointer :: NEIBPE (:)
28 integer(kind=kint ),
pointer :: STACK_IMPORT(:)
29 integer(kind=kint ),
pointer :: NOD_IMPORT (:)
30 integer(kind=kint ),
pointer :: STACK_EXPORT(:)
31 integer(kind=kint ),
pointer :: NOD_EXPORT (:)
32 integer(kind=kint ),
dimension(:) ,
intent(inout):: WS
33 integer(kind=kint ),
dimension(:) ,
intent(inout):: WR
34 integer(kind=kint ),
dimension(:) ,
intent(inout):: X
35 integer(kind=kint ) ,
intent(in) ::SOLVER_COMM
36 integer(kind=kint ) ,
intent(in) :: my_rank
39 integer(kind=kint ),
dimension(:,:),
allocatable :: sta1
40 integer(kind=kint ),
dimension(:,:),
allocatable :: sta2
41 integer(kind=kint ),
dimension(: ),
allocatable :: req1
42 integer(kind=kint ),
dimension(: ),
allocatable :: req2
44 integer(kind=kint ),
save :: NFLAG
47 integer(kind=kint ) :: neib, istart,inum,k,ierr,nreq1,nreq2
51 allocate (sta1(mpi_status_size,neibpetot))
52 allocate (sta2(mpi_status_size,neibpetot))
53 allocate (req1(neibpetot))
54 allocate (req2(neibpetot))
61 istart= stack_export(neib-1)
62 inum = stack_export(neib ) - istart
66 do k= istart+1, istart+inum
67 ws(k)= x(nod_export(k))
69 call mpi_isend (ws(istart+1), inum, mpi_integer, &
70 & neibpe(neib), 0, solver_comm, &
79 istart= stack_import(neib-1)
80 inum = stack_import(neib ) - istart
83 call mpi_irecv (wr(istart+1), inum, mpi_integer, &
84 & neibpe(neib), 0, solver_comm, &
88 call mpi_waitall (nreq2, req2, sta2, ierr)
91 istart= stack_import(neib-1)
92 inum = stack_import(neib ) - istart
93 do k= istart+1, istart+inum
94 x(nod_import(k))= wr(k)
98 call mpi_waitall (nreq1, req1, sta1, ierr)
100 deallocate (sta1, sta2, req1, req2)
subroutine hecmw_solve_send_recv_11i(N, NEIBPETOT, NEIBPE, STACK_IMPORT, NOD_IMPORT, STACK_EXPORT, NOD_EXPORT, WS, WR, X, SOLVER_COMM, my_rank)