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