Commit f1f26cc0 authored by Uliana Alekseeva's avatar Uliana Alekseeva

made PGI compiler happy again

parent bfa0f581
......@@ -63,7 +63,8 @@ CONTAINS
CALL vTot%resetPotDen()
CALL vCoul%resetPotDen()
CALL vx%resetPotDen()
ALLOCATE(vx%pw_w,vTot%pw_w,mold=vTot%pw)
ALLOCATE(vx%pw_w,mold=vTot%pw)
ALLOCATE(vTot%pw_w,mold=vTot%pw)
ALLOCATE(vCoul%pw_w(SIZE(den%pw,1),1))
CALL workDen%init(stars,atoms,sphhar,vacuum,input%jspins,noco%l_noco,0)
......
......@@ -112,7 +112,8 @@ contains
end if
rnm(:) = 1.0
PARALLEL_L_LOOP: do concurrent (l = 0: lmax)
!PARALLEL_L_LOOP: do concurrent (l = 0: lmax)
PARALLEL_L_LOOP: do l = 0,lmax
if ( abs( z ) >= l + 1.0 ) then
hl(l) = 0.0
nl(l) = 0.0
......
......@@ -115,7 +115,8 @@ contains
imax = atoms%jri(n)
lmax = sphhar%llh(sphhar%nlh(ns), ns)
if ( potdenType == POTDEN_TYPE_POTYUK ) then
do concurrent (i = 1:imax)
!do concurrent (i = 1:imax)
do i = 1,imax
call ModSphBessel( il(:, i), kl(:, i), input%preconditioning_param * atoms%rmsh(i, n), lmax )
end do
end if
......
......@@ -153,7 +153,8 @@ contains
imax = atoms%jri(n)
lmax = sphhar%llh(sphhar%nlh(nd), nd)
if ( potdenType == POTDEN_TYPE_POTYUK ) then
do concurrent (i = 1:imax)
!do concurrent (i = 1:imax)
do i = 1,imax
call ModSphBessel( il(0:,i), kl(0:,i), input%preconditioning_param * atoms%rmsh(i,n), lmax )
end do
end if
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment