diff options
Diffstat (limited to 'ML_CCZ4/src/ML_CCZ4_convertToADMBaseDtLapseShiftBoundary.cc')
| -rw-r--r-- | ML_CCZ4/src/ML_CCZ4_convertToADMBaseDtLapseShiftBoundary.cc | 43 |
1 files changed, 0 insertions, 43 deletions
diff --git a/ML_CCZ4/src/ML_CCZ4_convertToADMBaseDtLapseShiftBoundary.cc b/ML_CCZ4/src/ML_CCZ4_convertToADMBaseDtLapseShiftBoundary.cc index d36ef7f..0776724 100644 --- a/ML_CCZ4/src/ML_CCZ4_convertToADMBaseDtLapseShiftBoundary.cc +++ b/ML_CCZ4/src/ML_CCZ4_convertToADMBaseDtLapseShiftBoundary.cc @@ -234,7 +234,6 @@ static void ML_CCZ4_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]); @@ -287,10 +286,6 @@ static void ML_CCZ4_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(AL,ToReal(LapseACoeff),kmul(kmsub(ThetaL,ToReal(2),trKL),ToReal(-1 - + LapseACoeff)))))); - CCTK_REAL_VEC dtbetaxL CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC dtbetayL CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC dtbetazL CCTK_ATTRIBUTE_UNUSED; @@ -320,7 +315,6 @@ static void ML_CCZ4_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); @@ -344,43 +338,6 @@ extern "C" void ML_CCZ4_convertToADMBaseDtLapseShiftBoundary(CCTK_ARGUMENTS) return; } - const char* const groups[] = { - "ADMBase::dtlapse", - "ADMBase::dtshift", - "grid::coordinates", - "ML_CCZ4::ML_dtlapse", - "ML_CCZ4::ML_dtshift", - "ML_CCZ4::ML_Gamma", - "ML_CCZ4::ML_lapse", - "ML_CCZ4::ML_shift", - "ML_CCZ4::ML_Theta", - "ML_CCZ4::ML_trace_curv"}; - GenericFD_AssertGroupStorage(cctkGH, "ML_CCZ4_convertToADMBaseDtLapseShiftBoundary", 10, groups); - - switch (fdOrder) - { - case 2: - { - break; - } - - case 4: - { - break; - } - - case 6: - { - break; - } - - case 8: - { - break; - } - default: - CCTK_BUILTIN_UNREACHABLE(); - } GenericFD_LoopOverBoundaryWithGhosts(cctkGH, ML_CCZ4_convertToADMBaseDtLapseShiftBoundary_Body); |