#include <TenTusscher2006OdeSystem.hpp>
Public Member Functions | |
TenTusscher2006OdeSystem (AbstractIvpOdeSolver *pSolver, AbstractStimulusFunction *pIntracellularStimulus) | |
~TenTusscher2006OdeSystem () | |
void | EvaluateYDerivatives (double time, const std::vector< double > &rY, std::vector< double > &rDY) |
void | SetScaleFactorGks (double sfgks) |
void | SetScaleFactorIto (double sfito) |
void | SetScaleFactorGkr (double sfgkr) |
double | GetIIonic () |
Private Member Functions | |
void | VerifyStateVariables () |
Private Attributes | |
double | mScaleFactorGks |
double | mScaleFactorIto |
double | mScaleFactorGkr |
double | L_type_Ca_current_d_gate_alpha_d |
double | L_type_Ca_current_d_gate_beta_d |
double | L_type_Ca_current_d_gate_d_inf |
double | L_type_Ca_current_d_gate_gamma_d |
double | L_type_Ca_current_d_gate_tau_d |
double | L_type_Ca_current_f2_gate_f2_inf |
double | L_type_Ca_current_f2_gate_tau_f2 |
double | L_type_Ca_current_fCass_gate_fCass_inf |
double | L_type_Ca_current_fCass_gate_tau_fCass |
double | L_type_Ca_current_f_gate_f_inf |
double | L_type_Ca_current_f_gate_tau_f |
double | L_type_Ca_current_i_CaL |
double | calcium_background_current_i_b_Ca |
double | calcium_dynamics_Ca_i_bufc |
double | calcium_dynamics_Ca_sr_bufsr |
double | calcium_dynamics_Ca_ss_bufss |
double | calcium_dynamics_O |
double | calcium_dynamics_i_leak |
double | calcium_dynamics_i_rel |
double | calcium_dynamics_i_up |
double | calcium_dynamics_i_xfer |
double | calcium_dynamics_k1 |
double | calcium_dynamics_k2 |
double | calcium_dynamics_kcasr |
double | calcium_pump_current_i_p_Ca |
double | fast_sodium_current_h_gate_alpha_h |
double | fast_sodium_current_h_gate_beta_h |
double | fast_sodium_current_h_gate_h_inf |
double | fast_sodium_current_h_gate_tau_h |
double | fast_sodium_current_i_Na |
double | fast_sodium_current_j_gate_alpha_j |
double | fast_sodium_current_j_gate_beta_j |
double | fast_sodium_current_j_gate_j_inf |
double | fast_sodium_current_j_gate_tau_j |
double | fast_sodium_current_m_gate_alpha_m |
double | fast_sodium_current_m_gate_beta_m |
double | fast_sodium_current_m_gate_m_inf |
double | fast_sodium_current_m_gate_tau_m |
double | inward_rectifier_potassium_current_alpha_K1 |
double | inward_rectifier_potassium_current_beta_K1 |
double | inward_rectifier_potassium_current_i_K1 |
double | inward_rectifier_potassium_current_xK1_inf |
double | potassium_pump_current_i_p_K |
double | rapid_time_dependent_potassium_current_Xr1_gate_alpha_xr1 |
double | rapid_time_dependent_potassium_current_Xr1_gate_beta_xr1 |
double | rapid_time_dependent_potassium_current_Xr1_gate_tau_xr1 |
double | rapid_time_dependent_potassium_current_Xr1_gate_xr1_inf |
double | rapid_time_dependent_potassium_current_Xr2_gate_alpha_xr2 |
double | rapid_time_dependent_potassium_current_Xr2_gate_beta_xr2 |
double | rapid_time_dependent_potassium_current_Xr2_gate_tau_xr2 |
double | rapid_time_dependent_potassium_current_Xr2_gate_xr2_inf |
double | rapid_time_dependent_potassium_current_i_Kr |
double | reversal_potentials_E_Ca |
double | reversal_potentials_E_K |
double | reversal_potentials_E_Ks |
double | reversal_potentials_E_Na |
double | slow_time_dependent_potassium_current_Xs_gate_alpha_xs |
double | slow_time_dependent_potassium_current_Xs_gate_beta_xs |
double | slow_time_dependent_potassium_current_Xs_gate_tau_xs |
double | slow_time_dependent_potassium_current_Xs_gate_xs_inf |
double | slow_time_dependent_potassium_current_i_Ks |
double | sodium_background_current_i_b_Na |
double | sodium_calcium_exchanger_current_i_NaCa |
double | sodium_potassium_pump_current_i_NaK |
double | transient_outward_current_i_to |
double | transient_outward_current_r_gate_r_inf |
double | transient_outward_current_r_gate_tau_r |
double | transient_outward_current_s_gate_s_inf |
double | transient_outward_current_s_gate_tau_s |
Static Private Attributes | |
static const double | L_type_Ca_current_g_CaL = 0.0000398 |
static const double | calcium_background_current_g_bca = 0.000592 |
static const double | calcium_dynamics_Buf_c = 0.2 |
static const double | calcium_dynamics_Buf_sr = 10.0 |
static const double | calcium_dynamics_Buf_ss = 0.4 |
static const double | calcium_dynamics_Ca_o = 2.0 |
static const double | calcium_dynamics_EC = 1.5 |
static const double | calcium_dynamics_K_buf_c = 0.001 |
static const double | calcium_dynamics_K_buf_sr = 0.3 |
static const double | calcium_dynamics_K_buf_ss = 0.00025 |
static const double | calcium_dynamics_K_up = 0.00025 |
static const double | calcium_dynamics_V_leak = 0.00036 |
static const double | calcium_dynamics_V_rel = 0.102 |
static const double | calcium_dynamics_V_sr = 0.001094 |
static const double | calcium_dynamics_V_ss = 0.00005468 |
static const double | calcium_dynamics_V_xfer = 0.0038 |
static const double | calcium_dynamics_Vmax_up = 0.006375 |
static const double | calcium_dynamics_k1_prime = 0.15 |
static const double | calcium_dynamics_k2_prime = 0.045 |
static const double | calcium_dynamics_k3 = 0.06 |
static const double | calcium_dynamics_k4 = 0.005 |
static const double | calcium_dynamics_max_sr = 2.5 |
static const double | calcium_dynamics_min_sr = 1.0 |
static const double | calcium_pump_current_K_pCa = 0.0005 |
static const double | calcium_pump_current_g_pCa = 0.1238 |
static const double | fast_sodium_current_g_Na = 14.838 |
static const double | inward_rectifier_potassium_current_g_K1 = 5.405 |
static const double | membrane_Cm = 0.185 |
static const double | membrane_F = 96485.3415 |
static const double | membrane_R = 8314.472 |
static const double | membrane_T = 310.0 |
static const double | membrane_V_c = 0.016404 |
static const double | potassium_dynamics_K_o = 5.4 |
static const double | potassium_pump_current_g_pK = 0.0146 |
static const double | rapid_time_dependent_potassium_current_g_Kr = 0.153 |
static const double | reversal_potentials_P_kna = 0.03 |
static const double | slow_time_dependent_potassium_current_g_Ks = 0.392 |
static const double | sodium_background_current_g_bna = 0.00029 |
static const double | sodium_calcium_exchanger_current_K_NaCa = 1000.0 |
static const double | sodium_calcium_exchanger_current_K_sat = 0.1 |
static const double | sodium_calcium_exchanger_current_Km_Ca = 1.38 |
static const double | sodium_calcium_exchanger_current_Km_Nai = 87.5 |
static const double | sodium_calcium_exchanger_current_alpha = 2.5 |
static const double | sodium_calcium_exchanger_current_gamma = 0.35 |
static const double | sodium_dynamics_Na_o = 140.0 |
static const double | sodium_potassium_pump_current_K_mNa = 40.0 |
static const double | sodium_potassium_pump_current_K_mk = 1.0 |
static const double | sodium_potassium_pump_current_P_NaK = 2.724 |
static const double | transient_outward_current_g_to = 0.294 |
Definition at line 39 of file TenTusscher2006OdeSystem.hpp.
TenTusscher2006OdeSystem::TenTusscher2006OdeSystem | ( | AbstractIvpOdeSolver * | pSolver, | |
AbstractStimulusFunction * | pIntracellularStimulus | |||
) |
Constructor
Definition at line 36 of file TenTusscher2006OdeSystem.cpp.
References AbstractCardiacCell::Init(), OdeSystemInformation< ODE_SYSTEM >::Instance(), AbstractOdeSystem::mpSystemInfo, mScaleFactorGkr, mScaleFactorGks, and mScaleFactorIto.
TenTusscher2006OdeSystem::~TenTusscher2006OdeSystem | ( | void | ) |
Destructor
Definition at line 52 of file TenTusscher2006OdeSystem.cpp.
void TenTusscher2006OdeSystem::VerifyStateVariables | ( | ) | [private, virtual] |
Empty method which can be over-ridden in concrete cell class which should go through the current state vector and go range checking on the values (eg check that concentrations are positive and gating variables are between zero and one). This method is called in the ComputeExceptVoltage method.
Reimplemented from AbstractCardiacCell.
Definition at line 406 of file TenTusscher2006OdeSystem.cpp.
References AbstractOdeSystem::DumpState().
Referenced by EvaluateYDerivatives().
void TenTusscher2006OdeSystem::EvaluateYDerivatives | ( | double | time, | |
const std::vector< double > & | rY, | |||
std::vector< double > & | rDY | |||
) | [virtual] |
This method will compute the RHS of the TenTusscher model
Fill in a vector representing the RHS of the TenTusscher2006 system of Odes at each time step, y' = [y1' ... yn']. Some ODE solver will call this function repeatedly to solve for y = [y1 ... yn].
time | the current time, in milliseconds | |
rY | current values of the state variables | |
rDY | to be filled in with derivatives |
Implements AbstractOdeSystem.
Definition at line 77 of file TenTusscher2006OdeSystem.cpp.
References mScaleFactorGkr, mScaleFactorGks, mScaleFactorIto, and VerifyStateVariables().
void TenTusscher2006OdeSystem::SetScaleFactorGks | ( | double | sfgks | ) |
Set the scale factor for Gks in order to differentiate epi M and endo cells
Definition at line 56 of file TenTusscher2006OdeSystem.cpp.
References mScaleFactorGks.
void TenTusscher2006OdeSystem::SetScaleFactorIto | ( | double | sfito | ) |
Set the scale factor for Gks in order to differentiate epi M and endo cells
Definition at line 60 of file TenTusscher2006OdeSystem.cpp.
References mScaleFactorIto.
void TenTusscher2006OdeSystem::SetScaleFactorGkr | ( | double | sfgkr | ) |
Set the scale factor for Gkr (might be useful for simulating LQT2)
Definition at line 64 of file TenTusscher2006OdeSystem.cpp.
References mScaleFactorGkr.
double TenTusscher2006OdeSystem::GetIIonic | ( | ) | [virtual] |
Returns the ionic current
Implements AbstractCardiacCell.
Definition at line 252 of file TenTusscher2006OdeSystem.cpp.
References mScaleFactorGkr, mScaleFactorGks, and mScaleFactorIto.
double TenTusscher2006OdeSystem::mScaleFactorGks [private] |
Scale factor for Gks
Definition at line 44 of file TenTusscher2006OdeSystem.hpp.
Referenced by EvaluateYDerivatives(), GetIIonic(), SetScaleFactorGks(), and TenTusscher2006OdeSystem().
double TenTusscher2006OdeSystem::mScaleFactorIto [private] |
Scale factor for Gto
Definition at line 46 of file TenTusscher2006OdeSystem.hpp.
Referenced by EvaluateYDerivatives(), GetIIonic(), SetScaleFactorIto(), and TenTusscher2006OdeSystem().
double TenTusscher2006OdeSystem::mScaleFactorGkr [private] |
Scale factor for Gkr
Definition at line 48 of file TenTusscher2006OdeSystem.hpp.
Referenced by EvaluateYDerivatives(), GetIIonic(), SetScaleFactorGkr(), and TenTusscher2006OdeSystem().