@@ -57,7 +57,7 @@ void cmm_3d(ComMod& com_mod, const int eNoN, const double w, const Vector<double
5757
5858 int cEq = com_mod.cEq ;
5959 auto & eq = com_mod.eq [cEq];
60- int cDmn = com_mod.cDmn ;
60+ const auto cDmn = com_mod.cDmn ;
6161 auto & dmn = eq.dmn [cDmn];
6262 const double dt = com_mod.dt ;
6363
@@ -450,7 +450,7 @@ void cmm_mass(ComMod& com_mod, const double w, const Vector<double>& N, const Ar
450450 const int cEq = com_mod.cEq ;
451451 const auto & eq = com_mod.eq [cEq];
452452 const double dt = com_mod.dt ;
453- const int cDmn = com_mod.cDmn ;
453+ const auto cDmn = com_mod.cDmn ;
454454 #ifdef debug_cmm_mass
455455 dmsg << " nsd: " << nsd;
456456 dmsg << " dof: " << dof;
@@ -523,7 +523,7 @@ void cmm_stiffness(ComMod& com_mod, const Array<double>& Nxi, const Array<double
523523 const int cEq = com_mod.cEq ;
524524 const auto & eq = com_mod.eq [cEq];
525525 const double dt = com_mod.dt ;
526- const int cDmn = com_mod.cDmn ;
526+ const auto cDmn = com_mod.cDmn ;
527527
528528 double nu = eq.dmn [cDmn].prop .at (PhysicalProperyType::poisson_ratio);
529529 double ht, elM;
@@ -797,7 +797,7 @@ void construct_cmm(ComMod& com_mod, const mshType& lM, const Array<double>& Ag,
797797 const int dof = com_mod.dof ;
798798 const int cEq = com_mod.cEq ;
799799 const auto & eq = com_mod.eq [cEq];
800- auto cDmn = com_mod.cDmn ;
800+ auto & cDmn = com_mod.cDmn ;
801801 const int nsymd = com_mod.nsymd ;
802802 auto & pS0 = com_mod.pS0 ;
803803 auto & pSn = com_mod.pSn ;
@@ -827,7 +827,7 @@ void construct_cmm(ComMod& com_mod, const mshType& lM, const Array<double>& Ag,
827827
828828
829829 for (int e = 0 ; e < lM.nEl ; e++) {
830- // Update domain and proceed if domain phys and eqn phys match
830+ // Change the current domain which will be used in later function calls.
831831 cDmn = all_fun::domain (com_mod, lM, cEq, e);
832832 auto cPhys = eq.dmn [cDmn].phys ;
833833 if (cPhys != EquationType::phys_CMM) {
0 commit comments