Add test case

Author: kko27

Code/Source/solver/CepModTtp.cpp

@@ -151,14 +151,6 @@ void CepModTtp::getf(const int i, const int nX, const int nG, const Vector<doubl
 
   // I_Na: Fast sodium current
   double I_Na = G_Na * pow(m,3.0) * h * j * (V - E_Na);
-  
-  // // Debug print for G_Na usage
-  // static int debug_counter = 0;
-  // if (debug_counter % 50 == 0) {  // Only print every 50 times to avoid spam
-  //   std::cout << "[DEBUG] TTP getf: G_Na=" << G_Na << ", I_Na=" << I_Na 
-  //             << ", V=" << V << ", m=" << m << ", h=" << h << ", j=" << j << std::endl;
-  // }
-  // debug_counter++;
 
   // I_to: transient outward current
   double I_to = G_to[i-1] * r * s * (V - E_K);

Code/Source/solver/cep_ion.cpp

@@ -6,6 +6,7 @@
 #include "all_fun.h"
 #include "post.h"
 #include "utils.h"
+#include <iostream>
 #include <math.h>
 
 namespace cep_ion {
@@ -41,6 +42,21 @@ void cep_init(Simulation* simulation)
       continue;
     }
 
+    // Warn once on rank 0 if any TTP domain uses default initial conditions.
+    bool any_ttp_defaults = false;
+    for (int iDmn = 0; iDmn < eq.nDmn; iDmn++) {
+      if (eq.dmn[iDmn].cep.cepType == ElectrophysiologyModelType::TTP
+          && !eq.dmn[iDmn].cep.ttp_user_initial_state) {
+        any_ttp_defaults = true;
+        break;
+      }
+    }
+    static bool ttp_default_ic_warning_printed = false;
+    if (any_ttp_defaults && !ttp_default_ic_warning_printed && cm.idcm() == 0) {
+      std::cout << "[WARNING] Default initial conditions are being used for the TTP ionic model (zone-based: epi/endo/mid-myo). To set custom initial conditions, specify <TTP_initial_conditions> in the domain XML." << std::endl;
+      ttp_default_ic_warning_printed = true;
+    }
+
     if (com_mod.dmnId.size() != 0) {
       Vector<double> sA(tnNo); 
       Array<double> sF(nXion,tnNo);

Code/Source/solver/read_files.cpp

@@ -1002,12 +1002,7 @@ void read_cep_domain(Simulation* simulation, EquationParameters* eq_params, Doma
 
   // Set Ttp parameters.
   //
-  if (domain_params->G_Na.defined())  { 
-    lDmn.cep.ttp.G_Na = domain_params->G_Na.value(); 
-    // std::cout << "[DEBUG] G_Na assigned for domain " << domain_params->id.value() 
-    //           << " (myocardial zone: " << domain_params->myocardial_zone.value() 
-    //           << ", imyo: " << lDmn.cep.imyo << "): " << lDmn.cep.ttp.G_Na << std::endl;
-  }
+  if (domain_params->G_Na.defined())  { lDmn.cep.ttp.G_Na = domain_params->G_Na.value(); }
   if (domain_params->G_Kr.defined())  { lDmn.cep.ttp.G_Kr = domain_params->G_Kr.value(); }
   if (domain_params->G_Ks.defined())  { lDmn.cep.ttp.G_Ks[lDmn.cep.imyo - 1] = domain_params->G_Ks.value(); }
   if (domain_params->G_to.defined())  { lDmn.cep.ttp.G_to[lDmn.cep.imyo - 1] = domain_params->G_to.value(); }