From 8df3aae833eb6bf334b0decb57ebf9e1baf33ddf Mon Sep 17 00:00:00 2001 From: Ian Hinder Date: Tue, 10 Sep 2013 16:03:25 +0200 Subject: Regenerate code (Kranc commit 10daf4a83a783a812c7b5a3aba0df115f0eac98f) --- ML_ADMQuantities/src/ML_ADMQuantities.cc | 363 +++++++++++++++---------------- 1 file changed, 181 insertions(+), 182 deletions(-) (limited to 'ML_ADMQuantities') diff --git a/ML_ADMQuantities/src/ML_ADMQuantities.cc b/ML_ADMQuantities/src/ML_ADMQuantities.cc index bc63940..5b2268a 100644 --- a/ML_ADMQuantities/src/ML_ADMQuantities.cc +++ b/ML_ADMQuantities/src/ML_ADMQuantities.cc @@ -63,13 +63,15 @@ static void ML_ADMQuantities_Body(const cGH* restrict const cctkGH, const int di const CCTK_REAL dyi CCTK_ATTRIBUTE_UNUSED = INV(dy); const CCTK_REAL dzi CCTK_ATTRIBUTE_UNUSED = INV(dz); const CCTK_REAL khalf CCTK_ATTRIBUTE_UNUSED = 0.5; - const CCTK_REAL kthird CCTK_ATTRIBUTE_UNUSED = 1/3.0; - const CCTK_REAL ktwothird CCTK_ATTRIBUTE_UNUSED = 2.0/3.0; - const CCTK_REAL kfourthird CCTK_ATTRIBUTE_UNUSED = 4.0/3.0; - const CCTK_REAL keightthird CCTK_ATTRIBUTE_UNUSED = 8.0/3.0; - const CCTK_REAL hdxi CCTK_ATTRIBUTE_UNUSED = 0.5 * dxi; - const CCTK_REAL hdyi CCTK_ATTRIBUTE_UNUSED = 0.5 * dyi; - const CCTK_REAL hdzi CCTK_ATTRIBUTE_UNUSED = 0.5 * dzi; + const CCTK_REAL kthird CCTK_ATTRIBUTE_UNUSED = + 0.333333333333333333333333333333; + const CCTK_REAL ktwothird CCTK_ATTRIBUTE_UNUSED = + 0.666666666666666666666666666667; + const CCTK_REAL kfourthird CCTK_ATTRIBUTE_UNUSED = + 1.33333333333333333333333333333; + const CCTK_REAL hdxi CCTK_ATTRIBUTE_UNUSED = 0.5*dxi; + const CCTK_REAL hdyi CCTK_ATTRIBUTE_UNUSED = 0.5*dyi; + const CCTK_REAL hdzi CCTK_ATTRIBUTE_UNUSED = 0.5*dzi; /* Initialize predefined quantities */ const CCTK_REAL p1o12dx CCTK_ATTRIBUTE_UNUSED = 0.0833333333333333333333333333333*INV(dx); @@ -440,17 +442,17 @@ static void ML_ADMQuantities_Body(const cGH* restrict const cctkGH, const int di JacPDstandardNth1gt33 = J11L*PDstandardNth1gt33 + J21L*PDstandardNth2gt33 + J31L*PDstandardNth3gt33; - JacPDstandardNth1trK = J11L*PDstandardNth1trK + - J21L*PDstandardNth2trK + J31L*PDstandardNth3trK; + JacPDstandardNth1trK = J11L*PDstandardNth1trK + J21L*PDstandardNth2trK + + J31L*PDstandardNth3trK; - JacPDstandardNth1Xt1 = J11L*PDstandardNth1Xt1 + - J21L*PDstandardNth2Xt1 + J31L*PDstandardNth3Xt1; + JacPDstandardNth1Xt1 = J11L*PDstandardNth1Xt1 + J21L*PDstandardNth2Xt1 + + J31L*PDstandardNth3Xt1; - JacPDstandardNth1Xt2 = J11L*PDstandardNth1Xt2 + - J21L*PDstandardNth2Xt2 + J31L*PDstandardNth3Xt2; + JacPDstandardNth1Xt2 = J11L*PDstandardNth1Xt2 + J21L*PDstandardNth2Xt2 + + J31L*PDstandardNth3Xt2; - JacPDstandardNth1Xt3 = J11L*PDstandardNth1Xt3 + - J21L*PDstandardNth2Xt3 + J31L*PDstandardNth3Xt3; + JacPDstandardNth1Xt3 = J11L*PDstandardNth1Xt3 + J21L*PDstandardNth2Xt3 + + J31L*PDstandardNth3Xt3; JacPDstandardNth2gt11 = J12L*PDstandardNth1gt11 + J22L*PDstandardNth2gt11 + J32L*PDstandardNth3gt11; @@ -470,17 +472,17 @@ static void ML_ADMQuantities_Body(const cGH* restrict const cctkGH, const int di JacPDstandardNth2gt33 = J12L*PDstandardNth1gt33 + J22L*PDstandardNth2gt33 + J32L*PDstandardNth3gt33; - JacPDstandardNth2trK = J12L*PDstandardNth1trK + - J22L*PDstandardNth2trK + J32L*PDstandardNth3trK; + JacPDstandardNth2trK = J12L*PDstandardNth1trK + J22L*PDstandardNth2trK + + J32L*PDstandardNth3trK; - JacPDstandardNth2Xt1 = J12L*PDstandardNth1Xt1 + - J22L*PDstandardNth2Xt1 + J32L*PDstandardNth3Xt1; + JacPDstandardNth2Xt1 = J12L*PDstandardNth1Xt1 + J22L*PDstandardNth2Xt1 + + J32L*PDstandardNth3Xt1; - JacPDstandardNth2Xt2 = J12L*PDstandardNth1Xt2 + - J22L*PDstandardNth2Xt2 + J32L*PDstandardNth3Xt2; + JacPDstandardNth2Xt2 = J12L*PDstandardNth1Xt2 + J22L*PDstandardNth2Xt2 + + J32L*PDstandardNth3Xt2; - JacPDstandardNth2Xt3 = J12L*PDstandardNth1Xt3 + - J22L*PDstandardNth2Xt3 + J32L*PDstandardNth3Xt3; + JacPDstandardNth2Xt3 = J12L*PDstandardNth1Xt3 + J22L*PDstandardNth2Xt3 + + J32L*PDstandardNth3Xt3; JacPDstandardNth3gt11 = J13L*PDstandardNth1gt11 + J23L*PDstandardNth2gt11 + J33L*PDstandardNth3gt11; @@ -500,17 +502,17 @@ static void ML_ADMQuantities_Body(const cGH* restrict const cctkGH, const int di JacPDstandardNth3gt33 = J13L*PDstandardNth1gt33 + J23L*PDstandardNth2gt33 + J33L*PDstandardNth3gt33; - JacPDstandardNth3trK = J13L*PDstandardNth1trK + - J23L*PDstandardNth2trK + J33L*PDstandardNth3trK; + JacPDstandardNth3trK = J13L*PDstandardNth1trK + J23L*PDstandardNth2trK + + J33L*PDstandardNth3trK; - JacPDstandardNth3Xt1 = J13L*PDstandardNth1Xt1 + - J23L*PDstandardNth2Xt1 + J33L*PDstandardNth3Xt1; + JacPDstandardNth3Xt1 = J13L*PDstandardNth1Xt1 + J23L*PDstandardNth2Xt1 + + J33L*PDstandardNth3Xt1; - JacPDstandardNth3Xt2 = J13L*PDstandardNth1Xt2 + - J23L*PDstandardNth2Xt2 + J33L*PDstandardNth3Xt2; + JacPDstandardNth3Xt2 = J13L*PDstandardNth1Xt2 + J23L*PDstandardNth2Xt2 + + J33L*PDstandardNth3Xt2; - JacPDstandardNth3Xt3 = J13L*PDstandardNth1Xt3 + - J23L*PDstandardNth2Xt3 + J33L*PDstandardNth3Xt3; + JacPDstandardNth3Xt3 = J13L*PDstandardNth1Xt3 + J23L*PDstandardNth2Xt3 + + J33L*PDstandardNth3Xt3; JacPDstandardNth11gt11 = dJ111L*PDstandardNth1gt11 + 2*(J11L*(J21L*PDstandardNth12gt11 + J31L*PDstandardNth13gt11) + @@ -1079,324 +1081,324 @@ static void ML_ADMQuantities_Body(const cGH* restrict const cctkGH, const int di JacPDstandardNth32gt33 = PDstandardNth23gt33; } - CCTK_REAL CCTK_ATTRIBUTE_UNUSED detgt = 1; + CCTK_REAL detgt CCTK_ATTRIBUTE_UNUSED = 1; - CCTK_REAL CCTK_ATTRIBUTE_UNUSED gtu11 = INV(detgt)*(gt22L*gt33L - + CCTK_REAL gtu11 CCTK_ATTRIBUTE_UNUSED = INV(detgt)*(gt22L*gt33L - SQR(gt23L)); - CCTK_REAL CCTK_ATTRIBUTE_UNUSED gtu21 = (gt13L*gt23L - + CCTK_REAL gtu21 CCTK_ATTRIBUTE_UNUSED = (gt13L*gt23L - gt12L*gt33L)*INV(detgt); - CCTK_REAL CCTK_ATTRIBUTE_UNUSED gtu31 = (-(gt13L*gt22L) + + CCTK_REAL gtu31 CCTK_ATTRIBUTE_UNUSED = (-(gt13L*gt22L) + gt12L*gt23L)*INV(detgt); - CCTK_REAL CCTK_ATTRIBUTE_UNUSED gtu22 = INV(detgt)*(gt11L*gt33L - + CCTK_REAL gtu22 CCTK_ATTRIBUTE_UNUSED = INV(detgt)*(gt11L*gt33L - SQR(gt13L)); - CCTK_REAL CCTK_ATTRIBUTE_UNUSED gtu32 = (gt12L*gt13L - + CCTK_REAL gtu32 CCTK_ATTRIBUTE_UNUSED = (gt12L*gt13L - gt11L*gt23L)*INV(detgt); - CCTK_REAL CCTK_ATTRIBUTE_UNUSED gtu33 = INV(detgt)*(gt11L*gt22L - + CCTK_REAL gtu33 CCTK_ATTRIBUTE_UNUSED = INV(detgt)*(gt11L*gt22L - SQR(gt12L)); - CCTK_REAL CCTK_ATTRIBUTE_UNUSED dgtu111 = + CCTK_REAL dgtu111 CCTK_ATTRIBUTE_UNUSED = -2*(gtu11*(gtu21*JacPDstandardNth1gt12 + gtu31*JacPDstandardNth1gt13) + gtu21*gtu31*JacPDstandardNth1gt23) - JacPDstandardNth1gt11*SQR(gtu11) - JacPDstandardNth1gt22*SQR(gtu21) - JacPDstandardNth1gt33*SQR(gtu31); - CCTK_REAL CCTK_ATTRIBUTE_UNUSED dgtu211 = + CCTK_REAL dgtu211 CCTK_ATTRIBUTE_UNUSED = -(gtu21*(gtu11*JacPDstandardNth1gt11 + gtu21*JacPDstandardNth1gt12 + gtu31*JacPDstandardNth1gt13)) - gtu22*(gtu11*JacPDstandardNth1gt12 + gtu21*JacPDstandardNth1gt22 + gtu31*JacPDstandardNth1gt23) - gtu32*(gtu11*JacPDstandardNth1gt13 + gtu21*JacPDstandardNth1gt23 + gtu31*JacPDstandardNth1gt33); - CCTK_REAL CCTK_ATTRIBUTE_UNUSED dgtu311 = + CCTK_REAL dgtu311 CCTK_ATTRIBUTE_UNUSED = -(gtu11*(gtu31*JacPDstandardNth1gt11 + gtu32*JacPDstandardNth1gt12 + gtu33*JacPDstandardNth1gt13)) - gtu21*(gtu31*JacPDstandardNth1gt12 + gtu32*JacPDstandardNth1gt22 + gtu33*JacPDstandardNth1gt23) - gtu31*(gtu31*JacPDstandardNth1gt13 + gtu32*JacPDstandardNth1gt23 + gtu33*JacPDstandardNth1gt33); - CCTK_REAL CCTK_ATTRIBUTE_UNUSED dgtu221 = + CCTK_REAL dgtu221 CCTK_ATTRIBUTE_UNUSED = -2*(gtu21*(gtu22*JacPDstandardNth1gt12 + gtu32*JacPDstandardNth1gt13) + gtu22*gtu32*JacPDstandardNth1gt23) - JacPDstandardNth1gt11*SQR(gtu21) - JacPDstandardNth1gt22*SQR(gtu22) - JacPDstandardNth1gt33*SQR(gtu32); - CCTK_REAL CCTK_ATTRIBUTE_UNUSED dgtu321 = + CCTK_REAL dgtu321 CCTK_ATTRIBUTE_UNUSED = -(gtu21*(gtu31*JacPDstandardNth1gt11 + gtu32*JacPDstandardNth1gt12 + gtu33*JacPDstandardNth1gt13)) - gtu22*(gtu31*JacPDstandardNth1gt12 + gtu32*JacPDstandardNth1gt22 + gtu33*JacPDstandardNth1gt23) - gtu32*(gtu31*JacPDstandardNth1gt13 + gtu32*JacPDstandardNth1gt23 + gtu33*JacPDstandardNth1gt33); - CCTK_REAL CCTK_ATTRIBUTE_UNUSED dgtu331 = + CCTK_REAL dgtu331 CCTK_ATTRIBUTE_UNUSED = -2*(gtu31*(gtu32*JacPDstandardNth1gt12 + gtu33*JacPDstandardNth1gt13) + gtu32*gtu33*JacPDstandardNth1gt23) - JacPDstandardNth1gt11*SQR(gtu31) - JacPDstandardNth1gt22*SQR(gtu32) - JacPDstandardNth1gt33*SQR(gtu33); - CCTK_REAL CCTK_ATTRIBUTE_UNUSED dgtu112 = + CCTK_REAL dgtu112 CCTK_ATTRIBUTE_UNUSED = -2*(gtu11*(gtu21*JacPDstandardNth2gt12 + gtu31*JacPDstandardNth2gt13) + gtu21*gtu31*JacPDstandardNth2gt23) - JacPDstandardNth2gt11*SQR(gtu11) - JacPDstandardNth2gt22*SQR(gtu21) - JacPDstandardNth2gt33*SQR(gtu31); - CCTK_REAL CCTK_ATTRIBUTE_UNUSED dgtu212 = + CCTK_REAL dgtu212 CCTK_ATTRIBUTE_UNUSED = -(gtu21*(gtu11*JacPDstandardNth2gt11 + gtu21*JacPDstandardNth2gt12 + gtu31*JacPDstandardNth2gt13)) - gtu22*(gtu11*JacPDstandardNth2gt12 + gtu21*JacPDstandardNth2gt22 + gtu31*JacPDstandardNth2gt23) - gtu32*(gtu11*JacPDstandardNth2gt13 + gtu21*JacPDstandardNth2gt23 + gtu31*JacPDstandardNth2gt33); - CCTK_REAL CCTK_ATTRIBUTE_UNUSED dgtu312 = + CCTK_REAL dgtu312 CCTK_ATTRIBUTE_UNUSED = -(gtu11*(gtu31*JacPDstandardNth2gt11 + gtu32*JacPDstandardNth2gt12 + gtu33*JacPDstandardNth2gt13)) - gtu21*(gtu31*JacPDstandardNth2gt12 + gtu32*JacPDstandardNth2gt22 + gtu33*JacPDstandardNth2gt23) - gtu31*(gtu31*JacPDstandardNth2gt13 + gtu32*JacPDstandardNth2gt23 + gtu33*JacPDstandardNth2gt33); - CCTK_REAL CCTK_ATTRIBUTE_UNUSED dgtu222 = + CCTK_REAL dgtu222 CCTK_ATTRIBUTE_UNUSED = -2*(gtu21*(gtu22*JacPDstandardNth2gt12 + gtu32*JacPDstandardNth2gt13) + gtu22*gtu32*JacPDstandardNth2gt23) - JacPDstandardNth2gt11*SQR(gtu21) - JacPDstandardNth2gt22*SQR(gtu22) - JacPDstandardNth2gt33*SQR(gtu32); - CCTK_REAL CCTK_ATTRIBUTE_UNUSED dgtu322 = + CCTK_REAL dgtu322 CCTK_ATTRIBUTE_UNUSED = -(gtu21*(gtu31*JacPDstandardNth2gt11 + gtu32*JacPDstandardNth2gt12 + gtu33*JacPDstandardNth2gt13)) - gtu22*(gtu31*JacPDstandardNth2gt12 + gtu32*JacPDstandardNth2gt22 + gtu33*JacPDstandardNth2gt23) - gtu32*(gtu31*JacPDstandardNth2gt13 + gtu32*JacPDstandardNth2gt23 + gtu33*JacPDstandardNth2gt33); - CCTK_REAL CCTK_ATTRIBUTE_UNUSED dgtu332 = + CCTK_REAL dgtu332 CCTK_ATTRIBUTE_UNUSED = -2*(gtu31*(gtu32*JacPDstandardNth2gt12 + gtu33*JacPDstandardNth2gt13) + gtu32*gtu33*JacPDstandardNth2gt23) - JacPDstandardNth2gt11*SQR(gtu31) - JacPDstandardNth2gt22*SQR(gtu32) - JacPDstandardNth2gt33*SQR(gtu33); - CCTK_REAL CCTK_ATTRIBUTE_UNUSED dgtu113 = + CCTK_REAL dgtu113 CCTK_ATTRIBUTE_UNUSED = -2*(gtu11*(gtu21*JacPDstandardNth3gt12 + gtu31*JacPDstandardNth3gt13) + gtu21*gtu31*JacPDstandardNth3gt23) - JacPDstandardNth3gt11*SQR(gtu11) - JacPDstandardNth3gt22*SQR(gtu21) - JacPDstandardNth3gt33*SQR(gtu31); - CCTK_REAL CCTK_ATTRIBUTE_UNUSED dgtu213 = + CCTK_REAL dgtu213 CCTK_ATTRIBUTE_UNUSED = -(gtu21*(gtu11*JacPDstandardNth3gt11 + gtu21*JacPDstandardNth3gt12 + gtu31*JacPDstandardNth3gt13)) - gtu22*(gtu11*JacPDstandardNth3gt12 + gtu21*JacPDstandardNth3gt22 + gtu31*JacPDstandardNth3gt23) - gtu32*(gtu11*JacPDstandardNth3gt13 + gtu21*JacPDstandardNth3gt23 + gtu31*JacPDstandardNth3gt33); - CCTK_REAL CCTK_ATTRIBUTE_UNUSED dgtu313 = + CCTK_REAL dgtu313 CCTK_ATTRIBUTE_UNUSED = -(gtu11*(gtu31*JacPDstandardNth3gt11 + gtu32*JacPDstandardNth3gt12 + gtu33*JacPDstandardNth3gt13)) - gtu21*(gtu31*JacPDstandardNth3gt12 + gtu32*JacPDstandardNth3gt22 + gtu33*JacPDstandardNth3gt23) - gtu31*(gtu31*JacPDstandardNth3gt13 + gtu32*JacPDstandardNth3gt23 + gtu33*JacPDstandardNth3gt33); - CCTK_REAL CCTK_ATTRIBUTE_UNUSED dgtu223 = + CCTK_REAL dgtu223 CCTK_ATTRIBUTE_UNUSED = -2*(gtu21*(gtu22*JacPDstandardNth3gt12 + gtu32*JacPDstandardNth3gt13) + gtu22*gtu32*JacPDstandardNth3gt23) - JacPDstandardNth3gt11*SQR(gtu21) - JacPDstandardNth3gt22*SQR(gtu22) - JacPDstandardNth3gt33*SQR(gtu32); - CCTK_REAL CCTK_ATTRIBUTE_UNUSED dgtu323 = + CCTK_REAL dgtu323 CCTK_ATTRIBUTE_UNUSED = -(gtu21*(gtu31*JacPDstandardNth3gt11 + gtu32*JacPDstandardNth3gt12 + gtu33*JacPDstandardNth3gt13)) - gtu22*(gtu31*JacPDstandardNth3gt12 + gtu32*JacPDstandardNth3gt22 + gtu33*JacPDstandardNth3gt23) - gtu32*(gtu31*JacPDstandardNth3gt13 + gtu32*JacPDstandardNth3gt23 + gtu33*JacPDstandardNth3gt33); - CCTK_REAL CCTK_ATTRIBUTE_UNUSED dgtu333 = + CCTK_REAL dgtu333 CCTK_ATTRIBUTE_UNUSED = -2*(gtu31*(gtu32*JacPDstandardNth3gt12 + gtu33*JacPDstandardNth3gt13) + gtu32*gtu33*JacPDstandardNth3gt23) - JacPDstandardNth3gt11*SQR(gtu31) - JacPDstandardNth3gt22*SQR(gtu32) - JacPDstandardNth3gt33*SQR(gtu33); - CCTK_REAL CCTK_ATTRIBUTE_UNUSED Gtl111 = 0.5*JacPDstandardNth1gt11; + CCTK_REAL Gtl111 CCTK_ATTRIBUTE_UNUSED = 0.5*JacPDstandardNth1gt11; - CCTK_REAL CCTK_ATTRIBUTE_UNUSED Gtl112 = 0.5*JacPDstandardNth2gt11; + CCTK_REAL Gtl112 CCTK_ATTRIBUTE_UNUSED = 0.5*JacPDstandardNth2gt11; - CCTK_REAL CCTK_ATTRIBUTE_UNUSED Gtl113 = 0.5*JacPDstandardNth3gt11; + CCTK_REAL Gtl113 CCTK_ATTRIBUTE_UNUSED = 0.5*JacPDstandardNth3gt11; - CCTK_REAL CCTK_ATTRIBUTE_UNUSED Gtl122 = -0.5*JacPDstandardNth1gt22 + + CCTK_REAL Gtl122 CCTK_ATTRIBUTE_UNUSED = -0.5*JacPDstandardNth1gt22 + JacPDstandardNth2gt12; - CCTK_REAL CCTK_ATTRIBUTE_UNUSED Gtl123 = 0.5*(-JacPDstandardNth1gt23 + + CCTK_REAL Gtl123 CCTK_ATTRIBUTE_UNUSED = 0.5*(-JacPDstandardNth1gt23 + JacPDstandardNth2gt13 + JacPDstandardNth3gt12); - CCTK_REAL CCTK_ATTRIBUTE_UNUSED Gtl133 = -0.5*JacPDstandardNth1gt33 + + CCTK_REAL Gtl133 CCTK_ATTRIBUTE_UNUSED = -0.5*JacPDstandardNth1gt33 + JacPDstandardNth3gt13; - CCTK_REAL CCTK_ATTRIBUTE_UNUSED Gtl211 = JacPDstandardNth1gt12 - + CCTK_REAL Gtl211 CCTK_ATTRIBUTE_UNUSED = JacPDstandardNth1gt12 - 0.5*JacPDstandardNth2gt11; - CCTK_REAL CCTK_ATTRIBUTE_UNUSED Gtl212 = 0.5*JacPDstandardNth1gt22; + CCTK_REAL Gtl212 CCTK_ATTRIBUTE_UNUSED = 0.5*JacPDstandardNth1gt22; - CCTK_REAL CCTK_ATTRIBUTE_UNUSED Gtl213 = 0.5*(JacPDstandardNth1gt23 - + CCTK_REAL Gtl213 CCTK_ATTRIBUTE_UNUSED = 0.5*(JacPDstandardNth1gt23 - JacPDstandardNth2gt13 + JacPDstandardNth3gt12); - CCTK_REAL CCTK_ATTRIBUTE_UNUSED Gtl222 = 0.5*JacPDstandardNth2gt22; + CCTK_REAL Gtl222 CCTK_ATTRIBUTE_UNUSED = 0.5*JacPDstandardNth2gt22; - CCTK_REAL CCTK_ATTRIBUTE_UNUSED Gtl223 = 0.5*JacPDstandardNth3gt22; + CCTK_REAL Gtl223 CCTK_ATTRIBUTE_UNUSED = 0.5*JacPDstandardNth3gt22; - CCTK_REAL CCTK_ATTRIBUTE_UNUSED Gtl233 = -0.5*JacPDstandardNth2gt33 + + CCTK_REAL Gtl233 CCTK_ATTRIBUTE_UNUSED = -0.5*JacPDstandardNth2gt33 + JacPDstandardNth3gt23; - CCTK_REAL CCTK_ATTRIBUTE_UNUSED Gtl311 = JacPDstandardNth1gt13 - + CCTK_REAL Gtl311 CCTK_ATTRIBUTE_UNUSED = JacPDstandardNth1gt13 - 0.5*JacPDstandardNth3gt11; - CCTK_REAL CCTK_ATTRIBUTE_UNUSED Gtl312 = 0.5*(JacPDstandardNth1gt23 + + CCTK_REAL Gtl312 CCTK_ATTRIBUTE_UNUSED = 0.5*(JacPDstandardNth1gt23 + JacPDstandardNth2gt13 - JacPDstandardNth3gt12); - CCTK_REAL CCTK_ATTRIBUTE_UNUSED Gtl313 = 0.5*JacPDstandardNth1gt33; + CCTK_REAL Gtl313 CCTK_ATTRIBUTE_UNUSED = 0.5*JacPDstandardNth1gt33; - CCTK_REAL CCTK_ATTRIBUTE_UNUSED Gtl322 = JacPDstandardNth2gt23 - + CCTK_REAL Gtl322 CCTK_ATTRIBUTE_UNUSED = JacPDstandardNth2gt23 - 0.5*JacPDstandardNth3gt22; - CCTK_REAL CCTK_ATTRIBUTE_UNUSED Gtl323 = 0.5*JacPDstandardNth2gt33; + CCTK_REAL Gtl323 CCTK_ATTRIBUTE_UNUSED = 0.5*JacPDstandardNth2gt33; - CCTK_REAL CCTK_ATTRIBUTE_UNUSED Gtl333 = 0.5*JacPDstandardNth3gt33; + CCTK_REAL Gtl333 CCTK_ATTRIBUTE_UNUSED = 0.5*JacPDstandardNth3gt33; - CCTK_REAL CCTK_ATTRIBUTE_UNUSED Gtlu111 = Gtl111*gtu11 + Gtl112*gtu21 + CCTK_REAL Gtlu111 CCTK_ATTRIBUTE_UNUSED = Gtl111*gtu11 + Gtl112*gtu21 + Gtl113*gtu31; - CCTK_REAL CCTK_ATTRIBUTE_UNUSED Gtlu112 = Gtl111*gtu21 + Gtl112*gtu22 + CCTK_REAL Gtlu112 CCTK_ATTRIBUTE_UNUSED = Gtl111*gtu21 + Gtl112*gtu22 + Gtl113*gtu32; - CCTK_REAL CCTK_ATTRIBUTE_UNUSED Gtlu113 = Gtl111*gtu31 + Gtl112*gtu32 + CCTK_REAL Gtlu113 CCTK_ATTRIBUTE_UNUSED = Gtl111*gtu31 + Gtl112*gtu32 + Gtl113*gtu33; - CCTK_REAL CCTK_ATTRIBUTE_UNUSED Gtlu121 = Gtl112*gtu11 + Gtl122*gtu21 + CCTK_REAL Gtlu121 CCTK_ATTRIBUTE_UNUSED = Gtl112*gtu11 + Gtl122*gtu21 + Gtl123*gtu31; - CCTK_REAL CCTK_ATTRIBUTE_UNUSED Gtlu122 = Gtl112*gtu21 + Gtl122*gtu22 + CCTK_REAL Gtlu122 CCTK_ATTRIBUTE_UNUSED = Gtl112*gtu21 + Gtl122*gtu22 + Gtl123*gtu32; - CCTK_REAL CCTK_ATTRIBUTE_UNUSED Gtlu123 = Gtl112*gtu31 + Gtl122*gtu32 + CCTK_REAL Gtlu123 CCTK_ATTRIBUTE_UNUSED = Gtl112*gtu31 + Gtl122*gtu32 + Gtl123*gtu33; - CCTK_REAL CCTK_ATTRIBUTE_UNUSED Gtlu131 = Gtl113*gtu11 + Gtl123*gtu21 + CCTK_REAL Gtlu131 CCTK_ATTRIBUTE_UNUSED = Gtl113*gtu11 + Gtl123*gtu21 + Gtl133*gtu31; - CCTK_REAL CCTK_ATTRIBUTE_UNUSED Gtlu132 = Gtl113*gtu21 + Gtl123*gtu22 + CCTK_REAL Gtlu132 CCTK_ATTRIBUTE_UNUSED = Gtl113*gtu21 + Gtl123*gtu22 + Gtl133*gtu32; - CCTK_REAL CCTK_ATTRIBUTE_UNUSED Gtlu133 = Gtl113*gtu31 + Gtl123*gtu32 + CCTK_REAL Gtlu133 CCTK_ATTRIBUTE_UNUSED = Gtl113*gtu31 + Gtl123*gtu32 + Gtl133*gtu33; - CCTK_REAL CCTK_ATTRIBUTE_UNUSED Gtlu211 = Gtl211*gtu11 + Gtl212*gtu21 + CCTK_REAL Gtlu211 CCTK_ATTRIBUTE_UNUSED = Gtl211*gtu11 + Gtl212*gtu21 + Gtl213*gtu31; - CCTK_REAL CCTK_ATTRIBUTE_UNUSED Gtlu212 = Gtl211*gtu21 + Gtl212*gtu22 + CCTK_REAL Gtlu212 CCTK_ATTRIBUTE_UNUSED = Gtl211*gtu21 + Gtl212*gtu22 + Gtl213*gtu32; - CCTK_REAL CCTK_ATTRIBUTE_UNUSED Gtlu213 = Gtl211*gtu31 + Gtl212*gtu32 + CCTK_REAL Gtlu213 CCTK_ATTRIBUTE_UNUSED = Gtl211*gtu31 + Gtl212*gtu32 + Gtl213*gtu33; - CCTK_REAL CCTK_ATTRIBUTE_UNUSED Gtlu221 = Gtl212*gtu11 + Gtl222*gtu21 + CCTK_REAL Gtlu221 CCTK_ATTRIBUTE_UNUSED = Gtl212*gtu11 + Gtl222*gtu21 + Gtl223*gtu31; - CCTK_REAL CCTK_ATTRIBUTE_UNUSED Gtlu222 = Gtl212*gtu21 + Gtl222*gtu22 + CCTK_REAL Gtlu222 CCTK_ATTRIBUTE_UNUSED = Gtl212*gtu21 + Gtl222*gtu22 + Gtl223*gtu32; - CCTK_REAL CCTK_ATTRIBUTE_UNUSED Gtlu223 = Gtl212*gtu31 + Gtl222*gtu32 + CCTK_REAL Gtlu223 CCTK_ATTRIBUTE_UNUSED = Gtl212*gtu31 + Gtl222*gtu32 + Gtl223*gtu33; - CCTK_REAL CCTK_ATTRIBUTE_UNUSED Gtlu231 = Gtl213*gtu11 + Gtl223*gtu21 + CCTK_REAL Gtlu231 CCTK_ATTRIBUTE_UNUSED = Gtl213*gtu11 + Gtl223*gtu21 + Gtl233*gtu31; - CCTK_REAL CCTK_ATTRIBUTE_UNUSED Gtlu232 = Gtl213*gtu21 + Gtl223*gtu22 + CCTK_REAL Gtlu232 CCTK_ATTRIBUTE_UNUSED = Gtl213*gtu21 + Gtl223*gtu22 + Gtl233*gtu32; - CCTK_REAL CCTK_ATTRIBUTE_UNUSED Gtlu233 = Gtl213*gtu31 + Gtl223*gtu32 + CCTK_REAL Gtlu233 CCTK_ATTRIBUTE_UNUSED = Gtl213*gtu31 + Gtl223*gtu32 + Gtl233*gtu33; - CCTK_REAL CCTK_ATTRIBUTE_UNUSED Gtlu311 = Gtl311*gtu11 + Gtl312*gtu21 + CCTK_REAL Gtlu311 CCTK_ATTRIBUTE_UNUSED = Gtl311*gtu11 + Gtl312*gtu21 + Gtl313*gtu31; - CCTK_REAL CCTK_ATTRIBUTE_UNUSED Gtlu312 = Gtl311*gtu21 + Gtl312*gtu22 + CCTK_REAL Gtlu312 CCTK_ATTRIBUTE_UNUSED = Gtl311*gtu21 + Gtl312*gtu22 + Gtl313*gtu32; - CCTK_REAL CCTK_ATTRIBUTE_UNUSED Gtlu313 = Gtl311*gtu31 + Gtl312*gtu32 + CCTK_REAL Gtlu313 CCTK_ATTRIBUTE_UNUSED = Gtl311*gtu31 + Gtl312*gtu32 + Gtl313*gtu33; - CCTK_REAL CCTK_ATTRIBUTE_UNUSED Gtlu321 = Gtl312*gtu11 + Gtl322*gtu21 + CCTK_REAL Gtlu321 CCTK_ATTRIBUTE_UNUSED = Gtl312*gtu11 + Gtl322*gtu21 + Gtl323*gtu31; - CCTK_REAL CCTK_ATTRIBUTE_UNUSED Gtlu322 = Gtl312*gtu21 + Gtl322*gtu22 + CCTK_REAL Gtlu322 CCTK_ATTRIBUTE_UNUSED = Gtl312*gtu21 + Gtl322*gtu22 + Gtl323*gtu32; - CCTK_REAL CCTK_ATTRIBUTE_UNUSED Gtlu323 = Gtl312*gtu31 + Gtl322*gtu32 + CCTK_REAL Gtlu323 CCTK_ATTRIBUTE_UNUSED = Gtl312*gtu31 + Gtl322*gtu32 + Gtl323*gtu33; - CCTK_REAL CCTK_ATTRIBUTE_UNUSED Gtlu331 = Gtl313*gtu11 + Gtl323*gtu21 + CCTK_REAL Gtlu331 CCTK_ATTRIBUTE_UNUSED = Gtl313*gtu11 + Gtl323*gtu21 + Gtl333*gtu31; - CCTK_REAL CCTK_ATTRIBUTE_UNUSED Gtlu332 = Gtl313*gtu21 + Gtl323*gtu22 + CCTK_REAL Gtlu332 CCTK_ATTRIBUTE_UNUSED = Gtl313*gtu21 + Gtl323*gtu22 + Gtl333*gtu32; - CCTK_REAL CCTK_ATTRIBUTE_UNUSED Gtlu333 = Gtl313*gtu31 + Gtl323*gtu32 + CCTK_REAL Gtlu333 CCTK_ATTRIBUTE_UNUSED = Gtl313*gtu31 + Gtl323*gtu32 + Gtl333*gtu33; - CCTK_REAL CCTK_ATTRIBUTE_UNUSED Gt111 = Gtl111*gtu11 + Gtl211*gtu21 + + CCTK_REAL Gt111 CCTK_ATTRIBUTE_UNUSED = Gtl111*gtu11 + Gtl211*gtu21 + Gtl311*gtu31; - CCTK_REAL CCTK_ATTRIBUTE_UNUSED Gt211 = Gtl111*gtu21 + Gtl211*gtu22 + + CCTK_REAL Gt211 CCTK_ATTRIBUTE_UNUSED = Gtl111*gtu21 + Gtl211*gtu22 + Gtl311*gtu32; - CCTK_REAL CCTK_ATTRIBUTE_UNUSED Gt311 = Gtl111*gtu31 + Gtl211*gtu32 + + CCTK_REAL Gt311 CCTK_ATTRIBUTE_UNUSED = Gtl111*gtu31 + Gtl211*gtu32 + Gtl311*gtu33; - CCTK_REAL CCTK_ATTRIBUTE_UNUSED Gt112 = Gtl112*gtu11 + Gtl212*gtu21 + + CCTK_REAL Gt112 CCTK_ATTRIBUTE_UNUSED = Gtl112*gtu11 + Gtl212*gtu21 + Gtl312*gtu31; - CCTK_REAL CCTK_ATTRIBUTE_UNUSED Gt212 = Gtl112*gtu21 + Gtl212*gtu22 + + CCTK_REAL Gt212 CCTK_ATTRIBUTE_UNUSED = Gtl112*gtu21 + Gtl212*gtu22 + Gtl312*gtu32; - CCTK_REAL CCTK_ATTRIBUTE_UNUSED Gt312 = Gtl112*gtu31 + Gtl212*gtu32 + + CCTK_REAL Gt312 CCTK_ATTRIBUTE_UNUSED = Gtl112*gtu31 + Gtl212*gtu32 + Gtl312*gtu33; - CCTK_REAL CCTK_ATTRIBUTE_UNUSED Gt113 = Gtl113*gtu11 + Gtl213*gtu21 + + CCTK_REAL Gt113 CCTK_ATTRIBUTE_UNUSED = Gtl113*gtu11 + Gtl213*gtu21 + Gtl313*gtu31; - CCTK_REAL CCTK_ATTRIBUTE_UNUSED Gt213 = Gtl113*gtu21 + Gtl213*gtu22 + + CCTK_REAL Gt213 CCTK_ATTRIBUTE_UNUSED = Gtl113*gtu21 + Gtl213*gtu22 + Gtl313*gtu32; - CCTK_REAL CCTK_ATTRIBUTE_UNUSED Gt313 = Gtl113*gtu31 + Gtl213*gtu32 + + CCTK_REAL Gt313 CCTK_ATTRIBUTE_UNUSED = Gtl113*gtu31 + Gtl213*gtu32 + Gtl313*gtu33; - CCTK_REAL CCTK_ATTRIBUTE_UNUSED Gt122 = Gtl122*gtu11 + Gtl222*gtu21 + + CCTK_REAL Gt122 CCTK_ATTRIBUTE_UNUSED = Gtl122*gtu11 + Gtl222*gtu21 + Gtl322*gtu31; - CCTK_REAL CCTK_ATTRIBUTE_UNUSED Gt222 = Gtl122*gtu21 + Gtl222*gtu22 + + CCTK_REAL Gt222 CCTK_ATTRIBUTE_UNUSED = Gtl122*gtu21 + Gtl222*gtu22 + Gtl322*gtu32; - CCTK_REAL CCTK_ATTRIBUTE_UNUSED Gt322 = Gtl122*gtu31 + Gtl222*gtu32 + + CCTK_REAL Gt322 CCTK_ATTRIBUTE_UNUSED = Gtl122*gtu31 + Gtl222*gtu32 + Gtl322*gtu33; - CCTK_REAL CCTK_ATTRIBUTE_UNUSED Gt123 = Gtl123*gtu11 + Gtl223*gtu21 + + CCTK_REAL Gt123 CCTK_ATTRIBUTE_UNUSED = Gtl123*gtu11 + Gtl223*gtu21 + Gtl323*gtu31; - CCTK_REAL CCTK_ATTRIBUTE_UNUSED Gt223 = Gtl123*gtu21 + Gtl223*gtu22 + + CCTK_REAL Gt223 CCTK_ATTRIBUTE_UNUSED = Gtl123*gtu21 + Gtl223*gtu22 + Gtl323*gtu32; - CCTK_REAL CCTK_ATTRIBUTE_UNUSED Gt323 = Gtl123*gtu31 + Gtl223*gtu32 + + CCTK_REAL Gt323 CCTK_ATTRIBUTE_UNUSED = Gtl123*gtu31 + Gtl223*gtu32 + Gtl323*gtu33; - CCTK_REAL CCTK_ATTRIBUTE_UNUSED Gt133 = Gtl133*gtu11 + Gtl233*gtu21 + + CCTK_REAL Gt133 CCTK_ATTRIBUTE_UNUSED = Gtl133*gtu11 + Gtl233*gtu21 + Gtl333*gtu31; - CCTK_REAL CCTK_ATTRIBUTE_UNUSED Gt233 = Gtl133*gtu21 + Gtl233*gtu22 + + CCTK_REAL Gt233 CCTK_ATTRIBUTE_UNUSED = Gtl133*gtu21 + Gtl233*gtu22 + Gtl333*gtu32; - CCTK_REAL CCTK_ATTRIBUTE_UNUSED Gt333 = Gtl133*gtu31 + Gtl233*gtu32 + + CCTK_REAL Gt333 CCTK_ATTRIBUTE_UNUSED = Gtl133*gtu31 + Gtl233*gtu32 + Gtl333*gtu33; - CCTK_REAL CCTK_ATTRIBUTE_UNUSED Xtn1 = Gt111*gtu11 + Gt122*gtu22 + + CCTK_REAL Xtn1 CCTK_ATTRIBUTE_UNUSED = Gt111*gtu11 + Gt122*gtu22 + 2*(Gt112*gtu21 + Gt113*gtu31 + Gt123*gtu32) + Gt133*gtu33; - CCTK_REAL CCTK_ATTRIBUTE_UNUSED Xtn2 = Gt211*gtu11 + Gt222*gtu22 + + CCTK_REAL Xtn2 CCTK_ATTRIBUTE_UNUSED = Gt211*gtu11 + Gt222*gtu22 + 2*(Gt212*gtu21 + Gt213*gtu31 + Gt223*gtu32) + Gt233*gtu33; - CCTK_REAL CCTK_ATTRIBUTE_UNUSED Xtn3 = Gt311*gtu11 + Gt322*gtu22 + + CCTK_REAL Xtn3 CCTK_ATTRIBUTE_UNUSED = Gt311*gtu11 + Gt322*gtu22 + 2*(Gt312*gtu21 + Gt313*gtu31 + Gt323*gtu32) + Gt333*gtu33; - CCTK_REAL CCTK_ATTRIBUTE_UNUSED Rt11 = 3*(Gt111*Gtlu111 + + CCTK_REAL Rt11 CCTK_ATTRIBUTE_UNUSED = 3*(Gt111*Gtlu111 + Gt112*Gtlu112 + Gt113*Gtlu113) + 2*(Gt211*Gtlu121 + Gt212*Gtlu122 + Gt213*Gtlu123 + Gt311*Gtlu131 + Gt312*Gtlu132 + Gt313*Gtlu133) + Gt211*Gtlu211 + Gt212*Gtlu212 + Gt213*Gtlu213 + Gt311*Gtlu311 + @@ -1409,7 +1411,7 @@ static void ML_ADMQuantities_Body(const cGH* restrict const cctkGH, const int di gtu33*JacPDstandardNth33gt11) + Gtl111*Xtn1 + Gtl112*Xtn2 + Gtl113*Xtn3; - CCTK_REAL CCTK_ATTRIBUTE_UNUSED Rt12 = 0.5*(4*(Gt211*Gtlu221 + + CCTK_REAL Rt12 CCTK_ATTRIBUTE_UNUSED = 0.5*(4*(Gt211*Gtlu221 + Gt212*Gtlu222 + Gt213*Gtlu223) + 2*(Gt122*Gtlu112 + Gt123*Gtlu113 + Gt111*Gtlu121 + Gt212*Gtlu121 + Gt222*Gtlu122 + Gt113*Gtlu123 + Gt223*Gtlu123 + Gt312*Gtlu131 + Gt322*Gtlu132 + Gt323*Gtlu133 + @@ -1425,7 +1427,7 @@ static void ML_ADMQuantities_Body(const cGH* restrict const cctkGH, const int di JacPDstandardNth32gt12) - gtu33*JacPDstandardNth33gt12 + Gtl112*Xtn1 + Gtl211*Xtn1 + Gtl122*Xtn2 + Gtl212*Xtn2 + Gtl123*Xtn3 + Gtl213*Xtn3); - CCTK_REAL CCTK_ATTRIBUTE_UNUSED Rt13 = 0.5*(2*(Gt123*Gtlu112 + + CCTK_REAL Rt13 CCTK_ATTRIBUTE_UNUSED = 0.5*(2*(Gt123*Gtlu112 + Gt133*Gtlu113 + Gt213*Gtlu121 + Gt223*Gtlu122 + Gt233*Gtlu123 + Gt111*Gtlu131 + Gt313*Gtlu131 + Gt112*Gtlu132 + Gt323*Gtlu132 + Gt333*Gtlu133 + Gt211*Gtlu231 + Gt212*Gtlu232 + Gt213*Gtlu233 + @@ -1438,24 +1440,22 @@ static void ML_ADMQuantities_Body(const cGH* restrict const cctkGH, const int di gtu31*(-JacPDstandardNth13gt13 - JacPDstandardNth31gt13) + gtu32*(-JacPDstandardNth23gt13 - JacPDstandardNth32gt13) - gtu33*JacPDstandardNth33gt13 + gt11L*JacPDstandardNth3Xt1 + - gt12L*JacPDstandardNth3Xt2 + gt13L*JacPDstandardNth3Xt3 + - Gtl113*Xtn1 + Gtl311*Xtn1 + Gtl123*Xtn2 + Gtl312*Xtn2 + Gtl133*Xtn3 + - Gtl313*Xtn3); + gt12L*JacPDstandardNth3Xt2 + gt13L*JacPDstandardNth3Xt3 + Gtl113*Xtn1 + + Gtl311*Xtn1 + Gtl123*Xtn2 + Gtl312*Xtn2 + Gtl133*Xtn3 + Gtl313*Xtn3); - CCTK_REAL CCTK_ATTRIBUTE_UNUSED Rt22 = Gt112*(Gtlu121 + 2*Gtlu211) + + CCTK_REAL Rt22 CCTK_ATTRIBUTE_UNUSED = 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*JacPDstandardNth2Xt1 + - gt22L*JacPDstandardNth2Xt2 + gt23L*JacPDstandardNth2Xt3 + - 0.5*(-(gtu11*JacPDstandardNth11gt22) - gtu21*(JacPDstandardNth12gt22 + - JacPDstandardNth21gt22) - gtu22*JacPDstandardNth22gt22 + - gtu31*(-JacPDstandardNth13gt22 - JacPDstandardNth31gt22) + - gtu32*(-JacPDstandardNth23gt22 - JacPDstandardNth32gt22) - - gtu33*JacPDstandardNth33gt22) + Gtl212*Xtn1 + Gtl222*Xtn2 + - Gtl223*Xtn3; - - CCTK_REAL CCTK_ATTRIBUTE_UNUSED Rt23 = 0.5*(2*(Gt123*Gtlu133 + + Gt323*Gtlu323 + gt12L*JacPDstandardNth2Xt1 + gt22L*JacPDstandardNth2Xt2 + + gt23L*JacPDstandardNth2Xt3 + 0.5*(-(gtu11*JacPDstandardNth11gt22) - + gtu21*(JacPDstandardNth12gt22 + JacPDstandardNth21gt22) - + gtu22*JacPDstandardNth22gt22 + gtu31*(-JacPDstandardNth13gt22 - + JacPDstandardNth31gt22) + gtu32*(-JacPDstandardNth23gt22 - + JacPDstandardNth32gt22) - gtu33*JacPDstandardNth33gt22) + Gtl212*Xtn1 + + Gtl222*Xtn2 + Gtl223*Xtn3; + + CCTK_REAL Rt23 CCTK_ATTRIBUTE_UNUSED = 0.5*(2*(Gt123*Gtlu133 + Gt113*Gtlu211 + Gt123*Gtlu212 + Gt133*Gtlu213 + Gt213*Gtlu221 + Gt223*Gtlu222 + Gt233*Gtlu223 + Gt212*Gtlu231 + Gt313*Gtlu231 + Gt222*Gtlu232 + Gt323*Gtlu232 + Gt223*Gtlu233 + Gt333*Gtlu233 + @@ -1468,11 +1468,10 @@ static void ML_ADMQuantities_Body(const cGH* restrict const cctkGH, const int di gtu31*(-JacPDstandardNth13gt23 - JacPDstandardNth31gt23) + gtu32*(-JacPDstandardNth23gt23 - JacPDstandardNth32gt23) - gtu33*JacPDstandardNth33gt23 + gt12L*JacPDstandardNth3Xt1 + - gt22L*JacPDstandardNth3Xt2 + gt23L*JacPDstandardNth3Xt3 + - Gtl213*Xtn1 + Gtl312*Xtn1 + Gtl223*Xtn2 + Gtl322*Xtn2 + Gtl233*Xtn3 + - Gtl323*Xtn3); + gt22L*JacPDstandardNth3Xt2 + gt23L*JacPDstandardNth3Xt3 + Gtl213*Xtn1 + + Gtl312*Xtn1 + Gtl223*Xtn2 + Gtl322*Xtn2 + Gtl233*Xtn3 + Gtl323*Xtn3); - CCTK_REAL CCTK_ATTRIBUTE_UNUSED Rt33 = Gt113*(Gtlu131 + 2*Gtlu311) + + CCTK_REAL Rt33 CCTK_ATTRIBUTE_UNUSED = 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 + @@ -1484,54 +1483,54 @@ static void ML_ADMQuantities_Body(const cGH* restrict const cctkGH, const int di gt13L*JacPDstandardNth3Xt1 + gt23L*JacPDstandardNth3Xt2 + gt33L*JacPDstandardNth3Xt3 + Gtl313*Xtn1 + Gtl323*Xtn2 + Gtl333*Xtn3; - CCTK_REAL CCTK_ATTRIBUTE_UNUSED trRt = gtu11*Rt11 + gtu22*Rt22 + + CCTK_REAL trRt CCTK_ATTRIBUTE_UNUSED = gtu11*Rt11 + gtu22*Rt22 + 2*(gtu21*Rt12 + gtu31*Rt13 + gtu32*Rt23) + gtu33*Rt33; - CCTK_REAL CCTK_ATTRIBUTE_UNUSED ephi = + CCTK_REAL ephi CCTK_ATTRIBUTE_UNUSED = IfThen(conformalMethod,INV(sqrt(phiL)),exp(phiL)); - CCTK_REAL CCTK_ATTRIBUTE_UNUSED Atm11 = At11L*gtu11 + At12L*gtu21 - + At13L*gtu31; + CCTK_REAL Atm11 CCTK_ATTRIBUTE_UNUSED = At11L*gtu11 + At12L*gtu21 + + At13L*gtu31; - CCTK_REAL CCTK_ATTRIBUTE_UNUSED Atm21 = At11L*gtu21 + At12L*gtu22 - + At13L*gtu32; + CCTK_REAL Atm21 CCTK_ATTRIBUTE_UNUSED = At11L*gtu21 + At12L*gtu22 + + At13L*gtu32; - CCTK_REAL CCTK_ATTRIBUTE_UNUSED Atm31 = At11L*gtu31 + At12L*gtu32 - + At13L*gtu33; + CCTK_REAL Atm31 CCTK_ATTRIBUTE_UNUSED = At11L*gtu31 + At12L*gtu32 + + At13L*gtu33; - CCTK_REAL CCTK_ATTRIBUTE_UNUSED Atm12 = At12L*gtu11 + At22L*gtu21 - + At23L*gtu31; + CCTK_REAL Atm12 CCTK_ATTRIBUTE_UNUSED = At12L*gtu11 + At22L*gtu21 + + At23L*gtu31; - CCTK_REAL CCTK_ATTRIBUTE_UNUSED Atm22 = At12L*gtu21 + At22L*gtu22 - + At23L*gtu32; + CCTK_REAL Atm22 CCTK_ATTRIBUTE_UNUSED = At12L*gtu21 + At22L*gtu22 + + At23L*gtu32; - CCTK_REAL CCTK_ATTRIBUTE_UNUSED Atm32 = At12L*gtu31 + At22L*gtu32 - + At23L*gtu33; + CCTK_REAL Atm32 CCTK_ATTRIBUTE_UNUSED = At12L*gtu31 + At22L*gtu32 + + At23L*gtu33; - CCTK_REAL CCTK_ATTRIBUTE_UNUSED Atm13 = At13L*gtu11 + At23L*gtu21 - + At33L*gtu31; + CCTK_REAL Atm13 CCTK_ATTRIBUTE_UNUSED = At13L*gtu11 + At23L*gtu21 + + At33L*gtu31; - CCTK_REAL CCTK_ATTRIBUTE_UNUSED Atm23 = At13L*gtu21 + At23L*gtu22 - + At33L*gtu32; + CCTK_REAL Atm23 CCTK_ATTRIBUTE_UNUSED = At13L*gtu21 + At23L*gtu22 + + At33L*gtu32; - CCTK_REAL CCTK_ATTRIBUTE_UNUSED Atm33 = At13L*gtu31 + At23L*gtu32 - + At33L*gtu33; + CCTK_REAL Atm33 CCTK_ATTRIBUTE_UNUSED = At13L*gtu31 + At23L*gtu32 + + At33L*gtu33; - CCTK_REAL CCTK_ATTRIBUTE_UNUSED rho = INV(SQR(alphaL))*(eTttL - - 2*(beta2L*eTtyL + beta3L*eTtzL) + 2*(beta1L*(-eTtxL + - beta2L*eTxyL + beta3L*eTxzL) + beta2L*beta3L*eTyzL) + - eTxxL*SQR(beta1L) + eTyyL*SQR(beta2L) + eTzzL*SQR(beta3L)); + CCTK_REAL rho CCTK_ATTRIBUTE_UNUSED = INV(SQR(alphaL))*(eTttL - + 2*(beta2L*eTtyL + beta3L*eTtzL) + 2*(beta1L*(-eTtxL + beta2L*eTxyL + + beta3L*eTxzL) + beta2L*beta3L*eTyzL) + eTxxL*SQR(beta1L) + + eTyyL*SQR(beta2L) + eTzzL*SQR(beta3L)); - CCTK_REAL CCTK_ATTRIBUTE_UNUSED S1 = (-eTtxL + beta1L*eTxxL + + CCTK_REAL S1 CCTK_ATTRIBUTE_UNUSED = (-eTtxL + beta1L*eTxxL + beta2L*eTxyL + beta3L*eTxzL)*INV(alphaL); - CCTK_REAL CCTK_ATTRIBUTE_UNUSED S2 = (-eTtyL + beta1L*eTxyL + + CCTK_REAL S2 CCTK_ATTRIBUTE_UNUSED = (-eTtyL + beta1L*eTxyL + beta2L*eTyyL + beta3L*eTyzL)*INV(alphaL); - CCTK_REAL CCTK_ATTRIBUTE_UNUSED S3 = (-eTtzL + beta1L*eTxzL + + CCTK_REAL S3 CCTK_ATTRIBUTE_UNUSED = (-eTtzL + beta1L*eTxzL + beta2L*eTyzL + beta3L*eTzzL)*INV(alphaL); - CCTK_REAL CCTK_ATTRIBUTE_UNUSED MadmL = + CCTK_REAL MadmL CCTK_ATTRIBUTE_UNUSED = -0.0625*INV(Pi)*(Gtlu111*(Gt111*gtu11 + Gt112*gtu21 + Gt113*gtu31) + Gtlu112*(Gt112*gtu11 + Gt122*gtu21 + Gt123*gtu31) + Gtlu113*(Gt113*gtu11 + Gt123*gtu21 + Gt133*gtu31) + @@ -1560,30 +1559,30 @@ static void ML_ADMQuantities_Body(const cGH* restrict const cctkGH, const int di Gtlu332*(Gt312*gtu31 + Gt322*gtu32 + Gt323*gtu33) + Gtlu333*(Gt313*gtu31 + Gt323*gtu32 + Gt333*gtu33) + (-1 + ephi)*trRt - pow(ephi,5)*(2*(Atm12*Atm21 + Atm13*Atm31 + Atm23*Atm32) + 16*Pi*rho - - 0.666666666666666666666666666667*SQR(trKL) + SQR(Atm11) + SQR(Atm22) - + SQR(Atm33))); + 0.666666666666666666666666666667*SQR(trKL) + SQR(Atm11) + SQR(Atm22) + + SQR(Atm33))); - CCTK_REAL CCTK_ATTRIBUTE_UNUSED Jadm1L = 0.0625*(2*Atm23 - 2*Atm32 + - zL*(At11L*dgtu112 + At22L*dgtu222 + 2*(At12L*dgtu212 + - At13L*dgtu312 + At23L*dgtu322) + At33L*dgtu332 - + CCTK_REAL Jadm1L CCTK_ATTRIBUTE_UNUSED = 0.0625*(2*Atm23 - 2*Atm32 + + zL*(At11L*dgtu112 + At22L*dgtu222 + 2*(At12L*dgtu212 + At13L*dgtu312 + + At23L*dgtu322) + At33L*dgtu332 - 1.33333333333333333333333333333*JacPDstandardNth2trK - 16*Pi*S2) + yL*(-(At11L*dgtu113) - 2*At12L*dgtu213 - At22L*dgtu223 - 2*At13L*dgtu313 - 2*At23L*dgtu323 - At33L*dgtu333 + 1.33333333333333333333333333333*JacPDstandardNth3trK + 16*Pi*S3))*INV(Pi)*pow(ephi,6); - CCTK_REAL CCTK_ATTRIBUTE_UNUSED Jadm2L = 0.0625*(-2*Atm13 + 2*Atm31 + + CCTK_REAL Jadm2L CCTK_ATTRIBUTE_UNUSED = 0.0625*(-2*Atm13 + 2*Atm31 + zL*(-(At11L*dgtu111) - 2*At12L*dgtu211 - At22L*dgtu221 - 2*At13L*dgtu311 - 2*At23L*dgtu321 - At33L*dgtu331 + 1.33333333333333333333333333333*JacPDstandardNth1trK + 16*Pi*S1) + - xL*(At11L*dgtu113 + At22L*dgtu223 + 2*(At12L*dgtu213 + - At13L*dgtu313 + At23L*dgtu323) + At33L*dgtu333 - + xL*(At11L*dgtu113 + At22L*dgtu223 + 2*(At12L*dgtu213 + At13L*dgtu313 + + At23L*dgtu323) + At33L*dgtu333 - 1.33333333333333333333333333333*JacPDstandardNth3trK - 16*Pi*S3))*INV(Pi)*pow(ephi,6); - CCTK_REAL CCTK_ATTRIBUTE_UNUSED Jadm3L = 0.0625*(2*Atm12 - 2*Atm21 + - yL*(At11L*dgtu111 + At22L*dgtu221 + 2*(At12L*dgtu211 + - At13L*dgtu311 + At23L*dgtu321) + At33L*dgtu331 - + CCTK_REAL Jadm3L CCTK_ATTRIBUTE_UNUSED = 0.0625*(2*Atm12 - 2*Atm21 + + yL*(At11L*dgtu111 + At22L*dgtu221 + 2*(At12L*dgtu211 + At13L*dgtu311 + + At23L*dgtu321) + At33L*dgtu331 - 1.33333333333333333333333333333*JacPDstandardNth1trK - 16*Pi*S1) + xL*(-(At11L*dgtu112) - 2*At12L*dgtu212 - At22L*dgtu222 - 2*At13L*dgtu312 - 2*At23L*dgtu322 - At33L*dgtu332 + -- cgit v1.2.3