static int rates ( _threadargsprotocomma_ double _lv ) { a = Ra * qa * efun ( _threadargscomma_ - ( _lv - tha ) / qa ) ; b = Rb * qa * efun ( _threadargscomma_ ( _lv - tha ) / qa ) ; tadj = pow( q10 , ( ( celsius - temp ) / 10.0 ) ) ; ntau = 1.0 / tadj / ( a + b ) ; ninf = a / ( a + b ) ; return 0; }
double vtrap8 ( _threadargsprotocomma_ double _lx ) { double _lvtrap8; if ( fabs ( ( _lx + ahB ) / ahC ) < 1e-6 ) { _lvtrap8 = - ahA * ahC ; } else { _lvtrap8 = ( ahA * ( - ( _lx + ahB ) ) ) / ( 1.0 - Exp ( _threadargscomma_ ( _lx + ahB ) / ahC ) ) ; } return _lvtrap8; }
double vtrap7 ( _threadargsprotocomma_ double _lx ) { double _lvtrap7; if ( fabs ( ( _lx + bmB ) / bmC ) < 1e-6 ) { _lvtrap7 = - bmA * bmC ; } else { _lvtrap7 = ( bmA * ( - ( _lx + bmB ) ) ) / ( 1.0 - Exp ( _threadargscomma_ ( _lx + bmB ) / bmC ) ) ; } return _lvtrap7; }
double vtrap2 ( _threadargsprotocomma_ double _lx ) { double _lvtrap2; if ( fabs ( ( _lx + bmpB ) / bmpC ) < 1e-6 ) { _lvtrap2 = - bmpA * bmpC ; } else { _lvtrap2 = ( bmpA * ( - ( _lx + bmpB ) ) ) / ( 1.0 - Exp ( _threadargscomma_ ( _lx + bmpB ) / bmpC ) ) ; } return _lvtrap2; }
double vtrap ( _threadargsprotocomma_ double _lx ) { double _lvtrap; if ( _lx < - 50.0 ) { _lvtrap = 0.0 ; } else { _lvtrap = bsA / ( Exp ( _threadargscomma_ ( _lx + bsB ) / bsC ) + 1.0 ) ; } return _lvtrap; }
static int evaluate_fct ( _threadargsprotocomma_ double _lv ) { double _la , _lb , _lv2 ; _la = q10_1 * vtrap1 ( _threadargscomma_ _lv ) ; _lb = q10_1 * vtrap2 ( _threadargscomma_ _lv ) ; tau_mp = 1.0 / ( _la + _lb ) ; mp_inf = _la / ( _la + _lb ) ; _la = q10_1 * vtrap6 ( _threadargscomma_ _lv ) ; _lb = q10_1 * vtrap7 ( _threadargscomma_ _lv ) ; tau_m = 1.0 / ( _la + _lb ) ; m_inf = _la / ( _la + _lb ) ; _la = q10_2 * vtrap8 ( _threadargscomma_ _lv ) ; _lb = q10_2 * bhA / ( 1.0 + Exp ( _threadargscomma_ - ( _lv + bhB ) / bhC ) ) ; tau_h = 1.0 / ( _la + _lb ) ; h_inf = _la / ( _la + _lb ) ; _lv2 = _lv - vtraub ; _la = q10_3 * asA / ( Exp ( _threadargscomma_ ( _lv2 + asB ) / asC ) + 1.0 ) ; _lb = q10_3 * bsA / ( Exp ( _threadargscomma_ ( _lv2 + bsB ) / bsC ) + 1.0 ) ; tau_s = 1.0 / ( _la + _lb ) ; s_inf = _la / ( _la + _lb ) ; return 0; }
static int _f_rates ( _threadargsprotocomma_ double _lv ) { double _lalpha , _lbeta , _ltau , _linf , _lgamma , _lzeta , _ltemp_adj_m , _lA_alpha_m , _lB_alpha_m , _lVhalf_alpha_m , _lA_beta_m , _lB_beta_m , _lVhalf_beta_m , _ltemp_adj_h , _lA_alpha_h , _lB_alpha_h , _lVhalf_alpha_h , _lA_beta_h , _lB_beta_h , _lVhalf_beta_h , _lA_inf_h , _lB_inf_h , _lVhalf_inf_h ; _ltemp_adj_m = pow( 2.0 , ( ( celsius - 24.0 ) / 10.0 ) ) ; _ltemp_adj_h = pow( 2.0 , ( ( celsius - 24.0 ) / 10.0 ) ) ; _lA_alpha_m = 2.880000018 ; _lB_alpha_m = 7.2 ; _lVhalf_alpha_m = - 20.0 ; _lalpha = _lA_alpha_m * vtrap ( _threadargscomma_ ( _lv - _lVhalf_alpha_m ) , _lB_alpha_m ) ; _lA_beta_m = 0.892800005 ; _lB_beta_m = - 7.2 ; _lVhalf_beta_m = - 20.0 ; _lbeta = _lA_beta_m * vtrap ( _threadargscomma_ ( _lv - _lVhalf_beta_m ) , _lB_beta_m ) ; if ( 1.0 / ( ( _lalpha + _lbeta ) * _ltemp_adj_m ) < 0.02 ) { _ltau = ( 0.02 * _ltemp_adj_m ) ; } else { _ltau = 1.0 / ( _lalpha + _lbeta ) ; } mtau = _ltau / _ltemp_adj_m ; minf = _lalpha / ( _lalpha + _lbeta ) ; _lA_alpha_h = 0.045 ; _lB_alpha_h = 1.5 ; _lVhalf_alpha_h = - 35.0 ; _lalpha = _lA_alpha_h * vtrap ( _threadargscomma_ ( _lv - _lVhalf_alpha_h ) , _lB_alpha_h ) ; _lA_beta_h = 0.015 ; _lB_beta_h = - 1.5 ; _lVhalf_beta_h = - 35.0 ; _lbeta = _lA_beta_h * vtrap ( _threadargscomma_ ( _lv - _lVhalf_beta_h ) , _lB_beta_h ) ; if ( 1.0 / ( ( _lalpha + _lbeta ) * _ltemp_adj_h ) < 0.5 ) { _ltau = ( 0.5 * _ltemp_adj_h ) ; } else { _ltau = 1.0 / ( _lalpha + _lbeta ) ; } htau = _ltau / _ltemp_adj_h ; _lA_inf_h = 1.0 ; _lB_inf_h = 4.0 ; _lVhalf_inf_h = - 40.0 ; _linf = _lA_inf_h / ( exp ( ( _lv - _lVhalf_inf_h ) / _lB_inf_h ) + 1.0 ) ; hinf = _linf ; return 0; }