From 80b59ff6696a857928c1315a36bac2cb206384b5 Mon Sep 17 00:00:00 2001
From: Johanna <johanna@zijderveld.de>
Date: Thu, 11 Apr 2024 22:47:08 +0200
Subject: [PATCH] move where the densitymatrixTb is created very slightly
---
codes/model.py | 7 +++----
1 file changed, 3 insertions(+), 4 deletions(-)
diff --git a/codes/model.py b/codes/model.py
index cad2a5b..9606a2e 100644
--- a/codes/model.py
+++ b/codes/model.py
@@ -34,14 +34,13 @@ class Model:
def calculateEF(self):
self.EF = fermiOnGrid(self.kham, self.filling)
- def makeDensityMatrix(self, mf_model, nK=200):
+ def makeDensityMatrixTb(self, mf_model, nK=200):
self.kham = tb2kham(addTb(self.h_0, mf_model), nK=nK, ndim=self._ndim)
self.calculateEF()
- return densityMatrix(self.kham, self.EF)
+ return ifftn2tb(kdens2tbFFT(densityMatrix(self.kham, self.EF), self._ndim))
def mfield(self, mf_model, nK=200):
- densityMatrix = self.makeDensityMatrix(mf_model, nK=nK)
- densityMatrixTb = ifftn2tb(kdens2tbFFT(densityMatrix, self._ndim))
+ densityMatrixTb = self.makeDensityMatrixTb(mf_model, nK=nK)
return addTb(
meanField(densityMatrixTb, self.h_int, n=self._ndim),
{self._localKey: -self.EF * np.eye(self._size)},
--
GitLab