static void nrn_state(_NrnThread* _nt, _Memb_list* _ml, int _type){ double _break, _save; Node *_nd; double _v; int* _ni; int _iml, _cntml; #if CACHEVEC _ni = _ml->_nodeindices; #endif _cntml = _ml->_nodecount; for (_iml = 0; _iml < _cntml; ++_iml) { _p = _ml->_data[_iml]; _ppvar = _ml->_pdata[_iml]; _nd = _ml->_nodelist[_iml]; #if CACHEVEC if (use_cachevec) { _v = VEC_V(_ni[_iml]); }else #endif { _nd = _ml->_nodelist[_iml]; _v = NODEV(_nd); } _break = t + .5*dt; _save = t; v=_v; { ena = _ion_ena; ek = _ion_ek; { { for (; t < _break; t += dt) { error = states(); if(error){fprintf(stderr,"at line 66 in file hh2.mod:\n SOLVE states\n"); nrn_complain(_p); abort_run(error);} }} t = _save; } }} }
static void nrn_state(_NrnThread* _nt, _Memb_list* _ml, int _type){ double _break, _save; Node *_nd; double _v; int* _ni; int _iml, _cntml; #if CACHEVEC _ni = _ml->_nodeindices; #endif _cntml = _ml->_nodecount; for (_iml = 0; _iml < _cntml; ++_iml) { _p = _ml->_data[_iml]; _ppvar = _ml->_pdata[_iml]; _nd = _ml->_nodelist[_iml]; #if CACHEVEC if (use_cachevec) { _v = VEC_V(_ni[_iml]); }else #endif { _nd = _ml->_nodelist[_iml]; _v = NODEV(_nd); } _break = t + .5*dt; _save = t; v=_v; { ica = _ion_ica; cai = _ion_cai; cai = _ion_cai; { { for (; t < _break; t += dt) { error = euler(_ninits, 1, _slist1, _dlist1, _p, &t, dt, state, &_temp1); if(error){fprintf(stderr,"at line 87 in file capump.mod:\n SOLVE state METHOD euler\n"); nrn_complain(_p); abort_run(error);} }} t = _save; state(); } { } _ion_cai = cai; }} }
static void nrn_state(_NrnThread* _nt, _Memb_list* _ml, int _type){ double _break, _save; Node *_nd; double _v; int* _ni; int _iml, _cntml; #if CACHEVEC _ni = _ml->_nodeindices; #endif _cntml = _ml->_nodecount; for (_iml = 0; _iml < _cntml; ++_iml) { _p = _ml->_data[_iml]; _ppvar = _ml->_pdata[_iml]; _nd = _ml->_nodelist[_iml]; #if CACHEVEC if (use_cachevec) { _v = VEC_V(_ni[_iml]); }else #endif { _nd = _ml->_nodelist[_iml]; _v = NODEV(_nd); } _break = t + .5*dt; _save = t; v=_v; { ica = _ion_ica; cai = _ion_cai; cai = _ion_cai; { { for (; t < _break; t += dt) { error = _deriv1_advance = 1; derivimplicit(_ninits, 1, _slist1, _dlist1, _p, &t, dt, state, &_temp1); _deriv1_advance = 0; if(error){fprintf(stderr,"at line 72 in file cad_0403.mod:\n SOLVE state METHOD derivimplicit\n"); nrn_complain(_p); abort_run(error);} }} t = _save; } { /*VERBATIM*/ /* printf("CAD.MOD cai:%.10lf\t",cai);*/ } _ion_cai = cai; }} }