From 4b204f780a3169202600a35ab637d5fb2436ba39 Mon Sep 17 00:00:00 2001
From: Maximilian Schmidt <max.schmidt@fz-juelich.de>
Date: Mon, 11 Jun 2018 16:46:51 +0900
Subject: [PATCH] Fix index error in utils.py function

---
 figures/SchueckerSchmidt2017/utils.py | 6 +++---
 1 file changed, 3 insertions(+), 3 deletions(-)

diff --git a/figures/SchueckerSchmidt2017/utils.py b/figures/SchueckerSchmidt2017/utils.py
index ef8e03a..2fde31d 100644
--- a/figures/SchueckerSchmidt2017/utils.py
+++ b/figures/SchueckerSchmidt2017/utils.py
@@ -113,11 +113,11 @@ def compute_iteration(max_iter, fac_nu_ext_5E_list, theory_params, M_base, threa
         r = iteration_results(fac_nu_ext_5E_list, M_base, threads=threads)
         results = np.vstack((results, r))
         par_list = np.append(par_list, fac_nu_ext_5E_list)
-        i = np.argmax(np.diff(np.mean(r, axis=1)[:, -1]))
+        j = np.argmax(np.diff(np.mean(r, axis=1)[:, -1]))
         i += 1
-        fac_nu_ext_5E_list = np.arange(fac_nu_ext_5E_list[i],
+        fac_nu_ext_5E_list = np.arange(fac_nu_ext_5E_list[j],
                                        # to ensure that the array includes the last value, we add a small epsilon
-                                       fac_nu_ext_5E_list[i+1] + 1.e-10,
+                                       fac_nu_ext_5E_list[j+1] + 1.e-10,
                                        10**(-(i+2.)))
 
     ind = np.argsort(par_list)
-- 
GitLab