Check (or fix) the two Forrest models: PC2015aForrest & PC2015bForrest

[lungsi]

Although the models, both full 1088 compartment and reduced two compartment models run it does not appear to run correctly. For example,

from executive import ExecutiveControl
from models.cells.modelPC2015aForrest import PurkinjeCell
#from models.cells.modelPC2015bForrest import PurkinjeCell
from cerebunit.capabilities.cells.measurements import ProducesEphysMeasurement

runtimeparam = {"dt": 0.025, "celsius": 30, "tstop": 1000.0, "v_init": -80.}
stimparam = {"type": ["current", "IClamp"],
             "stimlist": [ {"amp": 0.006, "dur": 800.0, "delay": 100.0} ],
                            "tstop": runtimeparam["tstop"] }

ec = ExecutiveControl()
ec.chosenmodel = PurkinjeCell()

model = ec.launch_model( parameters = runtimeparam, stimparameters = stimparam,
                         stimloc = ec.chosenmodel.cell.soma, onmodel = ec.chosenmodel,
                         capabilities = {"model": "produce_restingVm",
                                         "vtest": ProducesEphysMeasurement},
                         mode = "capability" )

and then

>>> model.prediction
[array([nan]) * mV, array([nan]) * mV]

[lungsi]

If I use runtimeparam = {"dt": 0.025, "celsius": 30, "tstop": 1000.0, "v_init": -75.} for from models.cells.modelPC2015aForrest import PurkinjeCell I get

>>> model.prediction
[array([-49.03888712]) * mV, array([-29.59161421]) * mV]

and for from models.cells.modelPC2015bForrest import PurkinjeCell I get

>>> model.prediction
[array([-53.01861152]) * mV, array([-37.6118527]) * mV]

Although the values are not equal, the point is for "v_init" > 80. it no longer returns nan. I have checked this with "v_init": -70, "v_init": -65 and "v_init": -60.