@@ -127,7 +127,6 @@ subroutine SetInteriorFromEquation_MappedVector3D_t(this,geometry,time)
127
127
128
128
end subroutine SetInteriorFromEquation_MappedVector3D_t
129
129
130
-
131
130
subroutine MPIExchangeAsync_MappedVector3D_t (this ,mesh ,resetCount )
132
131
implicit none
133
132
class(MappedVector3D_t),intent (inout ) :: this
@@ -158,23 +157,23 @@ subroutine MPIExchangeAsync_MappedVector3D_t(this,mesh,resetCount)
158
157
159
158
s2 = mesh% sideInfo(4 ,s1,e1)/ 10
160
159
globalSideId = abs (mesh% sideInfo(2 ,s1,e1))
161
- tag = globalsideid+ mesh% nUniqueSides* (ivar-1 + this% nvar* (idir-1 ))
160
+ tag = globalsideid+ mesh% nUniqueSides* (ivar-1 + this% nvar* (idir-1 ))
162
161
163
162
msgCount = msgCount+1
164
163
call MPI_IRECV(this% extBoundary(:,:,s1,e1,ivar,idir), &
165
- (this% interp% N+1 )* (this% interp% N+1 ), &
166
- mesh% decomp% mpiPrec, &
167
- r2,globalSideId, &
168
- mesh% decomp% mpiComm, &
169
- mesh% decomp% requests(msgCount),iError)
164
+ (this% interp% N+1 )* (this% interp% N+1 ), &
165
+ mesh% decomp% mpiPrec, &
166
+ r2,globalSideId, &
167
+ mesh% decomp% mpiComm, &
168
+ mesh% decomp% requests(msgCount),iError)
170
169
171
170
msgCount = msgCount+1
172
171
call MPI_ISEND(this% boundary(:,:,s1,e1,ivar,idir), &
173
- (this% interp% N+1 )* (this% interp% N+1 ), &
174
- mesh% decomp% mpiPrec, &
175
- r2,globalSideId, &
176
- mesh% decomp% mpiComm, &
177
- mesh% decomp% requests(msgCount),iError)
172
+ (this% interp% N+1 )* (this% interp% N+1 ), &
173
+ mesh% decomp% mpiPrec, &
174
+ r2,globalSideId, &
175
+ mesh% decomp% mpiComm, &
176
+ mesh% decomp% requests(msgCount),iError)
178
177
endif
179
178
endif
180
179
@@ -199,7 +198,6 @@ subroutine ApplyFlip_MappedVector3D_t(this,mesh)
199
198
integer :: bcid
200
199
real (prec) :: extBuff(1 :this% interp% N+1 ,1 :this% interp% N+1 )
201
200
202
-
203
201
do idir = 1 ,3
204
202
do ivar = 1 ,this% nvar
205
203
do e1 = 1 ,this% nElem
@@ -208,7 +206,7 @@ subroutine ApplyFlip_MappedVector3D_t(this,mesh)
208
206
e2 = mesh% sideInfo(3 ,s1,e1) ! Neighbor Element
209
207
s2 = mesh% sideInfo(4 ,s1,e1)/ 10
210
208
bcid = mesh% sideInfo(5 ,s1,e1)
211
- if (e2 > 0 ) then ! Interior Element
209
+ if (e2 > 0 ) then ! Interior Element
212
210
r2 = mesh% decomp% elemToRank(e2) ! Neighbor Rank
213
211
214
212
if (r2 /= mesh% decomp% rankId) then
0 commit comments