Changeset 407

Timestamp:

07/03/08 03:17:50 (5 months ago)

Author:

pierre

Message:

Implement the Gutig rule and the van Rossum rule of plasticity in nest2, and change the old naming. Fix also the meanSpikeCount() function to have directly the spike count from the recorder, and not from the file

Files:

• trunk/src/common.py (modified) (2 diffs)
• trunk/src/nest2/__init__.py (modified) (2 diffs)
• trunk/src/nest2/synapses.py (modified) (2 diffs)

trunk/src/common.py

@@ -1513 +1513 @@
-    """
-    default_parameters = {
-:  2
-: 20
-
-
+  : 
+  : 1
+ 
+,
-    }
-    
@@ -1530 +1530 @@
-    """
-
-
+
+
+
+
+
+
+
+
-        pass
-
-    
-class GutigWeightDependence(STDPWeightDependence):
-
+
+
+
+
+
+
+
+
+
+
+
+
-
-# Not yet implemented for any module

trunk/src/nest2/__init__.py

@@ -765 +765 @@
-            tmp_list = record_from
-        else:
-
+
+
-            tmp_list = rng.permutation(numpy.reshape(self.cell, (self.cell.size,)))[0:n_rec]
-
@@ -863 +864 @@
-        Returns the mean number of spikes per neuron.
-        """
-len(self.recorders['spikes'].get(gather))
+nest.GetStatus(self.recorders['spikes']._device,'n_events')[0]
-        n_rec = len(self.recorders['spikes'].recorded)
-        return float(n_spikes)/n_rec

trunk/src/nest2/synapses.py

@@ -82 +82 @@
-    """
-    translations = common.build_translations(
-w
+W
-        ('w_min',     'w_min_always_zero_in_NEST'),
-aLTP
+lambda
-        ('A_minus',   'alpha', 'A_minus/A_plus', 'alpha*lambda'),
-    )
@@ -98 +98 @@
-        self.parameters['mu_plus'] = 1.0
-        self.parameters['mu_minus'] = 1.0
+
+class AdditivePotentiationMultiplicativeDepression(common.AdditivePotentiationMultiplicativeDepression):
+    """
+    The amplitude of the weight change depends on the current weight.
+    For depression, Dw propto w-w_min
+    For potentiation, Dw propto w_max-w
+    """
+    translations = common.build_translations(
+        ('w_max',     'Wmax',  1000.0), # unit conversion
+        ('w_min',     'w_min_always_zero_in_NEST'),
+        ('A_plus',    'lambda'),
+        ('A_minus',   'alpha', 'A_minus/A_plus', 'alpha*lambda'),
+    )
+    possible_models = set(['stdp_synapse',]) #'stdp_synapse_hom'
+        
+    def __init__(self, w_min=0.0, w_max=1.0, A_plus=0.01, A_minus=0.01):
+        if w_min != 0:
+            raise Exception("Non-zero minimum weight is not supported by NEST.")
+        common.AdditivePotentiationMultiplicativeDepression.__init__(self, w_min, w_max, A_plus, A_minus)
+        parameters = locals()
+        parameters.pop('self') 
+        self.parameters = self.translate(parameters)
+        self.parameters['mu_plus'] = 0.0
+        self.parameters['mu_minus'] = 1.0
+
+
+class GutigWeightDependence(common.GutigWeightDependence):
+    """
+    The amplitude of the weight change depends on the current weight.
+    For depression, Dw propto w-w_min
+    For potentiation, Dw propto w_max-w
+    """
+    translations = common.build_translations(
+        ('w_max',     'Wmax',  1000.0), # unit conversion
+        ('w_min',     'w_min_always_zero_in_NEST'),
+        ('A_plus',    'lambda'),
+        ('A_minus',   'alpha', 'A_minus/A_plus', 'alpha*lambda'),
+        ('mu_plus',   'mu_plus'),
+        ('mu_minus',  'mu_minus'),
+    )
+    possible_models = set(['stdp_synapse',]) #'stdp_synapse_hom'
+        
+    def __init__(self, w_min=0.0, w_max=1.0, A_plus=0.01, A_minus=0.01,mu_plus=0.5,mu_minus=0.5):
+        if w_min != 0:
+            raise Exception("Non-zero minimum weight is not supported by NEST.")
+        common.AdditivePotentiationMultiplicativeDepression.__init__(self, w_min, w_max, A_plus, A_minus)
+        parameters = locals()
+        parameters.pop('self') 
+        self.parameters = self.translate(parameters)
-
-