diff options
Diffstat (limited to 'ML_CCZ4/src/ML_CCZ4_convertToADMBaseDtLapseShift.cc')
-rw-r--r-- | ML_CCZ4/src/ML_CCZ4_convertToADMBaseDtLapseShift.cc | 50 |
1 files changed, 0 insertions, 50 deletions
diff --git a/ML_CCZ4/src/ML_CCZ4_convertToADMBaseDtLapseShift.cc b/ML_CCZ4/src/ML_CCZ4_convertToADMBaseDtLapseShift.cc index 4bdbdc9..c8b1cb8 100644 --- a/ML_CCZ4/src/ML_CCZ4_convertToADMBaseDtLapseShift.cc +++ b/ML_CCZ4/src/ML_CCZ4_convertToADMBaseDtLapseShift.cc @@ -234,7 +234,6 @@ static void ML_CCZ4_convertToADMBaseDtLapseShift_Body(const cGH* restrict const /* 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]); @@ -880,10 +879,6 @@ static void ML_CCZ4_convertToADMBaseDtLapseShift_Body(const cGH* restrict const CCTK_REAL_VEC Ddetgt3 CCTK_ATTRIBUTE_UNUSED = kmadd(gtu11,JacPDstandardNth3gt11,kmadd(gtu22,JacPDstandardNth3gt22,kmadd(gtu33,JacPDstandardNth3gt33,kmul(kmadd(gtu12,JacPDstandardNth3gt12,kmadd(gtu13,JacPDstandardNth3gt13,kmul(gtu23,JacPDstandardNth3gt23))),ToReal(2))))); - CCTK_REAL_VEC dtalpL CCTK_ATTRIBUTE_UNUSED = - kmsub(kmadd(beta1L,JacPDupwindNthAnti1alpha,kmadd(beta2L,JacPDupwindNthAnti2alpha,kmadd(beta3L,JacPDupwindNthAnti3alpha,kmadd(JacPDupwindNthSymm1alpha,kfabs(beta1L),kmadd(JacPDupwindNthSymm2alpha,kfabs(beta2L),kmul(JacPDupwindNthSymm3alpha,kfabs(beta3L))))))),ToReal(LapseAdvectionCoeff),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 = kmadd(kmadd(beta1L,JacPDupwindNthAnti1beta1,kmadd(beta2L,JacPDupwindNthAnti2beta1,kmadd(beta3L,JacPDupwindNthAnti3beta1,kmadd(JacPDupwindNthSymm1beta1,kfabs(beta1L),kmadd(JacPDupwindNthSymm2beta1,kfabs(beta2L),kmul(JacPDupwindNthSymm3beta1,kfabs(beta3L))))))),ToReal(ShiftAdvectionCoeff),IfThen(harmonicShift,kmul(alphaL,kmul(em4phi,kmul(ToReal(-0.5),kmadd(gtu11,kmadd(alphaL,kmadd(kmadd(gtu11,JacPDstandardNth1gt11,kmadd(gtu12,kadd(JacPDstandardNth1gt12,JacPDstandardNth2gt11),kmadd(gtu22,JacPDstandardNth2gt12,kmadd(gtu13,kadd(JacPDstandardNth1gt13,JacPDstandardNth3gt11),kmadd(gtu23,kadd(JacPDstandardNth2gt13,JacPDstandardNth3gt12),kmadd(gtu33,JacPDstandardNth3gt13,kmul(JacPDstandardNth1phi,IfThen(conformalMethod == @@ -916,7 +911,6 @@ static void ML_CCZ4_convertToADMBaseDtLapseShift_Body(const cGH* restrict const /* 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); @@ -940,50 +934,6 @@ extern "C" void ML_CCZ4_convertToADMBaseDtLapseShift(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_log_confac", - "ML_CCZ4::ML_metric", - "ML_CCZ4::ML_shift", - "ML_CCZ4::ML_Theta", - "ML_CCZ4::ML_trace_curv"}; - GenericFD_AssertGroupStorage(cctkGH, "ML_CCZ4_convertToADMBaseDtLapseShift", 12, groups); - - switch (fdOrder) - { - case 2: - { - GenericFD_EnsureStencilFits(cctkGH, "ML_CCZ4_convertToADMBaseDtLapseShift", 2, 2, 2); - break; - } - - case 4: - { - GenericFD_EnsureStencilFits(cctkGH, "ML_CCZ4_convertToADMBaseDtLapseShift", 3, 3, 3); - break; - } - - case 6: - { - GenericFD_EnsureStencilFits(cctkGH, "ML_CCZ4_convertToADMBaseDtLapseShift", 4, 4, 4); - break; - } - - case 8: - { - GenericFD_EnsureStencilFits(cctkGH, "ML_CCZ4_convertToADMBaseDtLapseShift", 5, 5, 5); - break; - } - default: - CCTK_BUILTIN_UNREACHABLE(); - } - GenericFD_LoopOverInterior(cctkGH, ML_CCZ4_convertToADMBaseDtLapseShift_Body); if (verbose > 1) |