scnc
/
quantum-espresso
mirror of https://gitlab.com/QEF/q-e.git
1
0
Fork 0
Code Issues Packages Projects Releases Wiki Activity

Cray compiler directive 

!$DIRS ...
enclosed in
#ifdef __CRAY
#endif


git-svn-id: http://qeforge.qe-forge.org/svn/q-e/trunk/espresso@560 c92efa57-630b-4861-b058-cf58834340f0
This commit is contained in:
cavazzon 2004-02-05 13:18:10 +00:00
parent 1a1a14859f
commit 7aca5c6a92
1 changed files with 29 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

29
Modules/ptoolkit.f90
View File

@ -725,13 +725,18 @@
! ...... COSTRUZIONE DEL VETTORE U
#if defined __T3E
!DIR$ UNROLL 8
#endif
DO k = 1,L
vtmp(k) = 0.0d0
END DO
k = ME + 1
#if defined __T3E
!DIR$ UNROLL 4
#endif
DO kl = 1,is(l)
vtmp(k) = A(kl,I)
UL(kl) = A(kl,I)
@ -757,18 +762,28 @@
! ...... COSTRUZIONE DEL VETTORE P
KAPPA = 0.0d0
DO J = 1,L
vtmp(j) = 0.0d0
#if defined __T3E
!DIR$ UNROLL 8
#endif
DO KL = 1, IS(J)
vtmp(J) = vtmp(J) + A(KL,J) * UL(KL)
END DO
IF(L.GT.J .AND. ME.EQ.RI(J)) then
#if defined __T3E
!DIR$ UNROLL 8
#endif
DO K = J+1,L
vtmp(J) = vtmp(J) + A(IS(J),K) * U(K)
END DO
END IF
vtmp(J) = vtmp(J) * ONE_OVER_H
KAPPA = KAPPA + vtmp(J) * U(J) * 0.5d0 * ONE_OVER_H
END DO
@ -792,7 +807,9 @@
CALL DAXPY(l, -kappa, u, 1, p, 1)
k = me + 1
#if defined __T3E
!DIR$ UNROLL 8
#endif
DO kl = 1,is(l)
PL(kl) = P(k)
k = k + NPROC
@ -801,7 +818,9 @@
DO J = 1,L
tu= U(J)
tp= P(J)
#if defined __T3E
!DIR$ UNROLL 12
#endif
DO KL = 1,is(l)
A(KL,j) = A(KL,j) - UL(KL) * tp
A(KL,j) = A(KL,j) - PL(KL) * tu
@ -844,7 +863,9 @@
! t1 = mclock()
DO J = 1,N
#if defined __T3E
!DIR$ UNROLL 12
#endif
DO I = 1,NRL
V(I,J) = 0.0d0
END DO
@ -1079,15 +1100,21 @@
end do
do i=m-1,l,-1
#if defined __T3E
!DIR$ UNROLL 12
#endif
do k=1,nrl
fv2(k) =z(k,i+1)
end do
#if defined __T3E
!DIR$ UNROLL 12
#endif
do k=1,nrl
fv1(k) =z(k,i)
end do
#if defined __T3E
!DIR$ UNROLL 12
#endif
do k=1,nrl
z(k,i+1) =sv(i)*fv1(k) + cv(i)*fv2(k)
z(k,i) =cv(i)*fv1(k) - sv(i)*fv2(k)
@ -1145,7 +1172,9 @@
!
! Exchange local elements of eigenvectors.
!
#if defined __T3E
!DIR$ UNROLL 12
#endif
do j=1,nrl
p=v(j,i)
v(j,i)=v(j,k)