Skip to content
Snippets Groups Projects
Commit b85839d8 authored by ESCOBAR Juan's avatar ESCOBAR Juan
Browse files

Juan 14/09/2020 : tridiag_thermo.f90 , OpenACC optimisation -> DO CONCURRENT +...

Juan 14/09/2020 : tridiag_thermo.f90 , OpenACC optimisation -> DO CONCURRENT +  loop gang vector collapse
parent a84e32c8
No related branches found
No related tags found
No related merge requests found
......@@ -371,19 +371,19 @@ END DO !CONCURRENT
!
!$acc loop seq
DO JK = IKB+KKL,IKE-KKL,KKL
!DO CONCURRENT ( JI=1:JIU,JJ=1:JJU)
!$acc loop collapse(2) independent
DO JJ=1,JJU
DO JI=1,JIU
!$acc loop gang, vector collapse(2)
DO CONCURRENT ( JI=1:JIU,JJ=1:JJU)
!DO JJ=1,JJU
! DO JI=1,JIU
ZGAM(JI,JJ,JK) = ZC(JI,JJ,JK-KKL) / ZBET(JI,JJ)
! gam(k) = c(k-1) / bet
ZBET(JI,JJ) = ZB(JI,JJ,JK) - ZA(JI,JJ,JK) * ZGAM(JI,JJ,JK)
! bet = b(k) - a(k)* gam(k)
PVARP(JI,JJ,JK)= ( ZY(JI,JJ,JK) - ZA(JI,JJ,JK) * PVARP(JI,JJ,JK-KKL) ) / ZBET(JI,JJ)
! res(k) = (y(k) -a(k)*res(k-1))/ bet
END DO
END DO
!END DO !CONCURRENT
! END DO
!END DO
END DO !CONCURRENT
END DO
! special treatment for the last level
DO CONCURRENT ( JI=1:JIU,JJ=1:JJU)
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment