- further OpenMP parallelization, no fortran lines have been changed,

so no new bug is expected for non OpenMP build.


git-svn-id: http://qeforge.qe-forge.org/svn/q-e/trunk/espresso@5718 c92efa57-630b-4861-b058-cf58834340f0
This commit is contained in:
ccavazzoni 2009-07-19 21:08:58 +00:00
parent 1cb5257e94
commit fd57cd39a5
1 changed files with 5 additions and 1 deletions

View File

@ -350,7 +350,8 @@ SUBROUTINE v_xc( rho, rho_core, rhog_core, etxc, vtxc, v )
!
! ... spin-unpolarized case
!
!$omp parallel do private( rhox, arhox, ex, ec, vx, vc ), reduction(+:etxc,vtxc,rhoneg)
!$omp parallel do private( rhox, arhox, ex, ec, vx, vc ), &
!$omp reduction(+:etxc,vtxc), reduction(-:rhoneg)
DO ir = 1, nrxx
!
rhox = rho%of_r(ir,1) + rho_core(ir)
@ -378,6 +379,8 @@ SUBROUTINE v_xc( rho, rho_core, rhog_core, etxc, vtxc, v )
!
! ... spin-polarized case
!
!$omp parallel do private( rhox, arhox, zeta, ex, ec, vx, vc ), &
!$omp reduction(+:etxc,vtxc), reduction(-:rhoneg)
DO ir = 1, nrxx
!
rhox = rho%of_r(ir,1) + rho%of_r(ir,2) + rho_core(ir)
@ -404,6 +407,7 @@ SUBROUTINE v_xc( rho, rho_core, rhog_core, etxc, vtxc, v )
END IF
!
END DO
!$omp end parallel do
!
ELSE IF ( nspin == 4 ) THEN
!