def mkgap(pc, sec, gid, secpos, sgid, dgid, w, gjlist):
myrank = int(pc.id())
seg = sec(secpos)
gj = h.ggap(seg)
gj.g = w
pc.source_var(seg._ref_v, sgid, sec=sec)
pc.target_var(gj, gj._ref_vgap, dgid)
if myrank == 0:
print('mkgap: gid %i: sec=%s sgid=%i dgid=%i w=%f' % (gid, str(sec), sgid, dgid, w))
gjlist.append(gj)
return gj
def mkgap(env, cell, gid, secpos, secidx, sgid, dgid, w):
"""
Create gap junctions
:param pc:
:param gjlist:
:param gid:
:param secidx:
:param sgid:
:param dgid:
:param w:
:return:
"""
sec = list(cell.sections)[secidx]
seg = sec(secpos)
gj = h.ggap(seg)
gj.g = w
env.pc.source_var(seg._ref_v, sgid, sec=sec)
env.pc.target_var(gj, gj._ref_vgap, dgid)
env.gjlist.append(gj)
return gj