Skip to content
Snippets Groups Projects
Commit 4b204f78 authored by Maximilian Schmidt's avatar Maximilian Schmidt
Browse files

Fix index error in utils.py function

parent 46143ca2
No related branches found
No related tags found
1 merge request!1Add all necessary files for the multi-area model
......@@ -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)
......
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment