diff options
Diffstat (limited to 'ML_ADMQuantities/src/ML_ADMQuantities.cc')
-rw-r--r-- | ML_ADMQuantities/src/ML_ADMQuantities.cc | 618 |
1 files changed, 313 insertions, 305 deletions
diff --git a/ML_ADMQuantities/src/ML_ADMQuantities.cc b/ML_ADMQuantities/src/ML_ADMQuantities.cc index 4879023..2479b7a 100644 --- a/ML_ADMQuantities/src/ML_ADMQuantities.cc +++ b/ML_ADMQuantities/src/ML_ADMQuantities.cc @@ -46,45 +46,45 @@ static void ML_ADMQuantities_Body(cGH const * restrict const cctkGH, int const d /* Include user-supplied include files */ /* Initialise finite differencing variables */ - ptrdiff_t const di CCTK_ATTRIBUTE_UNUSED = 1; - ptrdiff_t const dj CCTK_ATTRIBUTE_UNUSED = CCTK_GFINDEX3D(cctkGH,0,1,0) - CCTK_GFINDEX3D(cctkGH,0,0,0); - ptrdiff_t const dk CCTK_ATTRIBUTE_UNUSED = CCTK_GFINDEX3D(cctkGH,0,0,1) - CCTK_GFINDEX3D(cctkGH,0,0,0); - ptrdiff_t const cdi CCTK_ATTRIBUTE_UNUSED = sizeof(CCTK_REAL) * di; - ptrdiff_t const cdj CCTK_ATTRIBUTE_UNUSED = sizeof(CCTK_REAL) * dj; - ptrdiff_t const cdk CCTK_ATTRIBUTE_UNUSED = sizeof(CCTK_REAL) * dk; - CCTK_REAL const dx CCTK_ATTRIBUTE_UNUSED = ToReal(CCTK_DELTA_SPACE(0)); - CCTK_REAL const dy CCTK_ATTRIBUTE_UNUSED = ToReal(CCTK_DELTA_SPACE(1)); - CCTK_REAL const dz CCTK_ATTRIBUTE_UNUSED = ToReal(CCTK_DELTA_SPACE(2)); - CCTK_REAL const dt CCTK_ATTRIBUTE_UNUSED = ToReal(CCTK_DELTA_TIME); - CCTK_REAL const t CCTK_ATTRIBUTE_UNUSED = ToReal(cctk_time); - CCTK_REAL const dxi CCTK_ATTRIBUTE_UNUSED = INV(dx); - CCTK_REAL const dyi CCTK_ATTRIBUTE_UNUSED = INV(dy); - CCTK_REAL const dzi CCTK_ATTRIBUTE_UNUSED = INV(dz); - CCTK_REAL const khalf CCTK_ATTRIBUTE_UNUSED = 0.5; - CCTK_REAL const kthird CCTK_ATTRIBUTE_UNUSED = 1/3.0; - CCTK_REAL const ktwothird CCTK_ATTRIBUTE_UNUSED = 2.0/3.0; - CCTK_REAL const kfourthird CCTK_ATTRIBUTE_UNUSED = 4.0/3.0; - CCTK_REAL const keightthird CCTK_ATTRIBUTE_UNUSED = 8.0/3.0; - CCTK_REAL const hdxi CCTK_ATTRIBUTE_UNUSED = 0.5 * dxi; - CCTK_REAL const hdyi CCTK_ATTRIBUTE_UNUSED = 0.5 * dyi; - CCTK_REAL const hdzi CCTK_ATTRIBUTE_UNUSED = 0.5 * dzi; + ptrdiff_t /*const*/ di CCTK_ATTRIBUTE_UNUSED = 1; + ptrdiff_t /*const*/ dj CCTK_ATTRIBUTE_UNUSED = CCTK_GFINDEX3D(cctkGH,0,1,0) - CCTK_GFINDEX3D(cctkGH,0,0,0); + ptrdiff_t /*const*/ dk CCTK_ATTRIBUTE_UNUSED = CCTK_GFINDEX3D(cctkGH,0,0,1) - CCTK_GFINDEX3D(cctkGH,0,0,0); + ptrdiff_t /*const*/ cdi CCTK_ATTRIBUTE_UNUSED = sizeof(CCTK_REAL) * di; + ptrdiff_t /*const*/ cdj CCTK_ATTRIBUTE_UNUSED = sizeof(CCTK_REAL) * dj; + ptrdiff_t /*const*/ cdk CCTK_ATTRIBUTE_UNUSED = sizeof(CCTK_REAL) * dk; + CCTK_REAL /*const*/ dx CCTK_ATTRIBUTE_UNUSED = ToReal(CCTK_DELTA_SPACE(0)); + CCTK_REAL /*const*/ dy CCTK_ATTRIBUTE_UNUSED = ToReal(CCTK_DELTA_SPACE(1)); + CCTK_REAL /*const*/ dz CCTK_ATTRIBUTE_UNUSED = ToReal(CCTK_DELTA_SPACE(2)); + CCTK_REAL /*const*/ dt CCTK_ATTRIBUTE_UNUSED = ToReal(CCTK_DELTA_TIME); + CCTK_REAL /*const*/ t CCTK_ATTRIBUTE_UNUSED = ToReal(cctk_time); + CCTK_REAL /*const*/ dxi CCTK_ATTRIBUTE_UNUSED = INV(dx); + CCTK_REAL /*const*/ dyi CCTK_ATTRIBUTE_UNUSED = INV(dy); + CCTK_REAL /*const*/ dzi CCTK_ATTRIBUTE_UNUSED = INV(dz); + CCTK_REAL /*const*/ khalf CCTK_ATTRIBUTE_UNUSED = 0.5; + CCTK_REAL /*const*/ kthird CCTK_ATTRIBUTE_UNUSED = 1/3.0; + CCTK_REAL /*const*/ ktwothird CCTK_ATTRIBUTE_UNUSED = 2.0/3.0; + CCTK_REAL /*const*/ kfourthird CCTK_ATTRIBUTE_UNUSED = 4.0/3.0; + CCTK_REAL /*const*/ keightthird CCTK_ATTRIBUTE_UNUSED = 8.0/3.0; + CCTK_REAL /*const*/ hdxi CCTK_ATTRIBUTE_UNUSED = 0.5 * dxi; + CCTK_REAL /*const*/ hdyi CCTK_ATTRIBUTE_UNUSED = 0.5 * dyi; + CCTK_REAL /*const*/ hdzi CCTK_ATTRIBUTE_UNUSED = 0.5 * dzi; /* Initialize predefined quantities */ - CCTK_REAL const p1o12dx CCTK_ATTRIBUTE_UNUSED = 0.0833333333333333333333333333333*INV(dx); - CCTK_REAL const p1o12dy CCTK_ATTRIBUTE_UNUSED = 0.0833333333333333333333333333333*INV(dy); - CCTK_REAL const p1o12dz CCTK_ATTRIBUTE_UNUSED = 0.0833333333333333333333333333333*INV(dz); - CCTK_REAL const p1o144dxdy CCTK_ATTRIBUTE_UNUSED = 0.00694444444444444444444444444444*INV(dx*dy); - CCTK_REAL const p1o144dxdz CCTK_ATTRIBUTE_UNUSED = 0.00694444444444444444444444444444*INV(dx*dz); - CCTK_REAL const p1o144dydz CCTK_ATTRIBUTE_UNUSED = 0.00694444444444444444444444444444*INV(dy*dz); - CCTK_REAL const p1odx CCTK_ATTRIBUTE_UNUSED = INV(dx); - CCTK_REAL const p1ody CCTK_ATTRIBUTE_UNUSED = INV(dy); - CCTK_REAL const p1odz CCTK_ATTRIBUTE_UNUSED = INV(dz); - CCTK_REAL const pm1o12dx2 CCTK_ATTRIBUTE_UNUSED = -0.0833333333333333333333333333333*INV(SQR(dx)); - CCTK_REAL const pm1o12dy2 CCTK_ATTRIBUTE_UNUSED = -0.0833333333333333333333333333333*INV(SQR(dy)); - CCTK_REAL const pm1o12dz2 CCTK_ATTRIBUTE_UNUSED = -0.0833333333333333333333333333333*INV(SQR(dz)); + CCTK_REAL /*const*/ p1o12dx CCTK_ATTRIBUTE_UNUSED = 0.0833333333333333333333333333333*INV(dx); + CCTK_REAL /*const*/ p1o12dy CCTK_ATTRIBUTE_UNUSED = 0.0833333333333333333333333333333*INV(dy); + CCTK_REAL /*const*/ p1o12dz CCTK_ATTRIBUTE_UNUSED = 0.0833333333333333333333333333333*INV(dz); + CCTK_REAL /*const*/ p1o144dxdy CCTK_ATTRIBUTE_UNUSED = 0.00694444444444444444444444444444*INV(dx*dy); + CCTK_REAL /*const*/ p1o144dxdz CCTK_ATTRIBUTE_UNUSED = 0.00694444444444444444444444444444*INV(dx*dz); + CCTK_REAL /*const*/ p1o144dydz CCTK_ATTRIBUTE_UNUSED = 0.00694444444444444444444444444444*INV(dy*dz); + CCTK_REAL /*const*/ p1odx CCTK_ATTRIBUTE_UNUSED = INV(dx); + CCTK_REAL /*const*/ p1ody CCTK_ATTRIBUTE_UNUSED = INV(dy); + CCTK_REAL /*const*/ p1odz CCTK_ATTRIBUTE_UNUSED = INV(dz); + CCTK_REAL /*const*/ pm1o12dx2 CCTK_ATTRIBUTE_UNUSED = -0.0833333333333333333333333333333*INV(SQR(dx)); + CCTK_REAL /*const*/ pm1o12dy2 CCTK_ATTRIBUTE_UNUSED = -0.0833333333333333333333333333333*INV(SQR(dy)); + CCTK_REAL /*const*/ pm1o12dz2 CCTK_ATTRIBUTE_UNUSED = -0.0833333333333333333333333333333*INV(SQR(dz)); /* Jacobian variable pointers */ - bool const use_jacobian = (!CCTK_IsFunctionAliased("MultiPatch_GetMap") || MultiPatch_GetMap(cctkGH) != jacobian_identity_map) + bool const use_jacobian CCTK_ATTRIBUTE_UNUSED = (!CCTK_IsFunctionAliased("MultiPatch_GetMap") || MultiPatch_GetMap(cctkGH) != jacobian_identity_map) && strlen(jacobian_group) > 0; if (use_jacobian && strlen(jacobian_derivative_group) == 0) { @@ -142,7 +142,7 @@ static void ML_ADMQuantities_Body(cGH const * restrict const cctkGH, int const d i,j,k, imin[0],imin[1],imin[2], imax[0],imax[1],imax[2], cctk_ash[0],cctk_ash[1],cctk_ash[2]) { - ptrdiff_t const index CCTK_ATTRIBUTE_UNUSED = di*i + dj*j + dk*k; + ptrdiff_t /*const*/ index CCTK_ATTRIBUTE_UNUSED = di*i + dj*j + dk*k; /* Assign local copies of grid functions */ @@ -236,72 +236,72 @@ static void ML_ADMQuantities_Body(cGH const * restrict const cctkGH, int const d /* Include user supplied include files */ /* Precompute derivatives */ - CCTK_REAL const PDstandardNth1gt11 CCTK_ATTRIBUTE_UNUSED = PDstandardNth1(>11[index]); - CCTK_REAL const PDstandardNth2gt11 CCTK_ATTRIBUTE_UNUSED = PDstandardNth2(>11[index]); - CCTK_REAL const PDstandardNth3gt11 CCTK_ATTRIBUTE_UNUSED = PDstandardNth3(>11[index]); - CCTK_REAL const PDstandardNth11gt11 CCTK_ATTRIBUTE_UNUSED = PDstandardNth11(>11[index]); - CCTK_REAL const PDstandardNth22gt11 CCTK_ATTRIBUTE_UNUSED = PDstandardNth22(>11[index]); - CCTK_REAL const PDstandardNth33gt11 CCTK_ATTRIBUTE_UNUSED = PDstandardNth33(>11[index]); - CCTK_REAL const PDstandardNth12gt11 CCTK_ATTRIBUTE_UNUSED = PDstandardNth12(>11[index]); - CCTK_REAL const PDstandardNth13gt11 CCTK_ATTRIBUTE_UNUSED = PDstandardNth13(>11[index]); - CCTK_REAL const PDstandardNth23gt11 CCTK_ATTRIBUTE_UNUSED = PDstandardNth23(>11[index]); - CCTK_REAL const PDstandardNth1gt12 CCTK_ATTRIBUTE_UNUSED = PDstandardNth1(>12[index]); - CCTK_REAL const PDstandardNth2gt12 CCTK_ATTRIBUTE_UNUSED = PDstandardNth2(>12[index]); - CCTK_REAL const PDstandardNth3gt12 CCTK_ATTRIBUTE_UNUSED = PDstandardNth3(>12[index]); - CCTK_REAL const PDstandardNth11gt12 CCTK_ATTRIBUTE_UNUSED = PDstandardNth11(>12[index]); - CCTK_REAL const PDstandardNth22gt12 CCTK_ATTRIBUTE_UNUSED = PDstandardNth22(>12[index]); - CCTK_REAL const PDstandardNth33gt12 CCTK_ATTRIBUTE_UNUSED = PDstandardNth33(>12[index]); - CCTK_REAL const PDstandardNth12gt12 CCTK_ATTRIBUTE_UNUSED = PDstandardNth12(>12[index]); - CCTK_REAL const PDstandardNth13gt12 CCTK_ATTRIBUTE_UNUSED = PDstandardNth13(>12[index]); - CCTK_REAL const PDstandardNth23gt12 CCTK_ATTRIBUTE_UNUSED = PDstandardNth23(>12[index]); - CCTK_REAL const PDstandardNth1gt13 CCTK_ATTRIBUTE_UNUSED = PDstandardNth1(>13[index]); - CCTK_REAL const PDstandardNth2gt13 CCTK_ATTRIBUTE_UNUSED = PDstandardNth2(>13[index]); - CCTK_REAL const PDstandardNth3gt13 CCTK_ATTRIBUTE_UNUSED = PDstandardNth3(>13[index]); - CCTK_REAL const PDstandardNth11gt13 CCTK_ATTRIBUTE_UNUSED = PDstandardNth11(>13[index]); - CCTK_REAL const PDstandardNth22gt13 CCTK_ATTRIBUTE_UNUSED = PDstandardNth22(>13[index]); - CCTK_REAL const PDstandardNth33gt13 CCTK_ATTRIBUTE_UNUSED = PDstandardNth33(>13[index]); - CCTK_REAL const PDstandardNth12gt13 CCTK_ATTRIBUTE_UNUSED = PDstandardNth12(>13[index]); - CCTK_REAL const PDstandardNth13gt13 CCTK_ATTRIBUTE_UNUSED = PDstandardNth13(>13[index]); - CCTK_REAL const PDstandardNth23gt13 CCTK_ATTRIBUTE_UNUSED = PDstandardNth23(>13[index]); - CCTK_REAL const PDstandardNth1gt22 CCTK_ATTRIBUTE_UNUSED = PDstandardNth1(>22[index]); - CCTK_REAL const PDstandardNth2gt22 CCTK_ATTRIBUTE_UNUSED = PDstandardNth2(>22[index]); - CCTK_REAL const PDstandardNth3gt22 CCTK_ATTRIBUTE_UNUSED = PDstandardNth3(>22[index]); - CCTK_REAL const PDstandardNth11gt22 CCTK_ATTRIBUTE_UNUSED = PDstandardNth11(>22[index]); - CCTK_REAL const PDstandardNth22gt22 CCTK_ATTRIBUTE_UNUSED = PDstandardNth22(>22[index]); - CCTK_REAL const PDstandardNth33gt22 CCTK_ATTRIBUTE_UNUSED = PDstandardNth33(>22[index]); - CCTK_REAL const PDstandardNth12gt22 CCTK_ATTRIBUTE_UNUSED = PDstandardNth12(>22[index]); - CCTK_REAL const PDstandardNth13gt22 CCTK_ATTRIBUTE_UNUSED = PDstandardNth13(>22[index]); - CCTK_REAL const PDstandardNth23gt22 CCTK_ATTRIBUTE_UNUSED = PDstandardNth23(>22[index]); - CCTK_REAL const PDstandardNth1gt23 CCTK_ATTRIBUTE_UNUSED = PDstandardNth1(>23[index]); - CCTK_REAL const PDstandardNth2gt23 CCTK_ATTRIBUTE_UNUSED = PDstandardNth2(>23[index]); - CCTK_REAL const PDstandardNth3gt23 CCTK_ATTRIBUTE_UNUSED = PDstandardNth3(>23[index]); - CCTK_REAL const PDstandardNth11gt23 CCTK_ATTRIBUTE_UNUSED = PDstandardNth11(>23[index]); - CCTK_REAL const PDstandardNth22gt23 CCTK_ATTRIBUTE_UNUSED = PDstandardNth22(>23[index]); - CCTK_REAL const PDstandardNth33gt23 CCTK_ATTRIBUTE_UNUSED = PDstandardNth33(>23[index]); - CCTK_REAL const PDstandardNth12gt23 CCTK_ATTRIBUTE_UNUSED = PDstandardNth12(>23[index]); - CCTK_REAL const PDstandardNth13gt23 CCTK_ATTRIBUTE_UNUSED = PDstandardNth13(>23[index]); - CCTK_REAL const PDstandardNth23gt23 CCTK_ATTRIBUTE_UNUSED = PDstandardNth23(>23[index]); - CCTK_REAL const PDstandardNth1gt33 CCTK_ATTRIBUTE_UNUSED = PDstandardNth1(>33[index]); - CCTK_REAL const PDstandardNth2gt33 CCTK_ATTRIBUTE_UNUSED = PDstandardNth2(>33[index]); - CCTK_REAL const PDstandardNth3gt33 CCTK_ATTRIBUTE_UNUSED = PDstandardNth3(>33[index]); - CCTK_REAL const PDstandardNth11gt33 CCTK_ATTRIBUTE_UNUSED = PDstandardNth11(>33[index]); - CCTK_REAL const PDstandardNth22gt33 CCTK_ATTRIBUTE_UNUSED = PDstandardNth22(>33[index]); - CCTK_REAL const PDstandardNth33gt33 CCTK_ATTRIBUTE_UNUSED = PDstandardNth33(>33[index]); - CCTK_REAL const PDstandardNth12gt33 CCTK_ATTRIBUTE_UNUSED = PDstandardNth12(>33[index]); - CCTK_REAL const PDstandardNth13gt33 CCTK_ATTRIBUTE_UNUSED = PDstandardNth13(>33[index]); - CCTK_REAL const PDstandardNth23gt33 CCTK_ATTRIBUTE_UNUSED = PDstandardNth23(>33[index]); - CCTK_REAL const PDstandardNth1trK CCTK_ATTRIBUTE_UNUSED = PDstandardNth1(&trK[index]); - CCTK_REAL const PDstandardNth2trK CCTK_ATTRIBUTE_UNUSED = PDstandardNth2(&trK[index]); - CCTK_REAL const PDstandardNth3trK CCTK_ATTRIBUTE_UNUSED = PDstandardNth3(&trK[index]); - CCTK_REAL const PDstandardNth1Xt1 CCTK_ATTRIBUTE_UNUSED = PDstandardNth1(&Xt1[index]); - CCTK_REAL const PDstandardNth2Xt1 CCTK_ATTRIBUTE_UNUSED = PDstandardNth2(&Xt1[index]); - CCTK_REAL const PDstandardNth3Xt1 CCTK_ATTRIBUTE_UNUSED = PDstandardNth3(&Xt1[index]); - CCTK_REAL const PDstandardNth1Xt2 CCTK_ATTRIBUTE_UNUSED = PDstandardNth1(&Xt2[index]); - CCTK_REAL const PDstandardNth2Xt2 CCTK_ATTRIBUTE_UNUSED = PDstandardNth2(&Xt2[index]); - CCTK_REAL const PDstandardNth3Xt2 CCTK_ATTRIBUTE_UNUSED = PDstandardNth3(&Xt2[index]); - CCTK_REAL const PDstandardNth1Xt3 CCTK_ATTRIBUTE_UNUSED = PDstandardNth1(&Xt3[index]); - CCTK_REAL const PDstandardNth2Xt3 CCTK_ATTRIBUTE_UNUSED = PDstandardNth2(&Xt3[index]); - CCTK_REAL const PDstandardNth3Xt3 CCTK_ATTRIBUTE_UNUSED = PDstandardNth3(&Xt3[index]); + CCTK_REAL /*const*/ PDstandardNth1gt11 CCTK_ATTRIBUTE_UNUSED = PDstandardNth1(>11[index]); + CCTK_REAL /*const*/ PDstandardNth2gt11 CCTK_ATTRIBUTE_UNUSED = PDstandardNth2(>11[index]); + CCTK_REAL /*const*/ PDstandardNth3gt11 CCTK_ATTRIBUTE_UNUSED = PDstandardNth3(>11[index]); + CCTK_REAL /*const*/ PDstandardNth11gt11 CCTK_ATTRIBUTE_UNUSED = PDstandardNth11(>11[index]); + CCTK_REAL /*const*/ PDstandardNth22gt11 CCTK_ATTRIBUTE_UNUSED = PDstandardNth22(>11[index]); + CCTK_REAL /*const*/ PDstandardNth33gt11 CCTK_ATTRIBUTE_UNUSED = PDstandardNth33(>11[index]); + CCTK_REAL /*const*/ PDstandardNth12gt11 CCTK_ATTRIBUTE_UNUSED = PDstandardNth12(>11[index]); + CCTK_REAL /*const*/ PDstandardNth13gt11 CCTK_ATTRIBUTE_UNUSED = PDstandardNth13(>11[index]); + CCTK_REAL /*const*/ PDstandardNth23gt11 CCTK_ATTRIBUTE_UNUSED = PDstandardNth23(>11[index]); + CCTK_REAL /*const*/ PDstandardNth1gt12 CCTK_ATTRIBUTE_UNUSED = PDstandardNth1(>12[index]); + CCTK_REAL /*const*/ PDstandardNth2gt12 CCTK_ATTRIBUTE_UNUSED = PDstandardNth2(>12[index]); + CCTK_REAL /*const*/ PDstandardNth3gt12 CCTK_ATTRIBUTE_UNUSED = PDstandardNth3(>12[index]); + CCTK_REAL /*const*/ PDstandardNth11gt12 CCTK_ATTRIBUTE_UNUSED = PDstandardNth11(>12[index]); + CCTK_REAL /*const*/ PDstandardNth22gt12 CCTK_ATTRIBUTE_UNUSED = PDstandardNth22(>12[index]); + CCTK_REAL /*const*/ PDstandardNth33gt12 CCTK_ATTRIBUTE_UNUSED = PDstandardNth33(>12[index]); + CCTK_REAL /*const*/ PDstandardNth12gt12 CCTK_ATTRIBUTE_UNUSED = PDstandardNth12(>12[index]); + CCTK_REAL /*const*/ PDstandardNth13gt12 CCTK_ATTRIBUTE_UNUSED = PDstandardNth13(>12[index]); + CCTK_REAL /*const*/ PDstandardNth23gt12 CCTK_ATTRIBUTE_UNUSED = PDstandardNth23(>12[index]); + CCTK_REAL /*const*/ PDstandardNth1gt13 CCTK_ATTRIBUTE_UNUSED = PDstandardNth1(>13[index]); + CCTK_REAL /*const*/ PDstandardNth2gt13 CCTK_ATTRIBUTE_UNUSED = PDstandardNth2(>13[index]); + CCTK_REAL /*const*/ PDstandardNth3gt13 CCTK_ATTRIBUTE_UNUSED = PDstandardNth3(>13[index]); + CCTK_REAL /*const*/ PDstandardNth11gt13 CCTK_ATTRIBUTE_UNUSED = PDstandardNth11(>13[index]); + CCTK_REAL /*const*/ PDstandardNth22gt13 CCTK_ATTRIBUTE_UNUSED = PDstandardNth22(>13[index]); + CCTK_REAL /*const*/ PDstandardNth33gt13 CCTK_ATTRIBUTE_UNUSED = PDstandardNth33(>13[index]); + CCTK_REAL /*const*/ PDstandardNth12gt13 CCTK_ATTRIBUTE_UNUSED = PDstandardNth12(>13[index]); + CCTK_REAL /*const*/ PDstandardNth13gt13 CCTK_ATTRIBUTE_UNUSED = PDstandardNth13(>13[index]); + CCTK_REAL /*const*/ PDstandardNth23gt13 CCTK_ATTRIBUTE_UNUSED = PDstandardNth23(>13[index]); + CCTK_REAL /*const*/ PDstandardNth1gt22 CCTK_ATTRIBUTE_UNUSED = PDstandardNth1(>22[index]); + CCTK_REAL /*const*/ PDstandardNth2gt22 CCTK_ATTRIBUTE_UNUSED = PDstandardNth2(>22[index]); + CCTK_REAL /*const*/ PDstandardNth3gt22 CCTK_ATTRIBUTE_UNUSED = PDstandardNth3(>22[index]); + CCTK_REAL /*const*/ PDstandardNth11gt22 CCTK_ATTRIBUTE_UNUSED = PDstandardNth11(>22[index]); + CCTK_REAL /*const*/ PDstandardNth22gt22 CCTK_ATTRIBUTE_UNUSED = PDstandardNth22(>22[index]); + CCTK_REAL /*const*/ PDstandardNth33gt22 CCTK_ATTRIBUTE_UNUSED = PDstandardNth33(>22[index]); + CCTK_REAL /*const*/ PDstandardNth12gt22 CCTK_ATTRIBUTE_UNUSED = PDstandardNth12(>22[index]); + CCTK_REAL /*const*/ PDstandardNth13gt22 CCTK_ATTRIBUTE_UNUSED = PDstandardNth13(>22[index]); + CCTK_REAL /*const*/ PDstandardNth23gt22 CCTK_ATTRIBUTE_UNUSED = PDstandardNth23(>22[index]); + CCTK_REAL /*const*/ PDstandardNth1gt23 CCTK_ATTRIBUTE_UNUSED = PDstandardNth1(>23[index]); + CCTK_REAL /*const*/ PDstandardNth2gt23 CCTK_ATTRIBUTE_UNUSED = PDstandardNth2(>23[index]); + CCTK_REAL /*const*/ PDstandardNth3gt23 CCTK_ATTRIBUTE_UNUSED = PDstandardNth3(>23[index]); + CCTK_REAL /*const*/ PDstandardNth11gt23 CCTK_ATTRIBUTE_UNUSED = PDstandardNth11(>23[index]); + CCTK_REAL /*const*/ PDstandardNth22gt23 CCTK_ATTRIBUTE_UNUSED = PDstandardNth22(>23[index]); + CCTK_REAL /*const*/ PDstandardNth33gt23 CCTK_ATTRIBUTE_UNUSED = PDstandardNth33(>23[index]); + CCTK_REAL /*const*/ PDstandardNth12gt23 CCTK_ATTRIBUTE_UNUSED = PDstandardNth12(>23[index]); + CCTK_REAL /*const*/ PDstandardNth13gt23 CCTK_ATTRIBUTE_UNUSED = PDstandardNth13(>23[index]); + CCTK_REAL /*const*/ PDstandardNth23gt23 CCTK_ATTRIBUTE_UNUSED = PDstandardNth23(>23[index]); + CCTK_REAL /*const*/ PDstandardNth1gt33 CCTK_ATTRIBUTE_UNUSED = PDstandardNth1(>33[index]); + CCTK_REAL /*const*/ PDstandardNth2gt33 CCTK_ATTRIBUTE_UNUSED = PDstandardNth2(>33[index]); + CCTK_REAL /*const*/ PDstandardNth3gt33 CCTK_ATTRIBUTE_UNUSED = PDstandardNth3(>33[index]); + CCTK_REAL /*const*/ PDstandardNth11gt33 CCTK_ATTRIBUTE_UNUSED = PDstandardNth11(>33[index]); + CCTK_REAL /*const*/ PDstandardNth22gt33 CCTK_ATTRIBUTE_UNUSED = PDstandardNth22(>33[index]); + CCTK_REAL /*const*/ PDstandardNth33gt33 CCTK_ATTRIBUTE_UNUSED = PDstandardNth33(>33[index]); + CCTK_REAL /*const*/ PDstandardNth12gt33 CCTK_ATTRIBUTE_UNUSED = PDstandardNth12(>33[index]); + CCTK_REAL /*const*/ PDstandardNth13gt33 CCTK_ATTRIBUTE_UNUSED = PDstandardNth13(>33[index]); + CCTK_REAL /*const*/ PDstandardNth23gt33 CCTK_ATTRIBUTE_UNUSED = PDstandardNth23(>33[index]); + CCTK_REAL /*const*/ PDstandardNth1trK CCTK_ATTRIBUTE_UNUSED = PDstandardNth1(&trK[index]); + CCTK_REAL /*const*/ PDstandardNth2trK CCTK_ATTRIBUTE_UNUSED = PDstandardNth2(&trK[index]); + CCTK_REAL /*const*/ PDstandardNth3trK CCTK_ATTRIBUTE_UNUSED = PDstandardNth3(&trK[index]); + CCTK_REAL /*const*/ PDstandardNth1Xt1 CCTK_ATTRIBUTE_UNUSED = PDstandardNth1(&Xt1[index]); + CCTK_REAL /*const*/ PDstandardNth2Xt1 CCTK_ATTRIBUTE_UNUSED = PDstandardNth2(&Xt1[index]); + CCTK_REAL /*const*/ PDstandardNth3Xt1 CCTK_ATTRIBUTE_UNUSED = PDstandardNth3(&Xt1[index]); + CCTK_REAL /*const*/ PDstandardNth1Xt2 CCTK_ATTRIBUTE_UNUSED = PDstandardNth1(&Xt2[index]); + CCTK_REAL /*const*/ PDstandardNth2Xt2 CCTK_ATTRIBUTE_UNUSED = PDstandardNth2(&Xt2[index]); + CCTK_REAL /*const*/ PDstandardNth3Xt2 CCTK_ATTRIBUTE_UNUSED = PDstandardNth3(&Xt2[index]); + CCTK_REAL /*const*/ PDstandardNth1Xt3 CCTK_ATTRIBUTE_UNUSED = PDstandardNth1(&Xt3[index]); + CCTK_REAL /*const*/ PDstandardNth2Xt3 CCTK_ATTRIBUTE_UNUSED = PDstandardNth2(&Xt3[index]); + CCTK_REAL /*const*/ PDstandardNth3Xt3 CCTK_ATTRIBUTE_UNUSED = PDstandardNth3(&Xt3[index]); /* Calculate temporaries and grid functions */ CCTK_REAL JacPDstandardNth11gt11 CCTK_ATTRIBUTE_UNUSED ; @@ -482,109 +482,109 @@ static void ML_ADMQuantities_Body(cGH const * restrict const cctkGH, int const d J23L*PDstandardNth2Xt3 + J33L*PDstandardNth3Xt3; JacPDstandardNth11gt11 = dJ111L*PDstandardNth1gt11 + - 2*(J11L*(J21L*PDstandardNth12gt11 + J31L*PDstandardNth13gt11) + + 2.*(J11L*(J21L*PDstandardNth12gt11 + J31L*PDstandardNth13gt11) + J21L*J31L*PDstandardNth23gt11) + dJ211L*PDstandardNth2gt11 + dJ311L*PDstandardNth3gt11 + PDstandardNth11gt11*SQR(J11L) + PDstandardNth22gt11*SQR(J21L) + PDstandardNth33gt11*SQR(J31L); JacPDstandardNth11gt12 = dJ111L*PDstandardNth1gt12 + - 2*(J11L*(J21L*PDstandardNth12gt12 + J31L*PDstandardNth13gt12) + + 2.*(J11L*(J21L*PDstandardNth12gt12 + J31L*PDstandardNth13gt12) + J21L*J31L*PDstandardNth23gt12) + dJ211L*PDstandardNth2gt12 + dJ311L*PDstandardNth3gt12 + PDstandardNth11gt12*SQR(J11L) + PDstandardNth22gt12*SQR(J21L) + PDstandardNth33gt12*SQR(J31L); JacPDstandardNth11gt13 = dJ111L*PDstandardNth1gt13 + - 2*(J11L*(J21L*PDstandardNth12gt13 + J31L*PDstandardNth13gt13) + + 2.*(J11L*(J21L*PDstandardNth12gt13 + J31L*PDstandardNth13gt13) + J21L*J31L*PDstandardNth23gt13) + dJ211L*PDstandardNth2gt13 + dJ311L*PDstandardNth3gt13 + PDstandardNth11gt13*SQR(J11L) + PDstandardNth22gt13*SQR(J21L) + PDstandardNth33gt13*SQR(J31L); JacPDstandardNth11gt22 = dJ111L*PDstandardNth1gt22 + - 2*(J11L*(J21L*PDstandardNth12gt22 + J31L*PDstandardNth13gt22) + + 2.*(J11L*(J21L*PDstandardNth12gt22 + J31L*PDstandardNth13gt22) + J21L*J31L*PDstandardNth23gt22) + dJ211L*PDstandardNth2gt22 + dJ311L*PDstandardNth3gt22 + PDstandardNth11gt22*SQR(J11L) + PDstandardNth22gt22*SQR(J21L) + PDstandardNth33gt22*SQR(J31L); JacPDstandardNth11gt23 = dJ111L*PDstandardNth1gt23 + - 2*(J11L*(J21L*PDstandardNth12gt23 + J31L*PDstandardNth13gt23) + + 2.*(J11L*(J21L*PDstandardNth12gt23 + J31L*PDstandardNth13gt23) + J21L*J31L*PDstandardNth23gt23) + dJ211L*PDstandardNth2gt23 + dJ311L*PDstandardNth3gt23 + PDstandardNth11gt23*SQR(J11L) + PDstandardNth22gt23*SQR(J21L) + PDstandardNth33gt23*SQR(J31L); JacPDstandardNth11gt33 = dJ111L*PDstandardNth1gt33 + - 2*(J11L*(J21L*PDstandardNth12gt33 + J31L*PDstandardNth13gt33) + + 2.*(J11L*(J21L*PDstandardNth12gt33 + J31L*PDstandardNth13gt33) + J21L*J31L*PDstandardNth23gt33) + dJ211L*PDstandardNth2gt33 + dJ311L*PDstandardNth3gt33 + PDstandardNth11gt33*SQR(J11L) + PDstandardNth22gt33*SQR(J21L) + PDstandardNth33gt33*SQR(J31L); JacPDstandardNth22gt11 = dJ122L*PDstandardNth1gt11 + - 2*(J12L*(J22L*PDstandardNth12gt11 + J32L*PDstandardNth13gt11) + + 2.*(J12L*(J22L*PDstandardNth12gt11 + J32L*PDstandardNth13gt11) + J22L*J32L*PDstandardNth23gt11) + dJ222L*PDstandardNth2gt11 + dJ322L*PDstandardNth3gt11 + PDstandardNth11gt11*SQR(J12L) + PDstandardNth22gt11*SQR(J22L) + PDstandardNth33gt11*SQR(J32L); JacPDstandardNth22gt12 = dJ122L*PDstandardNth1gt12 + - 2*(J12L*(J22L*PDstandardNth12gt12 + J32L*PDstandardNth13gt12) + + 2.*(J12L*(J22L*PDstandardNth12gt12 + J32L*PDstandardNth13gt12) + J22L*J32L*PDstandardNth23gt12) + dJ222L*PDstandardNth2gt12 + dJ322L*PDstandardNth3gt12 + PDstandardNth11gt12*SQR(J12L) + PDstandardNth22gt12*SQR(J22L) + PDstandardNth33gt12*SQR(J32L); JacPDstandardNth22gt13 = dJ122L*PDstandardNth1gt13 + - 2*(J12L*(J22L*PDstandardNth12gt13 + J32L*PDstandardNth13gt13) + + 2.*(J12L*(J22L*PDstandardNth12gt13 + J32L*PDstandardNth13gt13) + J22L*J32L*PDstandardNth23gt13) + dJ222L*PDstandardNth2gt13 + dJ322L*PDstandardNth3gt13 + PDstandardNth11gt13*SQR(J12L) + PDstandardNth22gt13*SQR(J22L) + PDstandardNth33gt13*SQR(J32L); JacPDstandardNth22gt22 = dJ122L*PDstandardNth1gt22 + - 2*(J12L*(J22L*PDstandardNth12gt22 + J32L*PDstandardNth13gt22) + + 2.*(J12L*(J22L*PDstandardNth12gt22 + J32L*PDstandardNth13gt22) + J22L*J32L*PDstandardNth23gt22) + dJ222L*PDstandardNth2gt22 + dJ322L*PDstandardNth3gt22 + PDstandardNth11gt22*SQR(J12L) + PDstandardNth22gt22*SQR(J22L) + PDstandardNth33gt22*SQR(J32L); JacPDstandardNth22gt23 = dJ122L*PDstandardNth1gt23 + - 2*(J12L*(J22L*PDstandardNth12gt23 + J32L*PDstandardNth13gt23) + + 2.*(J12L*(J22L*PDstandardNth12gt23 + J32L*PDstandardNth13gt23) + J22L*J32L*PDstandardNth23gt23) + dJ222L*PDstandardNth2gt23 + dJ322L*PDstandardNth3gt23 + PDstandardNth11gt23*SQR(J12L) + PDstandardNth22gt23*SQR(J22L) + PDstandardNth33gt23*SQR(J32L); JacPDstandardNth22gt33 = dJ122L*PDstandardNth1gt33 + - 2*(J12L*(J22L*PDstandardNth12gt33 + J32L*PDstandardNth13gt33) + + 2.*(J12L*(J22L*PDstandardNth12gt33 + J32L*PDstandardNth13gt33) + J22L*J32L*PDstandardNth23gt33) + dJ222L*PDstandardNth2gt33 + dJ322L*PDstandardNth3gt33 + PDstandardNth11gt33*SQR(J12L) + PDstandardNth22gt33*SQR(J22L) + PDstandardNth33gt33*SQR(J32L); JacPDstandardNth33gt11 = dJ133L*PDstandardNth1gt11 + - 2*(J13L*(J23L*PDstandardNth12gt11 + J33L*PDstandardNth13gt11) + + 2.*(J13L*(J23L*PDstandardNth12gt11 + J33L*PDstandardNth13gt11) + J23L*J33L*PDstandardNth23gt11) + dJ233L*PDstandardNth2gt11 + dJ333L*PDstandardNth3gt11 + PDstandardNth11gt11*SQR(J13L) + PDstandardNth22gt11*SQR(J23L) + PDstandardNth33gt11*SQR(J33L); JacPDstandardNth33gt12 = dJ133L*PDstandardNth1gt12 + - 2*(J13L*(J23L*PDstandardNth12gt12 + J33L*PDstandardNth13gt12) + + 2.*(J13L*(J23L*PDstandardNth12gt12 + J33L*PDstandardNth13gt12) + J23L*J33L*PDstandardNth23gt12) + dJ233L*PDstandardNth2gt12 + dJ333L*PDstandardNth3gt12 + PDstandardNth11gt12*SQR(J13L) + PDstandardNth22gt12*SQR(J23L) + PDstandardNth33gt12*SQR(J33L); JacPDstandardNth33gt13 = dJ133L*PDstandardNth1gt13 + - 2*(J13L*(J23L*PDstandardNth12gt13 + J33L*PDstandardNth13gt13) + + 2.*(J13L*(J23L*PDstandardNth12gt13 + J33L*PDstandardNth13gt13) + J23L*J33L*PDstandardNth23gt13) + dJ233L*PDstandardNth2gt13 + dJ333L*PDstandardNth3gt13 + PDstandardNth11gt13*SQR(J13L) + PDstandardNth22gt13*SQR(J23L) + PDstandardNth33gt13*SQR(J33L); JacPDstandardNth33gt22 = dJ133L*PDstandardNth1gt22 + - 2*(J13L*(J23L*PDstandardNth12gt22 + J33L*PDstandardNth13gt22) + + 2.*(J13L*(J23L*PDstandardNth12gt22 + J33L*PDstandardNth13gt22) + J23L*J33L*PDstandardNth23gt22) + dJ233L*PDstandardNth2gt22 + dJ333L*PDstandardNth3gt22 + PDstandardNth11gt22*SQR(J13L) + PDstandardNth22gt22*SQR(J23L) + PDstandardNth33gt22*SQR(J33L); JacPDstandardNth33gt23 = dJ133L*PDstandardNth1gt23 + - 2*(J13L*(J23L*PDstandardNth12gt23 + J33L*PDstandardNth13gt23) + + 2.*(J13L*(J23L*PDstandardNth12gt23 + J33L*PDstandardNth13gt23) + J23L*J33L*PDstandardNth23gt23) + dJ233L*PDstandardNth2gt23 + dJ333L*PDstandardNth3gt23 + PDstandardNth11gt23*SQR(J13L) + PDstandardNth22gt23*SQR(J23L) + PDstandardNth33gt23*SQR(J33L); JacPDstandardNth33gt33 = dJ133L*PDstandardNth1gt33 + - 2*(J13L*(J23L*PDstandardNth12gt33 + J33L*PDstandardNth13gt33) + + 2.*(J13L*(J23L*PDstandardNth12gt33 + J33L*PDstandardNth13gt33) + J23L*J33L*PDstandardNth23gt33) + dJ233L*PDstandardNth2gt33 + dJ333L*PDstandardNth3gt33 + PDstandardNth11gt33*SQR(J13L) + PDstandardNth22gt33*SQR(J23L) + PDstandardNth33gt33*SQR(J33L); @@ -1048,133 +1048,142 @@ static void ML_ADMQuantities_Body(cGH const * restrict const cctkGH, int const d JacPDstandardNth32gt33 = PDstandardNth23gt33; } - CCTK_REAL CCTK_ATTRIBUTE_UNUSED detgt = 1; + CCTK_REAL CCTK_ATTRIBUTE_UNUSED detgt = 1.; CCTK_REAL CCTK_ATTRIBUTE_UNUSED gtu11 = INV(detgt)*(gt22L*gt33L - - SQR(gt23L)); + 1.*SQR(gt23L)); CCTK_REAL CCTK_ATTRIBUTE_UNUSED gtu21 = (gt13L*gt23L - - gt12L*gt33L)*INV(detgt); + 1.*gt12L*gt33L)*INV(detgt); - CCTK_REAL CCTK_ATTRIBUTE_UNUSED gtu31 = (-(gt13L*gt22L) + + CCTK_REAL CCTK_ATTRIBUTE_UNUSED gtu31 = (-1.*gt13L*gt22L + gt12L*gt23L)*INV(detgt); CCTK_REAL CCTK_ATTRIBUTE_UNUSED gtu22 = INV(detgt)*(gt11L*gt33L - - SQR(gt13L)); + 1.*SQR(gt13L)); CCTK_REAL CCTK_ATTRIBUTE_UNUSED gtu32 = (gt12L*gt13L - - gt11L*gt23L)*INV(detgt); + 1.*gt11L*gt23L)*INV(detgt); CCTK_REAL CCTK_ATTRIBUTE_UNUSED gtu33 = INV(detgt)*(gt11L*gt22L - - SQR(gt12L)); + 1.*SQR(gt12L)); CCTK_REAL CCTK_ATTRIBUTE_UNUSED dgtu111 = - -2*(gtu11*(gtu21*JacPDstandardNth1gt12 + gtu31*JacPDstandardNth1gt13) + - gtu21*gtu31*JacPDstandardNth1gt23) - JacPDstandardNth1gt11*SQR(gtu11) - - JacPDstandardNth1gt22*SQR(gtu21) - JacPDstandardNth1gt33*SQR(gtu31); + -2.*(gtu11*(gtu21*JacPDstandardNth1gt12 + gtu31*JacPDstandardNth1gt13) + + gtu21*gtu31*JacPDstandardNth1gt23) - + 1.*(JacPDstandardNth1gt11*SQR(gtu11) + JacPDstandardNth1gt22*SQR(gtu21) + + JacPDstandardNth1gt33*SQR(gtu31)); CCTK_REAL CCTK_ATTRIBUTE_UNUSED dgtu211 = - -(gtu21*(gtu11*JacPDstandardNth1gt11 + gtu21*JacPDstandardNth1gt12 + - gtu31*JacPDstandardNth1gt13)) - gtu22*(gtu11*JacPDstandardNth1gt12 + - gtu21*JacPDstandardNth1gt22 + gtu31*JacPDstandardNth1gt23) - + -1.*(gtu21*(gtu11*JacPDstandardNth1gt11 + gtu21*JacPDstandardNth1gt12 + + gtu31*JacPDstandardNth1gt13) + gtu22*(gtu11*JacPDstandardNth1gt12 + + gtu21*JacPDstandardNth1gt22 + gtu31*JacPDstandardNth1gt23) + gtu32*(gtu11*JacPDstandardNth1gt13 + gtu21*JacPDstandardNth1gt23 + - gtu31*JacPDstandardNth1gt33); + gtu31*JacPDstandardNth1gt33)); CCTK_REAL CCTK_ATTRIBUTE_UNUSED dgtu311 = - -(gtu11*(gtu31*JacPDstandardNth1gt11 + gtu32*JacPDstandardNth1gt12 + - gtu33*JacPDstandardNth1gt13)) - gtu21*(gtu31*JacPDstandardNth1gt12 + - gtu32*JacPDstandardNth1gt22 + gtu33*JacPDstandardNth1gt23) - + -1.*(gtu11*(gtu31*JacPDstandardNth1gt11 + gtu32*JacPDstandardNth1gt12 + + gtu33*JacPDstandardNth1gt13) + gtu21*(gtu31*JacPDstandardNth1gt12 + + gtu32*JacPDstandardNth1gt22 + gtu33*JacPDstandardNth1gt23) + gtu31*(gtu31*JacPDstandardNth1gt13 + gtu32*JacPDstandardNth1gt23 + - gtu33*JacPDstandardNth1gt33); + gtu33*JacPDstandardNth1gt33)); CCTK_REAL CCTK_ATTRIBUTE_UNUSED dgtu221 = - -2*(gtu21*(gtu22*JacPDstandardNth1gt12 + gtu32*JacPDstandardNth1gt13) + - gtu22*gtu32*JacPDstandardNth1gt23) - JacPDstandardNth1gt11*SQR(gtu21) - - JacPDstandardNth1gt22*SQR(gtu22) - JacPDstandardNth1gt33*SQR(gtu32); + -2.*(gtu21*(gtu22*JacPDstandardNth1gt12 + gtu32*JacPDstandardNth1gt13) + + gtu22*gtu32*JacPDstandardNth1gt23) - + 1.*(JacPDstandardNth1gt11*SQR(gtu21) + JacPDstandardNth1gt22*SQR(gtu22) + + JacPDstandardNth1gt33*SQR(gtu32)); CCTK_REAL CCTK_ATTRIBUTE_UNUSED dgtu321 = - -(gtu21*(gtu31*JacPDstandardNth1gt11 + gtu32*JacPDstandardNth1gt12 + - gtu33*JacPDstandardNth1gt13)) - gtu22*(gtu31*JacPDstandardNth1gt12 + - gtu32*JacPDstandardNth1gt22 + gtu33*JacPDstandardNth1gt23) - + -1.*(gtu21*(gtu31*JacPDstandardNth1gt11 + gtu32*JacPDstandardNth1gt12 + + gtu33*JacPDstandardNth1gt13) + gtu22*(gtu31*JacPDstandardNth1gt12 + + gtu32*JacPDstandardNth1gt22 + gtu33*JacPDstandardNth1gt23) + gtu32*(gtu31*JacPDstandardNth1gt13 + gtu32*JacPDstandardNth1gt23 + - gtu33*JacPDstandardNth1gt33); + gtu33*JacPDstandardNth1gt33)); CCTK_REAL CCTK_ATTRIBUTE_UNUSED dgtu331 = - -2*(gtu31*(gtu32*JacPDstandardNth1gt12 + gtu33*JacPDstandardNth1gt13) + - gtu32*gtu33*JacPDstandardNth1gt23) - JacPDstandardNth1gt11*SQR(gtu31) - - JacPDstandardNth1gt22*SQR(gtu32) - JacPDstandardNth1gt33*SQR(gtu33); + -2.*(gtu31*(gtu32*JacPDstandardNth1gt12 + gtu33*JacPDstandardNth1gt13) + + gtu32*gtu33*JacPDstandardNth1gt23) - + 1.*(JacPDstandardNth1gt11*SQR(gtu31) + JacPDstandardNth1gt22*SQR(gtu32) + + JacPDstandardNth1gt33*SQR(gtu33)); CCTK_REAL CCTK_ATTRIBUTE_UNUSED dgtu112 = - -2*(gtu11*(gtu21*JacPDstandardNth2gt12 + gtu31*JacPDstandardNth2gt13) + - gtu21*gtu31*JacPDstandardNth2gt23) - JacPDstandardNth2gt11*SQR(gtu11) - - JacPDstandardNth2gt22*SQR(gtu21) - JacPDstandardNth2gt33*SQR(gtu31); + -2.*(gtu11*(gtu21*JacPDstandardNth2gt12 + gtu31*JacPDstandardNth2gt13) + + gtu21*gtu31*JacPDstandardNth2gt23) - + 1.*(JacPDstandardNth2gt11*SQR(gtu11) + JacPDstandardNth2gt22*SQR(gtu21) + + JacPDstandardNth2gt33*SQR(gtu31)); CCTK_REAL CCTK_ATTRIBUTE_UNUSED dgtu212 = - -(gtu21*(gtu11*JacPDstandardNth2gt11 + gtu21*JacPDstandardNth2gt12 + - gtu31*JacPDstandardNth2gt13)) - gtu22*(gtu11*JacPDstandardNth2gt12 + - gtu21*JacPDstandardNth2gt22 + gtu31*JacPDstandardNth2gt23) - + -1.*(gtu21*(gtu11*JacPDstandardNth2gt11 + gtu21*JacPDstandardNth2gt12 + + gtu31*JacPDstandardNth2gt13) + gtu22*(gtu11*JacPDstandardNth2gt12 + + gtu21*JacPDstandardNth2gt22 + gtu31*JacPDstandardNth2gt23) + gtu32*(gtu11*JacPDstandardNth2gt13 + gtu21*JacPDstandardNth2gt23 + - gtu31*JacPDstandardNth2gt33); + gtu31*JacPDstandardNth2gt33)); CCTK_REAL CCTK_ATTRIBUTE_UNUSED dgtu312 = - -(gtu11*(gtu31*JacPDstandardNth2gt11 + gtu32*JacPDstandardNth2gt12 + - gtu33*JacPDstandardNth2gt13)) - gtu21*(gtu31*JacPDstandardNth2gt12 + - gtu32*JacPDstandardNth2gt22 + gtu33*JacPDstandardNth2gt23) - + -1.*(gtu11*(gtu31*JacPDstandardNth2gt11 + gtu32*JacPDstandardNth2gt12 + + gtu33*JacPDstandardNth2gt13) + gtu21*(gtu31*JacPDstandardNth2gt12 + + gtu32*JacPDstandardNth2gt22 + gtu33*JacPDstandardNth2gt23) + gtu31*(gtu31*JacPDstandardNth2gt13 + gtu32*JacPDstandardNth2gt23 + - gtu33*JacPDstandardNth2gt33); + gtu33*JacPDstandardNth2gt33)); CCTK_REAL CCTK_ATTRIBUTE_UNUSED dgtu222 = - -2*(gtu21*(gtu22*JacPDstandardNth2gt12 + gtu32*JacPDstandardNth2gt13) + - gtu22*gtu32*JacPDstandardNth2gt23) - JacPDstandardNth2gt11*SQR(gtu21) - - JacPDstandardNth2gt22*SQR(gtu22) - JacPDstandardNth2gt33*SQR(gtu32); + -2.*(gtu21*(gtu22*JacPDstandardNth2gt12 + gtu32*JacPDstandardNth2gt13) + + gtu22*gtu32*JacPDstandardNth2gt23) - + 1.*(JacPDstandardNth2gt11*SQR(gtu21) + JacPDstandardNth2gt22*SQR(gtu22) + + JacPDstandardNth2gt33*SQR(gtu32)); CCTK_REAL CCTK_ATTRIBUTE_UNUSED dgtu322 = - -(gtu21*(gtu31*JacPDstandardNth2gt11 + gtu32*JacPDstandardNth2gt12 + - gtu33*JacPDstandardNth2gt13)) - gtu22*(gtu31*JacPDstandardNth2gt12 + - gtu32*JacPDstandardNth2gt22 + gtu33*JacPDstandardNth2gt23) - + -1.*(gtu21*(gtu31*JacPDstandardNth2gt11 + gtu32*JacPDstandardNth2gt12 + + gtu33*JacPDstandardNth2gt13) + gtu22*(gtu31*JacPDstandardNth2gt12 + + gtu32*JacPDstandardNth2gt22 + gtu33*JacPDstandardNth2gt23) + gtu32*(gtu31*JacPDstandardNth2gt13 + gtu32*JacPDstandardNth2gt23 + - gtu33*JacPDstandardNth2gt33); + gtu33*JacPDstandardNth2gt33)); CCTK_REAL CCTK_ATTRIBUTE_UNUSED dgtu332 = - -2*(gtu31*(gtu32*JacPDstandardNth2gt12 + gtu33*JacPDstandardNth2gt13) + - gtu32*gtu33*JacPDstandardNth2gt23) - JacPDstandardNth2gt11*SQR(gtu31) - - JacPDstandardNth2gt22*SQR(gtu32) - JacPDstandardNth2gt33*SQR(gtu33); + -2.*(gtu31*(gtu32*JacPDstandardNth2gt12 + gtu33*JacPDstandardNth2gt13) + + gtu32*gtu33*JacPDstandardNth2gt23) - + 1.*(JacPDstandardNth2gt11*SQR(gtu31) + JacPDstandardNth2gt22*SQR(gtu32) + + JacPDstandardNth2gt33*SQR(gtu33)); CCTK_REAL CCTK_ATTRIBUTE_UNUSED dgtu113 = - -2*(gtu11*(gtu21*JacPDstandardNth3gt12 + gtu31*JacPDstandardNth3gt13) + - gtu21*gtu31*JacPDstandardNth3gt23) - JacPDstandardNth3gt11*SQR(gtu11) - - JacPDstandardNth3gt22*SQR(gtu21) - JacPDstandardNth3gt33*SQR(gtu31); + -2.*(gtu11*(gtu21*JacPDstandardNth3gt12 + gtu31*JacPDstandardNth3gt13) + + gtu21*gtu31*JacPDstandardNth3gt23) - + 1.*(JacPDstandardNth3gt11*SQR(gtu11) + JacPDstandardNth3gt22*SQR(gtu21) + + JacPDstandardNth3gt33*SQR(gtu31)); CCTK_REAL CCTK_ATTRIBUTE_UNUSED dgtu213 = - -(gtu21*(gtu11*JacPDstandardNth3gt11 + gtu21*JacPDstandardNth3gt12 + - gtu31*JacPDstandardNth3gt13)) - gtu22*(gtu11*JacPDstandardNth3gt12 + - gtu21*JacPDstandardNth3gt22 + gtu31*JacPDstandardNth3gt23) - + -1.*(gtu21*(gtu11*JacPDstandardNth3gt11 + gtu21*JacPDstandardNth3gt12 + + gtu31*JacPDstandardNth3gt13) + gtu22*(gtu11*JacPDstandardNth3gt12 + + gtu21*JacPDstandardNth3gt22 + gtu31*JacPDstandardNth3gt23) + gtu32*(gtu11*JacPDstandardNth3gt13 + gtu21*JacPDstandardNth3gt23 + - gtu31*JacPDstandardNth3gt33); + gtu31*JacPDstandardNth3gt33)); CCTK_REAL CCTK_ATTRIBUTE_UNUSED dgtu313 = - -(gtu11*(gtu31*JacPDstandardNth3gt11 + gtu32*JacPDstandardNth3gt12 + - gtu33*JacPDstandardNth3gt13)) - gtu21*(gtu31*JacPDstandardNth3gt12 + - gtu32*JacPDstandardNth3gt22 + gtu33*JacPDstandardNth3gt23) - + -1.*(gtu11*(gtu31*JacPDstandardNth3gt11 + gtu32*JacPDstandardNth3gt12 + + gtu33*JacPDstandardNth3gt13) + gtu21*(gtu31*JacPDstandardNth3gt12 + + gtu32*JacPDstandardNth3gt22 + gtu33*JacPDstandardNth3gt23) + gtu31*(gtu31*JacPDstandardNth3gt13 + gtu32*JacPDstandardNth3gt23 + - gtu33*JacPDstandardNth3gt33); + gtu33*JacPDstandardNth3gt33)); CCTK_REAL CCTK_ATTRIBUTE_UNUSED dgtu223 = - -2*(gtu21*(gtu22*JacPDstandardNth3gt12 + gtu32*JacPDstandardNth3gt13) + - gtu22*gtu32*JacPDstandardNth3gt23) - JacPDstandardNth3gt11*SQR(gtu21) - - JacPDstandardNth3gt22*SQR(gtu22) - JacPDstandardNth3gt33*SQR(gtu32); + -2.*(gtu21*(gtu22*JacPDstandardNth3gt12 + gtu32*JacPDstandardNth3gt13) + + gtu22*gtu32*JacPDstandardNth3gt23) - + 1.*(JacPDstandardNth3gt11*SQR(gtu21) + JacPDstandardNth3gt22*SQR(gtu22) + + JacPDstandardNth3gt33*SQR(gtu32)); CCTK_REAL CCTK_ATTRIBUTE_UNUSED dgtu323 = - -(gtu21*(gtu31*JacPDstandardNth3gt11 + gtu32*JacPDstandardNth3gt12 + - gtu33*JacPDstandardNth3gt13)) - gtu22*(gtu31*JacPDstandardNth3gt12 + - gtu32*JacPDstandardNth3gt22 + gtu33*JacPDstandardNth3gt23) - + -1.*(gtu21*(gtu31*JacPDstandardNth3gt11 + gtu32*JacPDstandardNth3gt12 + + gtu33*JacPDstandardNth3gt13) + gtu22*(gtu31*JacPDstandardNth3gt12 + + gtu32*JacPDstandardNth3gt22 + gtu33*JacPDstandardNth3gt23) + gtu32*(gtu31*JacPDstandardNth3gt13 + gtu32*JacPDstandardNth3gt23 + - gtu33*JacPDstandardNth3gt33); + gtu33*JacPDstandardNth3gt33)); CCTK_REAL CCTK_ATTRIBUTE_UNUSED dgtu333 = - -2*(gtu31*(gtu32*JacPDstandardNth3gt12 + gtu33*JacPDstandardNth3gt13) + - gtu32*gtu33*JacPDstandardNth3gt23) - JacPDstandardNth3gt11*SQR(gtu31) - - JacPDstandardNth3gt22*SQR(gtu32) - JacPDstandardNth3gt33*SQR(gtu33); + -2.*(gtu31*(gtu32*JacPDstandardNth3gt12 + gtu33*JacPDstandardNth3gt13) + + gtu32*gtu33*JacPDstandardNth3gt23) - + 1.*(JacPDstandardNth3gt11*SQR(gtu31) + JacPDstandardNth3gt22*SQR(gtu32) + + JacPDstandardNth3gt33*SQR(gtu33)); CCTK_REAL CCTK_ATTRIBUTE_UNUSED Gtl111 = 0.5*JacPDstandardNth1gt11; @@ -1185,8 +1194,9 @@ static void ML_ADMQuantities_Body(cGH const * restrict const cctkGH, int const d CCTK_REAL CCTK_ATTRIBUTE_UNUSED Gtl122 = -0.5*JacPDstandardNth1gt22 + JacPDstandardNth2gt12; - CCTK_REAL CCTK_ATTRIBUTE_UNUSED Gtl123 = 0.5*(-JacPDstandardNth1gt23 + - JacPDstandardNth2gt13 + JacPDstandardNth3gt12); + CCTK_REAL CCTK_ATTRIBUTE_UNUSED Gtl123 = + 0.5*(-1.*JacPDstandardNth1gt23 + JacPDstandardNth2gt13 + + JacPDstandardNth3gt12); CCTK_REAL CCTK_ATTRIBUTE_UNUSED Gtl133 = -0.5*JacPDstandardNth1gt33 + JacPDstandardNth3gt13; @@ -1197,7 +1207,7 @@ static void ML_ADMQuantities_Body(cGH const * restrict const cctkGH, int const d CCTK_REAL CCTK_ATTRIBUTE_UNUSED Gtl212 = 0.5*JacPDstandardNth1gt22; CCTK_REAL CCTK_ATTRIBUTE_UNUSED Gtl213 = 0.5*(JacPDstandardNth1gt23 - - JacPDstandardNth2gt13 + JacPDstandardNth3gt12); + 1.*JacPDstandardNth2gt13 + JacPDstandardNth3gt12); CCTK_REAL CCTK_ATTRIBUTE_UNUSED Gtl222 = 0.5*JacPDstandardNth2gt22; @@ -1210,7 +1220,7 @@ static void ML_ADMQuantities_Body(cGH const * restrict const cctkGH, int const d 0.5*JacPDstandardNth3gt11; CCTK_REAL CCTK_ATTRIBUTE_UNUSED Gtl312 = 0.5*(JacPDstandardNth1gt23 + - JacPDstandardNth2gt13 - JacPDstandardNth3gt12); + JacPDstandardNth2gt13 - 1.*JacPDstandardNth3gt12); CCTK_REAL CCTK_ATTRIBUTE_UNUSED Gtl313 = 0.5*JacPDstandardNth1gt33; @@ -1357,104 +1367,102 @@ static void ML_ADMQuantities_Body(cGH const * restrict const cctkGH, int const d Gtl333*gtu33; CCTK_REAL CCTK_ATTRIBUTE_UNUSED Xtn1 = Gt111*gtu11 + Gt122*gtu22 + - 2*(Gt112*gtu21 + Gt113*gtu31 + Gt123*gtu32) + Gt133*gtu33; + 2.*(Gt112*gtu21 + Gt113*gtu31 + Gt123*gtu32) + Gt133*gtu33; CCTK_REAL CCTK_ATTRIBUTE_UNUSED Xtn2 = Gt211*gtu11 + Gt222*gtu22 + - 2*(Gt212*gtu21 + Gt213*gtu31 + Gt223*gtu32) + Gt233*gtu33; + 2.*(Gt212*gtu21 + Gt213*gtu31 + Gt223*gtu32) + Gt233*gtu33; CCTK_REAL CCTK_ATTRIBUTE_UNUSED Xtn3 = Gt311*gtu11 + Gt322*gtu22 + - 2*(Gt312*gtu21 + Gt313*gtu31 + Gt323*gtu32) + Gt333*gtu33; + 2.*(Gt312*gtu21 + Gt313*gtu31 + Gt323*gtu32) + Gt333*gtu33; - CCTK_REAL CCTK_ATTRIBUTE_UNUSED Rt11 = 3*(Gt111*Gtlu111 + - Gt112*Gtlu112 + Gt113*Gtlu113) + 2*(Gt211*Gtlu121 + Gt212*Gtlu122 + + CCTK_REAL CCTK_ATTRIBUTE_UNUSED 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 + Gt213*Gtlu213 + Gt311*Gtlu311 + Gt312*Gtlu312 + Gt313*Gtlu313 + gt11L*JacPDstandardNth1Xt1 + - gt12L*JacPDstandardNth1Xt2 + gt13L*JacPDstandardNth1Xt3 + - 0.5*(-(gtu11*JacPDstandardNth11gt11) - gtu21*(JacPDstandardNth12gt11 + - JacPDstandardNth21gt11) - gtu22*JacPDstandardNth22gt11 + - gtu31*(-JacPDstandardNth13gt11 - JacPDstandardNth31gt11) + - gtu32*(-JacPDstandardNth23gt11 - JacPDstandardNth32gt11) - + gt12L*JacPDstandardNth1Xt2 + gt13L*JacPDstandardNth1Xt3 - + 0.5*(gtu11*JacPDstandardNth11gt11 + gtu21*(JacPDstandardNth12gt11 + + JacPDstandardNth21gt11) + gtu22*JacPDstandardNth22gt11 + + gtu31*(JacPDstandardNth13gt11 + JacPDstandardNth31gt11) + + gtu32*(JacPDstandardNth23gt11 + JacPDstandardNth32gt11) + gtu33*JacPDstandardNth33gt11) + Gtl111*Xtn1 + Gtl112*Xtn2 + Gtl113*Xtn3; - CCTK_REAL CCTK_ATTRIBUTE_UNUSED Rt12 = 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 + - Gt111*Gtlu211 + Gt112*(Gtlu111 + Gtlu122 + Gtlu212) + Gt113*Gtlu213 + - Gt311*Gtlu231 + Gt312*Gtlu232 + Gt313*Gtlu233 + Gt311*Gtlu321 + - Gt312*Gtlu322 + Gt313*Gtlu323) - gtu11*JacPDstandardNth11gt12 + - gt12L*JacPDstandardNth1Xt1 + gt22L*JacPDstandardNth1Xt2 + - gt23L*JacPDstandardNth1Xt3 + gtu21*(-JacPDstandardNth12gt12 - - JacPDstandardNth21gt12) - gtu22*JacPDstandardNth22gt12 + - gt11L*JacPDstandardNth2Xt1 + gt12L*JacPDstandardNth2Xt2 + - gt13L*JacPDstandardNth2Xt3 + gtu31*(-JacPDstandardNth13gt12 - - JacPDstandardNth31gt12) + gtu32*(-JacPDstandardNth23gt12 - - 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 + - Gt133*Gtlu113 + Gt213*Gtlu121 + Gt223*Gtlu122 + Gt233*Gtlu123 + - Gt111*Gtlu131 + Gt313*Gtlu131 + Gt112*Gtlu132 + Gt323*Gtlu132 + - Gt333*Gtlu133 + Gt211*Gtlu231 + Gt212*Gtlu232 + Gt213*Gtlu233 + - Gt111*Gtlu311 + Gt112*Gtlu312 + Gt113*(Gtlu111 + Gtlu133 + Gtlu313) + - Gt211*Gtlu321 + Gt212*Gtlu322 + Gt213*Gtlu323) + 4*(Gt311*Gtlu331 + - Gt312*Gtlu332 + Gt313*Gtlu333) - gtu11*JacPDstandardNth11gt13 + - gt13L*JacPDstandardNth1Xt1 + gt23L*JacPDstandardNth1Xt2 + - gt33L*JacPDstandardNth1Xt3 + gtu21*(-JacPDstandardNth12gt13 - - JacPDstandardNth21gt13) - gtu22*JacPDstandardNth22gt13 + - 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); - - CCTK_REAL CCTK_ATTRIBUTE_UNUSED 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 + + CCTK_REAL CCTK_ATTRIBUTE_UNUSED Rt12 = 0.5*(4.*(Gt211*Gtlu221 + + Gt212*Gtlu222 + Gt213*Gtlu223) + 2.*(Gt122*Gtlu112 + Gt123*Gtlu113 + + (Gt111 + Gt212)*Gtlu121 + Gt222*Gtlu122 + (Gt113 + Gt223)*Gtlu123 + + Gt322*Gtlu132 + Gt323*Gtlu133 + Gt111*Gtlu211 + Gt112*(Gtlu111 + + Gtlu122 + Gtlu212) + Gt113*Gtlu213 + Gt311*(Gtlu231 + Gtlu321) + + Gt312*(Gtlu131 + Gtlu232 + Gtlu322) + Gt313*(Gtlu233 + Gtlu323)) + + gt22L*JacPDstandardNth1Xt2 + gt23L*JacPDstandardNth1Xt3 + + gt11L*JacPDstandardNth2Xt1 + gt12L*(JacPDstandardNth1Xt1 + + JacPDstandardNth2Xt2) + gt13L*JacPDstandardNth2Xt3 - + 1.*(gtu11*JacPDstandardNth11gt12 + gtu21*(JacPDstandardNth12gt12 + + JacPDstandardNth21gt12) + gtu22*JacPDstandardNth22gt12 + + gtu31*(JacPDstandardNth13gt12 + JacPDstandardNth31gt12) + + gtu32*(JacPDstandardNth23gt12 + JacPDstandardNth32gt12) + + gtu33*JacPDstandardNth33gt12) + (Gtl112 + Gtl211)*Xtn1 + (Gtl122 + + Gtl212)*Xtn2 + (Gtl123 + Gtl213)*Xtn3); + + CCTK_REAL CCTK_ATTRIBUTE_UNUSED Rt13 = 0.5*(2.*(Gt123*Gtlu112 + + Gt133*Gtlu113 + Gt223*Gtlu122 + Gt233*Gtlu123 + (Gt111 + Gt313)*Gtlu131 + + (Gt112 + Gt323)*Gtlu132 + Gt333*Gtlu133 + Gt111*Gtlu311 + + Gt112*Gtlu312 + Gt113*(Gtlu111 + Gtlu133 + Gtlu313) + Gt211*(Gtlu231 + + Gtlu321) + Gt212*(Gtlu232 + Gtlu322) + Gt213*(Gtlu121 + Gtlu233 + + Gtlu323)) + 4.*(Gt311*Gtlu331 + Gt312*Gtlu332 + Gt313*Gtlu333) + + gt23L*JacPDstandardNth1Xt2 + gt33L*JacPDstandardNth1Xt3 - + 1.*(gtu11*JacPDstandardNth11gt13 + gtu21*(JacPDstandardNth12gt13 + + JacPDstandardNth21gt13) + gtu22*JacPDstandardNth22gt13 + + gtu31*(JacPDstandardNth13gt13 + JacPDstandardNth31gt13) + + gtu32*(JacPDstandardNth23gt13 + JacPDstandardNth32gt13) + + gtu33*JacPDstandardNth33gt13) + gt11L*JacPDstandardNth3Xt1 + + gt12L*JacPDstandardNth3Xt2 + gt13L*(JacPDstandardNth1Xt1 + + JacPDstandardNth3Xt3) + (Gtl113 + Gtl311)*Xtn1 + (Gtl123 + Gtl312)*Xtn2 + + (Gtl133 + Gtl313)*Xtn3); + + CCTK_REAL CCTK_ATTRIBUTE_UNUSED 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*JacPDstandardNth2Xt1 + - gt22L*JacPDstandardNth2Xt2 + gt23L*JacPDstandardNth2Xt3 + - 0.5*(-(gtu11*JacPDstandardNth11gt22) - gtu21*(JacPDstandardNth12gt22 + - JacPDstandardNth21gt22) - gtu22*JacPDstandardNth22gt22 + - gtu31*(-JacPDstandardNth13gt22 - JacPDstandardNth31gt22) + - gtu32*(-JacPDstandardNth23gt22 - JacPDstandardNth32gt22) - + 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 + - Gt113*Gtlu211 + Gt123*Gtlu212 + Gt133*Gtlu213 + Gt213*Gtlu221 + - Gt223*Gtlu222 + Gt233*Gtlu223 + Gt212*Gtlu231 + Gt313*Gtlu231 + - Gt222*Gtlu232 + Gt323*Gtlu232 + Gt223*Gtlu233 + Gt333*Gtlu233 + - Gt112*(Gtlu131 + Gtlu311) + Gt122*(Gtlu132 + Gtlu312) + Gt123*Gtlu313 + - Gt212*Gtlu321 + Gt222*Gtlu322 + Gt223*Gtlu323) + 4*(Gt312*Gtlu331 + - Gt322*Gtlu332 + Gt323*Gtlu333) - gtu11*JacPDstandardNth11gt23 + - gtu21*(-JacPDstandardNth12gt23 - JacPDstandardNth21gt23) - - gtu22*JacPDstandardNth22gt23 + gt13L*JacPDstandardNth2Xt1 + - gt23L*JacPDstandardNth2Xt2 + gt33L*JacPDstandardNth2Xt3 + - 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); - - CCTK_REAL CCTK_ATTRIBUTE_UNUSED 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*JacPDstandardNth11gt33) - - gtu21*(JacPDstandardNth12gt33 + JacPDstandardNth21gt33) - - gtu22*JacPDstandardNth22gt33 + gtu31*(-JacPDstandardNth13gt33 - - JacPDstandardNth31gt33) + gtu32*(-JacPDstandardNth23gt33 - - JacPDstandardNth32gt33) - gtu33*JacPDstandardNth33gt33) + + CCTK_REAL CCTK_ATTRIBUTE_UNUSED Rt23 = 0.5*(2.*(Gt113*Gtlu211 + + Gt133*Gtlu213 + Gt213*Gtlu221 + Gt233*Gtlu223 + (Gt212 + Gt313)*Gtlu231 + + (Gt222 + Gt323)*Gtlu232 + Gt333*Gtlu233 + Gt112*(Gtlu131 + Gtlu311) + + Gt122*(Gtlu132 + Gtlu312) + Gt123*(Gtlu133 + Gtlu212 + Gtlu313) + + Gt212*Gtlu321 + Gt222*Gtlu322 + Gt223*(Gtlu222 + Gtlu233 + Gtlu323)) + + 4.*(Gt312*Gtlu331 + Gt322*Gtlu332 + Gt323*Gtlu333) + + gt13L*JacPDstandardNth2Xt1 + gt33L*JacPDstandardNth2Xt3 - + 1.*(gtu11*JacPDstandardNth11gt23 + gtu21*(JacPDstandardNth12gt23 + + JacPDstandardNth21gt23) + gtu22*JacPDstandardNth22gt23 + + gtu31*(JacPDstandardNth13gt23 + JacPDstandardNth31gt23) + + gtu32*(JacPDstandardNth23gt23 + JacPDstandardNth32gt23) + + gtu33*JacPDstandardNth33gt23) + gt12L*JacPDstandardNth3Xt1 + + gt22L*JacPDstandardNth3Xt2 + gt23L*(JacPDstandardNth2Xt2 + + JacPDstandardNth3Xt3) + (Gtl213 + Gtl312)*Xtn1 + (Gtl223 + Gtl322)*Xtn2 + + (Gtl233 + Gtl323)*Xtn3); + + CCTK_REAL CCTK_ATTRIBUTE_UNUSED 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*JacPDstandardNth11gt33 + + gtu21*(JacPDstandardNth12gt33 + JacPDstandardNth21gt33) + + gtu22*JacPDstandardNth22gt33 + gtu31*(JacPDstandardNth13gt33 + + JacPDstandardNth31gt33) + gtu32*(JacPDstandardNth23gt33 + + JacPDstandardNth32gt33) + gtu33*JacPDstandardNth33gt33) + gt13L*JacPDstandardNth3Xt1 + gt23L*JacPDstandardNth3Xt2 + gt33L*JacPDstandardNth3Xt3 + Gtl313*Xtn1 + Gtl323*Xtn2 + Gtl333*Xtn3; CCTK_REAL CCTK_ATTRIBUTE_UNUSED trRt = gtu11*Rt11 + gtu22*Rt22 + - 2*(gtu21*Rt12 + gtu31*Rt13 + gtu32*Rt23) + gtu33*Rt33; + 2.*(gtu21*Rt12 + gtu31*Rt13 + gtu32*Rt23) + gtu33*Rt33; CCTK_REAL CCTK_ATTRIBUTE_UNUSED ephi = IfThen(conformalMethod,INV(sqrt(phiL)),exp(phiL)); @@ -1487,17 +1495,17 @@ static void ML_ADMQuantities_Body(cGH const * restrict const cctkGH, int const d + At33L*gtu33; CCTK_REAL CCTK_ATTRIBUTE_UNUSED rho = INV(SQR(alphaL))*(eTttL - - 2*(beta2L*eTtyL + beta3L*eTtzL) + 2*(beta1L*(-eTtxL + + 2.*(beta2L*eTtyL + beta3L*eTtzL) + 2.*(beta1L*(-1.*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 CCTK_ATTRIBUTE_UNUSED S1 = (-1.*eTtxL + beta1L*eTxxL + beta2L*eTxyL + beta3L*eTxzL)*INV(alphaL); - CCTK_REAL CCTK_ATTRIBUTE_UNUSED S2 = (-eTtyL + beta1L*eTxyL + + CCTK_REAL CCTK_ATTRIBUTE_UNUSED S2 = (-1.*eTtyL + beta1L*eTxyL + beta2L*eTyyL + beta3L*eTyzL)*INV(alphaL); - CCTK_REAL CCTK_ATTRIBUTE_UNUSED S3 = (-eTtzL + beta1L*eTxzL + + CCTK_REAL CCTK_ATTRIBUTE_UNUSED S3 = (-1.*eTtzL + beta1L*eTxzL + beta2L*eTyzL + beta3L*eTzzL)*INV(alphaL); CCTK_REAL CCTK_ATTRIBUTE_UNUSED MadmL = @@ -1527,37 +1535,37 @@ static void ML_ADMQuantities_Body(cGH const * restrict const cctkGH, int const d Gtlu323*(Gt213*gtu31 + Gt223*gtu32 + Gt233*gtu33) + Gtlu331*(Gt311*gtu31 + Gt312*gtu32 + Gt313*gtu33) + 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))); + Gtlu333*(Gt313*gtu31 + Gt323*gtu32 + Gt333*gtu33) + (-1. + ephi)*trRt - + 1.*pow(ephi,5.)*(2.*(Atm12*Atm21 + Atm13*Atm31 + Atm23*Atm32) + + 50.26548245743669*rho - 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 + + 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 - - 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 + - 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 - - 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 + + 1.33333333333333333333333333333*JacPDstandardNth2trK - + 50.26548245743669*S2) + yL*(-1.*(At11L*dgtu113 + At22L*dgtu223 + + 2.*(At12L*dgtu213 + At13L*dgtu313 + At23L*dgtu323) + + At33L*dgtu333) + 1.33333333333333333333333333333*JacPDstandardNth3trK + + 50.26548245743669*S3))*INV(Pi)*pow(ephi,6.); + + CCTK_REAL CCTK_ATTRIBUTE_UNUSED Jadm2L = 0.0625*(-2.*Atm13 + 2.*Atm31 + + zL*(-1.*(At11L*dgtu111 + At22L*dgtu221 + 2.*(At12L*dgtu211 + + At13L*dgtu311 + At23L*dgtu321) + At33L*dgtu331) + + 1.33333333333333333333333333333*JacPDstandardNth1trK + + 50.26548245743669*S1) + xL*(At11L*dgtu113 + At22L*dgtu223 + + 2.*(At12L*dgtu213 + At13L*dgtu313 + At23L*dgtu323) + + At33L*dgtu333 - 1.33333333333333333333333333333*JacPDstandardNth3trK + - 50.26548245743669*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 - - 1.33333333333333333333333333333*JacPDstandardNth1trK - 16*Pi*S1) + - xL*(-(At11L*dgtu112) - 2*At12L*dgtu212 - At22L*dgtu222 - - 2*At13L*dgtu312 - 2*At23L*dgtu322 - At33L*dgtu332 + - 1.33333333333333333333333333333*JacPDstandardNth2trK + - 16*Pi*S2))*INV(Pi)*pow(ephi,6); + 1.33333333333333333333333333333*JacPDstandardNth1trK - + 50.26548245743669*S1) + xL*(-1.*(At11L*dgtu112 + At22L*dgtu222 + + 2.*(At12L*dgtu212 + At13L*dgtu312 + At23L*dgtu322) + + At33L*dgtu332) + 1.33333333333333333333333333333*JacPDstandardNth2trK + + 50.26548245743669*S2))*INV(Pi)*pow(ephi,6.); /* Copy local copies back to grid functions */ Jadm1[index] = Jadm1L; |