ModelDB is moving. Check out our new site at https://modeldb.science. The corresponding page is https://modeldb.science/97917.

Cell splitting in neural networks extends strong scaling (Hines et al. 2008)

 Download zip file 
Help downloading and running models
Accession:97917
Neuron tree topology equations can be split into two subtrees and solved on different processors with no change in accuracy, stability, or computational effort; communication costs involve only sending and receiving two double precision values by each subtree at each time step. Application of the cell splitting method to two published network models exhibits good runtime scaling on twice as many processors as could be effectively used with whole-cell balancing.
Reference:
1 . Hines ML, Eichner H, Schürmann F (2008) Neuron splitting in compute-bound parallel network simulations enables runtime scaling with twice as many processors. J Comput Neurosci 25:203-10 [PubMed]
Model Information (Click on a link to find other models with that property)
Model Type: Realistic Network;
Brain Region(s)/Organism: Generic;
Cell Type(s):
Channel(s):
Gap Junctions:
Receptor(s):
Gene(s):
Transmitter(s):
Simulation Environment: NEURON;
Model Concept(s): Methods;
Implementer(s): Hines, Michael [Michael.Hines at Yale.edu];
/
splitcell
nrntraub
cells
hoc
mod
net
README
balcomp.hoc *
bgrunme
bgsmall.sh
bgsplit.sh
cell_templates.hoc *
clear.hoc *
finit.hoc *
fortmap.hoc *
gidcell.hoc
gidcell.ses *
init.hoc
manage_setup.hoc
metisbal.sh
mosinit_orig.hoc *
onecell.hoc
onecell.ses *
prcellstate.hoc *
prepare.sh
printcon.hoc *
spkplt.hoc *
vclampg.hoc *
vcompclamp.hoc *
vcompsim.hoc *
                            
objref type_vinit, vgaba, typename[14]
type_vinit = new Vector(14)
vgaba = new Vector(14)

//cell type id according to the mpirank in the fortran program
//cell.type() returns this value. Note that the table order is below is the
//same as the base order in groucho.hoc
for i=0, 13 {
	j = fscan()
	typename[j] = new String() getstr(typename[j].s, 1)
	type_vinit.x[j] = fscan()  vgaba.x[j] = fscan()
}

suppyrRS_template.hoc:  rank = 0  suppyrRS    -70 -81
suppyrFRB_template.hoc: rank = 1  suppyrFRB   -70 -81
supbask_template.hoc:   rank = 2  supbask     -65 -75
supaxax_template.hoc:   rank = 3  supaxax     -65 -75
supLTS_template.hoc:    rank = 4  supLTS      -65 -75
spinstell_template.hoc: rank = 5  spinstell   -65 -75
tuftIB_template.hoc:    rank = 6  tuftIB      -70 -75
tuftRS_template.hoc:    rank = 7  tuftRS      -70 -75
nontuftRS_template.hoc: rank = 8  nontuftRS   -70 -75
deepbask_template.hoc:  rank = 9  deepbask    -65 -75
deepaxax_template.hoc:  rank = 10 deepaxax    -65 -75
deepLTS_template.hoc:   rank = 11 deepLTS     -65 -75
TCR_template.hoc:       rank = 12 TCR         -70 -81
nRT_template.hoc:       rank = 13 nRT         -75 -75

//note that the TCR rank value is used explicitly in initialization since the
//states are initialized at a different voltage value than the membrane potential.

TCRtype = 12

Loading data, please wait...