From 9986ff09ddd0307d81076ebe11b8ea24fb20e31c Mon Sep 17 00:00:00 2001 From: Abdallah Ammar Date: Tue, 5 Nov 2024 07:56:10 +0100 Subject: [PATCH] q -> p in RGW_self_energy --- src/GW/RGW_self_energy.f90 | 8 ++++---- 1 file changed, 4 insertions(+), 4 deletions(-) diff --git a/src/GW/RGW_self_energy.f90 b/src/GW/RGW_self_energy.f90 index 0f0ce8a..b3d9c2c 100644 --- a/src/GW/RGW_self_energy.f90 +++ b/src/GW/RGW_self_energy.f90 @@ -93,7 +93,7 @@ subroutine RGW_self_energy(eta,nBas,nOrb,nC,nO,nV,nR,nS,e,Om,rho,EcGM,Sig,Z) !$OMP PARALLEL & !$OMP SHARED(Sig,Z,rho,eta,nS,nC,nO,nOrb,nR,e,Om) & -!$OMP PRIVATE(m,i,q,p,eps,num) & +!$OMP PRIVATE(m,i,p,eps,num) & !$OMP DEFAULT(NONE) !$OMP DO do p=nC+1,nOrb-nR @@ -101,7 +101,7 @@ subroutine RGW_self_energy(eta,nBas,nOrb,nC,nO,nV,nR,nS,e,Om,rho,EcGM,Sig,Z) do i=nC+1,nO eps = e(p) - e(i) + Om(m) - num = 2d0*rho(p,i,m)*rho(q,i,m) + num = 2d0*rho(p,i,m)*rho(p,i,m) Z(p) = Z(p) - num*(eps**2 - eta**2)/(eps**2 + eta**2)**2 end do @@ -114,7 +114,7 @@ subroutine RGW_self_energy(eta,nBas,nOrb,nC,nO,nV,nR,nS,e,Om,rho,EcGM,Sig,Z) !$OMP PARALLEL & !$OMP SHARED(Sig,Z,rho,eta,nS,nC,nO,nOrb,nR,e,Om) & - !$OMP PRIVATE(m,a,q,p,eps,num) & + !$OMP PRIVATE(m,a,p,eps,num) & !$OMP DEFAULT(NONE) !$OMP DO do p=nC+1,nOrb-nR @@ -122,7 +122,7 @@ subroutine RGW_self_energy(eta,nBas,nOrb,nC,nO,nV,nR,nS,e,Om,rho,EcGM,Sig,Z) do a=nO+1,nOrb-nR eps = e(p) - e(a) - Om(m) - num = 2d0*rho(p,a,m)*rho(q,a,m) + num = 2d0*rho(p,a,m)*rho(p,a,m) Z(p) = Z(p) - num*(eps**2 - eta**2)/(eps**2 + eta**2)**2 end do