diff options
Diffstat (limited to 'ML_CCZ4/src/ML_CCZ4_convertToADMBaseFakeDtLapseShift.cc')
| -rw-r--r-- | ML_CCZ4/src/ML_CCZ4_convertToADMBaseFakeDtLapseShift.cc | 6 |
1 files changed, 0 insertions, 6 deletions
diff --git a/ML_CCZ4/src/ML_CCZ4_convertToADMBaseFakeDtLapseShift.cc b/ML_CCZ4/src/ML_CCZ4_convertToADMBaseFakeDtLapseShift.cc index 6b4c608..d4cc7e0 100644 --- a/ML_CCZ4/src/ML_CCZ4_convertToADMBaseFakeDtLapseShift.cc +++ b/ML_CCZ4/src/ML_CCZ4_convertToADMBaseFakeDtLapseShift.cc @@ -217,7 +217,6 @@ static void ML_CCZ4_convertToADMBaseFakeDtLapseShift_Body(const cGH* restrict co /* 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]); @@ -270,10 +269,6 @@ static void ML_CCZ4_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(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; @@ -303,7 +298,6 @@ static void ML_CCZ4_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); |