From 7a14add2f7a0da25de5bd4de75fda1b2cf2666ae Mon Sep 17 00:00:00 2001
From: Natale Talarico <natale.talarico@utu.fi>
Date: Fri, 10 Jan 2020 17:28:06 +0200
Subject: [PATCH] Update lesser and greater components.

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

diff --git a/Libs_ness/welph.f90 b/Libs_ness/welph.f90
index 752549a..f9b6bb5 100644
--- a/Libs_ness/welph.f90
+++ b/Libs_ness/welph.f90
@@ -54,24 +54,19 @@ call fourierwt(context,nps,npw,wm,1,lnegg,desc,ltmpg)
 
 do lw=1, lldwc
 	liw=(lw-1)*nps+1
-!Polarization diagram
-	pol(:,:)=ltmpl(:,liw:liw+nps-1)*conjg(ltmpg(:,liw:liw+nps-1))
-	tmp=matmul(v,matmul(pol,v))
 !Self energy
-	lsl(:,liw:liw+nps-1)=lsl(:,liw:liw+nps-1)-ltmpl(:,liw:liw+nps-1)*tmp
+	lsl(:,liw:liw+nps-1)=lsl(:,liw:liw+nps-1)+ci*ltmpl(:,liw:liw+nps-1)*ld0l(:,liw:liw+nps-1)
 end do
 
 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!GREATER!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
 
 do lw=1, lldwc
 	liw=(lw-1)*nps+1
-!Polarization diagram
-	pol(:,:)=ltmpg(:,liw:liw+nps-1)*conjg(ltmpl(:,liw:liw+nps-1))
-	tmp=matmul(v,matmul(pol,v))
 !Self energy
-	lsg(:,liw:liw+nps-1)=lsg(:,liw:liw+nps-1)-ltmpg(:,liw:liw+nps-1)*tmp
+	lsg(:,liw:liw+nps-1)=lsg(:,liw:liw+nps-1)+ci*ltmpg(:,liw:liw+nps-1)*ld0g(:,liw:liw+nps-1)
 end do
 
+! WORK IN PROGRESS
 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!RETARDED!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
 ltmp=zero
 
@@ -89,6 +84,8 @@ do lw=1, lldwc
 	end if
 end do
 
+
+!  self energy in frequncy domain 
 ltmpl=zero
 call fourierwt(context,nps,npw,wm,-1,ltmp,desc,ltmpl)
 lsr = lsr + ltmpl
-- 
GitLab