aboutsummaryrefslogtreecommitdiff
path: root/ML_BSSN/src/ML_BSSN_RHS1.cc
diff options
context:
space:
mode:
Diffstat (limited to 'ML_BSSN/src/ML_BSSN_RHS1.cc')
-rw-r--r--ML_BSSN/src/ML_BSSN_RHS1.cc141
1 files changed, 67 insertions, 74 deletions
diff --git a/ML_BSSN/src/ML_BSSN_RHS1.cc b/ML_BSSN/src/ML_BSSN_RHS1.cc
index 0d11a97..f5456e3 100644
--- a/ML_BSSN/src/ML_BSSN_RHS1.cc
+++ b/ML_BSSN/src/ML_BSSN_RHS1.cc
@@ -534,7 +534,7 @@ static void ML_BSSN_RHS1_Body(cGH const * restrict const cctkGH, int const dir,
CCTK_REAL Xtn3 = Gt311*gtu11 + Gt322*gtu22 + 2*(Gt312*gtu12 +
Gt313*gtu13 + Gt323*gtu23) + Gt333*gtu33;
- CCTK_REAL fac1 = IfThen(conformalMethod,-0.5*INV(phiL),1);
+ CCTK_REAL fac1 = IfThen(ToReal(conformalMethod),-0.5*INV(phiL),1);
CCTK_REAL cdphi1 = fac1*PDstandardNth1phi;
@@ -572,7 +572,8 @@ static void ML_BSSN_RHS1_Body(cGH const * restrict const cctkGH, int const dir,
CCTK_REAL Atu33 = Atm31*gtu13 + Atm32*gtu23 + Atm33*gtu33;
- CCTK_REAL e4phi = IfThen(conformalMethod,INV(SQR(phiL)),exp(4*phiL));
+ CCTK_REAL e4phi =
+ IfThen(ToReal(conformalMethod),INV(SQR(phiL)),exp(4*phiL));
CCTK_REAL em4phi = INV(e4phi);
@@ -598,14 +599,11 @@ static void ML_BSSN_RHS1_Body(cGH const * restrict const cctkGH, int const dir,
beta1L*PDupwindNthAnti1phi + beta2L*PDupwindNthAnti2phi +
beta3L*PDupwindNthAnti3phi + PDupwindNthSymm1phi*Abs(beta1L) +
PDupwindNthSymm2phi*Abs(beta2L) + PDupwindNthSymm3phi*Abs(beta3L) +
- IfThen(conformalMethod,-0.333333333333333333333333333333*(PDstandardNth1beta1
- +
- PDstandardNth2beta2)*phiL,0.166666666666666666666666666667*(PDstandardNth1beta1
- + PDstandardNth2beta2)) +
- IfThen(conformalMethod,phiL*(-0.333333333333333333333333333333*PDstandardNth3beta3
- +
- 0.333333333333333333333333333333*alphaL*trKL),0.166666666666666666666666666667*PDstandardNth3beta3
- - 0.166666666666666666666666666667*alphaL*trKL);
+ IfThen(ToReal(conformalMethod),phiL*(-0.333333333333333333333333333333*(PDstandardNth1beta1
+ + PDstandardNth2beta2 + PDstandardNth3beta3) +
+ 0.333333333333333333333333333333*alphaL*trKL),0.166666666666666666666666666667*(PDstandardNth1beta1
+ + PDstandardNth2beta2 + PDstandardNth3beta3) -
+ 0.166666666666666666666666666667*alphaL*trKL);
CCTK_REAL gt11rhsL = -2*alphaL*At11L + epsdiss1*PDdissipationNth1gt11
+ epsdiss2*PDdissipationNth2gt11 + epsdiss3*PDdissipationNth3gt11 +
@@ -672,77 +670,72 @@ static void ML_BSSN_RHS1_Body(cGH const * restrict const cctkGH, int const dir,
CCTK_REAL dotXt1 =
0.333333333333333333333333333333*(7*(gtu12*PDstandardNth12beta1 +
- gtu13*PDstandardNth13beta1) + 6*gtu23*PDstandardNth23beta1 -
+ gtu13*PDstandardNth13beta1) + gtu11*(4*PDstandardNth11beta1 +
+ PDstandardNth12beta2 + PDstandardNth13beta3) +
+ gtu12*(PDstandardNth22beta2 + PDstandardNth23beta3) +
+ gtu13*(PDstandardNth23beta2 + PDstandardNth33beta3) -
6*(Atu11*PDstandardNth1alpha + Atu12*PDstandardNth2alpha +
- Atu13*PDstandardNth3alpha) + gtu11*(4*PDstandardNth11beta1 +
- PDstandardNth12beta2 + PDstandardNth13beta3 -
- 150.7964473723100754462068823974161384415*alphaL*S1) +
- gtu12*(PDstandardNth22beta2 + PDstandardNth23beta3 -
- 150.7964473723100754462068823974161384415*alphaL*S2) +
- gtu13*(PDstandardNth23beta2 + PDstandardNth33beta3 -
- 150.7964473723100754462068823974161384415*alphaL*S3) +
- (-PDstandardNth1beta1 + 2*PDstandardNth3beta3)*Xtn1 +
- 2*(alphaL*(18*(Atu11*cdphi1 + Atu12*cdphi2 + Atu13*cdphi3) +
- 6*(Atu12*Gt112 + Atu13*Gt113 + Atu23*Gt123) + 3*(Atu11*Gt111 +
- Atu22*Gt122 + Atu33*Gt133) - 2*(gtu11*PDstandardNth1trK +
- gtu12*PDstandardNth2trK + gtu13*PDstandardNth3trK)) +
- PDstandardNth2beta2*Xtn1) - 3*(PDstandardNth2beta1*Xtn2 +
- PDstandardNth3beta1*Xtn3) + 3*(epsdiss1*PDdissipationNth1Xt1 +
- epsdiss2*PDdissipationNth2Xt1 + epsdiss3*PDdissipationNth3Xt1 +
- gtu22*PDstandardNth22beta1 + gtu33*PDstandardNth33beta1 +
- beta1L*PDupwindNthAnti1Xt1 + beta2L*PDupwindNthAnti2Xt1 +
- beta3L*PDupwindNthAnti3Xt1 + PDupwindNthSymm1Xt1*Abs(beta1L) +
- PDupwindNthSymm2Xt1*Abs(beta2L) + PDupwindNthSymm3Xt1*Abs(beta3L)));
+ Atu13*PDstandardNth3alpha) + 6*(gtu23*PDstandardNth23beta1 +
+ alphaL*(6*(Atu11*cdphi1 + Atu12*cdphi2 + Atu13*cdphi3) + Atu11*Gt111 +
+ Atu22*Gt122 + 2*(Atu12*Gt112 + Atu13*Gt113 + Atu23*Gt123) + Atu33*Gt133
+ - 0.666666666666666666666666666667*(gtu11*PDstandardNth1trK +
+ gtu12*PDstandardNth2trK + gtu13*PDstandardNth3trK))) -
+ 150.7964473723100754462068823974161384415*alphaL*(gtu11*S1 + gtu12*S2 +
+ gtu13*S3) + (-3*PDstandardNth1beta1 + 2*(PDstandardNth1beta1 +
+ PDstandardNth2beta2 + PDstandardNth3beta3))*Xtn1 -
+ 3*(PDstandardNth2beta1*Xtn2 + PDstandardNth3beta1*Xtn3) +
+ 3*(epsdiss1*PDdissipationNth1Xt1 + epsdiss2*PDdissipationNth2Xt1 +
+ epsdiss3*PDdissipationNth3Xt1 + gtu22*PDstandardNth22beta1 +
+ gtu33*PDstandardNth33beta1 + beta1L*PDupwindNthAnti1Xt1 +
+ beta2L*PDupwindNthAnti2Xt1 + beta3L*PDupwindNthAnti3Xt1 +
+ PDupwindNthSymm1Xt1*Abs(beta1L) + PDupwindNthSymm2Xt1*Abs(beta2L) +
+ PDupwindNthSymm3Xt1*Abs(beta3L)));
CCTK_REAL dotXt2 =
- 0.333333333333333333333333333333*(6*gtu13*PDstandardNth13beta2 -
- 6*(Atu12*PDstandardNth1alpha + Atu22*PDstandardNth2alpha +
- Atu23*PDstandardNth3alpha) + gtu12*(PDstandardNth11beta1 +
- 7*PDstandardNth12beta2 + PDstandardNth13beta3 -
- 150.7964473723100754462068823974161384415*alphaL*S1) +
+ 0.333333333333333333333333333333*(gtu12*(PDstandardNth11beta1 +
+ 7*PDstandardNth12beta2 + PDstandardNth13beta3) +
gtu22*(PDstandardNth12beta1 + 4*PDstandardNth22beta2 +
- PDstandardNth23beta3 -
- 150.7964473723100754462068823974161384415*alphaL*S2) +
- gtu23*(PDstandardNth13beta1 + 7*PDstandardNth23beta2 +
- PDstandardNth33beta3 -
- 150.7964473723100754462068823974161384415*alphaL*S3) +
- (-PDstandardNth2beta2 + 2*PDstandardNth3beta3)*Xtn2 +
- 2*(alphaL*(18*(Atu12*cdphi1 + Atu22*cdphi2 + Atu23*cdphi3) +
- 6*(Atu12*Gt212 + Atu13*Gt213 + Atu23*Gt223) + 3*(Atu11*Gt211 +
- Atu22*Gt222 + Atu33*Gt233) - 2*(gtu12*PDstandardNth1trK +
- gtu22*PDstandardNth2trK + gtu23*PDstandardNth3trK)) +
- PDstandardNth1beta1*Xtn2) - 3*(PDstandardNth1beta2*Xtn1 +
- PDstandardNth3beta2*Xtn3) + 3*(epsdiss1*PDdissipationNth1Xt2 +
- epsdiss2*PDdissipationNth2Xt2 + epsdiss3*PDdissipationNth3Xt2 +
- gtu11*PDstandardNth11beta2 + gtu33*PDstandardNth33beta2 +
- beta1L*PDupwindNthAnti1Xt2 + beta2L*PDupwindNthAnti2Xt2 +
- beta3L*PDupwindNthAnti3Xt2 + PDupwindNthSymm1Xt2*Abs(beta1L) +
- PDupwindNthSymm2Xt2*Abs(beta2L) + PDupwindNthSymm3Xt2*Abs(beta3L)));
+ PDstandardNth23beta3) + gtu23*(PDstandardNth13beta1 +
+ 7*PDstandardNth23beta2 + PDstandardNth33beta3) -
+ 6*(Atu12*PDstandardNth1alpha + Atu22*PDstandardNth2alpha +
+ Atu23*PDstandardNth3alpha) + 6*(gtu13*PDstandardNth13beta2 +
+ alphaL*(6*(Atu12*cdphi1 + Atu22*cdphi2 + Atu23*cdphi3) + Atu11*Gt211 +
+ Atu22*Gt222 + 2*(Atu12*Gt212 + Atu13*Gt213 + Atu23*Gt223) + Atu33*Gt233
+ - 0.666666666666666666666666666667*(gtu12*PDstandardNth1trK +
+ gtu22*PDstandardNth2trK + gtu23*PDstandardNth3trK))) -
+ 150.7964473723100754462068823974161384415*alphaL*(gtu12*S1 + gtu22*S2 +
+ gtu23*S3) + 2*(PDstandardNth1beta1 + PDstandardNth2beta2 +
+ PDstandardNth3beta3)*Xtn2 - 3*(PDstandardNth1beta2*Xtn1 +
+ PDstandardNth2beta2*Xtn2 + PDstandardNth3beta2*Xtn3) +
+ 3*(epsdiss1*PDdissipationNth1Xt2 + epsdiss2*PDdissipationNth2Xt2 +
+ epsdiss3*PDdissipationNth3Xt2 + gtu11*PDstandardNth11beta2 +
+ gtu33*PDstandardNth33beta2 + beta1L*PDupwindNthAnti1Xt2 +
+ beta2L*PDupwindNthAnti2Xt2 + beta3L*PDupwindNthAnti3Xt2 +
+ PDupwindNthSymm1Xt2*Abs(beta1L) + PDupwindNthSymm2Xt2*Abs(beta2L) +
+ PDupwindNthSymm3Xt2*Abs(beta3L)));
CCTK_REAL dotXt3 =
- 0.333333333333333333333333333333*(6*gtu12*PDstandardNth12beta3 -
- 6*(Atu13*PDstandardNth1alpha + Atu23*PDstandardNth2alpha +
- Atu33*PDstandardNth3alpha) + gtu13*(PDstandardNth11beta1 +
- PDstandardNth12beta2 + 7*PDstandardNth13beta3 -
- 150.7964473723100754462068823974161384415*alphaL*S1) +
+ 0.333333333333333333333333333333*(gtu13*(PDstandardNth11beta1 +
+ PDstandardNth12beta2 + 7*PDstandardNth13beta3) +
gtu23*(PDstandardNth12beta1 + PDstandardNth22beta2 +
- 7*PDstandardNth23beta3 -
- 150.7964473723100754462068823974161384415*alphaL*S2) +
- gtu33*(PDstandardNth13beta1 + PDstandardNth23beta2 +
- 4*PDstandardNth33beta3 -
- 150.7964473723100754462068823974161384415*alphaL*S3) -
- 3*(PDstandardNth1beta3*Xtn1 + PDstandardNth2beta3*Xtn2) +
- (2*PDstandardNth2beta2 - PDstandardNth3beta3)*Xtn3 +
- 2*(alphaL*(18*(Atu13*cdphi1 + Atu23*cdphi2 + Atu33*cdphi3) +
- 6*(Atu12*Gt312 + Atu13*Gt313 + Atu23*Gt323) + 3*(Atu11*Gt311 +
- Atu22*Gt322 + Atu33*Gt333) - 2*(gtu13*PDstandardNth1trK +
- gtu23*PDstandardNth2trK + gtu33*PDstandardNth3trK)) +
- PDstandardNth1beta1*Xtn3) + 3*(epsdiss1*PDdissipationNth1Xt3 +
- epsdiss2*PDdissipationNth2Xt3 + epsdiss3*PDdissipationNth3Xt3 +
- gtu11*PDstandardNth11beta3 + gtu22*PDstandardNth22beta3 +
- beta1L*PDupwindNthAnti1Xt3 + beta2L*PDupwindNthAnti2Xt3 +
- beta3L*PDupwindNthAnti3Xt3 + PDupwindNthSymm1Xt3*Abs(beta1L) +
- PDupwindNthSymm2Xt3*Abs(beta2L) + PDupwindNthSymm3Xt3*Abs(beta3L)));
+ 7*PDstandardNth23beta3) + gtu33*(PDstandardNth13beta1 +
+ PDstandardNth23beta2 + 4*PDstandardNth33beta3) -
+ 6*(Atu13*PDstandardNth1alpha + Atu23*PDstandardNth2alpha +
+ Atu33*PDstandardNth3alpha) + 6*(gtu12*PDstandardNth12beta3 +
+ alphaL*(6*(Atu13*cdphi1 + Atu23*cdphi2 + Atu33*cdphi3) + Atu11*Gt311 +
+ Atu22*Gt322 + 2*(Atu12*Gt312 + Atu13*Gt313 + Atu23*Gt323) + Atu33*Gt333
+ - 0.666666666666666666666666666667*(gtu13*PDstandardNth1trK +
+ gtu23*PDstandardNth2trK + gtu33*PDstandardNth3trK))) -
+ 150.7964473723100754462068823974161384415*alphaL*(gtu13*S1 + gtu23*S2 +
+ gtu33*S3) + 2*(PDstandardNth1beta1 + PDstandardNth2beta2 +
+ PDstandardNth3beta3)*Xtn3 - 3*(PDstandardNth1beta3*Xtn1 +
+ PDstandardNth2beta3*Xtn2 + PDstandardNth3beta3*Xtn3) +
+ 3*(epsdiss1*PDdissipationNth1Xt3 + epsdiss2*PDdissipationNth2Xt3 +
+ epsdiss3*PDdissipationNth3Xt3 + gtu11*PDstandardNth11beta3 +
+ gtu22*PDstandardNth22beta3 + beta1L*PDupwindNthAnti1Xt3 +
+ beta2L*PDupwindNthAnti2Xt3 + beta3L*PDupwindNthAnti3Xt3 +
+ PDupwindNthSymm1Xt3*Abs(beta1L) + PDupwindNthSymm2Xt3*Abs(beta2L) +
+ PDupwindNthSymm3Xt3*Abs(beta3L)));
CCTK_REAL Xt1rhsL = dotXt1;
generated by cgit v1.2.3 (git 2.39.1) at 2024-06-01 04:31:21 +0000