diff options
Diffstat (limited to 'ML_BSSN/src/ML_BSSN_Dissipation.cc')
| -rw-r--r-- | ML_BSSN/src/ML_BSSN_Dissipation.cc | 90 |
1 files changed, 1 insertions, 89 deletions
diff --git a/ML_BSSN/src/ML_BSSN_Dissipation.cc b/ML_BSSN/src/ML_BSSN_Dissipation.cc index 2f06644..ac22c87 100644 --- a/ML_BSSN/src/ML_BSSN_Dissipation.cc +++ b/ML_BSSN/src/ML_BSSN_Dissipation.cc @@ -38,18 +38,12 @@ extern "C" void ML_BSSN_Dissipation_SelectBCs(CCTK_ARGUMENTS) ierr = Boundary_SelectGroupForBC(cctkGH, CCTK_ALL_FACES, GenericFD_GetBoundaryWidth(cctkGH), -1 /* no table */, "ML_BSSN::ML_curvrhs","flat"); if (ierr < 0) CCTK_WARN(1, "Failed to register flat BC for ML_BSSN::ML_curvrhs."); - ierr = Boundary_SelectGroupForBC(cctkGH, CCTK_ALL_FACES, GenericFD_GetBoundaryWidth(cctkGH), -1 /* no table */, "ML_BSSN::ML_dtlapserhs","flat"); - if (ierr < 0) - CCTK_WARN(1, "Failed to register flat BC for ML_BSSN::ML_dtlapserhs."); ierr = Boundary_SelectGroupForBC(cctkGH, CCTK_ALL_FACES, GenericFD_GetBoundaryWidth(cctkGH), -1 /* no table */, "ML_BSSN::ML_dtshiftrhs","flat"); if (ierr < 0) CCTK_WARN(1, "Failed to register flat BC for ML_BSSN::ML_dtshiftrhs."); ierr = Boundary_SelectGroupForBC(cctkGH, CCTK_ALL_FACES, GenericFD_GetBoundaryWidth(cctkGH), -1 /* no table */, "ML_BSSN::ML_Gammarhs","flat"); if (ierr < 0) CCTK_WARN(1, "Failed to register flat BC for ML_BSSN::ML_Gammarhs."); - ierr = Boundary_SelectGroupForBC(cctkGH, CCTK_ALL_FACES, GenericFD_GetBoundaryWidth(cctkGH), -1 /* no table */, "ML_BSSN::ML_lapserhs","flat"); - if (ierr < 0) - CCTK_WARN(1, "Failed to register flat BC for ML_BSSN::ML_lapserhs."); ierr = Boundary_SelectGroupForBC(cctkGH, CCTK_ALL_FACES, GenericFD_GetBoundaryWidth(cctkGH), -1 /* no table */, "ML_BSSN::ML_log_confacrhs","flat"); if (ierr < 0) CCTK_WARN(1, "Failed to register flat BC for ML_BSSN::ML_log_confacrhs."); @@ -255,10 +249,6 @@ static void ML_BSSN_Dissipation_Body(const cGH* restrict const cctkGH, const int /* 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 alpharhsL CCTK_ATTRIBUTE_UNUSED = vec_load(alpharhs[index]); - CCTK_REAL_VEC ArhsL CCTK_ATTRIBUTE_UNUSED = vec_load(Arhs[index]); CCTK_REAL_VEC At11L CCTK_ATTRIBUTE_UNUSED = vec_load(At11[index]); CCTK_REAL_VEC At11rhsL CCTK_ATTRIBUTE_UNUSED = vec_load(At11rhs[index]); CCTK_REAL_VEC At12L CCTK_ATTRIBUTE_UNUSED = vec_load(At12[index]); @@ -325,12 +315,6 @@ static void ML_BSSN_Dissipation_Body(const cGH* restrict const cctkGH, const int /* Include user supplied include files */ /* Precompute derivatives */ - CCTK_REAL_VEC PDdissipationNth1A CCTK_ATTRIBUTE_UNUSED; - CCTK_REAL_VEC PDdissipationNth2A CCTK_ATTRIBUTE_UNUSED; - CCTK_REAL_VEC PDdissipationNth3A CCTK_ATTRIBUTE_UNUSED; - CCTK_REAL_VEC PDdissipationNth1alpha CCTK_ATTRIBUTE_UNUSED; - CCTK_REAL_VEC PDdissipationNth2alpha CCTK_ATTRIBUTE_UNUSED; - CCTK_REAL_VEC PDdissipationNth3alpha CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC PDdissipationNth1At11 CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC PDdissipationNth2At11 CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC PDdissipationNth3At11 CCTK_ATTRIBUTE_UNUSED; @@ -405,12 +389,6 @@ static void ML_BSSN_Dissipation_Body(const cGH* restrict const cctkGH, const int { case 2: { - PDdissipationNth1A = PDdissipationNthfdOrder21(&A[index]); - PDdissipationNth2A = PDdissipationNthfdOrder22(&A[index]); - PDdissipationNth3A = PDdissipationNthfdOrder23(&A[index]); - PDdissipationNth1alpha = PDdissipationNthfdOrder21(&alpha[index]); - PDdissipationNth2alpha = PDdissipationNthfdOrder22(&alpha[index]); - PDdissipationNth3alpha = PDdissipationNthfdOrder23(&alpha[index]); PDdissipationNth1At11 = PDdissipationNthfdOrder21(&At11[index]); PDdissipationNth2At11 = PDdissipationNthfdOrder22(&At11[index]); PDdissipationNth3At11 = PDdissipationNthfdOrder23(&At11[index]); @@ -485,12 +463,6 @@ static void ML_BSSN_Dissipation_Body(const cGH* restrict const cctkGH, const int case 4: { - PDdissipationNth1A = PDdissipationNthfdOrder41(&A[index]); - PDdissipationNth2A = PDdissipationNthfdOrder42(&A[index]); - PDdissipationNth3A = PDdissipationNthfdOrder43(&A[index]); - PDdissipationNth1alpha = PDdissipationNthfdOrder41(&alpha[index]); - PDdissipationNth2alpha = PDdissipationNthfdOrder42(&alpha[index]); - PDdissipationNth3alpha = PDdissipationNthfdOrder43(&alpha[index]); PDdissipationNth1At11 = PDdissipationNthfdOrder41(&At11[index]); PDdissipationNth2At11 = PDdissipationNthfdOrder42(&At11[index]); PDdissipationNth3At11 = PDdissipationNthfdOrder43(&At11[index]); @@ -565,12 +537,6 @@ static void ML_BSSN_Dissipation_Body(const cGH* restrict const cctkGH, const int case 6: { - PDdissipationNth1A = PDdissipationNthfdOrder61(&A[index]); - PDdissipationNth2A = PDdissipationNthfdOrder62(&A[index]); - PDdissipationNth3A = PDdissipationNthfdOrder63(&A[index]); - PDdissipationNth1alpha = PDdissipationNthfdOrder61(&alpha[index]); - PDdissipationNth2alpha = PDdissipationNthfdOrder62(&alpha[index]); - PDdissipationNth3alpha = PDdissipationNthfdOrder63(&alpha[index]); PDdissipationNth1At11 = PDdissipationNthfdOrder61(&At11[index]); PDdissipationNth2At11 = PDdissipationNthfdOrder62(&At11[index]); PDdissipationNth3At11 = PDdissipationNthfdOrder63(&At11[index]); @@ -645,12 +611,6 @@ static void ML_BSSN_Dissipation_Body(const cGH* restrict const cctkGH, const int case 8: { - PDdissipationNth1A = PDdissipationNthfdOrder81(&A[index]); - PDdissipationNth2A = PDdissipationNthfdOrder82(&A[index]); - PDdissipationNth3A = PDdissipationNthfdOrder83(&A[index]); - PDdissipationNth1alpha = PDdissipationNthfdOrder81(&alpha[index]); - PDdissipationNth2alpha = PDdissipationNthfdOrder82(&alpha[index]); - PDdissipationNth3alpha = PDdissipationNthfdOrder83(&alpha[index]); PDdissipationNth1At11 = PDdissipationNthfdOrder81(&At11[index]); PDdissipationNth2At11 = PDdissipationNthfdOrder82(&At11[index]); PDdissipationNth3At11 = PDdissipationNthfdOrder83(&At11[index]); @@ -727,8 +687,6 @@ static void ML_BSSN_Dissipation_Body(const cGH* restrict const cctkGH, const int } /* Calculate temporaries and grid functions */ - CCTK_REAL_VEC JacPDdissipationNth1A CCTK_ATTRIBUTE_UNUSED; - CCTK_REAL_VEC JacPDdissipationNth1alpha CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC JacPDdissipationNth1At11 CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC JacPDdissipationNth1At12 CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC JacPDdissipationNth1At13 CCTK_ATTRIBUTE_UNUSED; @@ -752,8 +710,6 @@ static void ML_BSSN_Dissipation_Body(const cGH* restrict const cctkGH, const int CCTK_REAL_VEC JacPDdissipationNth1Xt1 CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC JacPDdissipationNth1Xt2 CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC JacPDdissipationNth1Xt3 CCTK_ATTRIBUTE_UNUSED; - CCTK_REAL_VEC JacPDdissipationNth2A CCTK_ATTRIBUTE_UNUSED; - CCTK_REAL_VEC JacPDdissipationNth2alpha CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC JacPDdissipationNth2At11 CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC JacPDdissipationNth2At12 CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC JacPDdissipationNth2At13 CCTK_ATTRIBUTE_UNUSED; @@ -777,8 +733,6 @@ static void ML_BSSN_Dissipation_Body(const cGH* restrict const cctkGH, const int CCTK_REAL_VEC JacPDdissipationNth2Xt1 CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC JacPDdissipationNth2Xt2 CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC JacPDdissipationNth2Xt3 CCTK_ATTRIBUTE_UNUSED; - CCTK_REAL_VEC JacPDdissipationNth3A CCTK_ATTRIBUTE_UNUSED; - CCTK_REAL_VEC JacPDdissipationNth3alpha CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC JacPDdissipationNth3At11 CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC JacPDdissipationNth3At12 CCTK_ATTRIBUTE_UNUSED; CCTK_REAL_VEC JacPDdissipationNth3At13 CCTK_ATTRIBUTE_UNUSED; @@ -805,12 +759,6 @@ static void ML_BSSN_Dissipation_Body(const cGH* restrict const cctkGH, const int if (use_jacobian) { - JacPDdissipationNth1A = - kmadd(J11L,PDdissipationNth1A,kmadd(J21L,PDdissipationNth2A,kmul(J31L,PDdissipationNth3A))); - - JacPDdissipationNth1alpha = - kmadd(J11L,PDdissipationNth1alpha,kmadd(J21L,PDdissipationNth2alpha,kmul(J31L,PDdissipationNth3alpha))); - JacPDdissipationNth1At11 = kmadd(J11L,PDdissipationNth1At11,kmadd(J21L,PDdissipationNth2At11,kmul(J31L,PDdissipationNth3At11))); @@ -880,12 +828,6 @@ static void ML_BSSN_Dissipation_Body(const cGH* restrict const cctkGH, const int JacPDdissipationNth1Xt3 = kmadd(J11L,PDdissipationNth1Xt3,kmadd(J21L,PDdissipationNth2Xt3,kmul(J31L,PDdissipationNth3Xt3))); - JacPDdissipationNth2A = - kmadd(J12L,PDdissipationNth1A,kmadd(J22L,PDdissipationNth2A,kmul(J32L,PDdissipationNth3A))); - - JacPDdissipationNth2alpha = - kmadd(J12L,PDdissipationNth1alpha,kmadd(J22L,PDdissipationNth2alpha,kmul(J32L,PDdissipationNth3alpha))); - JacPDdissipationNth2At11 = kmadd(J12L,PDdissipationNth1At11,kmadd(J22L,PDdissipationNth2At11,kmul(J32L,PDdissipationNth3At11))); @@ -955,12 +897,6 @@ static void ML_BSSN_Dissipation_Body(const cGH* restrict const cctkGH, const int JacPDdissipationNth2Xt3 = kmadd(J12L,PDdissipationNth1Xt3,kmadd(J22L,PDdissipationNth2Xt3,kmul(J32L,PDdissipationNth3Xt3))); - JacPDdissipationNth3A = - kmadd(J13L,PDdissipationNth1A,kmadd(J23L,PDdissipationNth2A,kmul(J33L,PDdissipationNth3A))); - - JacPDdissipationNth3alpha = - kmadd(J13L,PDdissipationNth1alpha,kmadd(J23L,PDdissipationNth2alpha,kmul(J33L,PDdissipationNth3alpha))); - JacPDdissipationNth3At11 = kmadd(J13L,PDdissipationNth1At11,kmadd(J23L,PDdissipationNth2At11,kmul(J33L,PDdissipationNth3At11))); @@ -1032,10 +968,6 @@ static void ML_BSSN_Dissipation_Body(const cGH* restrict const cctkGH, const int } else { - JacPDdissipationNth1A = PDdissipationNth1A; - - JacPDdissipationNth1alpha = PDdissipationNth1alpha; - JacPDdissipationNth1At11 = PDdissipationNth1At11; JacPDdissipationNth1At12 = PDdissipationNth1At12; @@ -1082,10 +1014,6 @@ static void ML_BSSN_Dissipation_Body(const cGH* restrict const cctkGH, const int JacPDdissipationNth1Xt3 = PDdissipationNth1Xt3; - JacPDdissipationNth2A = PDdissipationNth2A; - - JacPDdissipationNth2alpha = PDdissipationNth2alpha; - JacPDdissipationNth2At11 = PDdissipationNth2At11; JacPDdissipationNth2At12 = PDdissipationNth2At12; @@ -1132,10 +1060,6 @@ static void ML_BSSN_Dissipation_Body(const cGH* restrict const cctkGH, const int JacPDdissipationNth2Xt3 = PDdissipationNth2Xt3; - JacPDdissipationNth3A = PDdissipationNth3A; - - JacPDdissipationNth3alpha = PDdissipationNth3alpha; - JacPDdissipationNth3At11 = PDdissipationNth3At11; JacPDdissipationNth3At12 = PDdissipationNth3At12; @@ -1240,12 +1164,6 @@ static void ML_BSSN_Dissipation_Body(const cGH* restrict const cctkGH, const int At33rhsL = kadd(At33rhsL,kmadd(epsdiss1,JacPDdissipationNth1At33,kmadd(epsdiss2,JacPDdissipationNth2At33,kmul(epsdiss3,JacPDdissipationNth3At33)))); - alpharhsL = - kadd(alpharhsL,kmadd(epsdiss1,JacPDdissipationNth1alpha,kmadd(epsdiss2,JacPDdissipationNth2alpha,kmul(epsdiss3,JacPDdissipationNth3alpha)))); - - ArhsL = - kadd(ArhsL,kmadd(epsdiss1,JacPDdissipationNth1A,kmadd(epsdiss2,JacPDdissipationNth2A,kmul(epsdiss3,JacPDdissipationNth3A)))); - beta1rhsL = kadd(beta1rhsL,kmadd(epsdiss1,JacPDdissipationNth1beta1,kmadd(epsdiss2,JacPDdissipationNth2beta1,kmul(epsdiss3,JacPDdissipationNth3beta1)))); @@ -1266,8 +1184,6 @@ static void ML_BSSN_Dissipation_Body(const cGH* restrict const cctkGH, const int /* Copy local copies back to grid functions */ vec_store_partial_prepare(i,vecimin,vecimax); - vec_store_nta_partial(alpharhs[index],alpharhsL); - vec_store_nta_partial(Arhs[index],ArhsL); vec_store_nta_partial(At11rhs[index],At11rhsL); vec_store_nta_partial(At12rhs[index],At12rhsL); vec_store_nta_partial(At13rhs[index],At13rhsL); @@ -1314,14 +1230,10 @@ extern "C" void ML_BSSN_Dissipation(CCTK_ARGUMENTS) const char* const groups[] = { "ML_BSSN::ML_curv", "ML_BSSN::ML_curvrhs", - "ML_BSSN::ML_dtlapse", - "ML_BSSN::ML_dtlapserhs", "ML_BSSN::ML_dtshift", "ML_BSSN::ML_dtshiftrhs", "ML_BSSN::ML_Gamma", "ML_BSSN::ML_Gammarhs", - "ML_BSSN::ML_lapse", - "ML_BSSN::ML_lapserhs", "ML_BSSN::ML_log_confac", "ML_BSSN::ML_log_confacrhs", "ML_BSSN::ML_metric", @@ -1330,7 +1242,7 @@ extern "C" void ML_BSSN_Dissipation(CCTK_ARGUMENTS) "ML_BSSN::ML_shiftrhs", "ML_BSSN::ML_trace_curv", "ML_BSSN::ML_trace_curvrhs"}; - GenericFD_AssertGroupStorage(cctkGH, "ML_BSSN_Dissipation", 18, groups); + GenericFD_AssertGroupStorage(cctkGH, "ML_BSSN_Dissipation", 14, groups); switch (fdOrder) { |