PTEROSOR/QuAcK
10
1
Fork 0
mirror of https://github.com/pfloos/quack synced 2024-12-23 04:43:42 +01:00
Code Releases Activity

kappa in reg

Browse Source
This commit is contained in:
Pierre-Francois Loos 2023-08-01 17:18:10 +02:00
parent f3f2135074
commit f3f97a1daf
3 changed files with 6 additions and 6 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

4
src/GT/GTeh_regularization.f90
View File

@ -37,14 +37,14 @@ subroutine GTeh_regularization(nBas,nC,nO,nR,nS,e,Om,rhoL,rhoR)
do i=nC+1,nO
Dpim = e(p) - e(i) + Om(m)
kappa = 1d0 - exp(-Dpim*Dpim*s)
kappa = exp(-Dpim*Dpim*s)
rhoL(i,p,m) = kappa*rhoL(i,p,m)
rhoR(i,p,m) = kappa*rhoR(i,p,m)
enddo
do a=nO+1,nBas-nR
Dpam = e(p) - e(a) - Om(m)
kappa = 1d0 - exp(-Dpam*Dpam*s)
kappa = exp(-Dpam*Dpam*s)
rhoL(p,a,m) = kappa*rhoL(p,a,m)
rhoR(p,a,m) = kappa*rhoR(p,a,m)
enddo

4
src/GT/GTpp_regularization.f90
View File

@ -40,7 +40,7 @@ subroutine GTpp_regularization(nBas,nC,nO,nR,nOO,nVV,e,Om1,rho1,Om2,rho2)
do m=1,nVV
do i=nC+1,nO
Dpim = e(p) + e(i) - Om1(m)
kappa = 1d0 - exp(-Dpim*Dpim*s)
kappa = exp(-Dpim*Dpim*s)
rho1(p,i,m) = kappa*rho1(p,i,m)
end do
end do
@ -48,7 +48,7 @@ subroutine GTpp_regularization(nBas,nC,nO,nR,nOO,nVV,e,Om1,rho1,Om2,rho2)
do m=1,nOO
do a=nO+1,nBas-nR
Dpam = e(p) + e(a) - Om2(m)
kappa = 1d0 - exp(-Dpam*Dpam*s)
kappa = exp(-Dpam*Dpam*s)
rho2(p,a,m) = kappa*rho2(p,a,m)
end do
end do

4
src/GW/GW_regularization.f90
View File

@ -36,13 +36,13 @@ subroutine GW_regularization(nBas,nC,nO,nR,nS,e,Om,rho)
do i=nC+1,nO
Dpim = e(p) - e(i) + Om(m)
kappa = 1d0 - exp(-Dpim*Dpim*s)
kappa = exp(-Dpim*Dpim*s)
rho(p,i,m) = kappa*rho(p,i,m)
enddo
do a=nO+1,nBas-nR
Dpam = e(p) - e(a) - Om(m)
kappa = 1d0 - exp(-Dpam*Dpam*s)
kappa = exp(-Dpam*Dpam*s)
rho(p,a,m) = kappa*rho(p,a,m)
enddo