Changeset 5

Timestamp:

05/02/07 10:41:50 (5 years ago)

Author:

JensKremkow

Message:

 

Files:

• trunk/sandbox.py (modified) (2 diffs)

trunk/sandbox.py

@@ -49 +49 @@
     return tmpdict
 
+
+
+def get_connectivity(params):
+    """
+    
+    """
+    a = params['a']
+    radius = params['radius']
+    radius_normalized = radius/a # when a is  
+    population = params['population']
+
+    center = (int(population.dim[0]/2.),int(population.dim[1]/2.))
+    offset = (int(round(population.dim[0]*radius_normalized)),int(round(population.dim[1]*radius_normalized)))
+    
+    targets={}
+    targets_gid={}
+    for n in range(center[0]-offset[0],center[0]+offset[0]+1):
+        for m in range(center[1]-offset[1],center[1]+offset[1]+1):
+            #print 'n: ',n, ' m: ',m
+            gid = population[n,m]
+            targets_tmp = pynest.getDict([gid])[0]['targets']
+            targets_n_m = []
+            for tgid in targets_tmp:
+                targets_n_m.append(population.locate(tgid))
+            
+            targets[(n,m)]=targets_n_m
+            targets_gid[gid]=targets_tmp
+    return targets, targets_n_m
+    
+    
+    
+
+
+
 """
 SpikeLists functions
 
 """
-
 def tmpfile2spikelist(filename,dt):
     """
@@ -82 +115 @@
     return spike_list
 
+# we need that for the hardware stuff, since printspike does not yet exist in the low-level api. and therefore the spiketrains are in the old format.
+def tmpfile2spikelist_old(filename):
+    """ returns a spike list from the tmp file saved by PyNN-NEST
+    the spike list is according to latest standards a list of tuples (relative time
+    of spike, neuron_id); both are integers the time is obtained"""
+    from numpy import *
+    from pylab import load
+    import os
+    
+
+    spike_list = []
+
+    if os.path.getsize(filename) > 0: # check that the size is not zero which raises an error in load
+        #DATA = load(filename) # TODO: replace with scipy function (fragile if spike list is empty)
+        neuron_id, spike_time = load(filename, usecols=(1,2), unpack=True)
+
+        spike_list=[]
+
+        if size(neuron_id) >1:
+            # TODO : does not work with one spike. therefore this stupid if
+            neuron_id = array([int(k) for k in neuron_id])
+            spike_time = array([int(k) for k in spike_time])
+            rel_spiketime = concatenate(([spike_time[0]] , spike_time[1:]-spike_time[:-1]))
+
+            for index in range(len(neuron_id)):
+                spike_list.append( (rel_spiketime[index], neuron_id[index]) )
+
+        else:
+            spike_list.append( (spike_time, neuron_id) )
 
 
+    return spike_list