/usr/share/pyshared/pyNN/neuron/electrodes.py is in python-pynn 0.7.4-1.
This file is owned by root:root, with mode 0o644.
The actual contents of the file can be viewed below.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 | """
Current source classes for the neuron module.
Classes:
DCSource -- a single pulse of current of constant amplitude.
StepCurrentSource -- a step-wise time-varying current.
:copyright: Copyright 2006-2011 by the PyNN team, see AUTHORS.
:license: CeCILL, see LICENSE for details.
$Id: electrodes.py 957 2011-05-03 13:44:15Z apdavison $
"""
from neuron import h
class CurrentSource(object):
"""Base class for a source of current to be injected into a neuron."""
pass
class DCSource(CurrentSource):
"""Source producing a single pulse of current of constant amplitude."""
def __init__(self, amplitude=1.0, start=0.0, stop=None):
"""Construct the current source.
Arguments:
start -- onset time of pulse in ms
stop -- end of pulse in ms
amplitude -- pulse amplitude in nA
"""
self.amplitude = amplitude
self.start = start
self.stop = stop or 1e12
self._devices = []
def inject_into(self, cell_list):
"""Inject this current source into some cells."""
for id in cell_list:
if id.local:
iclamp = h.IClamp(0.5, sec=id._cell.source_section)
iclamp.delay = self.start
iclamp.dur = self.stop-self.start
iclamp.amp = self.amplitude
self._devices.append(iclamp)
class StepCurrentSource(CurrentSource):
"""A step-wise time-varying current source."""
def __init__(self, times, amplitudes):
"""Construct the current source.
Arguments:
times -- list/array of times at which the injected current changes.
amplitudes -- list/array of current amplitudes to be injected at the
times specified in `times`.
The injected current will be zero up until the first time in `times`. The
current will continue at the final value in `amplitudes` until the end
of the simulation.
"""
self.times = h.Vector(times)
self.amplitudes = h.Vector(amplitudes)
self._devices = []
def inject_into(self, cell_list):
"""Inject this current source into some cells."""
for id in cell_list:
if id.local:
if not id.celltype.injectable:
raise TypeError("Can't inject current into a spike source.")
iclamp = h.IClamp(0.5, sec=id._cell.source_section)
iclamp.delay = 0.0
iclamp.dur = 1e12
iclamp.amp = 0.0
self._devices.append(iclamp)
self.amplitudes.play(iclamp._ref_amp, self.times)
|