From f0835dd902980fb14395f772177a90cc061b30e3 Mon Sep 17 00:00:00 2001
From: Joseph Weston <joseph@weston.cloud>
Date: Fri, 14 Dec 2018 18:56:30 +0100
Subject: [PATCH] use shift-invert in other instances of 'sla.eigsh' in
tutorial
---
doc/source/code/figure/closed_system.py.diff | 4 ++--
1 file changed, 2 insertions(+), 2 deletions(-)
diff --git a/doc/source/code/figure/closed_system.py.diff b/doc/source/code/figure/closed_system.py.diff
index 1eb0080c..6b5d2b83 100644
--- a/doc/source/code/figure/closed_system.py.diff
+++ b/doc/source/code/figure/closed_system.py.diff
@@ -105,7 +105,7 @@
+
# Calculate the wave functions in the system.
ham_mat = syst.hamiltonian_submatrix(sparse=True, args=[B])
- evals, evecs = sorted_eigs(sla.eigsh(ham_mat, k=20, which='SM'))
+ evals, evecs = sorted_eigs(sla.eigsh(ham_mat, k=20, sigma=0))
# Plot the probability density of the 10th eigenmode.
- kwant.plotter.map(syst, np.abs(evecs[:, 9])**2,
@@ -124,7 +124,7 @@
+
# Calculate the wave functions in the system.
ham_mat = syst.hamiltonian_submatrix(sparse=True, args=[B])
- evals, evecs = sorted_eigs(sla.eigsh(ham_mat, k=20, which='SM'))
+ evals, evecs = sorted_eigs(sla.eigsh(ham_mat, k=20, sigma=0))
# Calculate and plot the local current of the 10th eigenmode.
J = kwant.operator.Current(syst)
--
GitLab