aboutsummaryrefslogtreecommitdiff
path: root/ML_BSSN/src/ML_BSSN_Dissipation.cc
diff options
context:
space:
mode:
Diffstat (limited to 'ML_BSSN/src/ML_BSSN_Dissipation.cc')
-rw-r--r--ML_BSSN/src/ML_BSSN_Dissipation.cc46
1 files changed, 1 insertions, 45 deletions
diff --git a/ML_BSSN/src/ML_BSSN_Dissipation.cc b/ML_BSSN/src/ML_BSSN_Dissipation.cc
index 40ff98c..98378a0 100644
--- a/ML_BSSN/src/ML_BSSN_Dissipation.cc
+++ b/ML_BSSN/src/ML_BSSN_Dissipation.cc
@@ -53,9 +53,6 @@ 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_shiftrhs","flat");
if (ierr < 0)
CCTK_WARN(1, "Failed to register flat BC for ML_BSSN::ML_shiftrhs.");
- ierr = Boundary_SelectGroupForBC(cctkGH, CCTK_ALL_FACES, GenericFD_GetBoundaryWidth(cctkGH), -1 /* no table */, "ML_BSSN::ML_Thetarhs","flat");
- if (ierr < 0)
- CCTK_WARN(1, "Failed to register flat BC for ML_BSSN::ML_Thetarhs.");
ierr = Boundary_SelectGroupForBC(cctkGH, CCTK_ALL_FACES, GenericFD_GetBoundaryWidth(cctkGH), -1 /* no table */, "ML_BSSN::ML_trace_curvrhs","flat");
if (ierr < 0)
CCTK_WARN(1, "Failed to register flat BC for ML_BSSN::ML_trace_curvrhs.");
@@ -280,8 +277,6 @@ static void ML_BSSN_Dissipation_Body(cGH const * restrict const cctkGH, int cons
CCTK_REAL_VEC gt33rhsL = vec_load(gt33rhs[index]);
CCTK_REAL_VEC phiL = vec_load(phi[index]);
CCTK_REAL_VEC phirhsL = vec_load(phirhs[index]);
- CCTK_REAL_VEC ThetaL = vec_load(Theta[index]);
- CCTK_REAL_VEC ThetarhsL = vec_load(Thetarhs[index]);
CCTK_REAL_VEC trKL = vec_load(trK[index]);
CCTK_REAL_VEC trKrhsL = vec_load(trKrhs[index]);
CCTK_REAL_VEC Xt1L = vec_load(Xt1[index]);
@@ -373,9 +368,6 @@ static void ML_BSSN_Dissipation_Body(cGH const * restrict const cctkGH, int cons
CCTK_REAL_VEC PDdissipationNth1phi;
CCTK_REAL_VEC PDdissipationNth2phi;
CCTK_REAL_VEC PDdissipationNth3phi;
- CCTK_REAL_VEC PDdissipationNth1Theta;
- CCTK_REAL_VEC PDdissipationNth2Theta;
- CCTK_REAL_VEC PDdissipationNth3Theta;
CCTK_REAL_VEC PDdissipationNth1trK;
CCTK_REAL_VEC PDdissipationNth2trK;
CCTK_REAL_VEC PDdissipationNth3trK;
@@ -455,9 +447,6 @@ static void ML_BSSN_Dissipation_Body(cGH const * restrict const cctkGH, int cons
PDdissipationNth1phi = PDdissipationNthfdOrder21(&phi[index]);
PDdissipationNth2phi = PDdissipationNthfdOrder22(&phi[index]);
PDdissipationNth3phi = PDdissipationNthfdOrder23(&phi[index]);
- PDdissipationNth1Theta = PDdissipationNthfdOrder21(&Theta[index]);
- PDdissipationNth2Theta = PDdissipationNthfdOrder22(&Theta[index]);
- PDdissipationNth3Theta = PDdissipationNthfdOrder23(&Theta[index]);
PDdissipationNth1trK = PDdissipationNthfdOrder21(&trK[index]);
PDdissipationNth2trK = PDdissipationNthfdOrder22(&trK[index]);
PDdissipationNth3trK = PDdissipationNthfdOrder23(&trK[index]);
@@ -536,9 +525,6 @@ static void ML_BSSN_Dissipation_Body(cGH const * restrict const cctkGH, int cons
PDdissipationNth1phi = PDdissipationNthfdOrder41(&phi[index]);
PDdissipationNth2phi = PDdissipationNthfdOrder42(&phi[index]);
PDdissipationNth3phi = PDdissipationNthfdOrder43(&phi[index]);
- PDdissipationNth1Theta = PDdissipationNthfdOrder41(&Theta[index]);
- PDdissipationNth2Theta = PDdissipationNthfdOrder42(&Theta[index]);
- PDdissipationNth3Theta = PDdissipationNthfdOrder43(&Theta[index]);
PDdissipationNth1trK = PDdissipationNthfdOrder41(&trK[index]);
PDdissipationNth2trK = PDdissipationNthfdOrder42(&trK[index]);
PDdissipationNth3trK = PDdissipationNthfdOrder43(&trK[index]);
@@ -617,9 +603,6 @@ static void ML_BSSN_Dissipation_Body(cGH const * restrict const cctkGH, int cons
PDdissipationNth1phi = PDdissipationNthfdOrder61(&phi[index]);
PDdissipationNth2phi = PDdissipationNthfdOrder62(&phi[index]);
PDdissipationNth3phi = PDdissipationNthfdOrder63(&phi[index]);
- PDdissipationNth1Theta = PDdissipationNthfdOrder61(&Theta[index]);
- PDdissipationNth2Theta = PDdissipationNthfdOrder62(&Theta[index]);
- PDdissipationNth3Theta = PDdissipationNthfdOrder63(&Theta[index]);
PDdissipationNth1trK = PDdissipationNthfdOrder61(&trK[index]);
PDdissipationNth2trK = PDdissipationNthfdOrder62(&trK[index]);
PDdissipationNth3trK = PDdissipationNthfdOrder63(&trK[index]);
@@ -698,9 +681,6 @@ static void ML_BSSN_Dissipation_Body(cGH const * restrict const cctkGH, int cons
PDdissipationNth1phi = PDdissipationNthfdOrder81(&phi[index]);
PDdissipationNth2phi = PDdissipationNthfdOrder82(&phi[index]);
PDdissipationNth3phi = PDdissipationNthfdOrder83(&phi[index]);
- PDdissipationNth1Theta = PDdissipationNthfdOrder81(&Theta[index]);
- PDdissipationNth2Theta = PDdissipationNthfdOrder82(&Theta[index]);
- PDdissipationNth3Theta = PDdissipationNthfdOrder83(&Theta[index]);
PDdissipationNth1trK = PDdissipationNthfdOrder81(&trK[index]);
PDdissipationNth2trK = PDdissipationNthfdOrder82(&trK[index]);
PDdissipationNth3trK = PDdissipationNthfdOrder83(&trK[index]);
@@ -738,7 +718,6 @@ static void ML_BSSN_Dissipation_Body(cGH const * restrict const cctkGH, int cons
CCTK_REAL_VEC JacPDdissipationNth1gt23;
CCTK_REAL_VEC JacPDdissipationNth1gt33;
CCTK_REAL_VEC JacPDdissipationNth1phi;
- CCTK_REAL_VEC JacPDdissipationNth1Theta;
CCTK_REAL_VEC JacPDdissipationNth1trK;
CCTK_REAL_VEC JacPDdissipationNth1Xt1;
CCTK_REAL_VEC JacPDdissipationNth1Xt2;
@@ -764,7 +743,6 @@ static void ML_BSSN_Dissipation_Body(cGH const * restrict const cctkGH, int cons
CCTK_REAL_VEC JacPDdissipationNth2gt23;
CCTK_REAL_VEC JacPDdissipationNth2gt33;
CCTK_REAL_VEC JacPDdissipationNth2phi;
- CCTK_REAL_VEC JacPDdissipationNth2Theta;
CCTK_REAL_VEC JacPDdissipationNth2trK;
CCTK_REAL_VEC JacPDdissipationNth2Xt1;
CCTK_REAL_VEC JacPDdissipationNth2Xt2;
@@ -790,7 +768,6 @@ static void ML_BSSN_Dissipation_Body(cGH const * restrict const cctkGH, int cons
CCTK_REAL_VEC JacPDdissipationNth3gt23;
CCTK_REAL_VEC JacPDdissipationNth3gt33;
CCTK_REAL_VEC JacPDdissipationNth3phi;
- CCTK_REAL_VEC JacPDdissipationNth3Theta;
CCTK_REAL_VEC JacPDdissipationNth3trK;
CCTK_REAL_VEC JacPDdissipationNth3Xt1;
CCTK_REAL_VEC JacPDdissipationNth3Xt2;
@@ -861,9 +838,6 @@ static void ML_BSSN_Dissipation_Body(cGH const * restrict const cctkGH, int cons
JacPDdissipationNth1phi =
kmadd(J11L,PDdissipationNth1phi,kmadd(J21L,PDdissipationNth2phi,kmul(J31L,PDdissipationNth3phi)));
- JacPDdissipationNth1Theta =
- kmadd(J11L,PDdissipationNth1Theta,kmadd(J21L,PDdissipationNth2Theta,kmul(J31L,PDdissipationNth3Theta)));
-
JacPDdissipationNth1trK =
kmadd(J11L,PDdissipationNth1trK,kmadd(J21L,PDdissipationNth2trK,kmul(J31L,PDdissipationNth3trK)));
@@ -939,9 +913,6 @@ static void ML_BSSN_Dissipation_Body(cGH const * restrict const cctkGH, int cons
JacPDdissipationNth2phi =
kmadd(J12L,PDdissipationNth1phi,kmadd(J22L,PDdissipationNth2phi,kmul(J32L,PDdissipationNth3phi)));
- JacPDdissipationNth2Theta =
- kmadd(J12L,PDdissipationNth1Theta,kmadd(J22L,PDdissipationNth2Theta,kmul(J32L,PDdissipationNth3Theta)));
-
JacPDdissipationNth2trK =
kmadd(J12L,PDdissipationNth1trK,kmadd(J22L,PDdissipationNth2trK,kmul(J32L,PDdissipationNth3trK)));
@@ -1017,9 +988,6 @@ static void ML_BSSN_Dissipation_Body(cGH const * restrict const cctkGH, int cons
JacPDdissipationNth3phi =
kmadd(J13L,PDdissipationNth1phi,kmadd(J23L,PDdissipationNth2phi,kmul(J33L,PDdissipationNth3phi)));
- JacPDdissipationNth3Theta =
- kmadd(J13L,PDdissipationNth1Theta,kmadd(J23L,PDdissipationNth2Theta,kmul(J33L,PDdissipationNth3Theta)));
-
JacPDdissipationNth3trK =
kmadd(J13L,PDdissipationNth1trK,kmadd(J23L,PDdissipationNth2trK,kmul(J33L,PDdissipationNth3trK)));
@@ -1076,8 +1044,6 @@ static void ML_BSSN_Dissipation_Body(cGH const * restrict const cctkGH, int cons
JacPDdissipationNth1phi = PDdissipationNth1phi;
- JacPDdissipationNth1Theta = PDdissipationNth1Theta;
-
JacPDdissipationNth1trK = PDdissipationNth1trK;
JacPDdissipationNth1Xt1 = PDdissipationNth1Xt1;
@@ -1128,8 +1094,6 @@ static void ML_BSSN_Dissipation_Body(cGH const * restrict const cctkGH, int cons
JacPDdissipationNth2phi = PDdissipationNth2phi;
- JacPDdissipationNth2Theta = PDdissipationNth2Theta;
-
JacPDdissipationNth2trK = PDdissipationNth2trK;
JacPDdissipationNth2Xt1 = PDdissipationNth2Xt1;
@@ -1180,8 +1144,6 @@ static void ML_BSSN_Dissipation_Body(cGH const * restrict const cctkGH, int cons
JacPDdissipationNth3phi = PDdissipationNth3phi;
- JacPDdissipationNth3Theta = PDdissipationNth3Theta;
-
JacPDdissipationNth3trK = PDdissipationNth3trK;
JacPDdissipationNth3Xt1 = PDdissipationNth3Xt1;
@@ -1227,9 +1189,6 @@ static void ML_BSSN_Dissipation_Body(cGH const * restrict const cctkGH, int cons
Xt3rhsL =
kadd(Xt3rhsL,kmadd(epsdiss1,JacPDdissipationNth1Xt3,kmadd(epsdiss2,JacPDdissipationNth2Xt3,kmul(epsdiss3,JacPDdissipationNth3Xt3))));
- ThetarhsL =
- kadd(ThetarhsL,kmadd(epsdiss1,JacPDdissipationNth1Theta,kmadd(epsdiss2,JacPDdissipationNth2Theta,kmul(epsdiss3,JacPDdissipationNth3Theta))));
-
trKrhsL =
kadd(trKrhsL,kmadd(epsdiss1,JacPDdissipationNth1trK,kmadd(epsdiss2,JacPDdissipationNth2trK,kmul(epsdiss3,JacPDdissipationNth3trK))));
@@ -1298,7 +1257,6 @@ static void ML_BSSN_Dissipation_Body(cGH const * restrict const cctkGH, int cons
vec_store_nta_partial(gt23rhs[index],gt23rhsL);
vec_store_nta_partial(gt33rhs[index],gt33rhsL);
vec_store_nta_partial(phirhs[index],phirhsL);
- vec_store_nta_partial(Thetarhs[index],ThetarhsL);
vec_store_nta_partial(trKrhs[index],trKrhsL);
vec_store_nta_partial(Xt1rhs[index],Xt1rhsL);
vec_store_nta_partial(Xt2rhs[index],Xt2rhsL);
@@ -1340,11 +1298,9 @@ extern "C" void ML_BSSN_Dissipation(CCTK_ARGUMENTS)
"ML_BSSN::ML_metricrhs",
"ML_BSSN::ML_shift",
"ML_BSSN::ML_shiftrhs",
- "ML_BSSN::ML_Theta",
- "ML_BSSN::ML_Thetarhs",
"ML_BSSN::ML_trace_curv",
"ML_BSSN::ML_trace_curvrhs"};
- GenericFD_AssertGroupStorage(cctkGH, "ML_BSSN_Dissipation", 20, groups);
+ GenericFD_AssertGroupStorage(cctkGH, "ML_BSSN_Dissipation", 18, groups);
switch(fdOrder)
{
generated by cgit v1.2.3 (git 2.39.1) at 2024-06-01 01:32:50 +0000