// nFit
// (c) Charles CH Cohen, 2014-present
// this software is released to the public under a Creative Commons Attribution-NonCommercial-NoDerivatives 4.0
// International license (CC BY-NC-ND 4.0, in English).
// for any questions, please email c.cohen@gmx.com
// -------------------------------Variables----------------------------------------
objref recsecxvec, recsecyvec
// --------------------------------------------------------------------------------
// -----------------------------writecontrol---------------------------------------
proc writecontrol() {local chk
if (numarg() == 1 && argtype(1) == 2) {
if (!strcmp($s1, "help")) {
print "Enter:"
print "1) mode"
print "output: default runcontrol panel in given mode"
stop
}
}
chk = chkcwd(cwd)
if (chk) cwd = getcwd()
chdir(ses)
writefile.wopen(runcontrolfilestr)
writefile.printf("%s", "objectvar save_window_, rvp_\nobjectvar scene_vector_[5]\nobjectvar ocbox_, ocbox_list_, scene_, scene_list_\n{ocbox_list_ = new List() scene_list_ = new List()}")
writefile.printf("%s", "\n{\nxpanel(\"RunControl\", 0)\nv_init = ")
getRMP()
writefile.printf(vform, v_init)
writefile.printf("%s", "\nxvalue(\"Init\",\"v_init\", 1,\"stdinit()\", 1, 1 )\nxbutton(\"Init & Run\",\"run()\")")
writefile.printf("%s", "\nxbutton(\"Stop\",\"stoprun=1\")\nrunStopAt = ")
writefile.printf(tform, gettstop())
writefile.printf("%s", "\nxvalue(\"Continue til\",\"runStopAt\", 1,\"{continuerun(runStopAt) stoprun=1}\", 1, 1 )\nrunStopIn = 0.1")
writefile.printf("%s", "\nxvalue(\"Continue for\",\"runStopIn\", 1,\"{continuerun(t + runStopIn) stoprun=1}\", 1, 1 )\nxbutton(\"Single Step\",\"steprun()\")\nt = 0")
writefile.printf("%s", "\nxvalue(\"t\",\"t\", 2 )\ntstop = ")
writefile.printf(tform, gettstop())
writefile.printf("%s", "\nxvalue(\"Tstop\",\"tstop\", 1,\"tstop_changed()\", 0, 1 )\ndt = ")
setdt()
dt = setnsig(dt, nsigt)
writefile.printf(tform, dt)
writefile.printf("%s", "\nxvalue(\"dt\",\"dt\", 1,\"setdt()\", 0, 1 )\nsteps_per_ms = ")
writefile.printf(vform, getsteps())
writefile.printf("%s", "\nxvalue(\"Points plotted/ms\",\"steps_per_ms\", 1,\"setdt()\", 0, 1 )\nscreen_update_invl = 0.05\nxvalue(\"Scrn update invl\",\"screen_update_invl\", 1,\"\", 0, 1 )\nrealtime = 0")
writefile.printf("%s", "\nxvalue(\"Real Time\",\"realtime\", 0,\"\", 0, 1 )\nxpanel(54,110)\n}")
writefile.close()
chdir(cwd)
if (chk) cwd = ""
}
// --------------------------------------------------------------------------------
// --------------------------------------------------------------------------------
proc writevtplot() {local chk, ev, rec, iamp
if (numarg() == 1 && argtype(1) == 2) {
if (!strcmp($s1, "help")) {
print "Enter:"
print "1) optional: iiamp associated with desired vtplot. Default: current iiamp"
print "output: default vtplot according to given mode and iiamp saved in vtplot-*mode*.ses"
stop
}
}
chk = chkcwd(cwd)
if (chk) cwd = getcwd()
if (numarg() == 1 && argtype(1) == 0) {
iamp = $1
} else {
iamp = getvarp(ses, iampfilestr)
}
chdir(ses)
writefile.wopen(vtplotfilestr)
writefile.printf("%s", "objectvar save_window_, rvp_\nobjectvar scene_vector_[5]\nobjectvar ocbox_, ocbox_list_, scene_, scene_list_\n{ocbox_list_ = new List() scene_list_ = new List()}")
tstop = gettstop()
crvvarlist(recloclist, vvarlist)
for rec = 0, nrec-1 {
writefile.printf("%s", "\n{")
writefile.printf("%s", "\nsave_window_ = new Graph(0)\nsave_window_.size(-2, ")
writefile.printf(tform, tstop)
writefile.printf("%s", ", ")
if (mode == 0) {
ev = getev(iamp, recseclist.o(rec).label, getsubloc(recloclist.o(rec).label))
if (ev < 0) {
ev -= 1
writefile.printf(vform, ev)
writefile.printf("%s", ", -0.15")
} else {
ev += 1
writefile.printf("%s", "-0.15, ")
writefile.printf(vform, ev)
}
} else if (mode == 1) {
writefile.printf("%s", "-80, 40")
}
writefile.printf("%s", ")")
writefile.printf("%s", "\nscene_vector_[")
writefile.printf("%d", 2+rec)
writefile.printf("%s", "] = save_window_")
writefile.printf("%s", "\n{save_window_.view(-2, ")
if (mode == 0) {
if (ev < 0) {
writefile.printf(vform, ev)
writefile.printf("%s", ", ")
writefile.printf(tform, tstop+2)
writefile.printf("%s", ",")
writefile.printf(vform, abs(ev-0.05))
} else {
writefile.printf("%s", "-0.15, ")
writefile.printf(tform, tstop+2)
writefile.printf("%s", ", ")
writefile.printf(vform, ev)
}
} else if (mode == 1) {
writefile.printf("%s", "-80, ")
writefile.printf(tform, tstop+2)
writefile.printf("%s", ", ")
writefile.printf("%s", "120")
}
writefile.printf("%s", ", 341, ")
writefile.printf("%g", 110+rec*(374-115))
writefile.printf("%s", ", 302.4, 193.6)}")
writefile.printf("%s", "\ngraphList[0].append(save_window_)\nsave_window_.save_name(\"graphList[0].\")\nsave_window_.addexpr(\"")
writefile.printf("%s", vvarlist.o(rec).label)
writefile.printf("%s", "\", 1, 1, 0.660318, ")
if (mode == 0) {
if (ev < 0) {
writefile.printf("%s", "0.141322")
} else {
writefile.printf("%s", "0.966103")
}
} else if (mode == 1) {
writefile.printf("%s", "0.141322")
}
writefile.printf("%s", ", 2)\n}")
}
writefile.printf("%s", "\nobjectvar scene_vector_[1]\n{doNotify()}\n")
writefile.close()
chdir(cwd)
if (chk) cwd = ""
}
// --------------------------------------------------------------------------------
// ------------------------------writevxplot---------------------------------------
proc writevxplot() {local chk, ev, dist, distsoma, ystart, ext
if (numarg() == 1 && argtype(1) == 2) {
if (!strcmp($s1, "help")) {
print "Enter:"
print "1) optional: if passive, enter the amplitude of choice (nA), to help set the voltage axis."
print "output: default vxplot according to mode (vxplot-*mode*.ses)"
stop
}
}
chk = chkcwd(cwd)
if (chk) cwd = getcwd()
if (mode == 0) {
if (numarg() == 1 && argtype(2) == 0) {
ev = getmaxev($1)
} else {
ev = getvarp(ses, iampfilestr)
}
} else if (mode == 1) {
ev = 40
}
chdir(ses)
writefile.wopen(vxplotfilestr)
writefile.printf("%s", "objectvar save_window_, rvp_")
writefile.printf("%s", "\nobjectvar scene_vector_[5]")
writefile.printf("%s", "\nobjectvar ocbox_, ocbox_list_, scene_, scene_list_")
writefile.printf("%s", "\n{ocbox_list_ = new List() scene_list_ = new List()}")
writefile.printf("%s", "\n{")
writefile.printf("%s", "\nsave_window_ = new Graph(0)")
writefile.printf("%s", "\nsave_window_.size(")
dist = getdist(recseclist.o(0).label, 0.5, recseclist.o(recseclist.count-1).label, 0.5) + 20
distsoma = getdist("soma", 0.5, recseclist.o(0).label, getsubloc(recloclist.o(0).label))
setmaxform()
writefile.printf(maxform, distsoma)
writefile.printf("%s", ", ")
writefile.printf(maxform, dist)
writefile.printf("%s", ", ")
if (mode == 0) {
if (ev < 0) {
ev -= 1
writefile.printf(maxform, ev)
writefile.printf("%s", ", ")
ystart = 1
writefile.printf("%g", ystart)
} else {
ystart = -1
writefile.printf("%g", ystart)
writefile.printf("%s", ", ")
ev += 1
writefile.printf(maxform, ev)
}
} else if (mode == 1) {
ystart = -80
writefile.printf("%g", ystart)
writefile.printf("%s", ", ")
writefile.printf(maxform, ev)
}
writefile.printf("%s", ")")
writefile.printf("%s", "\nscene_vector_[4] = save_window_")
writefile.printf("%s", "\n{save_window_.view(")
writefile.printf(maxform, distsoma)
writefile.printf("%s", ", ")
if (ev < 0) {
writefile.printf(maxform, ev)
} else {
writefile.printf("%g", ystart)
}
writefile.printf("%s", ", ")
writefile.printf(maxform, dist)
writefile.printf("%s", ", ")
writefile.printf(maxform, abs(ev) + abs(ystart))
writefile.printf("%s", ", 666, 110, 586.56, 233.92)}")
writefile.printf("%s", "\nflush_list.append(save_window_)")
writefile.printf("%s", "\nsave_window_.save_name(\"flush_list.\")")
writefile.printf("%s", "\nobjectvar rvp_")
writefile.printf("%s", "\nrvp_ = new RangeVarPlot(\"v\")")
writefile.printf("\n%s", recseclist.o(0).label)
writefile.printf("%s", " rvp_.begin(0.5)")
writefile.printf("\n%s", recseclist.o(recseclist.count-1).label)
writefile.printf("%s", " rvp_.end(0.5)")
writefile.printf("%s", "\nrvp_.origin(0)")
writefile.printf("%s", "\nsave_window_.addobject(rvp_, 2, 1, 0.146864, 0.629207)")
ext = getaxontype(0)
if (ext > 2) {
writefile.printf("%s", "\nobjectvar rvp_")
writefile.printf("%s", "\nrvp_ = new RangeVarPlot(\"vext\")")
writefile.printf("\n%s", recseclist.o(0).label)
writefile.printf("%s", " rvp_.begin(0.5)")
writefile.printf("\n%s", recseclist.o(recseclist.count-1).label)
writefile.printf("%s", " rvp_.end(0.5)")
writefile.printf("%s", "\nrvp_.origin(0)")
writefile.printf("%s", "\nsave_window_.addobject(rvp_, 3, 1, 0.146864, 0.594483)")
writefile.printf("%s", "\nobjectvar rvp_")
writefile.printf("%s", "\nrvp_ = new RangeVarPlot(\"v($1)+vext($1)\")")
writefile.printf("\n%s", recseclist.o(0).label)
writefile.printf("%s", " rvp_.begin(0.5)")
writefile.printf("\n%s", recseclist.o(recseclist.count-1).label)
writefile.printf("%s", " rvp_.end(0.5)")
writefile.printf("%s", "\nrvp_.origin(0)")
writefile.printf("%s", "\nsave_window_.addobject(rvp_, 9, 1, 0.146864, 0.553377)")
}
writefile.printf("%s", "\n}")
writefile.printf("%s", "\nobjectvar scene_vector_[1]")
writefile.printf("%s", "\n{doNotify()}")
writefile.close()
chdir(cwd)
if (chk) cwd = ""
}
// --------------------------------------------------------------------------------
// --------------------------------------------------------------------------------
proc writeshape() {local chk, k, mleft, mright, mbottom, mtop, sleft, screentop, swidth, sheight, maxx, minx, maxy, miny, minmax
if (numarg() == 1 && argtype(1) == 2) {
if (!strcmp($s1, "help")) {
print "Enter:"
print "1) optional: if passive, enter the amplitude of choice (nA), to help set voltage colour range."
print "output: default shapeplot for Vm saved (shapeplot-*mode*.ses)"
stop
}
}
chk = chkcwd(cwd)
if (chk) cwd = getcwd()
chdir(ses)
writefile.wopen(shapeplotfilestr)
writefile.printf("%s", "\nobjectvar save_window_, rvp_")
writefile.printf("%s", "\nobjectvar scene_vector_[7]")
writefile.printf("%s", "\nobjectvar ocbox_, ocbox_list_, scene_, scene_list_")
writefile.printf("%s", "\n{ocbox_list_ = new List() scene_list_ = new List()}")
writefile.printf("%s", "\n{")
writefile.printf("%s", "\nsave_window_ = new PlotShape(0)")
writefile.printf("%s", "\nsave_window_.size(")
recsecxvec = new Vector()
for k = 0, recseclist.count-1 recsecxvec.append(getsecx(recseclist.o(k).label))
recsecyvec = new Vector()
for k = 0, recseclist.count-1 recsecyvec.append(getsecy(recseclist.o(k).label))
maxx = recsecxvec.max
minx = recsecxvec.min
maxy = recsecyvec.max
miny = recsecyvec.min
minmax = abs(maxx) + abs(minx) + abs(maxy) + abs(miny)
sleft = 666
screentop = 413
swidth = 586.56
sheight = 235.84
if (minmax > 300) {
if ((abs(maxy) + abs(miny)) > (abs(maxx) + abs(minx))) {
mwidth = (swidth/sheight) * (abs(maxy)+abs(miny))
mleft = -0.5*mwidth
mright = 0.5*mwidth
mheight = (sheight/swidth)*mwidth
if (miny < 0 && maxy >= 0) {
mtop = ((maxy/(abs(miny)+maxy))+0.05)*mheight
mbottom = -((abs(miny)/(abs(miny)+maxy))-0.05)*mheight
} else if (miny >= 0 && maxy > 0) {
mtop = ((maxy/(abs(miny)+maxy))+0.05)*mheight
mbottom = ((abs(miny)/(abs(miny)+maxy))-0.05)*mheight
} else if (miny < 0 && maxy < 0) {
mtop = ((maxy/(abs(miny)+maxy))+0.05)*mheight
mbottom = -((abs(miny)/(abs(miny)+maxy))-0.05)*mheight
}
} else {
mwidth = (swidth/sheight) * (abs(maxx)+abs(minx))
mheight = (sheight/swidth)*mwidth
mtop = 0.5*mheight
mbottom = -0.5*mheight
if (minx < 0 && maxx >= 0) {
mright = ((maxx/(abs(minx)+maxx))+0.05)*mheight
mleft = -((abs(minx)/(abs(minx)+maxx))-0.05)*mheight
} else if (minx >= 0 && maxx > 0) {
mright = ((maxx/(abs(minx)+maxx))+0.05)*mheight
mleft = ((abs(minx)/(abs(minx)+maxx))-0.05)*mheight
} else if (minx < 0 && maxx < 0) {
mright = ((maxx/(abs(minx)+maxx))+0.05)*mheight
mleft = -((abs(minx)/(abs(minx)+maxx))-0.05)*mheight
}
}
} else {
mleft = -100
mright = 100
mbottom = -75
mtop = 25
mwidth = 200
mheight = 100
}
writefile.printf(maxform, mleft)
writefile.printf("%s", ", ")
writefile.printf(maxform, mright)
writefile.printf("%s", ", ")
writefile.printf(maxform, mbottom)
writefile.printf("%s", ", ")
writefile.printf(maxform, mtop)
writefile.printf("%s", ")")
writefile.printf("%s", "\nsave_window_.variable(\"v\")")
writefile.printf("%s", "\nscene_vector_[6] = save_window_")
writefile.printf("%s", "\n{save_window_.view(")
writefile.printf(maxform, mleft)
writefile.printf("%s", ", ")
writefile.printf(maxform, mbottom)
writefile.printf("%s", ", ")
writefile.printf(maxform, mwidth)
writefile.printf("%s", ", ")
writefile.printf(maxform, mheight)
writefile.printf("%s", ", ")
writefile.printf(maxform, sleft)
writefile.printf("%s", ", ")
writefile.printf(maxform, screentop)
writefile.printf("%s", ", ")
writefile.printf(maxform, swidth)
writefile.printf("%s", ", ")
writefile.printf(maxform, sheight)
writefile.printf("%s", ")}")
// writefile.printf("%s", "\nsave_window_.exec_menu(\"View = plot\")")
writefile.printf("%s", "\nfast_flush_list.append(save_window_)")
writefile.printf("%s", "\nsave_window_.save_name(\"fast_flush_list.\")")
writefile.printf("%s", "\nsave_window_.exec_menu(\"Show Diam\")")
if (mode == 0) {
ev = getmaxev($1)
} else if (mode == 1) {
ev = 40
}
if (ev < 0) {
// 25 point colour map (purple-yellow)
// writefile.printf("%s", "\nsave_window_.colormap(25, 0)")
// writefile.printf("%s", "\nsave_window_.colormap(0, 255, 255, 255)")
// writefile.printf("%s", "\nsave_window_.colormap(1, 255, 255, 0)")
// writefile.printf("%s", "\nsave_window_.colormap(2, 255, 247, 0)")
// writefile.printf("%s", "\nsave_window_.colormap(3, 255, 239, 0)")
// writefile.printf("%s", "\nsave_window_.colormap(4, 255, 223, 0)")
// writefile.printf("%s", "\nsave_window_.colormap(5, 255, 207, 0)")
// writefile.printf("%s", "\nsave_window_.colormap(6, 255, 191, 0)")
// writefile.printf("%s", "\nsave_window_.colormap(7, 255, 175, 0)")
// writefile.printf("%s", "\nsave_window_.colormap(8, 255, 159, 0)")
// writefile.printf("%s", "\nsave_window_.colormap(9, 255, 143, 0)")
// writefile.printf("%s", "\nsave_window_.colormap(10, 255, 127, 0)")
// writefile.printf("%s", "\nsave_window_.colormap(11, 255, 111, 0)")
// writefile.printf("%s", "\nsave_window_.colormap(12, 255, 95, 7)")
// writefile.printf("%s", "\nsave_window_.colormap(13, 255, 79, 15)")
// writefile.printf("%s", "\nsave_window_.colormap(14, 239, 63, 31)")
// writefile.printf("%s", "\nsave_window_.colormap(15, 223, 47, 47)")
// writefile.printf("%s", "\nsave_window_.colormap(16, 207, 31, 47)")
// writefile.printf("%s", "\nsave_window_.colormap(17, 207, 0, 63)")
// writefile.printf("%s", "\nsave_window_.colormap(18, 191, 0, 79)")
// writefile.printf("%s", "\nsave_window_.colormap(19, 175, 0, 95)")
// writefile.printf("%s", "\nsave_window_.colormap(20, 159, 0, 111)")
// writefile.printf("%s", "\nsave_window_.colormap(21, 143, 0, 127)")
// writefile.printf("%s", "\nsave_window_.colormap(22, 127, 0, 143)")
// writefile.printf("%s", "\nsave_window_.colormap(23, 111, 0, 111)")
// writefile.printf("%s", "\nsave_window_.colormap(24, 95, 0, 95)")
// 13 point colour map
// writefile.printf("%s", "\nsave_window_.colormap(13, 0)")
// writefile.printf("%s", "\nsave_window_.colormap(0, 255, 255, 200)")
// writefile.printf("%s", "\nsave_window_.colormap(1, 255, 255, 0)")
// writefile.printf("%s", "\nsave_window_.colormap(2, 255, 239, 0)")
// writefile.printf("%s", "\nsave_window_.colormap(3, 255, 207, 0)")
// writefile.printf("%s", "\nsave_window_.colormap(4, 255, 175, 0)")
// writefile.printf("%s", "\nsave_window_.colormap(5, 255, 143, 0)")
// writefile.printf("%s", "\nsave_window_.colormap(6, 255, 111, 0)")
// writefile.printf("%s", "\nsave_window_.colormap(7, 255, 79, 15)")
// writefile.printf("%s", "\nsave_window_.colormap(8, 223, 47, 47)")
// writefile.printf("%s", "\nsave_window_.colormap(9, 207, 0, 63)")
// writefile.printf("%s", "\nsave_window_.colormap(10, 175, 0, 95)")
// writefile.printf("%s", "\nsave_window_.colormap(11, 143, 0, 127)")
// writefile.printf("%s", "\nsave_window_.colormap(12, 111, 0, 111)")
// 7 point colour map (purple-yellow)
// writefile.printf("%s", "\nsave_window_.colormap(7, 0)")
// writefile.printf("%s", "\nsave_window_.colormap(0, 255, 255, 200)")
// writefile.printf("%s", "\nsave_window_.colormap(1, 255, 239, 0)")
// writefile.printf("%s", "\nsave_window_.colormap(2, 255, 175, 0)")
// writefile.printf("%s", "\nsave_window_.colormap(3, 255, 111, 0)")
// writefile.printf("%s", "\nsave_window_.colormap(4, 223, 47, 47)")
// writefile.printf("%s", "\nsave_window_.colormap(5, 175, 0, 95)")
// writefile.printf("%s", "\nsave_window_.colormap(6, 111, 0, 111)")
if (mode == 0) {
// 7 point colour map (blue-red)
writefile.printf("%s", "\nsave_window_.colormap(7, 0)")
writefile.printf("%s", "\nsave_window_.colormap(0, 233, 34, 16)")
writefile.printf("%s", "\nsave_window_.colormap(1, 225, 157, 18)")
writefile.printf("%s", "\nsave_window_.colormap(2, 174, 219, 20)")
writefile.printf("%s", "\nsave_window_.colormap(3, 23, 208, 44)")
writefile.printf("%s", "\nsave_window_.colormap(4, 25, 198, 195)")
writefile.printf("%s", "\nsave_window_.colormap(5, 23, 88, 188)")
writefile.printf("%s", "\nsave_window_.colormap(6, 7, 0, 149)")
} else if (mode == 1) {
// 13 point colour map (blue-red)
writefile.printf("%s", "\nsave_window_.colormap(13, 0)")
writefile.printf("%s", "\nsave_window_.colormap(0, 233, 34, 16)")
writefile.printf("%s", "\nsave_window_.colormap(1, 229, 97, 17)")
writefile.printf("%s", "\nsave_window_.colormap(2, 225, 157, 18)")
writefile.printf("%s", "\nsave_window_.colormap(3, 222, 214, 19)")
writefile.printf("%s", "\nsave_window_.colormap(4, 174, 219, 20)")
writefile.printf("%s", "\nsave_window_.colormap(5, 101, 215, 21)")
writefile.printf("%s", "\nsave_window_.colormap(6, 23, 208, 44)")
writefile.printf("%s", "\nsave_window_.colormap(7, 24, 202, 138)")
writefile.printf("%s", "\nsave_window_.colormap(8, 25, 198, 195)")
writefile.printf("%s", "\nsave_window_.colormap(9, 26, 143, 194)")
writefile.printf("%s", "\nsave_window_.colormap(10, 23, 88, 188)")
writefile.printf("%s", "\nsave_window_.colormap(11, 6, 36, 169)")
writefile.printf("%s", "\nsave_window_.colormap(12, 7, 0, 149)")
}
writefile.printf("%s", "\nsave_window_.scale(")
writefile.printf(vform, ev)
writefile.printf("%s", ", ")
writefile.printf("%d", 0)
} else {
// 25 point colour map (purple-yellow)
// writefile.printf("%s", "\nsave_window_.colormap(25, 0)")
// writefile.printf("%s", "\nsave_window_.colormap(0, 95, 0, 95)")
// writefile.printf("%s", "\nsave_window_.colormap(1, 111, 0, 111)")
// writefile.printf("%s", "\nsave_window_.colormap(2, 127, 0, 143)")
// writefile.printf("%s", "\nsave_window_.colormap(3, 143, 0, 127)")
// writefile.printf("%s", "\nsave_window_.colormap(4, 159, 0, 111)")
// writefile.printf("%s", "\nsave_window_.colormap(5, 175, 0, 95)")
// writefile.printf("%s", "\nsave_window_.colormap(6, 191, 0, 79)")
// writefile.printf("%s", "\nsave_window_.colormap(7, 207, 0, 63)")
// writefile.printf("%s", "\nsave_window_.colormap(8, 207, 31, 47)")
// writefile.printf("%s", "\nsave_window_.colormap(9, 223, 47, 47)")
// writefile.printf("%s", "\nsave_window_.colormap(10, 239, 63, 31)")
// writefile.printf("%s", "\nsave_window_.colormap(11, 255, 79, 15)")
// writefile.printf("%s", "\nsave_window_.colormap(12, 255, 95, 7)")
// writefile.printf("%s", "\nsave_window_.colormap(13, 255, 111, 0)")
// writefile.printf("%s", "\nsave_window_.colormap(14, 255, 127, 0)")
// writefile.printf("%s", "\nsave_window_.colormap(15, 255, 143, 0)")
// writefile.printf("%s", "\nsave_window_.colormap(16, 255, 159, 0)")
// writefile.printf("%s", "\nsave_window_.colormap(17, 255, 175, 0)")
// writefile.printf("%s", "\nsave_window_.colormap(18, 255, 191, 0)")
// writefile.printf("%s", "\nsave_window_.colormap(19, 255, 207, 0)")
// writefile.printf("%s", "\nsave_window_.colormap(20, 255, 223, 0)")
// writefile.printf("%s", "\nsave_window_.colormap(21, 255, 239, 0)")
// writefile.printf("%s", "\nsave_window_.colormap(22, 255, 247, 0)")
// writefile.printf("%s", "\nsave_window_.colormap(23, 255, 255, 0)")
// writefile.printf("%s", "\nsave_window_.colormap(24, 255, 255, 255)")
// 13 point colour map (purple-yellow)
// writefile.printf("%s", "\nsave_window_.colormap(13, 0)")
// writefile.printf("%s", "\nsave_window_.colormap(0, 111, 0, 111)")
// writefile.printf("%s", "\nsave_window_.colormap(1, 143, 0, 127)")
// writefile.printf("%s", "\nsave_window_.colormap(2, 175, 0, 95)")
// writefile.printf("%s", "\nsave_window_.colormap(3, 207, 0, 63)")
// writefile.printf("%s", "\nsave_window_.colormap(4, 223, 47, 47)")
// writefile.printf("%s", "\nsave_window_.colormap(5, 255, 79, 15)")
// writefile.printf("%s", "\nsave_window_.colormap(6, 255, 111, 0)")
// writefile.printf("%s", "\nsave_window_.colormap(7, 255, 143, 0)")
// writefile.printf("%s", "\nsave_window_.colormap(8, 255, 175, 0)")
// writefile.printf("%s", "\nsave_window_.colormap(9, 255, 207, 0)")
// writefile.printf("%s", "\nsave_window_.colormap(10, 255, 239, 0)")
// writefile.printf("%s", "\nsave_window_.colormap(11, 255, 255, 0)")
// writefile.printf("%s", "\nsave_window_.colormap(12, 255, 255, 200)")
// 7 point colour map (purple-yellow)
// writefile.printf("%s", "\nsave_window_.colormap(7, 0)")
// writefile.printf("%s", "\nsave_window_.colormap(0, 111, 0, 111)")
// writefile.printf("%s", "\nsave_window_.colormap(1, 175, 0, 95)")
// writefile.printf("%s", "\nsave_window_.colormap(2, 223, 47, 47)")
// writefile.printf("%s", "\nsave_window_.colormap(3, 255, 111, 0)")
// writefile.printf("%s", "\nsave_window_.colormap(4, 255, 175, 0)")
// writefile.printf("%s", "\nsave_window_.colormap(5, 255, 239, 0)")
// writefile.printf("%s", "\nsave_window_.colormap(6, 255, 255, 200)")
if (mode == 0) {
// 7 point colour map (blue-red)
writefile.printf("%s", "\nsave_window_.colormap(7, 0)")
writefile.printf("%s", "\nsave_window_.colormap(0, 7, 0, 149)")
writefile.printf("%s", "\nsave_window_.colormap(1, 23, 88, 188)")
writefile.printf("%s", "\nsave_window_.colormap(2, 25, 198, 195)")
writefile.printf("%s", "\nsave_window_.colormap(3, 23, 208, 44)")
writefile.printf("%s", "\nsave_window_.colormap(4, 174, 219, 20)")
writefile.printf("%s", "\nsave_window_.colormap(5, 225, 157, 18)")
writefile.printf("%s", "\nsave_window_.colormap(6, 233, 34, 16)")
} else if (mode == 1) {
// 13 point colour map (blue-red)
writefile.printf("%s", "\nsave_window_.colormap(13, 0)")
writefile.printf("%s", "\nsave_window_.colormap(0, 7, 0, 149)")
writefile.printf("%s", "\nsave_window_.colormap(1, 6, 36, 169)")
writefile.printf("%s", "\nsave_window_.colormap(2, 23, 88, 188)")
writefile.printf("%s", "\nsave_window_.colormap(3, 26, 143, 194)")
writefile.printf("%s", "\nsave_window_.colormap(4, 25, 198, 195)")
writefile.printf("%s", "\nsave_window_.colormap(5, 24, 202, 138)")
writefile.printf("%s", "\nsave_window_.colormap(6, 23, 208, 44)")
writefile.printf("%s", "\nsave_window_.colormap(7, 101, 215, 21)")
writefile.printf("%s", "\nsave_window_.colormap(8, 174, 219, 20)")
writefile.printf("%s", "\nsave_window_.colormap(9, 222, 214, 19)")
writefile.printf("%s", "\nsave_window_.colormap(10, 225, 157, 18)")
writefile.printf("%s", "\nsave_window_.colormap(11, 229, 97, 17)")
writefile.printf("%s", "\nsave_window_.colormap(12, 233, 34, 16)")
}
writefile.printf("%s", "\nsave_window_.scale(")
if (mode == 0) {
writefile.printf("%d", 0)
} else if (mode == 1) {
writefile.printf("%d", -80)
}
writefile.printf("%s", ", ")
writefile.printf(vform, ev)
}
writefile.printf("%s", ")")
writefile.printf("%s", "\nsave_window_.exec_menu(\"Shape Plot\")")
writefile.printf("%s", "\n}")
writefile.printf("%s", "\nobjectvar scene_vector_[1]")
writefile.printf("%s", "\n{doNotify()}")
writefile.close()
chdir(cwd)
if (chk) cwd = ""
}
// --------------------------------------------------------------------------------
// --------------------------------------------------------------------------------
proc writemulfit() {local chk
if (numarg() == 1 && argtype(1) == 2) {
if (!strcmp($s1, "help")) {
print "takes no arguments. to initialize mulrunfitter. used wih writefitdata and writefitparam"
stop
}
}
chk = chkcwd(cwd)
if (chk) cwd = getcwd()
chdir(ses)
writefile.wopen(savedatafilestr)
writefile.printf("%s", "objectvar save_window_, rvp_")
writefile.printf("%s", "\nobjectvar scene_vector_[8]")
writefile.printf("%s", "\nobjectvar ocbox_, ocbox_list_, scene_, scene_list_")
writefile.printf("%s", "\n{ocbox_list_ = new List() scene_list_ = new List()}")
writefile.printf("%s", "\n\n//Begin MulRunFitter[0]")
writefile.printf("%s", "\n{")
writefile.printf("%s", "\nload_file(\"mulfit.hoc\", \"MulRunFitter\")")
writefile.printf("%s", "\n}")
writefile.printf("%s", "\n{")
writefile.printf("%s", "\nocbox_ = new MulRunFitter(1)")
writefile.printf("%s", "\n}")
writefile.printf("%s", "\n{object_push(ocbox_)}")
writefile.printf("%s", "\n{")
writefile.printf("%s", "\nversion(6)")
writefile.printf("%s", "\nranfac = 2")
writefile.printf("%s", "\nfspec = new File(\"")
writefile.printf("%s", savedataparamfilestr)
writefile.printf("%s", "\")")
writefile.printf("%s", "\nfdat = new File(\"")
writefile.printf("%s", savedatadatafilestr)
writefile.printf("%s", "\")")
writefile.printf("%s", "\nread_data()")
writefile.printf("%s", "\nbuild()")
writefile.printf("%s", "\n}")
writefile.printf("%s", "\nopt.set_optimizer(\"MulfitPraxWrap\")")
writefile.printf("%s", "\n{object_push(opt.optimizer)}")
writefile.printf("%s", "\n{")
writefile.printf("%s", "\nnstep = 0")
writefile.printf("%s", "\n}")
writefile.printf("%s", "\n{object_pop()}")
writefile.printf("%s", "\n{p.gengui(0, 1234, 114, 608.64, 553.92)}")
// writefile.printf("%s", "\n{p.showargs(1328, 729, 269.76, 257.28)}")
// writefile.printf("%s", "\n{p.showdomain(1581, 729, 272.64, 213.12)}")
writefile.printf("%s", "\n{optrestore(995, 729, 312.96, 338.88)}")
writefile.printf("%s", "\n{object_pop()}")
writefile.printf("%s", "\n{")
writefile.printf("%s", "\nocbox_.map(\"MulRunFitter[0]\", 612, 729, 360.96, 199.68)")
writefile.printf("%s", "\n}")
writefile.printf("%s", "\nobjref ocbox_")
writefile.printf("%s", "\n//End MulRunFitter[0]")
writefile.printf("%s", "\nobjectvar scene_vector_[1]")
writefile.printf("%s", "\n{doNotify()}")
writefile.close()
chdir(cwd)
if (chk) cwd = ""
}
// --------------------------------------------------------------------------------
// --------------------------------------------------------------------------------
proc writefitdata() {local chk, iiamp, lines, f, row
if (numarg() == 1 && argtype(1) == 2) {
if (!strcmp($s1, "help")) {
print "Enter:"
print "1) optional: iiamp index within iampvec. Default: current iiamp (saved in iiamp-*mode*.dat)"
print "output: saved vdat associated with current or given iiamp in given mode (savedata-*mode*.ses.fd1)"
stop
}
}
chk = chkcwd(cwd)
if (chk) cwd = getcwd()
if (numarg() == 1) {
iiamp = $1
} else {
iiamp = 0
}
chdir(ses)
writefile.wopen(savedatadatafilestr)
lines = 2*nrows + 8
f = 0
while (f < nf) {
writefile.printf("\n%s%d%s", "RegionFitness xdat ydat boundary weight (lines=", lines, ") 1")
writefile.printf("\n%s", "||")
writefile.printf("\n%d", nrows)
sprint(tempstr, "%s%s", "\n", tform)
for row = 0, nrows-1 writefile.printf(tempstr, vdatmat.x[row][0])
sprint(tempstr, "%s%s", "\n", vform)
for row = 0, nrows-1 writefile.printf(tempstr, vdatmat.x[row][iiamp*nf+f+1])
sprint(tempstr, "%s%s%s%s%s", "\n\n%d\n", tform, "\n", tform, "\n\n%d\n%d\n")
writefile.printf(tempstr, 2, vdatmat.x[0][0], vdatmat.x[nrows-1][0], 0, 1)
f += 1
}
writefile.printf("\n%g\n", 1)
writefile.close()
tempstr = ""
chdir(cwd)
if (chk) cwd = ""
}
// --------------------------------------------------------------------------------
// --------------------------------------------------------------------------------
proc writefitparam() {local chk, k, f
if (numarg() == 1 && argtype(1) == 2) {
if (!strcmp($s1, "help")) {
print "saves default parameters associated with given mode (savedata-*mode*.ses.ft1)"
stop
}
}
chk = chkcwd(cwd)
if (chk) cwd = getcwd()
chdir(ses)
writefile.wopen(savedataparamfilestr)
writefile.printf("%s", "ParmFitness: nFit")
if (mode == 0) {
writefile.printf("%s", "\n\tFitnessGenerator: nFit Passive")
} else if (mode == 1) {
writefile.printf("%s", "\n\tFitnessGenerator: nFit Active")
}
for f = 0, nf-1 {
writefile.printf("%s", "\n\t\tRegionFitness:\t")
writefile.printf("%s", fitvarlist.o(f).label)
}
writefile.printf("%s", "\n\n\tParameters:")
for k = 0, pnamelist.count-1 {
sprint(tempstr, "%s%s%s%s", "\"", pnamelist.o(k).label, "_norm", "\"")
writefile.printf("\n\t\t%s", tempstr)
writefile.printf("%s", ", ")
writefile.printf(maxform, pdefvec.x[k]/plowvec.x[k])
writefile.printf("%s", ", ")
writefile.printf("%d", 1)
writefile.printf("%s", ", ")
writefile.printf(maxform, phivec.x[k]/plowvec.x[k])
writefile.printf("%s", ", ")
writefile.printf("%d", doargvec.x[k])
writefile.printf("%s", ", ")
writefile.printf("%d", uselogvec.x[k])
}
writefile.printf("%s", "\nEnd ParmFitness")
writefile.close()
tempstr = ""
chdir(cwd)
if (chk) cwd = ""
}
// --------------------------------------------------------------------------------
// --------------------------------------------------------------------------------
proc writenagating() {local chk
if (numarg() == 1 && argtype(1) == 2) {
if (!strcmp($s1, "help")) {
print "No arguments required."
print "output: saves default Na gating parameters to session file (na_gating.ses)"
stop
}
}
chk = chkcwd(cwd)
if (chk) cwd = getcwd()
chdir(ses)
writefile.wopen("na_gating.ses")
writefile.printf("%s", "objectvar save_window_, rvp_")
writefile.printf("%s", "\nobjectvar scene_vector_[8]")
writefile.printf("%s", "\nobjectvar ocbox_, ocbox_list_, scene_, scene_list_")
writefile.printf("%s", "\n{ocbox_list_ = new List() scene_list_ = new List()}")
writefile.printf("%s", "\n{")
writefile.printf("%s", "\nxpanel(\"nais (Globals)\", 0)")
writefile.printf("%s", "\nvShift_nais = ")
writefile.printf("%g", vShift_nais)
writefile.printf("%s", "\nxvalue(\"vShift_nais\",\"vShift_nais\", 1,\"\", 0, 0 )")
writefile.printf("%s", "\nvShift_inact_nais = ")
writefile.printf("%g", vShift_inact_nais)
writefile.printf("%s", "\nxvalue(\"vShift_inact_nais\",\"vShift_inact_nais\", 1,\"\", 0, 0 )")
writefile.printf("%s", "\nmaxrate_nais = ")
writefile.printf("%g", maxrate_nais)
writefile.printf("%s", "\nxvalue(\"maxrate_nais\",\"maxrate_nais\", 1,\"\", 0, 0 )")
writefile.printf("%s", "\ntemp_nais = ")
writefile.printf("%g", temp_nais)
writefile.printf("%s", "\nxvalue(\"temp_nais\",\"temp_nais\", 1,\"\", 0, 0 )")
writefile.printf("%s", "\nq10_nais = ")
writefile.printf("%g", q10_nais)
writefile.printf("%s", "\nxvalue(\"q10_nais\",\"q10_nais\", 1,\"\", 0, 0 )")
writefile.printf("%s", "\nq10h_nais = ")
writefile.printf("%g", q10h_nais)
writefile.printf("%s", "\nxvalue(\"q10h_nais\",\"q10h_nais\", 1,\"\", 0, 0 )")
writefile.printf("%s", "\nxpanel(294,636)")
writefile.printf("%s", "\n}")
writefile.printf("%s", "\n{")
writefile.printf("%s", "\nxpanel(\"nax (Globals)\", 0)")
writefile.printf("%s", "\nvShift_nax = ")
writefile.printf("%g", vShift_nax)
writefile.printf("%s", "\nxvalue(\"vShift_nax\",\"vShift_nax\", 1,\"\", 0, 0 )")
writefile.printf("%s", "\nvShift_inact_nax = ")
writefile.printf("%g", vShift_inact_nax)
writefile.printf("%s", "\nxvalue(\"vShift_inact_nax\",\"vShift_inact_nax\", 1,\"\", 0, 0 )")
writefile.printf("%s", "\nmaxrate_nax = ")
writefile.printf("%g", maxrate_nax)
writefile.printf("%s", "\nxvalue(\"maxrate_nax\",\"maxrate_nax\", 1,\"\", 0, 0 )")
writefile.printf("%s", "\ntemp_nax = ")
writefile.printf("%g", temp_nax)
writefile.printf("%s", "\nxvalue(\"temp_nax\",\"temp_nax\", 1,\"\", 0, 0 )")
writefile.printf("%s", "\nq10_nax = ")
writefile.printf("%g", q10_nax)
writefile.printf("%s", "\nxvalue(\"q10_nax\",\"q10_nax\", 1,\"\", 0, 0 )")
writefile.printf("%s", "\nq10h_nax = ")
writefile.printf("%g", q10h_nax)
writefile.printf("%s", "\nxvalue(\"q10h_nax\",\"q10h_nax\", 1,\"\", 0, 0 )")
writefile.printf("%s", "\nxpanel(294,866)")
writefile.printf("%s", "\n}")
writefile.printf("%s", "\n{")
writefile.printf("%s", "\nxpanel(\"na (Globals)\", 0)")
writefile.printf("%s", "\nvShift_na = ")
writefile.printf("%g", vShift_na)
writefile.printf("%s", "\nxvalue(\"vShift_na\",\"vShift_na\", 1,\"\", 0, 0 )")
writefile.printf("%s", "\nvShift_inact_na = ")
writefile.printf("%g", vShift_inact_na)
writefile.printf("%s", "\nxvalue(\"vShift_inact_na\",\"vShift_inact_na\", 1,\"\", 0, 0 )")
writefile.printf("%s", "\nmaxrate_na = ")
writefile.printf("%g", maxrate_na)
writefile.printf("%s", "\nxvalue(\"maxrate_na\",\"maxrate_na\", 1,\"\", 0, 0 )")
writefile.printf("%s", "\ntemp_na = ")
writefile.printf("%g", temp_na)
writefile.printf("%s", "\nxvalue(\"temp_na\",\"temp_na\", 1,\"\", 0, 0 )")
writefile.printf("%s", "\nq10_na = ")
writefile.printf("%g", q10_na)
writefile.printf("%s", "\nxvalue(\"q10_na\",\"q10_na\", 1,\"\", 0, 0 )")
writefile.printf("%s", "\nq10h_na = ")
writefile.printf("%g", q10h_na)
writefile.printf("%s", "\nxvalue(\"q10h_na\",\"q10h_na\", 1,\"\", 0, 0 )")
writefile.printf("%s", "\nxpanel(294,1096)")
writefile.printf("%s", "\n}")
writefile.printf("%s", "\nobjectvar scene_vector_[1]")
writefile.printf("%s", "\n{doNotify()}")
writefile.close()
chdir(cwd)
if (chk) cwd = ""
}
// --------------------------------------------------------------------------------
// ----------------------------writemiscparam--------------------------------------
proc writemiscparam() {local chk
if (numarg() == 1 && argtype(1) == 2) {
if (!strcmp($s1, "help")) {
print "No arguments required."
print "output: save new default misc param panel window in session file (miscparam.ses)"
stop
}
}
chk = chkcwd(cwd)
if (chk) cwd = getcwd()
chdir(ses)
writefile.wopen("miscparam.ses")
writefile.printf("%s", "objectvar save_window_, rvp_")
writefile.printf("%s", "\nobjectvar scene_vector_[8]")
writefile.printf("%s", "\nobjectvar ocbox_, ocbox_list_, scene_, scene_list_")
writefile.printf("%s", "\n{ocbox_list_ = new List() scene_list_ = new List()}")
writefile.printf("%s", "\n{")
writefile.printf("%s", "\nxpanel(\"Misc Parameters\", 0)")
writefile.printf("%s", "\nHoff = ")
writefile.printf("%g", Hoff)
writefile.printf("%s", "\nxvalue(\"Hoff\",\"Hoff\", 1,\"\", 0, 1 )")
writefile.printf("%s", "\ntauh = ")
writefile.printf("%g", tauh)
writefile.printf("%s", "\nxvalue(\"tauh\",\"tauh\", 1,\"\", 0, 1 )")
writefile.printf("%s", "\nmaxdist_Na_apic = ")
writefile.printf("%g", maxdist_Na_apic)
writefile.printf("%s", "\nxvalue(\"maxdist_Na_apic\",\"maxdist_Na_apic\", 1,\"\", 0, 1 )")
writefile.printf("%s", "\nmaxdist_Na_bas = ")
writefile.printf("%g", maxdist_Na_bas)
writefile.printf("%s", "\nxvalue(\"maxdist_Na_bas\",\"maxdist_Na_bas\", 1,\"\", 0, 1 )")
writefile.printf("%s", "\nlc_Kv_Kv1 = ")
writefile.printf("%g", lc_Kv_Kv1)
writefile.printf("%s", "\nxvalue(\"lc_Kv_Kv1\",\"lc_Kv_Kv1\", 1,\"\", 0, 1 )")
writefile.printf("%s", "\nxpanel(0,791)")
writefile.printf("%s", "\n}")
writefile.printf("%s", "\nobjectvar scene_vector_[1]")
writefile.printf("%s", "\n{doNotify()}")
writefile.close()
chdir(cwd)
if (chk) cwd = ""
}
// --------------------------------------------------------------------------------
