/*END CVODE*/ static int castate (double* _p, Datum* _ppvar, Datum* _thread, _NrnThread* _nt) { { evaluate_fct ( _threadargscomma_ v ) ; m = m + (1. - exp(dt*(( ( ( - 1.0 ) ) ) / tau_m)))*(- ( ( ( m_inf ) ) / tau_m ) / ( ( ( ( - 1.0) ) ) / tau_m ) - m) ; h = h + (1. - exp(dt*(( ( ( - 1.0 ) ) ) / tau_h)))*(- ( ( ( h_inf ) ) / tau_h ) / ( ( ( ( - 1.0) ) ) / tau_h ) - h) ; } return 0; }
/*CVODE*/ static int _ode_spec1 (double* _p, Datum* _ppvar, Datum* _thread, _NrnThread* _nt) {int _reset = 0; { evaluate_fct ( _threadargscomma_ v ) ; Dm = ( m_inf - m ) / tau_m ; Dh = ( h_inf - h ) / tau_h ; } return _reset; }
/*END CVODE*/ static int states () {_reset=0; { evaluate_fct ( _threadargscomma_ v ) ; m = m + (1. - exp(dt*(( ( ( - 1.0 ) ) ) / taum)))*(- ( ( ( m_inf ) ) / taum ) / ( ( ( ( - 1.0) ) ) / taum ) - m) ; } return 0; }
/*CVODE*/ static int _ode_spec1 () {_reset=0; { evaluate_fct ( _threadargscomma_ v ) ; Dm = ( m_inf - m ) / taum ; } return _reset; }
static int _ode_matsol1 (double* _p, Datum* _ppvar, Datum* _thread, _NrnThread* _nt) { evaluate_fct ( _threadargscomma_ v ) ; Dmp = Dmp / (1. - dt*( ( ( ( - 1.0 ) ) ) / tau_mp )) ; Dm = Dm / (1. - dt*( ( ( ( - 1.0 ) ) ) / tau_m )) ; Dh = Dh / (1. - dt*( ( ( ( - 1.0 ) ) ) / tau_h )) ; Ds = Ds / (1. - dt*( ( ( ( - 1.0 ) ) ) / tau_s )) ; return 0; }
/*END CVODE*/ static int castate () {_reset=0; { evaluate_fct ( _threadargscomma_ v ) ; m = m + (1. - exp(dt*(( ( ( - 1.0 ) ) ) / tau_m)))*(- ( ( ( m_inf ) ) / tau_m ) / ( ( ( ( - 1.0) ) ) / tau_m ) - m) ; h = h + (1. - exp(dt*(( ( ( - 1.0 ) ) ) / tau_h)))*(- ( ( ( h_inf ) ) / tau_h ) / ( ( ( ( - 1.0) ) ) / tau_h ) - h) ; } return 0; }
static int _ode_matsol1 () { evaluate_fct ( _threadargscomma_ v ) ; DW = DW / (1. - dt*( ( ( ( - 1.0 ) ) ) / tau_W )) ; DX = DX / (1. - dt*( ( ( ( - 1.0 ) ) ) / tau_X )) ; Dcai = Dcai / (1. - dt*( (caGain)*(( ( - (kca)*(1.0) ) )) )) ; DCsk = DCsk / (1. - dt*( ( - ( ( 1.0 ) ) / tau_sk ) )) ; DB = DB / (1. - dt*( ( ( ( - 1.0 ) ) ) / tau_B )) ; return 0; }
/*END CVODE*/ static int states (double* _p, Datum* _ppvar, Datum* _thread, _NrnThread* _nt) { { evaluate_fct ( _threadargscomma_ v ) ; mp = mp + (1. - exp(dt*(( ( ( - 1.0 ) ) ) / tau_mp)))*(- ( ( ( mp_inf ) ) / tau_mp ) / ( ( ( ( - 1.0) ) ) / tau_mp ) - mp) ; m = m + (1. - exp(dt*(( ( ( - 1.0 ) ) ) / tau_m)))*(- ( ( ( m_inf ) ) / tau_m ) / ( ( ( ( - 1.0) ) ) / tau_m ) - m) ; h = h + (1. - exp(dt*(( ( ( - 1.0 ) ) ) / tau_h)))*(- ( ( ( h_inf ) ) / tau_h ) / ( ( ( ( - 1.0) ) ) / tau_h ) - h) ; s = s + (1. - exp(dt*(( ( ( - 1.0 ) ) ) / tau_s)))*(- ( ( ( s_inf ) ) / tau_s ) / ( ( ( ( - 1.0) ) ) / tau_s ) - s) ; } return 0; }
static int states ( ) { evaluate_fct ( _threadargscomma_ v ) ; m = m + m_exp * ( m_inf - m ) ; h = h + h_exp * ( h_inf - h ) ; n = n + n_exp * ( n_inf - n ) ; /*VERBATIM*/ return 0; return 0; }
static void _hoc_evaluate_fct(void) { double _r; double* _p; Datum* _ppvar; Datum* _thread; _NrnThread* _nt; if (_extcall_prop) {_p = _extcall_prop->param; _ppvar = _extcall_prop->dparam;}else{ _p = (double*)0; _ppvar = (Datum*)0; } _thread = _extcall_thread; _nt = nrn_threads; _r = 1.; evaluate_fct ( _p, _ppvar, _thread, _nt, *getarg(1) ); hoc_retpushx(_r); }
/*CVODE*/ static int _ode_spec1 () {_reset=0; { evaluate_fct ( _threadargscomma_ v ) ; DW = ( W_inf - W ) / tau_W ; DX = ( X_inf - X ) / tau_X ; Dcai = caGain * ( - ica * area * 1e-11 / ( 2.0 * FARADAY * vol ) - kca * cai ) ; DCsk = - a_sk * ica - ( Csk - C_gamma ) / tau_sk ; DB = ( B_inf - B ) / tau_B ; } return _reset; }
/*END CVODE*/ static int states () {_reset=0; { evaluate_fct ( _threadargscomma_ v ) ; W = W + (1. - exp(dt*(( ( ( - 1.0 ) ) ) / tau_W)))*(- ( ( ( W_inf ) ) / tau_W ) / ( ( ( ( - 1.0) ) ) / tau_W ) - W) ; X = X + (1. - exp(dt*(( ( ( - 1.0 ) ) ) / tau_X)))*(- ( ( ( X_inf ) ) / tau_X ) / ( ( ( ( - 1.0) ) ) / tau_X ) - X) ; cai = cai + (1. - exp(dt*((caGain)*(( ( - (kca)*(1.0) ) )))))*(- ( (caGain)*(( ( ((- ica)*(area))*(1e-11) ) / ( 2.0 * FARADAY * vol ) )) ) / ( (caGain)*(( ( - (kca)*(1.0)) )) ) - cai) ; Csk = Csk + (1. - exp(dt*(( - ( ( 1.0 ) ) / tau_sk ))))*(- ( (- a_sk)*(ica) - ( ( ( - C_gamma ) ) ) / tau_sk ) / ( ( - ( ( 1.0 ) ) / tau_sk) ) - Csk) ; B = B + (1. - exp(dt*(( ( ( - 1.0 ) ) ) / tau_B)))*(- ( ( ( B_inf ) ) / tau_B ) / ( ( ( ( - 1.0) ) ) / tau_B ) - B) ; } return 0; }
static void initmodel(double* _p, Datum* _ppvar, Datum* _thread, _NrnThread* _nt) { int _i; double _save;{ h = h0; m = m0; { phi_m = pow( 5.0 , ( ( celsius - 24.0 ) / 10.0 ) ) ; phi_h = pow( 3.0 , ( ( celsius - 24.0 ) / 10.0 ) ) ; evaluate_fct ( _threadargscomma_ v ) ; m = m_inf ; h = h_inf ; } } }
static void initmodel() { int _i; double _save;_ninits++; _save = t; t = 0.0; { m = m0; { tadj = pow( 3.0 , ( ( celsius - 21.0 ) / 10.0 ) ) ; evaluate_fct ( _threadargscomma_ v ) ; m = m_inf ; } _sav_indep = t; t = _save; } }
static void initmodel(double* _p, Datum* _ppvar, Datum* _thread, _NrnThread* _nt) { int _i; double _save;{ h = h0; m = m0; mp = mp0; s = s0; { q10_1 = pow( 2.2 , ( ( celsius - 20.0 ) / 10.0 ) ) ; q10_2 = pow( 2.9 , ( ( celsius - 20.0 ) / 10.0 ) ) ; q10_3 = pow( 3.0 , ( ( celsius - 36.0 ) / 10.0 ) ) ; evaluate_fct ( _threadargscomma_ v ) ; mp = mp_inf ; m = m_inf ; h = h_inf ; s = s_inf ; } } }
static void initmodel() { int _i; double _save;_ninits++; _save = t; t = 0.0; { B = B0; Csk = Csk0; W = W0; X = X0; { evaluate_fct ( _threadargscomma_ v ) ; W = W_inf ; X = X_inf ; Csk = C_gamma + 0.00001 ; B = B_inf ; } _sav_indep = t; t = _save; } }
static void initmodel() { int _i; double _save;_ninits++; _save = t; t = 0.0; { h = h0; m = m0; { /*VERBATIM*/ Cai = _ion_Cai; Cao = _ion_Cao; phi_m = pow( mx , ( ( celsius - 23.5 ) / 10.0 ) ) ; phi_h = pow( hx , ( ( celsius - 23.5 ) / 10.0 ) ) ; evaluate_fct ( _threadargscomma_ v ) ; m = m_inf ; h = h_inf ; } _sav_indep = t; t = _save; } }
static int _hoc_evaluate_fct() { double _r; _r = 1.; evaluate_fct ( *getarg(1) ); ret(_r); }
static int _ode_matsol1 () { evaluate_fct ( _threadargscomma_ v ) ; Dm = Dm / (1. - dt*( ( ( ( - 1.0 ) ) ) / tau_m )) ; Dh = Dh / (1. - dt*( ( ( ( - 1.0 ) ) ) / tau_h )) ; }
static void _hoc_evaluate_fct(void) { double _r; _r = 1.; evaluate_fct ( *getarg(1) ); hoc_retpushx(_r); }