diff options
Diffstat (limited to 'ML_BSSN/src/ML_BSSN_enforce.cc')
-rw-r--r-- | ML_BSSN/src/ML_BSSN_enforce.cc | 10 |
1 files changed, 4 insertions, 6 deletions
diff --git a/ML_BSSN/src/ML_BSSN_enforce.cc b/ML_BSSN/src/ML_BSSN_enforce.cc index f1483bf..30f7538 100644 --- a/ML_BSSN/src/ML_BSSN_enforce.cc +++ b/ML_BSSN/src/ML_BSSN_enforce.cc @@ -217,7 +217,6 @@ static void ML_BSSN_enforce_Body(const cGH* restrict const cctkGH, const int dir /* Assign local copies of grid functions */ - CCTK_REAL_VEC alphaL CCTK_ATTRIBUTE_UNUSED = vec_load(alpha[index]); CCTK_REAL_VEC At11L CCTK_ATTRIBUTE_UNUSED = vec_load(At11[index]); CCTK_REAL_VEC At12L CCTK_ATTRIBUTE_UNUSED = vec_load(At12[index]); CCTK_REAL_VEC At13L CCTK_ATTRIBUTE_UNUSED = vec_load(At13[index]); @@ -303,12 +302,12 @@ static void ML_BSSN_enforce_Body(const cGH* restrict const cctkGH, const int dir At33L = kmadd(gt33L,kmul(trAt,ToReal(-0.333333333333333333333333333333)),At33L); - - alphaL = kfmax(alphaL,ToReal(MinimumLapse)); + + if (alpha[index] < MinimumLapse) + alpha[index] = MinimumLapse; /* Copy local copies back to grid functions */ vec_store_partial_prepare(i,vecimin,vecimax); - vec_store_nta_partial(alpha[index],alphaL); vec_store_nta_partial(At11[index],At11L); vec_store_nta_partial(At12[index],At12L); vec_store_nta_partial(At13[index],At13L); @@ -337,9 +336,8 @@ extern "C" void ML_BSSN_enforce(CCTK_ARGUMENTS) const char* const groups[] = { "ML_BSSN::ML_curv", - "ML_BSSN::ML_lapse", "ML_BSSN::ML_metric"}; - GenericFD_AssertGroupStorage(cctkGH, "ML_BSSN_enforce", 3, groups); + GenericFD_AssertGroupStorage(cctkGH, "ML_BSSN_enforce", 2, groups); switch (fdOrder) { |