Example #1
0
 /*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;
}
Example #2
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;
}
Example #3
0
 /*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;
}
Example #4
0
/*CVODE*/
 static int _ode_spec1 () {_reset=0;
 {
   evaluate_fct ( _threadargscomma_ v ) ;
   Dm = ( m_inf - m ) / taum ;
   }
 return _reset;
}
Example #5
0
 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;
}
Example #6
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;
}
Example #8
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;
}
Example #9
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; }
Example #10
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;
}
Example #13
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 ;
   }
 
}
}
Example #14
0
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;

}
}
Example #15
0
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;

}
}
Example #17
0
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;

}
}
Example #18
0
static int _hoc_evaluate_fct() {
  double _r;
   _r = 1.;
 evaluate_fct (  *getarg(1) );
 ret(_r);
}
Example #19
0
 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);
}