Fix external siegert neurons

multiarea_model/default_params.py

@@ -281,7 +281,7 @@ theory_params = {'neuron_params': neuron_params,
                  # The simulation time of the mean-field theory integration
                  'T': 50.,
                  # The time step of the mean-field theory integration
-                 'dt': 0.1,
+                 'dt': 0.01,
                  # Time interval for recording the trajectory of the mean-field calcuation
                  # If None, then the interval is set to dt
                  'rec_interval': None}

multiarea_model/theory.py

@@ -78,13 +78,14 @@ class Theory:
                               'use_wfr': False,
                               'print_time': False,
                               'overwrite_files': True})
+
+        nest.SetDefaults('siegert_neuron', self.NP)
         # create neurons for external drive
         drive = nest.Create(
             'siegert_neuron', 1, params={'rate': rate_ext, 'mean': rate_ext})
 
         # create neurons representing populations
-        neurons = nest.Create(
-            'siegert_neuron', dim, params=self.NP)
+        neurons = nest.Create('siegert_neuron', dim)
         # external drive
         syn_dict = {'drift_factor': tau * np.array([K[:, -1] * J[:, -1]]).transpose(),
                     'diffusion_factor': tau * np.array([K[:, -1] * J[:, -1]**2]).transpose(),
@@ -132,7 +133,7 @@ class Theory:
         interval = self.params['rec_interval']
         if interval is None:
             interval = dt
-            
+
         multimeter = nest.Create('multimeter', params={'record_from': ['rate'],
                                                        'interval': interval,
                                                        'to_screen': False,