Regenerate code

2 files changed, 10 insertions, 10 deletions

diff --git a/ML_CCZ4/param.ccl b/ML_CCZ4/param.ccl
index f3cd91b..2e148b9 100644
--- a/ML_CCZ4/param.ccl
+++ b/ML_CCZ4/param.ccl
@@ -54,6 +54,12 @@ CCTK_INT verbose "verbose" STEERABLE=ALWAYS
 } 0
 
 restricted:
+CCTK_REAL GammaShift "Covariant shift term in Gamma"
+{
+  "*:*" :: ""
+} 0.5
+
+restricted:
 CCTK_REAL dampk1 "CCZ4 damping term 1 for Theta and Z"
 {
   "*:*" :: ""
@@ -138,12 +144,6 @@ CCTK_REAL EpsDiss "Dissipation strength"
 } 0
 
 restricted:
-CCTK_INT GammaShift "Covariant shift term in Gamma"
-{
-  *:* :: ""
-} 0
-
-restricted:
 CCTK_INT harmonicN "d/dt alpha = - f alpha^n K  (harmonic=2, 1+log=1)"
 {
   *:* :: ""
diff --git a/ML_CCZ4/src/ML_CCZ4_RHS1.cc b/ML_CCZ4/src/ML_CCZ4_RHS1.cc
index 7a5d675..ff38e40 100644
--- a/ML_CCZ4/src/ML_CCZ4_RHS1.cc
+++ b/ML_CCZ4/src/ML_CCZ4_RHS1.cc
@@ -2216,13 +2216,13 @@ static void ML_CCZ4_RHS1_Body(cGH const * restrict const cctkGH, int const dir,
       kmul(ToReal(0.666666666666666666666666666667),knmsub(gt33L,kadd(JacPDstandardNth1beta1,kmadd(JacPDstandardNth3beta3,ToReal(-2),JacPDstandardNth2beta2)),kmadd(alphaL,kmadd(At33L,kmadd(gt33L,gtu33,ToReal(-3)),kmul(gt33L,kmadd(At11L,gtu11,kmadd(At22L,gtu22,kmul(kmadd(At12L,gtu12,kmadd(At13L,gtu13,kmul(At23L,gtu23))),ToReal(2)))))),kmul(kmadd(gt13L,JacPDstandardNth3beta1,kmul(gt23L,JacPDstandardNth3beta2)),ToReal(3)))));
     
     CCTK_REAL_VEC dotXt1 = 
-      kmul(ToReal(0.333333333333333333333333333333),kmadd(kmadd(Atu11,JacPDstandardNth1alpha,kmadd(Atu12,JacPDstandardNth2alpha,kmul(Atu13,JacPDstandardNth3alpha))),ToReal(-6),kmadd(kmadd(JacPDstandardNth2beta1,Xtn2,kmul(JacPDstandardNth3beta1,Xtn3)),ToReal(-3),kmadd(Xtn1,kmsub(kadd(JacPDstandardNth2beta2,JacPDstandardNth3beta3),ToReal(2),JacPDstandardNth1beta1),kmadd(kmadd(gtu12,JacPDstandardNth12beta1,kmadd(gtu13,JacPDstandardNth13beta1,kmadd(gtu22,JacPDstandardNth22beta1,kmadd(gtu23,kadd(JacPDstandardNth23beta1,JacPDstandardNth32beta1),kmadd(gtu33,JacPDstandardNth33beta1,IfThen(GammaShift,kmul(e4phi,kmul(ToReal(-0.666666666666666666666666666667),kmadd(Z1,kmadd(kadd(JacPDstandardNth2beta2,JacPDstandardNth3beta3),ToReal(-2),JacPDstandardNth1beta1),kmul(kmadd(JacPDstandardNth2beta1,Z2,kmul(JacPDstandardNth3beta1,Z3)),ToReal(3))))),ToReal(0))))))),ToReal(3),kmadd(gtu11,kadd(JacPDstandardNth12beta2,kadd(JacPDstandardNth13beta3,kmadd(ThetaL,kmul(JacPDstandardNth1alpha,ToReal(-6)),kmul(JacPDstandardNth11beta1,ToReal(4))))),kmadd(gtu12,kadd(JacPDstandardNth22beta2,kadd(JacPDstandardNth23beta3,kmadd(ThetaL,kmul(JacPDstandardNth2alpha,ToReal(-6)),kmul(JacPDstandardNth21beta1,ToReal(4))))),kmadd(gtu13,kadd(JacPDstandardNth32beta2,kadd(JacPDstandardNth33beta3,kmadd(ThetaL,kmul(JacPDstandardNth3alpha,ToReal(-6)),kmul(JacPDstandardNth31beta1,ToReal(4))))),kmul(alphaL,kmul(ToReal(2),kmadd(kmadd(gtu11,JacPDstandardNth1trK,kmadd(gtu12,JacPDstandardNth2trK,kmadd(gtu13,JacPDstandardNth3trK,kmul(trKL,kmul(e4phi,Z1))))),ToReal(-2),kmadd(kmadd(Atu23,Gt123,kmul(Atu12,kmadd(cdphi2,ToReal(3),Gt112))),ToReal(6),kmadd(ToReal(3),kmadd(Atu22,Gt122,kmadd(Atu33,Gt133,kmadd(gtu11,JacPDstandardNth1Theta,kmadd(gtu12,JacPDstandardNth2Theta,kmadd(gtu13,JacPDstandardNth3Theta,kmul(Atu11,kmadd(cdphi1,ToReal(6),Gt111))))))),kmadd(Atu13,kmadd(Gt113,ToReal(6),kmul(cdphi3,ToReal(18))),kmadd(e4phi,kmul(Z1,kmul(ToReal(-3),ToReal(dampk1))),kmul(kmadd(gtu11,S1,kmadd(gtu12,S2,kmul(gtu13,S3))),kmul(ToReal(-24),ToReal(Pi))))))))))))))))));
+      kmul(ToReal(0.333333333333333333333333333333),kmadd(kmadd(JacPDstandardNth2beta1,Xtn2,kmul(JacPDstandardNth3beta1,Xtn3)),ToReal(-3),kmadd(Xtn1,kmsub(kadd(JacPDstandardNth2beta2,JacPDstandardNth3beta3),ToReal(2),JacPDstandardNth1beta1),kmadd(kmadd(gtu12,JacPDstandardNth12beta1,kmadd(gtu13,JacPDstandardNth13beta1,kmadd(gtu22,JacPDstandardNth22beta1,kmadd(gtu23,kadd(JacPDstandardNth23beta1,JacPDstandardNth32beta1),kmul(gtu33,JacPDstandardNth33beta1))))),ToReal(3),kmadd(gtu11,kadd(JacPDstandardNth12beta2,kadd(JacPDstandardNth13beta3,kmadd(ThetaL,kmul(JacPDstandardNth1alpha,ToReal(-6)),kmul(JacPDstandardNth11beta1,ToReal(4))))),kmadd(gtu12,kadd(JacPDstandardNth22beta2,kadd(JacPDstandardNth23beta3,kmadd(ThetaL,kmul(JacPDstandardNth2alpha,ToReal(-6)),kmul(JacPDstandardNth21beta1,ToReal(4))))),kmadd(gtu13,kadd(JacPDstandardNth32beta2,kadd(JacPDstandardNth33beta3,kmadd(ThetaL,kmul(JacPDstandardNth3alpha,ToReal(-6)),kmul(JacPDstandardNth31beta1,ToReal(4))))),kmadd(e4phi,kmul(kmadd(JacPDstandardNth3beta1,kmul(Z3,ToReal(-6)),kmul(Z1,kmadd(JacPDstandardNth1beta1,ToReal(-2),kmul(kadd(JacPDstandardNth2beta2,JacPDstandardNth3beta3),ToReal(4))))),ToReal(GammaShift)),kmadd(ToReal(-6),kmadd(Atu11,JacPDstandardNth1alpha,kmadd(Atu12,JacPDstandardNth2alpha,kmadd(Atu13,JacPDstandardNth3alpha,kmul(e4phi,kmul(JacPDstandardNth2beta1,kmul(Z2,ToReal(GammaShift))))))),kmul(alphaL,kmul(ToReal(2),kmadd(kmadd(gtu11,JacPDstandardNth1trK,kmadd(gtu12,JacPDstandardNth2trK,kmadd(gtu13,JacPDstandardNth3trK,kmul(trKL,kmul(e4phi,Z1))))),ToReal(-2),kmadd(kmadd(Atu23,Gt123,kmul(Atu12,kmadd(cdphi2,ToReal(3),Gt112))),ToReal(6),kmadd(ToReal(3),kmadd(Atu22,Gt122,kmadd(Atu33,Gt133,kmadd(gtu11,JacPDstandardNth1Theta,kmadd(gtu12,JacPDstandardNth2Theta,kmadd(gtu13,JacPDstandardNth3Theta,kmul(Atu11,kmadd(cdphi1,ToReal(6),Gt111))))))),kmadd(Atu13,kmadd(Gt113,ToReal(6),kmul(cdphi3,ToReal(18))),kmadd(e4phi,kmul(Z1,kmul(ToReal(-3),ToReal(dampk1))),kmul(kmadd(gtu11,S1,kmadd(gtu12,S2,kmul(gtu13,S3))),kmul(ToReal(-24),ToReal(Pi)))))))))))))))))));
     
     CCTK_REAL_VEC dotXt2 = 
-      kmul(ToReal(0.333333333333333333333333333333),kmadd(kmadd(Atu12,JacPDstandardNth1alpha,kmadd(Atu22,JacPDstandardNth2alpha,kmul(Atu23,JacPDstandardNth3alpha))),ToReal(-6),kmadd(kmadd(JacPDstandardNth1beta2,Xtn1,kmul(JacPDstandardNth3beta2,Xtn3)),ToReal(-3),kmadd(Xtn2,kmsub(kadd(JacPDstandardNth1beta1,JacPDstandardNth3beta3),ToReal(2),JacPDstandardNth2beta2),kmadd(kmadd(gtu11,JacPDstandardNth11beta2,kmadd(gtu12,JacPDstandardNth21beta2,kmadd(gtu23,JacPDstandardNth23beta2,kmadd(gtu13,kadd(JacPDstandardNth13beta2,JacPDstandardNth31beta2),kmadd(gtu33,JacPDstandardNth33beta2,IfThen(GammaShift,kmul(e4phi,kmul(ToReal(-0.666666666666666666666666666667),kmadd(Z2,kmadd(kadd(JacPDstandardNth1beta1,JacPDstandardNth3beta3),ToReal(-2),JacPDstandardNth2beta2),kmul(kmadd(JacPDstandardNth1beta2,Z1,kmul(JacPDstandardNth3beta2,Z3)),ToReal(3))))),ToReal(0))))))),ToReal(3),kmadd(gtu12,kadd(JacPDstandardNth11beta1,kadd(JacPDstandardNth13beta3,kmadd(ThetaL,kmul(JacPDstandardNth1alpha,ToReal(-6)),kmul(JacPDstandardNth12beta2,ToReal(4))))),kmadd(gtu22,kadd(JacPDstandardNth21beta1,kadd(JacPDstandardNth23beta3,kmadd(ThetaL,kmul(JacPDstandardNth2alpha,ToReal(-6)),kmul(JacPDstandardNth22beta2,ToReal(4))))),kmadd(gtu23,kadd(JacPDstandardNth31beta1,kadd(JacPDstandardNth33beta3,kmadd(ThetaL,kmul(JacPDstandardNth3alpha,ToReal(-6)),kmul(JacPDstandardNth32beta2,ToReal(4))))),kmul(alphaL,kmul(ToReal(2),kmadd(kmadd(gtu12,JacPDstandardNth1trK,kmadd(gtu22,JacPDstandardNth2trK,kmadd(gtu23,JacPDstandardNth3trK,kmul(trKL,kmul(e4phi,Z2))))),ToReal(-2),kmadd(kmadd(Atu13,Gt213,kmul(Atu12,kmadd(cdphi1,ToReal(3),Gt212))),ToReal(6),kmadd(ToReal(3),kmadd(Atu11,Gt211,kmadd(Atu33,Gt233,kmadd(gtu12,JacPDstandardNth1Theta,kmadd(gtu22,JacPDstandardNth2Theta,kmadd(gtu23,JacPDstandardNth3Theta,kmul(Atu22,kmadd(cdphi2,ToReal(6),Gt222))))))),kmadd(Atu23,kmadd(Gt223,ToReal(6),kmul(cdphi3,ToReal(18))),kmadd(e4phi,kmul(Z2,kmul(ToReal(-3),ToReal(dampk1))),kmul(kmadd(gtu12,S1,kmadd(gtu22,S2,kmul(gtu23,S3))),kmul(ToReal(-24),ToReal(Pi))))))))))))))))));
+      kmul(ToReal(0.333333333333333333333333333333),kmadd(kmadd(JacPDstandardNth1beta2,Xtn1,kmul(JacPDstandardNth3beta2,Xtn3)),ToReal(-3),kmadd(Xtn2,kmsub(kadd(JacPDstandardNth1beta1,JacPDstandardNth3beta3),ToReal(2),JacPDstandardNth2beta2),kmadd(kmadd(gtu11,JacPDstandardNth11beta2,kmadd(gtu23,JacPDstandardNth23beta2,kmadd(gtu13,kadd(JacPDstandardNth13beta2,JacPDstandardNth31beta2),kmul(gtu33,JacPDstandardNth33beta2)))),ToReal(3),kmadd(gtu12,kadd(JacPDstandardNth11beta1,kadd(JacPDstandardNth13beta3,kmadd(ThetaL,kmul(JacPDstandardNth1alpha,ToReal(-6)),kmadd(JacPDstandardNth21beta2,ToReal(3),kmul(JacPDstandardNth12beta2,ToReal(4)))))),kmadd(gtu22,kadd(JacPDstandardNth21beta1,kadd(JacPDstandardNth23beta3,kmadd(ThetaL,kmul(JacPDstandardNth2alpha,ToReal(-6)),kmul(JacPDstandardNth22beta2,ToReal(4))))),kmadd(gtu23,kadd(JacPDstandardNth31beta1,kadd(JacPDstandardNth33beta3,kmadd(ThetaL,kmul(JacPDstandardNth3alpha,ToReal(-6)),kmul(JacPDstandardNth32beta2,ToReal(4))))),kmadd(e4phi,kmul(kmadd(JacPDstandardNth3beta2,kmul(Z3,ToReal(-6)),kmul(Z2,kmadd(JacPDstandardNth2beta2,ToReal(-2),kmul(kadd(JacPDstandardNth1beta1,JacPDstandardNth3beta3),ToReal(4))))),ToReal(GammaShift)),kmadd(ToReal(-6),kmadd(Atu12,JacPDstandardNth1alpha,kmadd(Atu22,JacPDstandardNth2alpha,kmadd(Atu23,JacPDstandardNth3alpha,kmul(e4phi,kmul(JacPDstandardNth1beta2,kmul(Z1,ToReal(GammaShift))))))),kmul(alphaL,kmul(ToReal(2),kmadd(kmadd(gtu12,JacPDstandardNth1trK,kmadd(gtu22,JacPDstandardNth2trK,kmadd(gtu23,JacPDstandardNth3trK,kmul(trKL,kmul(e4phi,Z2))))),ToReal(-2),kmadd(kmadd(Atu13,Gt213,kmul(Atu12,kmadd(cdphi1,ToReal(3),Gt212))),ToReal(6),kmadd(ToReal(3),kmadd(Atu11,Gt211,kmadd(Atu33,Gt233,kmadd(gtu12,JacPDstandardNth1Theta,kmadd(gtu22,JacPDstandardNth2Theta,kmadd(gtu23,JacPDstandardNth3Theta,kmul(Atu22,kmadd(cdphi2,ToReal(6),Gt222))))))),kmadd(Atu23,kmadd(Gt223,ToReal(6),kmul(cdphi3,ToReal(18))),kmadd(e4phi,kmul(Z2,kmul(ToReal(-3),ToReal(dampk1))),kmul(kmadd(gtu12,S1,kmadd(gtu22,S2,kmul(gtu23,S3))),kmul(ToReal(-24),ToReal(Pi)))))))))))))))))));
     
     CCTK_REAL_VEC dotXt3 = 
-      kmul(ToReal(0.333333333333333333333333333333),kmadd(kmadd(Atu13,JacPDstandardNth1alpha,kmadd(Atu23,JacPDstandardNth2alpha,kmul(Atu33,JacPDstandardNth3alpha))),ToReal(-6),kmadd(kmadd(JacPDstandardNth1beta3,Xtn1,kmul(JacPDstandardNth2beta3,Xtn2)),ToReal(-3),knmsub(JacPDstandardNth3beta3,Xtn3,kmadd(kmadd(gtu11,JacPDstandardNth11beta3,kmadd(gtu12,kadd(JacPDstandardNth12beta3,JacPDstandardNth21beta3),kmadd(gtu22,JacPDstandardNth22beta3,kmadd(gtu23,JacPDstandardNth32beta3,IfThen(GammaShift,kmul(e4phi,kmul(ToReal(-0.666666666666666666666666666667),kmadd(Z3,kmadd(kadd(JacPDstandardNth1beta1,JacPDstandardNth2beta2),ToReal(-2),JacPDstandardNth3beta3),kmul(kmadd(JacPDstandardNth1beta3,Z1,kmul(JacPDstandardNth2beta3,Z2)),ToReal(3))))),ToReal(0)))))),ToReal(3),kmadd(gtu13,kadd(JacPDstandardNth11beta1,kadd(JacPDstandardNth12beta2,kmadd(ThetaL,kmul(JacPDstandardNth1alpha,ToReal(-6)),kmadd(JacPDstandardNth31beta3,ToReal(3),kmul(JacPDstandardNth13beta3,ToReal(4)))))),kmadd(gtu23,kadd(JacPDstandardNth21beta1,kadd(JacPDstandardNth22beta2,kmadd(ThetaL,kmul(JacPDstandardNth2alpha,ToReal(-6)),kmul(JacPDstandardNth23beta3,ToReal(4))))),kmadd(gtu33,kadd(JacPDstandardNth31beta1,kadd(JacPDstandardNth32beta2,kmadd(ThetaL,kmul(JacPDstandardNth3alpha,ToReal(-6)),kmul(JacPDstandardNth33beta3,ToReal(4))))),kmul(ToReal(2),kmadd(kadd(JacPDstandardNth1beta1,JacPDstandardNth2beta2),Xtn3,kmul(alphaL,kmadd(kmadd(gtu13,JacPDstandardNth1trK,kmadd(gtu23,JacPDstandardNth2trK,kmadd(gtu33,JacPDstandardNth3trK,kmul(trKL,kmul(e4phi,Z3))))),ToReal(-2),kmadd(kmadd(Atu11,Gt311,kmadd(Atu22,Gt322,kmadd(gtu13,JacPDstandardNth1Theta,kmadd(gtu23,JacPDstandardNth2Theta,kmul(gtu33,JacPDstandardNth3Theta))))),ToReal(3),kmadd(kmadd(Atu12,Gt312,kmadd(Atu13,kmadd(cdphi1,ToReal(3),Gt313),kmul(Atu23,kmadd(cdphi2,ToReal(3),Gt323)))),ToReal(6),kmadd(Atu33,kmadd(Gt333,ToReal(3),kmul(cdphi3,ToReal(18))),kmadd(e4phi,kmul(Z3,kmul(ToReal(-3),ToReal(dampk1))),kmul(kmadd(gtu13,S1,kmadd(gtu23,S2,kmul(gtu33,S3))),kmul(ToReal(-24),ToReal(Pi)))))))))))))))))));
+      kmul(ToReal(0.333333333333333333333333333333),kmadd(kmadd(JacPDstandardNth1beta3,Xtn1,kmul(JacPDstandardNth2beta3,Xtn2)),ToReal(-3),kmadd(kmadd(gtu11,JacPDstandardNth11beta3,kmadd(gtu12,kadd(JacPDstandardNth12beta3,JacPDstandardNth21beta3),kmadd(gtu22,JacPDstandardNth22beta3,kmul(gtu23,JacPDstandardNth32beta3)))),ToReal(3),kmadd(gtu13,kadd(JacPDstandardNth11beta1,kadd(JacPDstandardNth12beta2,kmadd(ThetaL,kmul(JacPDstandardNth1alpha,ToReal(-6)),kmadd(JacPDstandardNth31beta3,ToReal(3),kmul(JacPDstandardNth13beta3,ToReal(4)))))),kmadd(gtu23,kadd(JacPDstandardNth21beta1,kadd(JacPDstandardNth22beta2,kmadd(ThetaL,kmul(JacPDstandardNth2alpha,ToReal(-6)),kmul(JacPDstandardNth23beta3,ToReal(4))))),kmadd(gtu33,kadd(JacPDstandardNth31beta1,kadd(JacPDstandardNth32beta2,kmadd(ThetaL,kmul(JacPDstandardNth3alpha,ToReal(-6)),kmul(JacPDstandardNth33beta3,ToReal(4))))),kmadd(e4phi,kmul(kmadd(JacPDstandardNth2beta3,kmul(Z2,ToReal(-6)),kmul(kadd(JacPDstandardNth1beta1,JacPDstandardNth2beta2),kmul(Z3,ToReal(4)))),ToReal(GammaShift)),kmadd(ToReal(-6),kmadd(Atu13,JacPDstandardNth1alpha,kmadd(Atu23,JacPDstandardNth2alpha,kmadd(Atu33,JacPDstandardNth3alpha,kmul(e4phi,kmul(JacPDstandardNth1beta3,kmul(Z1,ToReal(GammaShift))))))),kmadd(JacPDstandardNth3beta3,kmsub(e4phi,kmul(Z3,kmul(ToReal(-2),ToReal(GammaShift))),Xtn3),kmul(ToReal(2),kmadd(kadd(JacPDstandardNth1beta1,JacPDstandardNth2beta2),Xtn3,kmul(alphaL,kmadd(kmadd(gtu13,JacPDstandardNth1trK,kmadd(gtu23,JacPDstandardNth2trK,kmadd(gtu33,JacPDstandardNth3trK,kmul(trKL,kmul(e4phi,Z3))))),ToReal(-2),kmadd(kmadd(Atu11,Gt311,kmadd(Atu22,Gt322,kmadd(gtu13,JacPDstandardNth1Theta,kmadd(gtu23,JacPDstandardNth2Theta,kmul(gtu33,JacPDstandardNth3Theta))))),ToReal(3),kmadd(kmadd(Atu12,Gt312,kmadd(Atu13,kmadd(cdphi1,ToReal(3),Gt313),kmul(Atu23,kmadd(cdphi2,ToReal(3),Gt323)))),ToReal(6),kmadd(Atu33,kmadd(Gt333,ToReal(3),kmul(cdphi3,ToReal(18))),kmadd(e4phi,kmul(Z3,kmul(ToReal(-3),ToReal(dampk1))),kmul(kmadd(gtu13,S1,kmadd(gtu23,S2,kmul(gtu33,S3))),kmul(ToReal(-24),ToReal(Pi))))))))))))))))))));
     
     CCTK_REAL_VEC Xt1rhsL = dotXt1;
     
@@ -2231,7 +2231,7 @@ static void ML_CCZ4_RHS1_Body(cGH const * restrict const cctkGH, int const dir,
     CCTK_REAL_VEC Xt3rhsL = dotXt3;
     
     CCTK_REAL_VEC dotTheta = 
-      knmsub(JacPDstandardNth1alpha,Z1,knmsub(JacPDstandardNth2alpha,Z2,kmsub(alphaL,kmsub(ToReal(0.5),kmadd(gu11,R11,kmadd(gu22,R22,kmadd(gu33,R33,kmadd(ThetaL,kmul(trKL,ToReal(-2)),kmadd(Atm12,kmul(Atm21,ToReal(-2)),kmadd(Atm13,kmul(Atm31,ToReal(-2)),kmadd(Atm23,kmul(Atm32,ToReal(-2)),ksub(ksub(ksub(kmadd(SQR(trKL),ToReal(0.666666666666666666666666666667),kmadd(gu12,kmul(R12,ToReal(2)),kmadd(gu13,kmul(R13,ToReal(2)),kmul(gu23,kmul(R23,ToReal(2)))))),SQR(Atm33)),SQR(Atm22)),SQR(Atm11))))))))),kmul(ThetaL,kmul(ToReal(dampk1),kadd(ToReal(2),ToReal(dampk2))))),kmul(JacPDstandardNth3alpha,Z3))));
+      knmsub(JacPDstandardNth1alpha,Z1,knmsub(JacPDstandardNth2alpha,Z2,kmsub(alphaL,kmadd(ToReal(0.5),kmadd(gu11,R11,kmadd(gu22,R22,kmadd(gu33,R33,kmadd(ThetaL,kmul(trKL,ToReal(-2)),kmadd(Atm12,kmul(Atm21,ToReal(-2)),kmadd(Atm13,kmul(Atm31,ToReal(-2)),kmadd(Atm23,kmul(Atm32,ToReal(-2)),ksub(ksub(ksub(kmadd(SQR(trKL),ToReal(0.666666666666666666666666666667),kmadd(gu12,kmul(R12,ToReal(2)),kmadd(gu13,kmul(R13,ToReal(2)),kmul(gu23,kmul(R23,ToReal(2)))))),SQR(Atm33)),SQR(Atm22)),SQR(Atm11))))))))),kmsub(rho,kmul(ToReal(-8),ToReal(Pi)),kmul(ThetaL,kmul(ToReal(dampk1),kadd(ToReal(2),ToReal(dampk2)))))),kmul(JacPDstandardNth3alpha,Z3))));
     
     CCTK_REAL_VEC ThetarhsL = dotTheta;