diff options
Diffstat (limited to 'ML_BSSN/src/ML_BSSN_convertToADMBaseDtLapseShiftBoundary.cc')
-rw-r--r-- | ML_BSSN/src/ML_BSSN_convertToADMBaseDtLapseShiftBoundary.cc | 18 |
1 files changed, 2 insertions, 16 deletions
diff --git a/ML_BSSN/src/ML_BSSN_convertToADMBaseDtLapseShiftBoundary.cc b/ML_BSSN/src/ML_BSSN_convertToADMBaseDtLapseShiftBoundary.cc index 57469de..9c0b3e5 100644 --- a/ML_BSSN/src/ML_BSSN_convertToADMBaseDtLapseShiftBoundary.cc +++ b/ML_BSSN/src/ML_BSSN_convertToADMBaseDtLapseShiftBoundary.cc @@ -35,9 +35,6 @@ extern "C" void ML_BSSN_convertToADMBaseDtLapseShiftBoundary_SelectBCs(CCTK_ARGU if (cctk_iteration % ML_BSSN_convertToADMBaseDtLapseShiftBoundary_calc_every != ML_BSSN_convertToADMBaseDtLapseShiftBoundary_calc_offset) return; CCTK_INT ierr CCTK_ATTRIBUTE_UNUSED = 0; - ierr = Boundary_SelectGroupForBC(cctkGH, CCTK_ALL_FACES, GenericFD_GetBoundaryWidth(cctkGH), -1 /* no table */, "ADMBase::dtlapse","flat"); - if (ierr < 0) - CCTK_WARN(1, "Failed to register flat BC for ADMBase::dtlapse."); ierr = Boundary_SelectGroupForBC(cctkGH, CCTK_ALL_FACES, GenericFD_GetBoundaryWidth(cctkGH), -1 /* no table */, "ADMBase::dtshift","flat"); if (ierr < 0) CCTK_WARN(1, "Failed to register flat BC for ADMBase::dtshift."); @@ -234,8 +231,6 @@ static void ML_BSSN_convertToADMBaseDtLapseShiftBoundary_Body(const cGH* restric /* 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]); @@ -243,7 +238,6 @@ static void ML_BSSN_convertToADMBaseDtLapseShiftBoundary_Body(const cGH* restric 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]); @@ -286,9 +280,6 @@ static void ML_BSSN_convertToADMBaseDtLapseShiftBoundary_Body(const cGH* restric 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; @@ -318,7 +309,6 @@ static void ML_BSSN_convertToADMBaseDtLapseShiftBoundary_Body(const cGH* restric /* 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); @@ -343,16 +333,12 @@ extern "C" void ML_BSSN_convertToADMBaseDtLapseShiftBoundary(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_convertToADMBaseDtLapseShiftBoundary", 9, groups); + "ML_BSSN::ML_shift"}; + GenericFD_AssertGroupStorage(cctkGH, "ML_BSSN_convertToADMBaseDtLapseShiftBoundary", 5, groups); switch (fdOrder) { |