diff options
Diffstat (limited to 'ML_BSSN/src')
-rw-r--r-- | ML_BSSN/src/ML_BSSN_RHS1.cc | 96 | ||||
-rw-r--r-- | ML_BSSN/src/ML_BSSN_RHS2.cc | 55 | ||||
-rw-r--r-- | ML_BSSN/src/ML_BSSN_constraints1.cc | 55 |
3 files changed, 110 insertions, 96 deletions
diff --git a/ML_BSSN/src/ML_BSSN_RHS1.cc b/ML_BSSN/src/ML_BSSN_RHS1.cc index d988bba..1c2e051 100644 --- a/ML_BSSN/src/ML_BSSN_RHS1.cc +++ b/ML_BSSN/src/ML_BSSN_RHS1.cc @@ -458,57 +458,65 @@ static void ML_BSSN_RHS1_Body(cGH const * restrict const cctkGH, int const dir, CCTK_REAL dotXt1 = 0.333333333333333333333333333333*(7*(gtu12*PDstandardNth12beta1 + - gtu13*PDstandardNth13beta1) + gtu11*(4*PDstandardNth11beta1 + - PDstandardNth12beta2 + PDstandardNth13beta3) + - gtu12*(PDstandardNth22beta2 + PDstandardNth23beta3) + - 3*(gtu22*PDstandardNth22beta1 + gtu33*PDstandardNth33beta1) + - gtu13*(PDstandardNth23beta2 + PDstandardNth33beta3) - + gtu13*PDstandardNth13beta1) + 6*gtu23*PDstandardNth23beta1 + + 3*(gtu22*PDstandardNth22beta1 + gtu33*PDstandardNth33beta1) - 6*(Atu11*PDstandardNth1alpha + Atu12*PDstandardNth2alpha + - 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)); + 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)); CCTK_REAL dotXt2 = - 0.333333333333333333333333333333*(gtu12*(PDstandardNth11beta1 + - 7*PDstandardNth12beta2 + PDstandardNth13beta3) + - gtu22*(PDstandardNth12beta1 + 4*PDstandardNth22beta2 + - PDstandardNth23beta3) + 3*(gtu11*PDstandardNth11beta2 + - gtu33*PDstandardNth33beta2) + gtu23*(PDstandardNth13beta1 + - 7*PDstandardNth23beta2 + PDstandardNth33beta3) - + 0.333333333333333333333333333333*(6*gtu13*PDstandardNth13beta2 + + 3*(gtu11*PDstandardNth11beta2 + gtu33*PDstandardNth33beta2) - 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)); + Atu23*PDstandardNth3alpha) + gtu12*(PDstandardNth11beta1 + + 7*PDstandardNth12beta2 + PDstandardNth13beta3 - + 150.7964473723100754462068823974161384415*alphaL*S1) + + 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)); CCTK_REAL dotXt3 = - 0.333333333333333333333333333333*(gtu13*(PDstandardNth11beta1 + - PDstandardNth12beta2 + 7*PDstandardNth13beta3) + - 3*(gtu11*PDstandardNth11beta3 + gtu22*PDstandardNth22beta3) + - gtu23*(PDstandardNth12beta1 + PDstandardNth22beta2 + - 7*PDstandardNth23beta3) + gtu33*(PDstandardNth13beta1 + - PDstandardNth23beta2 + 4*PDstandardNth33beta3) - + 0.333333333333333333333333333333*(6*gtu12*PDstandardNth12beta3 + + 3*(gtu11*PDstandardNth11beta3 + gtu22*PDstandardNth22beta3) - 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)); + Atu33*PDstandardNth3alpha) + gtu13*(PDstandardNth11beta1 + + PDstandardNth12beta2 + 7*PDstandardNth13beta3 - + 150.7964473723100754462068823974161384415*alphaL*S1) + + 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)); CCTK_REAL Xt1rhsL = dotXt1; diff --git a/ML_BSSN/src/ML_BSSN_RHS2.cc b/ML_BSSN/src/ML_BSSN_RHS2.cc index a54dc92..e78e16c 100644 --- a/ML_BSSN/src/ML_BSSN_RHS2.cc +++ b/ML_BSSN/src/ML_BSSN_RHS2.cc @@ -407,15 +407,16 @@ static void ML_BSSN_RHS2_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 Rt11 = 3*(Gt111*Gtlu111 + Gt112*Gtlu112 + Gt113*Gtlu113) + - 2*(Gt211*Gtlu121 + Gt212*Gtlu122 + Gt213*Gtlu123 + Gt311*Gtlu131 + - Gt312*Gtlu132 + Gt313*Gtlu133) + Gt211*Gtlu211 + Gt212*Gtlu212 + + CCTK_REAL Rt11 = 0.5*(6*(Gt111*Gtlu111 + Gt112*Gtlu112 + + Gt113*Gtlu113) + 4*(Gt211*Gtlu121 + Gt212*Gtlu122 + Gt213*Gtlu123 + + Gt311*Gtlu131 + Gt312*Gtlu132 + Gt313*Gtlu133) - + gtu11*PDstandardNth11gt11 - 2*gtu12*PDstandardNth12gt11 - + 2*gtu13*PDstandardNth13gt11 + 2*(Gt211*Gtlu211 + Gt212*Gtlu212 + Gt213*Gtlu213 + Gt311*Gtlu311 + Gt312*Gtlu312 + Gt313*Gtlu313 + - gt11L*PDstandardNth1Xt1 + gt12L*PDstandardNth1Xt2 + - gt13L*PDstandardNth1Xt3 + 0.5*(-(gtu11*PDstandardNth11gt11) - - 2*gtu12*PDstandardNth12gt11 - 2*gtu13*PDstandardNth13gt11 - - gtu22*PDstandardNth22gt11 - 2*gtu23*PDstandardNth23gt11 - - gtu33*PDstandardNth33gt11) + Gtl111*Xtn1 + Gtl112*Xtn2 + Gtl113*Xtn3; + gt11L*PDstandardNth1Xt1) + 2*gt12L*PDstandardNth1Xt2 + + 2*gt13L*PDstandardNth1Xt3 - gtu22*PDstandardNth22gt11 - + 2*gtu23*PDstandardNth23gt11 - gtu33*PDstandardNth33gt11 + 2*Gtl111*Xtn1 + + 2*Gtl112*Xtn2 + 2*Gtl113*Xtn3); CCTK_REAL Rt12 = 0.5*(4*(Gt211*Gtlu221 + Gt212*Gtlu222 + Gt213*Gtlu223) + 2*(Gt112*Gtlu111 + Gt122*Gtlu112 + Gt123*Gtlu113 + @@ -447,15 +448,16 @@ static void ML_BSSN_RHS2_Body(cGH const * restrict const cctkGH, int const dir, gt13L*PDstandardNth3Xt3 + Gtl113*Xtn1 + Gtl311*Xtn1 + Gtl123*Xtn2 + Gtl312*Xtn2 + Gtl133*Xtn3 + Gtl313*Xtn3); - CCTK_REAL Rt22 = Gt112*(Gtlu121 + 2*Gtlu211) + Gt122*(Gtlu122 + - 2*Gtlu212) + Gt123*(Gtlu123 + 2*Gtlu213) + 3*(Gt212*Gtlu221 + - Gt222*Gtlu222 + Gt223*Gtlu223) + 2*(Gt312*Gtlu231 + Gt322*Gtlu232 + - Gt323*Gtlu233) + Gt312*Gtlu321 + Gt322*Gtlu322 + Gt323*Gtlu323 + - gt12L*PDstandardNth2Xt1 + gt22L*PDstandardNth2Xt2 + - gt23L*PDstandardNth2Xt3 + 0.5*(-(gtu11*PDstandardNth11gt22) - - 2*gtu12*PDstandardNth12gt22 - 2*gtu13*PDstandardNth13gt22 - - gtu22*PDstandardNth22gt22 - 2*gtu23*PDstandardNth23gt22 - - gtu33*PDstandardNth33gt22) + Gtl212*Xtn1 + Gtl222*Xtn2 + Gtl223*Xtn3; + CCTK_REAL Rt22 = 0.5*(6*(Gt212*Gtlu221 + Gt222*Gtlu222 + + Gt223*Gtlu223) + 4*(Gt112*Gtlu211 + Gt122*Gtlu212 + Gt123*Gtlu213 + + Gt312*Gtlu231 + Gt322*Gtlu232 + Gt323*Gtlu233) - + gtu11*PDstandardNth11gt22 - 2*gtu12*PDstandardNth12gt22 - + 2*gtu13*PDstandardNth13gt22 - gtu22*PDstandardNth22gt22 - + 2*gtu23*PDstandardNth23gt22 + 2*(Gt112*Gtlu121 + Gt122*Gtlu122 + + Gt123*Gtlu123 + Gt312*Gtlu321 + Gt322*Gtlu322 + Gt323*Gtlu323 + + gt12L*PDstandardNth2Xt1) + 2*gt22L*PDstandardNth2Xt2 + + 2*gt23L*PDstandardNth2Xt3 - gtu33*PDstandardNth33gt22 + 2*Gtl212*Xtn1 + + 2*Gtl222*Xtn2 + 2*Gtl223*Xtn3); CCTK_REAL Rt23 = 0.5*(2*(Gt112*Gtlu131 + Gt122*Gtlu132 + Gt123*Gtlu133 + Gt113*Gtlu211 + Gt123*Gtlu212 + Gt133*Gtlu213 + Gt213*Gtlu221 + @@ -472,15 +474,16 @@ static void ML_BSSN_RHS2_Body(cGH const * restrict const cctkGH, int const dir, gt23L*PDstandardNth3Xt3 + Gtl213*Xtn1 + Gtl312*Xtn1 + Gtl223*Xtn2 + Gtl322*Xtn2 + Gtl233*Xtn3 + Gtl323*Xtn3); - CCTK_REAL Rt33 = Gt113*(Gtlu131 + 2*Gtlu311) + Gt123*(Gtlu132 + - 2*Gtlu312) + Gt133*(Gtlu133 + 2*Gtlu313) + Gt213*(Gtlu231 + 2*Gtlu321) - + Gt223*(Gtlu232 + 2*Gtlu322) + Gt233*(Gtlu233 + 2*Gtlu323) + - 3*(Gt313*Gtlu331 + Gt323*Gtlu332 + Gt333*Gtlu333) + - 0.5*(-(gtu11*PDstandardNth11gt33) - 2*gtu12*PDstandardNth12gt33 - - 2*gtu13*PDstandardNth13gt33 - gtu22*PDstandardNth22gt33 - - 2*gtu23*PDstandardNth23gt33 - gtu33*PDstandardNth33gt33) + - gt13L*PDstandardNth3Xt1 + gt23L*PDstandardNth3Xt2 + - gt33L*PDstandardNth3Xt3 + Gtl313*Xtn1 + Gtl323*Xtn2 + Gtl333*Xtn3; + CCTK_REAL Rt33 = 0.5*(4*(Gt113*Gtlu311 + Gt123*Gtlu312 + Gt133*Gtlu313 + + Gt213*Gtlu321 + Gt223*Gtlu322 + Gt233*Gtlu323) + 6*(Gt313*Gtlu331 + + Gt323*Gtlu332 + Gt333*Gtlu333) - gtu11*PDstandardNth11gt33 - + 2*gtu12*PDstandardNth12gt33 - 2*gtu13*PDstandardNth13gt33 - + gtu22*PDstandardNth22gt33 - 2*gtu23*PDstandardNth23gt33 - + gtu33*PDstandardNth33gt33 + 2*(Gt113*Gtlu131 + Gt123*Gtlu132 + + Gt133*Gtlu133 + Gt213*Gtlu231 + Gt223*Gtlu232 + Gt233*Gtlu233 + + gt13L*PDstandardNth3Xt1) + 2*gt23L*PDstandardNth3Xt2 + + 2*gt33L*PDstandardNth3Xt3 + 2*Gtl313*Xtn1 + 2*Gtl323*Xtn2 + + 2*Gtl333*Xtn3); CCTK_REAL fac1 = IfThen(conformalMethod,-0.5*INV(phiL),1); diff --git a/ML_BSSN/src/ML_BSSN_constraints1.cc b/ML_BSSN/src/ML_BSSN_constraints1.cc index 2814517..77156fc 100644 --- a/ML_BSSN/src/ML_BSSN_constraints1.cc +++ b/ML_BSSN/src/ML_BSSN_constraints1.cc @@ -391,15 +391,16 @@ static void ML_BSSN_constraints1_Body(cGH const * restrict const cctkGH, int con CCTK_REAL Xtn3 = Gt311*gtu11 + Gt322*gtu22 + 2*(Gt312*gtu12 + Gt313*gtu13 + Gt323*gtu23) + Gt333*gtu33; - CCTK_REAL Rt11 = 3*(Gt111*Gtlu111 + Gt112*Gtlu112 + Gt113*Gtlu113) + - 2*(Gt211*Gtlu121 + Gt212*Gtlu122 + Gt213*Gtlu123 + Gt311*Gtlu131 + - Gt312*Gtlu132 + Gt313*Gtlu133) + Gt211*Gtlu211 + Gt212*Gtlu212 + + CCTK_REAL Rt11 = 0.5*(6*(Gt111*Gtlu111 + Gt112*Gtlu112 + + Gt113*Gtlu113) + 4*(Gt211*Gtlu121 + Gt212*Gtlu122 + Gt213*Gtlu123 + + Gt311*Gtlu131 + Gt312*Gtlu132 + Gt313*Gtlu133) - + gtu11*PDstandardNth11gt11 - 2*gtu12*PDstandardNth12gt11 - + 2*gtu13*PDstandardNth13gt11 + 2*(Gt211*Gtlu211 + Gt212*Gtlu212 + Gt213*Gtlu213 + Gt311*Gtlu311 + Gt312*Gtlu312 + Gt313*Gtlu313 + - gt11L*PDstandardNth1Xt1 + gt12L*PDstandardNth1Xt2 + - gt13L*PDstandardNth1Xt3 + 0.5*(-(gtu11*PDstandardNth11gt11) - - 2*gtu12*PDstandardNth12gt11 - 2*gtu13*PDstandardNth13gt11 - - gtu22*PDstandardNth22gt11 - 2*gtu23*PDstandardNth23gt11 - - gtu33*PDstandardNth33gt11) + Gtl111*Xtn1 + Gtl112*Xtn2 + Gtl113*Xtn3; + gt11L*PDstandardNth1Xt1) + 2*gt12L*PDstandardNth1Xt2 + + 2*gt13L*PDstandardNth1Xt3 - gtu22*PDstandardNth22gt11 - + 2*gtu23*PDstandardNth23gt11 - gtu33*PDstandardNth33gt11 + 2*Gtl111*Xtn1 + + 2*Gtl112*Xtn2 + 2*Gtl113*Xtn3); CCTK_REAL Rt12 = 0.5*(4*(Gt211*Gtlu221 + Gt212*Gtlu222 + Gt213*Gtlu223) + 2*(Gt112*Gtlu111 + Gt122*Gtlu112 + Gt123*Gtlu113 + @@ -431,15 +432,16 @@ static void ML_BSSN_constraints1_Body(cGH const * restrict const cctkGH, int con gt13L*PDstandardNth3Xt3 + Gtl113*Xtn1 + Gtl311*Xtn1 + Gtl123*Xtn2 + Gtl312*Xtn2 + Gtl133*Xtn3 + Gtl313*Xtn3); - CCTK_REAL Rt22 = Gt112*(Gtlu121 + 2*Gtlu211) + Gt122*(Gtlu122 + - 2*Gtlu212) + Gt123*(Gtlu123 + 2*Gtlu213) + 3*(Gt212*Gtlu221 + - Gt222*Gtlu222 + Gt223*Gtlu223) + 2*(Gt312*Gtlu231 + Gt322*Gtlu232 + - Gt323*Gtlu233) + Gt312*Gtlu321 + Gt322*Gtlu322 + Gt323*Gtlu323 + - gt12L*PDstandardNth2Xt1 + gt22L*PDstandardNth2Xt2 + - gt23L*PDstandardNth2Xt3 + 0.5*(-(gtu11*PDstandardNth11gt22) - - 2*gtu12*PDstandardNth12gt22 - 2*gtu13*PDstandardNth13gt22 - - gtu22*PDstandardNth22gt22 - 2*gtu23*PDstandardNth23gt22 - - gtu33*PDstandardNth33gt22) + Gtl212*Xtn1 + Gtl222*Xtn2 + Gtl223*Xtn3; + CCTK_REAL Rt22 = 0.5*(6*(Gt212*Gtlu221 + Gt222*Gtlu222 + + Gt223*Gtlu223) + 4*(Gt112*Gtlu211 + Gt122*Gtlu212 + Gt123*Gtlu213 + + Gt312*Gtlu231 + Gt322*Gtlu232 + Gt323*Gtlu233) - + gtu11*PDstandardNth11gt22 - 2*gtu12*PDstandardNth12gt22 - + 2*gtu13*PDstandardNth13gt22 - gtu22*PDstandardNth22gt22 - + 2*gtu23*PDstandardNth23gt22 + 2*(Gt112*Gtlu121 + Gt122*Gtlu122 + + Gt123*Gtlu123 + Gt312*Gtlu321 + Gt322*Gtlu322 + Gt323*Gtlu323 + + gt12L*PDstandardNth2Xt1) + 2*gt22L*PDstandardNth2Xt2 + + 2*gt23L*PDstandardNth2Xt3 - gtu33*PDstandardNth33gt22 + 2*Gtl212*Xtn1 + + 2*Gtl222*Xtn2 + 2*Gtl223*Xtn3); CCTK_REAL Rt23 = 0.5*(2*(Gt112*Gtlu131 + Gt122*Gtlu132 + Gt123*Gtlu133 + Gt113*Gtlu211 + Gt123*Gtlu212 + Gt133*Gtlu213 + Gt213*Gtlu221 + @@ -456,15 +458,16 @@ static void ML_BSSN_constraints1_Body(cGH const * restrict const cctkGH, int con gt23L*PDstandardNth3Xt3 + Gtl213*Xtn1 + Gtl312*Xtn1 + Gtl223*Xtn2 + Gtl322*Xtn2 + Gtl233*Xtn3 + Gtl323*Xtn3); - CCTK_REAL Rt33 = Gt113*(Gtlu131 + 2*Gtlu311) + Gt123*(Gtlu132 + - 2*Gtlu312) + Gt133*(Gtlu133 + 2*Gtlu313) + Gt213*(Gtlu231 + 2*Gtlu321) - + Gt223*(Gtlu232 + 2*Gtlu322) + Gt233*(Gtlu233 + 2*Gtlu323) + - 3*(Gt313*Gtlu331 + Gt323*Gtlu332 + Gt333*Gtlu333) + - 0.5*(-(gtu11*PDstandardNth11gt33) - 2*gtu12*PDstandardNth12gt33 - - 2*gtu13*PDstandardNth13gt33 - gtu22*PDstandardNth22gt33 - - 2*gtu23*PDstandardNth23gt33 - gtu33*PDstandardNth33gt33) + - gt13L*PDstandardNth3Xt1 + gt23L*PDstandardNth3Xt2 + - gt33L*PDstandardNth3Xt3 + Gtl313*Xtn1 + Gtl323*Xtn2 + Gtl333*Xtn3; + CCTK_REAL Rt33 = 0.5*(4*(Gt113*Gtlu311 + Gt123*Gtlu312 + Gt133*Gtlu313 + + Gt213*Gtlu321 + Gt223*Gtlu322 + Gt233*Gtlu323) + 6*(Gt313*Gtlu331 + + Gt323*Gtlu332 + Gt333*Gtlu333) - gtu11*PDstandardNth11gt33 - + 2*gtu12*PDstandardNth12gt33 - 2*gtu13*PDstandardNth13gt33 - + gtu22*PDstandardNth22gt33 - 2*gtu23*PDstandardNth23gt33 - + gtu33*PDstandardNth33gt33 + 2*(Gt113*Gtlu131 + Gt123*Gtlu132 + + Gt133*Gtlu133 + Gt213*Gtlu231 + Gt223*Gtlu232 + Gt233*Gtlu233 + + gt13L*PDstandardNth3Xt1) + 2*gt23L*PDstandardNth3Xt2 + + 2*gt33L*PDstandardNth3Xt3 + 2*Gtl313*Xtn1 + 2*Gtl323*Xtn2 + + 2*Gtl333*Xtn3); CCTK_REAL fac1 = IfThen(conformalMethod,-0.5*INV(phiL),1); |