Commit 2061ad19 authored by Henning Janssen's avatar Henning Janssen

Small change in get_gf/set_gf calls

parent d50adb78
......@@ -50,7 +50,7 @@ MODULE m_gfDOS
re = 0.0
DO ipm = 1, 2 !Sum over G^+ and G^-
!Get the full gf matrix at the energy point
CALL g%get_gf(gmat,atoms,input,iz,l,nType,ipm.EQ.2)
CALL g%get(gmat,atoms,input,iz,l,nType,ipm.EQ.2)
!Convert to eV^-1
gmat%data_c = gmat%data_c/hartree_to_eV_const
!Calculate up/down dos
......
......@@ -42,7 +42,7 @@ MODULE m_hybridization
!--------------------------------------------------
! Get the full Greens function matrix for the current energy point
!--------------------------------------------------
CALL gf%get_gf(gmat,atoms,input,iz,l,nType,ipm.EQ.2)
CALL gf%get(gmat,atoms,input,iz,l,nType,ipm.EQ.2)
!--------------------------------------------------
! Invert the matrix using the routines in types_mat
!--------------------------------------------------
......@@ -54,7 +54,7 @@ MODULE m_hybridization
ENDDO
Delta(iz) = -1/(tpi_const*gmat%matsize1) * AIMAG(tr)
WRITE(1337,"(2f14.8)") REAL(gf%e(iz)-ef)*hartree_to_ev_const, Delta(iz)
!Free up the gmat matrix (it is initialized in gf%get_gf)
!Free up the gmat matrix (it is initialized in gf%get)
CALL gmat%free()
ENDDO
......
......@@ -107,7 +107,7 @@ MODULE m_j0
delta%data_c = 0.0
DO ispin = 1, input%jspins
DO ipm = 1, 2
CALL g0%get_gf(calc,atoms,input,iz,l,nType,ipm.EQ.2,spin=ispin)
CALL g0%get(calc,atoms,input,iz,l,nType,ipm.EQ.2,spin=ispin)
CALL calc%inverse()
delta%data_c = delta%data_c + 1/2.0 * (-1)**(ispin-1) * calc%data_c
ENDDO
......@@ -123,8 +123,8 @@ MODULE m_j0
sumdwn = 0.0
sumupdwn = 0.0
DO ipm = 1, 2
CALL g0%get_gf(calcup,atoms,input,iz,l,nType,ipm.EQ.2,spin=1)
CALL g0%get_gf(calcdwn,atoms,input,iz,l,nType,ipm.EQ.2,spin=2)
CALL g0%get(calcup,atoms,input,iz,l,nType,ipm.EQ.2,spin=1)
CALL g0%get(calcdwn,atoms,input,iz,l,nType,ipm.EQ.2,spin=2)
calcup%data_c = matmul(delta%data_c,calcup%data_c)
calcdwn%data_c = matmul(delta%data_c,calcdwn%data_c)
calc%data_c = matmul(calcup%data_c,calcdwn%data_c)
......
......@@ -64,7 +64,7 @@ CONTAINS
!Integrate over the contour:
DO iz = 1, g%nz
!get the corresponding gf-matrix
CALL g%get_gf(gmat,atoms,input,iz,l,nType,ipm.EQ.2,spin=ispin,lp=lp,nTypep=nTypep)
CALL g%get(gmat,atoms,input,iz,l,nType,ipm.EQ.2,spin=ispin,lp=lp,nTypep=nTypep)
ind1 = 0
DO m = -l, l
ind1 = ind1 + 1
......@@ -80,7 +80,7 @@ CONTAINS
!For the contour 3 (real Axis just shifted with sigma) we can add the tails on both ends
IF(g%mode.EQ.3.AND.input%gf_anacont) THEN
!left tail
CALL g%get_gf(gmat,atoms,input,1,l,nType,ipm.EQ.2,spin=ispin,lp=lp,nTypep=nTypep)
CALL g%get(gmat,atoms,input,1,l,nType,ipm.EQ.2,spin=ispin,lp=lp,nTypep=nTypep)
ind1 = 0
DO m = -l, l
ind1 = ind1 + 1
......@@ -93,7 +93,7 @@ CONTAINS
ENDDO
CALL gmat%free()
!right tail
CALL g%get_gf(gmat,atoms,input,g%nz,l,nType,ipm.EQ.2,spin=ispin,lp=lp,nTypep=nTypep)
CALL g%get(gmat,atoms,input,g%nz,l,nType,ipm.EQ.2,spin=ispin,lp=lp,nTypep=nTypep)
ind1 = 0
DO m = -l, l
ind1 = ind1 + 1
......
......@@ -86,15 +86,15 @@ MODULE m_add_selfen
vmat%data_c(ns+1:2*ns,1:ns) = 0.0
ENDIF
IF(l_match_both_spins) THEN
CALL g%get_gf(gmat,atoms,input,iz,l,nType,ipm.EQ.2)
CALL g%get(gmat,atoms,input,iz,l,nType,ipm.EQ.2)
ELSE
CALL g%get_gf(gmat,atoms,input,iz,l,nType,ipm.EQ.2,spin=i_match)
CALL g%get(gmat,atoms,input,iz,l,nType,ipm.EQ.2,spin=i_match)
ENDIF
CALL add_pot(gmat,vmat,mu)
IF(l_match_both_spins) THEN
CALL gp%set_gf(gmat,atoms,input,iz,l,nType,ipm.EQ.2)
CALL gp%set(gmat,atoms,input,iz,l,nType,ipm.EQ.2)
ELSE
CALL gp%set_gf(gmat,atoms,input,iz,l,nType,ipm.EQ.2,spin=i_match)
CALL gp%set(gmat,atoms,input,iz,l,nType,ipm.EQ.2,spin=i_match)
ENDIF
CALL gmat%free()
ENDDO
......@@ -142,15 +142,15 @@ MODULE m_add_selfen
vmat%data_c(ns+1:2*ns,1:ns) = 0.0
ENDIF
IF(l_match_both_spins) THEN
CALL g%get_gf(gmat,atoms,input,iz,l,nType,ipm.EQ.2)
CALL g%get(gmat,atoms,input,iz,l,nType,ipm.EQ.2)
ELSE
CALL g%get_gf(gmat,atoms,input,iz,l,nType,ipm.EQ.2,spin=i_match)
CALL g%get(gmat,atoms,input,iz,l,nType,ipm.EQ.2,spin=i_match)
ENDIF
CALL add_pot(gmat,vmat,mu)
IF(l_match_both_spins) THEN
CALL gp%set_gf(gmat,atoms,input,iz,l,nType,ipm.EQ.2)
CALL gp%set(gmat,atoms,input,iz,l,nType,ipm.EQ.2)
ELSE
CALL gp%set_gf(gmat,atoms,input,iz,l,nType,ipm.EQ.2,spin=i_match)
CALL gp%set(gmat,atoms,input,iz,l,nType,ipm.EQ.2,spin=i_match)
ENDIF
CALL gmat%free()
ENDDO
......
......@@ -46,11 +46,11 @@ MODULE m_types_greensf
COMPLEX, ALLOCATABLE :: ud(:,:,:,:,:,:)
CONTAINS
PROCEDURE, PASS :: init => greensf_init
PROCEDURE, PASS :: init => greensf_init
PROCEDURE :: get => get_gf
PROCEDURE :: set => set_gf
PROCEDURE :: reset => reset_gf
PROCEDURE :: e_contour
PROCEDURE :: get_gf
PROCEDURE :: set_gf
PROCEDURE :: reset => reset_gf
END TYPE t_greensf
......
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