From 1a6dfdb013919262b85972377e9ddb24ffcf7869 Mon Sep 17 00:00:00 2001
From: Maximilian Schmidt <max.schmidt@fz-juelich.de>
Date: Tue, 29 May 2018 17:31:14 +0900
Subject: [PATCH] Fix model initialization in Fig9 script
---
figures/SchueckerSchmidt2017/Fig9_simulation.eps | 2 +-
figures/SchueckerSchmidt2017/Fig9_simulation.py | 7 ++-----
2 files changed, 3 insertions(+), 6 deletions(-)
diff --git a/figures/SchueckerSchmidt2017/Fig9_simulation.eps b/figures/SchueckerSchmidt2017/Fig9_simulation.eps
index 6ae2d4d..eed822c 100644
--- a/figures/SchueckerSchmidt2017/Fig9_simulation.eps
+++ b/figures/SchueckerSchmidt2017/Fig9_simulation.eps
@@ -3,7 +3,7 @@
%%HiResBoundingBox: 0.000000 0.000000 374.400000 231.391925
%%Creator: GPL Ghostscript 922 (ps2write)
%%LanguageLevel: 2
-%%CreationDate: D:20180529165002+09'00'
+%%CreationDate: D:20180529173051+09'00'
%%EndComments
%%BeginProlog
save
diff --git a/figures/SchueckerSchmidt2017/Fig9_simulation.py b/figures/SchueckerSchmidt2017/Fig9_simulation.py
index b907d88..38245aa 100644
--- a/figures/SchueckerSchmidt2017/Fig9_simulation.py
+++ b/figures/SchueckerSchmidt2017/Fig9_simulation.py
@@ -50,12 +50,13 @@ network_params = {'N_scaling': 1.,
ana_spec = {'load_areas': []}
M = MultiAreaModel(network_params,
analysis=True,
+ ana_spec=ana_spec,
simulation=True,
sim_spec=sim_params,
theory=True,
theory_spec=theory_params)
p, r_stab = M.theory.integrate_siegert()
-M.analysis.create_pop_rates()
+M.analysis.create_pop_rates(areas=M.area_list)
"""
@@ -83,8 +84,6 @@ for i, area in enumerate(area_list):
matrix = np.transpose(matrix)
rate_matrix_plot(panel_factory.figure, ax, matrix, position='left')
-print("Meanfield: Area-averaged population rates: ", np.mean(matrix, axis=1))
-
"""
Panel B
@@ -111,8 +110,6 @@ for i, area in enumerate(M.area_list):
matrix_sim = np.transpose(matrix_sim)
rate_matrix_plot(panel_factory.figure, ax, matrix, position='right')
-print("Area-averaged population rates: ", np.mean(matrix, axis=0))
-
"""
Panel C
--
GitLab