From 7f532132fd82347ed5675d5a8af9b91678a0e40a Mon Sep 17 00:00:00 2001
From: Anton Akhmerov <anton.akhmerov@gmail.com>
Date: Wed, 22 May 2013 16:24:27 -0400
Subject: [PATCH] use reality of h and t in modes
---
kwant/physics/leads.py | 5 +++++
1 file changed, 5 insertions(+)
diff --git a/kwant/physics/leads.py b/kwant/physics/leads.py
index d417fc4c..f2516028 100644
--- a/kwant/physics/leads.py
+++ b/kwant/physics/leads.py
@@ -71,6 +71,11 @@ def setup_linsys(h_onslice, h_hop, tol=1e6, algorithm=None):
# False if h_hop is purely imaginary
raise ValueError("Inter-slice hopping is exactly zero.")
+ # If both h and t are real, it may be possible to use the real eigenproblem.
+ if (not np.any(h_hop.imag)) and (not np.any(h_onslice.imag)):
+ h_hop = h_hop.real
+ h_onslice = h_onslice.real
+
eps = np.finfo(np.common_type(h_onslice, h_hop)).eps * tol
# First check if the hopping matrix has singular values close to 0.
--
GitLab