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; }
