The model consist on a detailed morphology coded by a long list of pt3dadd calls.
Does anybody have an idea about how to connect the bipolar inputs to the corresponding dendrite?
For clarity, I'll provide the code I use in the simplified model.
Code: Select all
//Initial values for the bipolar inputs
xmin = -ldend
xmax = ldend
ymin = -ldend
ymax = ldend
y = ymin
fila = 0
n_is = 0
//Synaptic inputs location
while (y<=ymax) {
x = xmin + (1+(-1)^fila)*d_is/4
while (x<=xmax) {
a = r1.repick()
b = r2.repick()
if (sqrt((x+a*cos(b))^2+(y+a*sin(b))^2)<ldend) { //Circular array
xaux.append(x+a*cos(b))
yaux.append(y+a*sin(b))
n_is += 1
}
x += d_is
}
y += d_is*sin(PI/3)
fila += 1
}
//Synaptic input to dendrite connection
objref numdend, raux
numdend = new Vector(n_is) //Dendrite number
raux = new Vector(n_is) //Distance from the soma
for k = 0, n_is-1 {
angle = atan2(yaux.x[k], xaux.x[k]) //Determines the angle
if (angle < 0) {
angle += 2*PI
}
numdend.x[k] = int(angle*ndend/(2*PI)+0.5) //Determines the corresponding dendrite
if (numdend.x[k]==ndend) {
numdend.x[k]=0
}
raux.x[k]=sqrt(xaux.x[k]^2+yaux.x[k]^2)/ldend //Normalised distance to the soma
}
However, this approach doesn't work with this more sophisticated model.
I've thought about getting the pt3 coordinates somehow, but the documentation hasn't been helpful.
Any ideas will be really appreciated.
Thanks!