-
Notifications
You must be signed in to change notification settings - Fork 0
/
double_test.hoc
131 lines (109 loc) · 3.08 KB
/
double_test.hoc
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
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
/*==========================================
This file is used to test independence of two
synaptic objects
RTH Midnight Project by
Ruben A. Tikidji-Hamburyan <[email protected]>
May 2015
============================================*/
TPOST = 101 // time of postsynaptic spike
load_file("nrngui.hoc")
///// model cell /////
create soma1
access soma1
soma1.L = 20.
soma1.diam=20/PI
soma1.nseg=1
soma1.cm=1
insert hh
soma1.ek=-77
soma1.ena=50
objref edrv
edrv = new IClamp(0.5)
edrv.del = 0
edrv.amp = 0.024
edrv.dur = 500
objref syn1, netcon1, netcon2, netstim1, netstim2
//Synapse
syn1 = new exp2synGBstdp(0.5)
/*
//===== Use default parameters to reproduce classical STDP curveon Fig.2 top left ====//
*/
//Presynaptic input for STDP
netstim1 = new NetStim(0.5)
netstim1.interval = 10
netstim1.number = 0
netstim1.start = 10
netstim1.noise = 0
//Connections
netcon1 = new NetCon(netstim1,syn1)
netcon1.weight[0] = 1e-3 // synaptic conductance
netcon1.weight[1] = 0.5 // rho - synaptic efficacy
netcon1.weight[2] = 0.0 // calcium concentration
netcon1.delay = 0
netstim2 = new NetStim(0.5)
netstim2.interval = 10
netstim2.number = 0
netstim2.start = 10
netstim2.noise = 1
//Presynaptic input for STDP
netcon2 = new NetCon(netstim2,syn1)
netcon2.weight[0] = 1e-3 // synaptic conductance
netcon2.weight[1] = 0.5 // rho - synaptic efficacy
netcon2.weight[2] = 0.0 // calcium concentration
netcon2.delay = 0
create soma2
access soma2
soma2.L = 20.
soma2.diam=20/PI
soma2.nseg=1
soma2.cm=1
insert hh
soma2.ek=-77
soma2.ena=50
objref syn2, netcon3, netcon4, netstim3, netstim4
syn2 = new exp2synGBstdp(0.5)
//Presynaptic input for STDP
netstim3 = new NetStim(0.5)
netstim3.interval = 10
netstim3.number = 0
netstim3.start = 10
netstim3.noise = 0
//Connections
netcon3 = new NetCon(netstim3,syn2)
netcon3.weight[0] = 1e-3 // synaptic conductance
netcon3.weight[1] = 0.5 // rho - synaptic efficacy
netcon3.weight[2] = 0.0 // calcium concentration
netcon3.delay = 0
netstim4 = new NetStim(0.5)
netstim4.interval = 100
netstim4.number = 10
netstim4.start = 10
netstim4.noise = 1
//Presynaptic input for STDP
netcon4 = new NetCon(netstim4,syn2)
netcon4.weight[0] = 1e-3 // synaptic conductance
netcon4.weight[1] = 0.5 // rho - synaptic efficacy
netcon4.weight[2] = 0.0 // calcium concentration
netcon4.delay = 0
proc init() {
netcon1.weight[0] = 1e-3 // synaptic conductance
netcon1.weight[1] = 0.5 // rho - synaptic efficacy
netcon1.weight[2] = 0.0 // calcium concentration
netcon1.delay = 0
netcon2.weight[0] = 1e-3 // synaptic conductance
netcon2.weight[1] = 0.5 // rho - synaptic efficacy
netcon2.weight[2] = 0.0 // calcium concentration
netcon2.delay = 0
netcon3.weight[0] = 1e-3 // synaptic conductance
netcon3.weight[1] = 0.5 // rho - synaptic efficacy
netcon3.weight[2] = 0.0 // calcium concentration
netcon3.delay = 0
netcon4.weight[0] = 1e-3 // synaptic conductance
netcon4.weight[1] = 0.5 // rho - synaptic efficacy
netcon4.weight[2] = 0.0 // calcium concentration
netcon4.delay = 0
finitialize()
fcurrent()
frecord_init()
}
load_file("double_test.ses")