Optimize average2: use DO CONCURRENT loop form
(cherry picked from commit 1a518cd3f6)
This commit is contained in:
@@ -1832,13 +1832,9 @@ deallocate(f_flat)
|
|||||||
real*8,parameter::C1=3.d0/8.d0,C2=3.d0/4.d0,C3=-1.d0/8.d0
|
real*8,parameter::C1=3.d0/8.d0,C2=3.d0/4.d0,C3=-1.d0/8.d0
|
||||||
integer :: i,j,k
|
integer :: i,j,k
|
||||||
|
|
||||||
do k=1,ext(3)
|
do concurrent (k=1:ext(3), j=1:ext(2), i=1:ext(1))
|
||||||
do j=1,ext(2)
|
|
||||||
do i=1,ext(1)
|
|
||||||
fout(i,j,k) = C1*f1(i,j,k)+C2*f2(i,j,k)+C3*f3(i,j,k)
|
fout(i,j,k) = C1*f1(i,j,k)+C2*f2(i,j,k)+C3*f3(i,j,k)
|
||||||
end do
|
end do
|
||||||
enddo
|
|
||||||
enddo
|
|
||||||
|
|
||||||
return
|
return
|
||||||
|
|
||||||
|
|||||||
Reference in New Issue
Block a user