diff options
Diffstat (limited to 'ML_BSSN/src/ML_BSSN_convertToADMBaseFakeDtLapseShift.cc')
-rw-r--r-- | ML_BSSN/src/ML_BSSN_convertToADMBaseFakeDtLapseShift.cc | 17 |
1 files changed, 4 insertions, 13 deletions
diff --git a/ML_BSSN/src/ML_BSSN_convertToADMBaseFakeDtLapseShift.cc b/ML_BSSN/src/ML_BSSN_convertToADMBaseFakeDtLapseShift.cc index f8b7487..2d78427 100644 --- a/ML_BSSN/src/ML_BSSN_convertToADMBaseFakeDtLapseShift.cc +++ b/ML_BSSN/src/ML_BSSN_convertToADMBaseFakeDtLapseShift.cc @@ -214,11 +214,11 @@ static void ML_BSSN_convertToADMBaseFakeDtLapseShift_Body(const cGH* restrict co { const ptrdiff_t index CCTK_ATTRIBUTE_UNUSED = di*i + dj*j + dk*k; // vec_iter_counter+=CCTK_REAL_VEC_SIZE; + if (fabs(y[index]) > 1e-8) + continue; /* Assign local copies of grid functions */ - CCTK_REAL_VEC AL CCTK_ATTRIBUTE_UNUSED = vec_load(A[index]); - CCTK_REAL_VEC alphaL CCTK_ATTRIBUTE_UNUSED = vec_load(alpha[index]); CCTK_REAL_VEC B1L CCTK_ATTRIBUTE_UNUSED = vec_load(B1[index]); CCTK_REAL_VEC B2L CCTK_ATTRIBUTE_UNUSED = vec_load(B2[index]); CCTK_REAL_VEC B3L CCTK_ATTRIBUTE_UNUSED = vec_load(B3[index]); @@ -226,7 +226,6 @@ static void ML_BSSN_convertToADMBaseFakeDtLapseShift_Body(const cGH* restrict co CCTK_REAL_VEC beta2L CCTK_ATTRIBUTE_UNUSED = vec_load(beta2[index]); CCTK_REAL_VEC beta3L CCTK_ATTRIBUTE_UNUSED = vec_load(beta3[index]); CCTK_REAL_VEC rL CCTK_ATTRIBUTE_UNUSED = vec_load(r[index]); - CCTK_REAL_VEC trKL CCTK_ATTRIBUTE_UNUSED = vec_load(trK[index]); CCTK_REAL_VEC Xt1L CCTK_ATTRIBUTE_UNUSED = vec_load(Xt1[index]); CCTK_REAL_VEC Xt2L CCTK_ATTRIBUTE_UNUSED = vec_load(Xt2[index]); CCTK_REAL_VEC Xt3L CCTK_ATTRIBUTE_UNUSED = vec_load(Xt3[index]); @@ -269,9 +268,6 @@ static void ML_BSSN_convertToADMBaseFakeDtLapseShift_Body(const cGH* restrict co CCTK_REAL_VEC theta CCTK_ATTRIBUTE_UNUSED = kfmin(ToReal(1),kexp(knmsub(rL,ToReal(ScalarINV(SpatialShiftGammaCoeffRadius)),ToReal(1)))); - CCTK_REAL_VEC dtalpL CCTK_ATTRIBUTE_UNUSED = - kneg(kmul(kpow(alphaL,harmonicN),kmul(ToReal(harmonicF),kmadd(ksub(AL,trKL),ToReal(LapseACoeff),trKL)))); - CCTK_REAL_VEC dtbetaxL CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC dtbetayL CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC dtbetazL CCTK_ATTRIBUTE_UNUSED; @@ -301,7 +297,6 @@ static void ML_BSSN_convertToADMBaseFakeDtLapseShift_Body(const cGH* restrict co /* Copy local copies back to grid functions */ vec_store_partial_prepare(i,vecimin,vecimax); - vec_store_nta_partial(dtalp[index],dtalpL); vec_store_nta_partial(dtbetax[index],dtbetaxL); vec_store_nta_partial(dtbetay[index],dtbetayL); vec_store_nta_partial(dtbetaz[index],dtbetazL); @@ -326,16 +321,12 @@ extern "C" void ML_BSSN_convertToADMBaseFakeDtLapseShift(CCTK_ARGUMENTS) } const char* const groups[] = { - "ADMBase::dtlapse", "ADMBase::dtshift", "grid::coordinates", - "ML_BSSN::ML_dtlapse", "ML_BSSN::ML_dtshift", "ML_BSSN::ML_Gamma", - "ML_BSSN::ML_lapse", - "ML_BSSN::ML_shift", - "ML_BSSN::ML_trace_curv"}; - GenericFD_AssertGroupStorage(cctkGH, "ML_BSSN_convertToADMBaseFakeDtLapseShift", 9, groups); + "ML_BSSN::ML_shift"}; + GenericFD_AssertGroupStorage(cctkGH, "ML_BSSN_convertToADMBaseFakeDtLapseShift", 5, groups); switch (fdOrder) { |