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
|
: Axial current and derivative with respect to x
: INCLUDE in density and point process
NEURON {
THREADSAFE
RANGE ri, ia, im
POINTER pv, pia, pim
}
UNITS {
(mv) = (millivolt)
(nA) = (nanoamp)
(MOhm) = (megohm)
}
PARAMETER {
ri (MOhm) : center of compartment to center of parent compartment
}
ASSIGNED {
v (mV)
ia (nA)
im (nA) : im - i_membrane_ = electode_currents
pv (mV) : parent v
pia (nA) : parent ia
pim (nA) : parent im
}
AFTER INITIAL {
im = 0 : pia not ready yet
if (ri > 0) {
ia = (pv - v) / ri
foo()
}
}
AFTER SOLVE {
if (ri > 0) {
ia = (pv - v) / ri
im = ia : contribution from parent side
}
}
BEFORE STEP {
if (ri > 0) {
pim = pim - ia : child contributions
}
}
PROCEDURE foo() { : neuronsimulator/nrn#1782 is fixed
if (ri > 0) {
: just for compile test. Although this overwrites im from BEFORE STEP
: the recording of im takes place before the AFTER INITIAL
pim = pim + 1
}
}
|
