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

clean up RGWC

Browse Source
This commit is contained in:
Pierre-Francois Loos 2024-02-08 16:42:44 +01:00
parent 959cc22535
commit 0e974dcfdb
2 changed files with 4 additions and 4 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

4
src/GW/GW_plot_self_energy.f90
View File

@ -43,8 +43,8 @@ subroutine GW_plot_self_energy(nBas,eta,nC,nO,nV,nR,nS,eHF,eGW,Om,rho)
! Minimum and maximum frequency values
wmin = -5d0
wmax = +5d0
wmin = -10d0
wmax = +0d0
dw = (wmax - wmin)/dble(ngrid)
do g=1,nGrid

4
src/GW/RGWC.f90
View File

@ -29,7 +29,7 @@ subroutine RGWC(dotest,eta,nBas,nC,nO,nV,nR,nS,Om,rho,eHF,e,eGW,Z)
integer :: p,q,i,a,m
integer :: iSat
double precision :: num,eps
double precision,parameter :: cutoff = 1d-2
double precision,parameter :: cutoff = 0d-3
logical,parameter :: do_hole_branch = .true.
logical,parameter :: do_electron_branch = .false.
@ -190,7 +190,7 @@ subroutine RGWC(dotest,eta,nBas,nC,nO,nV,nR,nS,Om,rho,eHF,e,eGW,Z)
! Construct grid
nGrid = 1000
nGrid = 5000
allocate(w(nGrid),AGWC(nBas,nGrid))
! Minimum and maximum frequency values