From c9e8da2990f50a41c7dc26232dea90521a822737 Mon Sep 17 00:00:00 2001
From: Natale Talarico <natale.talarico@utu.fi>
Date: Mon, 13 Jan 2020 17:36:05 +0200
Subject: [PATCH] Update welph.f90

---
 Libs_ness/welph.f90 | 19 ++++++++-----------
 1 file changed, 8 insertions(+), 11 deletions(-)

diff --git a/Libs_ness/welph.f90 b/Libs_ness/welph.f90
index f9b6bb5..040e47b 100644
--- a/Libs_ness/welph.f90
+++ b/Libs_ness/welph.f90
@@ -1,7 +1,7 @@
 !^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
 !	electron-phonon self-energy (Start)
 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
-subroutine elph(context,nps,npw,wm,vw,v,lnegl,lnegg,d0l,d0g,d0r,desc,lsl,lsg,lsr)
+subroutine elph(context,nps,npw,wm,vw,v,g,lnegl,lnegg,d0l,d0g,d0r,desc,lsl,lsg,lsr)
 	use mpi
 	implicit none
 !Useful numerical parameters
@@ -11,7 +11,7 @@ subroutine elph(context,nps,npw,wm,vw,v,lnegl,lnegg,d0l,d0g,d0r,desc,lsl,lsg,lsr
 	double complex, parameter :: zero = (0.0d0,0.0d0)		! Complex zero
 !Input variables
 	integer :: nps,npw
-	double precision :: wm
+	double precision :: wm, g
 	integer, dimension (:) :: vw,desc
 	double precision, dimension (:,:) :: v
 	double complex, dimension(:,:) :: lnegl
@@ -33,6 +33,7 @@ subroutine elph(context,nps,npw,wm,vw,v,lnegl,lnegg,d0l,d0g,d0r,desc,lsl,lsg,lsr
 	double complex, dimension(nps,nps) :: intg, tmp,pol
 	double complex, dimension(:,:), allocatable :: ltmpl, ltmpg, ltmp
 
+
 call blacs_gridinfo(context,nprow,npcol,myrow,mycol)
 nprocs=nprow*npcol
 
@@ -51,27 +52,23 @@ call fourierwt(context,nps,npw,wm,1,lnegl,desc,ltmpl)
 call fourierwt(context,nps,npw,wm,1,lnegg,desc,ltmpg)
 
 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!LESSER!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
-
+lsl=zero
 do lw=1, lldwc
 	liw=(lw-1)*nps+1
 !Self energy
-	lsl(:,liw:liw+nps-1)=lsl(:,liw:liw+nps-1)+ci*ltmpl(:,liw:liw+nps-1)*ld0l(:,liw:liw+nps-1)
+	lsl(:,liw:liw+nps-1)=ci*ltmpl(:,liw:liw+nps-1)*ld0l(:,liw:liw+nps-1)*(g**2)
 end do
 
 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!GREATER!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
-
+lsg=zero
 do lw=1, lldwc
 	liw=(lw-1)*nps+1
 !Self energy
-	lsg(:,liw:liw+nps-1)=lsg(:,liw:liw+nps-1)+ci*ltmpg(:,liw:liw+nps-1)*ld0g(:,liw:liw+nps-1)
+	lsg(:,liw:liw+nps-1)=ci*ltmpg(:,liw:liw+nps-1)*ld0g(:,liw:liw+nps-1)*(g**2)
 end do
 
-! WORK IN PROGRESS
-!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!RETARDED!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
-ltmp=zero
 
-tm=pi*(npw-1)/(2.d0*wm)
-dt=pi/wm
+!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!RETARDED!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
 
 do lw=1, lldwc
 	liw=(lw-1)*nps+1
-- 
GitLab