diff options
| author | Barry Wardell <barry.wardell@gmail.com> | 2012-12-18 12:45:42 +0000 |
|---|---|---|
| committer | Barry Wardell <barry.wardell@gmail.com> | 2012-12-18 12:45:42 +0000 |
| commit | 962eb452acb41cdcb83824a412cdabf26481f0de (patch) | |
| tree | e4f5a2fca947a183184d54e56e01a727792a6657 | |
| parent | a25f808933aee6548579ee20c8b5cafa183ce6de (diff) | |
Regenerate thorns which were missed in 45bcf7046474c39dc6fb7c3e7729cdd1143a70fe.
| -rw-r--r-- | ML_ADM/schedule.ccl | 84 | ||||
| -rw-r--r-- | ML_ADM/src/ML_ADM_RHS.cc | 12 | ||||
| -rw-r--r-- | ML_ADM/src/ML_ADM_constraints.cc | 12 | ||||
| -rw-r--r-- | ML_ADMConstraints/schedule.ccl | 30 | ||||
| -rw-r--r-- | ML_ADMConstraints/src/ML_ADMConstraints.cc | 304 | ||||
| -rw-r--r-- | ML_ADMQuantities/schedule.ccl | 36 | ||||
| -rw-r--r-- | ML_ADMQuantities/src/ML_ADMQuantities.cc | 192 | ||||
| -rw-r--r-- | ML_WaveToy/src/WT_Standing.cc | 2 | ||||
| -rw-r--r-- | ML_WaveToy_CL/src/WT_CL_Standing.cc | 2 |
9 files changed, 348 insertions, 326 deletions
diff --git a/ML_ADM/schedule.ccl b/ML_ADM/schedule.ccl index ce81eaf..b39fb2a 100644 --- a/ML_ADM/schedule.ccl +++ b/ML_ADM/schedule.ccl @@ -158,22 +158,22 @@ if (CCTK_EQUALS(my_initial_data, "ADMBase")) schedule ML_ADM_convertFromADMBase AT initial AFTER ADMBase_PostInitial { LANG: C + READS: ADMBase::kxx + READS: ADMBase::kxy + READS: ADMBase::kxz + READS: ADMBase::kyy + READS: ADMBase::kyz + READS: ADMBase::kzz READS: ADMBase::alp - READS: ADMBase::betax - READS: ADMBase::betay - READS: ADMBase::betaz READS: ADMBase::gxx READS: ADMBase::gxy READS: ADMBase::gxz READS: ADMBase::gyy READS: ADMBase::gyz READS: ADMBase::gzz - READS: ADMBase::kxx - READS: ADMBase::kxy - READS: ADMBase::kxz - READS: ADMBase::kyy - READS: ADMBase::kyz - READS: ADMBase::kzz + READS: ADMBase::betax + READS: ADMBase::betay + READS: ADMBase::betaz WRITES: ML_ADM::alpha WRITES: ML_ADM::beta1 WRITES: ML_ADM::beta2 @@ -196,22 +196,22 @@ if (CCTK_EQUALS(my_initial_data, "ADMBase")) schedule ML_ADM_RHS IN MoL_CalcRHS { LANG: C + READS: ML_ADM::K11 + READS: ML_ADM::K12 + READS: ML_ADM::K13 + READS: ML_ADM::K22 + READS: ML_ADM::K23 + READS: ML_ADM::K33 READS: ML_ADM::alpha - READS: ML_ADM::beta1 - READS: ML_ADM::beta2 - READS: ML_ADM::beta3 READS: ML_ADM::g11 READS: ML_ADM::g12 READS: ML_ADM::g13 READS: ML_ADM::g22 READS: ML_ADM::g23 READS: ML_ADM::g33 - READS: ML_ADM::K11 - READS: ML_ADM::K12 - READS: ML_ADM::K13 - READS: ML_ADM::K22 - READS: ML_ADM::K23 - READS: ML_ADM::K33 + READS: ML_ADM::beta1 + READS: ML_ADM::beta2 + READS: ML_ADM::beta3 WRITES: ML_ADM::alpharhs WRITES: ML_ADM::beta1rhs WRITES: ML_ADM::beta2rhs @@ -237,22 +237,22 @@ schedule ML_ADM_RHS AT analysis SYNC: ML_lapserhs SYNC: ML_metricrhs SYNC: ML_shiftrhs + READS: ML_ADM::K11 + READS: ML_ADM::K12 + READS: ML_ADM::K13 + READS: ML_ADM::K22 + READS: ML_ADM::K23 + READS: ML_ADM::K33 READS: ML_ADM::alpha - READS: ML_ADM::beta1 - READS: ML_ADM::beta2 - READS: ML_ADM::beta3 READS: ML_ADM::g11 READS: ML_ADM::g12 READS: ML_ADM::g13 READS: ML_ADM::g22 READS: ML_ADM::g23 READS: ML_ADM::g33 - READS: ML_ADM::K11 - READS: ML_ADM::K12 - READS: ML_ADM::K13 - READS: ML_ADM::K22 - READS: ML_ADM::K23 - READS: ML_ADM::K33 + READS: ML_ADM::beta1 + READS: ML_ADM::beta2 + READS: ML_ADM::beta3 WRITES: ML_ADM::alpharhs WRITES: ML_ADM::beta1rhs WRITES: ML_ADM::beta2rhs @@ -299,22 +299,22 @@ if (CCTK_EQUALS(my_boundary_condition, "Minkowski")) schedule ML_ADM_convertToADMBase IN MoL_PostStep AFTER ML_ADM_ApplyBCs { LANG: C + READS: ML_ADM::K11 + READS: ML_ADM::K12 + READS: ML_ADM::K13 + READS: ML_ADM::K22 + READS: ML_ADM::K23 + READS: ML_ADM::K33 READS: ML_ADM::alpha - READS: ML_ADM::beta1 - READS: ML_ADM::beta2 - READS: ML_ADM::beta3 READS: ML_ADM::g11 READS: ML_ADM::g12 READS: ML_ADM::g13 READS: ML_ADM::g22 READS: ML_ADM::g23 READS: ML_ADM::g33 - READS: ML_ADM::K11 - READS: ML_ADM::K12 - READS: ML_ADM::K13 - READS: ML_ADM::K22 - READS: ML_ADM::K23 - READS: ML_ADM::K33 + READS: ML_ADM::beta1 + READS: ML_ADM::beta2 + READS: ML_ADM::beta3 WRITES: ADMBase::alp WRITES: ADMBase::betax WRITES: ADMBase::betay @@ -342,18 +342,18 @@ schedule ML_ADM_constraints AT analysis LANG: C SYNC: ML_Ham SYNC: ML_mom - READS: ML_ADM::g11 - READS: ML_ADM::g12 - READS: ML_ADM::g13 - READS: ML_ADM::g22 - READS: ML_ADM::g23 - READS: ML_ADM::g33 READS: ML_ADM::K11 READS: ML_ADM::K12 READS: ML_ADM::K13 READS: ML_ADM::K22 READS: ML_ADM::K23 READS: ML_ADM::K33 + READS: ML_ADM::g11 + READS: ML_ADM::g12 + READS: ML_ADM::g13 + READS: ML_ADM::g22 + READS: ML_ADM::g23 + READS: ML_ADM::g33 WRITES: ML_ADM::H WRITES: ML_ADM::M1 WRITES: ML_ADM::M2 diff --git a/ML_ADM/src/ML_ADM_RHS.cc b/ML_ADM/src/ML_ADM_RHS.cc index 16b162d..3049ef0 100644 --- a/ML_ADM/src/ML_ADM_RHS.cc +++ b/ML_ADM/src/ML_ADM_RHS.cc @@ -1388,7 +1388,7 @@ static void ML_ADM_RHS_Body(cGH const * restrict const cctkGH, int const dir, in kmul(ToReal(0.5),kmadd(gu33,JacPDstandardNth3g33,kmadd(gu31,kmsub(JacPDstandardNth3g13,ToReal(2),JacPDstandardNth1g33),kmul(gu32,kmsub(JacPDstandardNth3g23,ToReal(2),JacPDstandardNth2g33))))); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED R11 = - kmul(ToReal(0.5),kmadd(gu23,JacPDstandardNth31g12,kmadd(gu32,JacPDstandardNth31g12,kmadd(G111,kmul(kadd(G212,G313),ToReal(-2)),knmsub(JacPDstandardNth11g23,kadd(gu32,gu23),kmadd(gu12,ksub(JacPDstandardNth21g11,JacPDstandardNth12g11),kmadd(gu13,ksub(JacPDstandardNth31g11,JacPDstandardNth13g11),kmadd(gu23,ksub(JacPDstandardNth21g13,JacPDstandardNth23g11),kmadd(gu32,ksub(JacPDstandardNth21g13,JacPDstandardNth32g11),kmadd(kmul(G212,G212),ToReal(2),kmadd(kmul(G313,G313),ToReal(2),kmadd(G211,kmadd(G222,ToReal(-2),kmadd(G323,ToReal(-2),kmul(G112,ToReal(2)))),kmadd(G311,kmadd(G223,ToReal(-2),kmadd(G333,ToReal(-2),kmul(G113,ToReal(2)))),kmadd(gu22,ksub(kmsub(JacPDstandardNth21g12,ToReal(2),JacPDstandardNth22g11),JacPDstandardNth11g22),kmadd(gu33,ksub(kmsub(JacPDstandardNth31g13,ToReal(2),JacPDstandardNth33g11),JacPDstandardNth11g33),kmul(G213,kmul(G312,ToReal(4)))))))))))))))))); + kmadd(G111,G111,kmadd(G212,G212,kmadd(G313,G313,knmsub(G111,kadd(G111,kadd(G212,G313)),knmsub(G211,kadd(G112,kadd(G222,G323)),knmsub(G311,kadd(G113,kadd(G223,G333)),kmadd(knmsub(gu22,kadd(JacPDstandardNth11g22,kmadd(JacPDstandardNth21g12,ToReal(-2),JacPDstandardNth22g11)),knmsub(gu33,kadd(JacPDstandardNth11g33,kmadd(JacPDstandardNth31g13,ToReal(-2),JacPDstandardNth33g11)),kmadd(gu12,ksub(JacPDstandardNth21g11,JacPDstandardNth12g11),kmadd(gu13,ksub(JacPDstandardNth31g11,JacPDstandardNth13g11),kmadd(gu23,kadd(JacPDstandardNth21g13,ksub(JacPDstandardNth31g12,kadd(JacPDstandardNth23g11,JacPDstandardNth11g23))),kmul(gu32,kadd(JacPDstandardNth21g13,ksub(JacPDstandardNth31g12,kadd(JacPDstandardNth32g11,JacPDstandardNth11g23))))))))),ToReal(0.5),kmul(kmadd(G112,G211,kmadd(G113,G311,kmul(G213,G312))),ToReal(2))))))))); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED R12 = kmul(ToReal(0.5),kmadd(gu12,kadd(JacPDstandardNth11g22,JacPDstandardNth22g11),kmadd(gu32,JacPDstandardNth31g22,kmadd(gu13,JacPDstandardNth32g11,kmadd(gu23,JacPDstandardNth32g12,kmadd(gu33,JacPDstandardNth32g13,kmadd(kmadd(G112,kadd(G212,G313),kmadd(G212,G323,kmadd(G312,G333,kmul(gu12,JacPDstandardNth12g12)))),ToReal(-2),knmsub(JacPDstandardNth12g23,kadd(gu32,gu23),kmadd(gu22,ksub(JacPDstandardNth21g22,JacPDstandardNth12g22),kmadd(gu13,ksub(JacPDstandardNth11g23,kadd(JacPDstandardNth13g12,JacPDstandardNth12g13)),kmadd(gu23,ksub(JacPDstandardNth21g23,JacPDstandardNth23g12),kmadd(gu32,ksub(JacPDstandardNth22g13,JacPDstandardNth32g12),kmadd(gu33,ksub(JacPDstandardNth31g23,kadd(JacPDstandardNth33g12,JacPDstandardNth12g33)),kmul(kmadd(G122,G211,kmadd(G123,G311,kmadd(G213,G322,kmul(G313,G323)))),ToReal(2))))))))))))))); @@ -1397,13 +1397,13 @@ static void ML_ADM_RHS_Body(cGH const * restrict const cctkGH, int const dir, in kmul(ToReal(0.5),kmadd(gu22,JacPDstandardNth23g12,kmadd(gu32,JacPDstandardNth31g23,kmadd(gu13,kadd(JacPDstandardNth11g33,JacPDstandardNth33g11),kmadd(gu23,JacPDstandardNth33g12,kmadd(kmadd(G213,G222,kmadd(G223,G313,kmadd(G113,kadd(G212,G313),kmul(gu13,JacPDstandardNth13g13)))),ToReal(-2),knmsub(JacPDstandardNth13g23,kadd(gu32,gu23),kmadd(gu12,kadd(JacPDstandardNth11g23,ksub(JacPDstandardNth23g11,kadd(JacPDstandardNth13g12,JacPDstandardNth12g13))),kmadd(gu33,ksub(JacPDstandardNth31g33,JacPDstandardNth13g33),kmadd(gu22,ksub(JacPDstandardNth21g23,kadd(JacPDstandardNth22g13,JacPDstandardNth13g22)),kmadd(gu23,ksub(JacPDstandardNth21g33,JacPDstandardNth23g13),kmadd(gu32,ksub(JacPDstandardNth23g13,JacPDstandardNth32g13),kmul(kmadd(G123,G211,kmadd(G212,G223,kmadd(G133,G311,kmul(G233,G312)))),ToReal(2)))))))))))))); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED R22 = - kmul(ToReal(0.5),kmadd(kadd(gu13,gu31),JacPDstandardNth12g23,kmadd(gu13,JacPDstandardNth32g12,kmadd(gu31,JacPDstandardNth32g12,kmadd(kmadd(G112,G222,kmadd(G113,G322,kmul(G222,G323))),ToReal(-2),kmadd(gu21,ksub(JacPDstandardNth12g22,JacPDstandardNth21g22),knmsub(gu13,kadd(JacPDstandardNth22g13,JacPDstandardNth13g22),kmadd(gu23,ksub(JacPDstandardNth32g22,JacPDstandardNth23g22),knmsub(gu31,kadd(JacPDstandardNth31g22,JacPDstandardNth22g13),kmadd(kmul(G112,G112),ToReal(2),kmadd(kmul(G323,G323),ToReal(2),kmadd(G122,kmadd(G111,ToReal(-2),kmadd(G313,ToReal(-2),kmul(G212,ToReal(2)))),kmadd(G322,kmadd(G333,ToReal(-2),kmul(G223,ToReal(2))),kmadd(gu11,ksub(kmsub(JacPDstandardNth12g12,ToReal(2),JacPDstandardNth22g11),JacPDstandardNth11g22),kmadd(gu33,ksub(kmsub(JacPDstandardNth32g23,ToReal(2),JacPDstandardNth33g22),JacPDstandardNth22g33),kmul(G123,kmul(G312,ToReal(4)))))))))))))))))); + kmadd(G112,G112,kmadd(G222,G222,kmadd(G323,G323,knmsub(G122,kadd(G111,kadd(G212,G313)),knmsub(G222,kadd(G112,kadd(G222,G323)),knmsub(G322,kadd(G113,kadd(G223,G333)),kmadd(knmsub(gu11,kadd(JacPDstandardNth11g22,kmadd(JacPDstandardNth12g12,ToReal(-2),JacPDstandardNth22g11)),knmsub(gu33,kadd(JacPDstandardNth22g33,kmadd(JacPDstandardNth32g23,ToReal(-2),JacPDstandardNth33g22)),kmadd(gu21,ksub(JacPDstandardNth12g22,JacPDstandardNth21g22),kmadd(gu13,kadd(JacPDstandardNth12g23,ksub(JacPDstandardNth32g12,kadd(JacPDstandardNth22g13,JacPDstandardNth13g22))),kmadd(gu23,ksub(JacPDstandardNth32g22,JacPDstandardNth23g22),kmul(gu31,kadd(JacPDstandardNth12g23,ksub(JacPDstandardNth32g12,kadd(JacPDstandardNth31g22,JacPDstandardNth22g13))))))))),ToReal(0.5),kmul(kmadd(G122,G212,kmadd(G123,G312,kmul(G223,G322))),ToReal(2))))))))); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED R23 = kmul(ToReal(0.5),kmadd(gu23,kadd(JacPDstandardNth22g33,JacPDstandardNth33g22),kmadd(kmadd(G111,G123,kmadd(kadd(G113,G223),G323,kmul(gu23,JacPDstandardNth23g23))),ToReal(-2),kmadd(gu11,kadd(JacPDstandardNth12g13,ksub(JacPDstandardNth13g12,kadd(JacPDstandardNth23g11,JacPDstandardNth11g23))),kmadd(gu21,kadd(JacPDstandardNth13g22,ksub(JacPDstandardNth22g13,kadd(JacPDstandardNth23g12,JacPDstandardNth21g23))),kmadd(gu13,kadd(JacPDstandardNth12g33,ksub(JacPDstandardNth33g12,kadd(JacPDstandardNth23g13,JacPDstandardNth13g23))),kmadd(gu33,ksub(JacPDstandardNth32g33,JacPDstandardNth23g33),kmadd(gu31,kadd(JacPDstandardNth13g23,ksub(JacPDstandardNth32g13,kadd(JacPDstandardNth31g23,JacPDstandardNth23g13))),kmul(kmadd(G122,G213,kmadd(G133,G312,kmadd(G233,G322,kmul(G112,ksub(G113,G223))))),ToReal(2)))))))))); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED R33 = - kmul(ToReal(0.5),kmadd(gu12,JacPDstandardNth23g13,kmadd(gu21,JacPDstandardNth23g13,kmadd(kmadd(G111,G133,kmadd(G133,G212,kmadd(G112,G233,kmadd(G222,G233,kmadd(G113,G333,kmul(G223,G333)))))),ToReal(-2),kmadd(gu31,ksub(JacPDstandardNth13g33,JacPDstandardNth31g33),kmadd(gu32,ksub(JacPDstandardNth23g33,JacPDstandardNth32g33),knmsub(gu11,kadd(JacPDstandardNth33g11,JacPDstandardNth11g33),kmadd(gu12,ksub(JacPDstandardNth13g23,kadd(JacPDstandardNth33g12,JacPDstandardNth12g33)),kmadd(gu21,ksub(JacPDstandardNth13g23,kadd(JacPDstandardNth33g12,JacPDstandardNth21g33)),kmadd(kmadd(G133,G313,kmadd(G233,G323,kmul(gu11,JacPDstandardNth13g13))),ToReal(2),kmadd(kmul(G113,G113),ToReal(2),kmadd(kmul(G223,G223),ToReal(2),kmadd(gu22,ksub(kmsub(JacPDstandardNth23g23,ToReal(2),JacPDstandardNth33g22),JacPDstandardNth22g33),kmul(G123,kmul(G213,ToReal(4)))))))))))))))); + kmadd(G113,G113,kmadd(G223,G223,kmadd(G333,G333,knmsub(G133,kadd(G111,kadd(G212,G313)),knmsub(G333,kadd(G113,kadd(G223,G333)),kmadd(G233,ksub(G323,kadd(G222,G112)),kmadd(knmsub(gu11,kadd(JacPDstandardNth11g33,kmadd(JacPDstandardNth13g13,ToReal(-2),JacPDstandardNth33g11)),knmsub(gu22,kadd(JacPDstandardNth22g33,kmadd(JacPDstandardNth23g23,ToReal(-2),JacPDstandardNth33g22)),kmadd(gu31,ksub(JacPDstandardNth13g33,JacPDstandardNth31g33),kmadd(gu32,ksub(JacPDstandardNth23g33,JacPDstandardNth32g33),kmadd(gu12,kadd(JacPDstandardNth13g23,ksub(JacPDstandardNth23g13,kadd(JacPDstandardNth33g12,JacPDstandardNth12g33))),kmul(gu21,kadd(JacPDstandardNth13g23,ksub(JacPDstandardNth23g13,kadd(JacPDstandardNth33g12,JacPDstandardNth21g33))))))))),ToReal(0.5),kmul(kmadd(G123,G213,kmul(G133,G313)),ToReal(2))))))))); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED Km11 = kmadd(K11L,gu11,kmadd(K12L,gu12,kmul(K13L,gu13))); @@ -1453,16 +1453,16 @@ static void ML_ADM_RHS_Body(cGH const * restrict const cctkGH, int const dir, in kmadd(beta1L,JacPDstandardNth1g33,kmadd(beta2L,JacPDstandardNth2g33,kmadd(beta3L,JacPDstandardNth3g33,kmadd(alphaL,kmul(K33L,ToReal(-2)),kmul(kmadd(g13L,JacPDstandardNth3beta1,kmadd(g23L,JacPDstandardNth3beta2,kmul(g33L,JacPDstandardNth3beta3))),ToReal(2)))))); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED K11rhsL = - kmadd(G111,JacPDstandardNth1alpha,kmadd(beta1L,JacPDstandardNth1K11,kmadd(G211,JacPDstandardNth2alpha,kmadd(beta2L,JacPDstandardNth2K11,kmadd(G311,JacPDstandardNth3alpha,kmadd(beta3L,JacPDstandardNth3K11,kmadd(alphaL,kadd(R11,kmadd(kmadd(K12L,Km21,kmul(K13L,Km31)),ToReal(-2),kmul(K11L,kmadd(Km11,ToReal(-2),trK)))),kmsub(kmadd(K11L,JacPDstandardNth1beta1,kmadd(K12L,JacPDstandardNth1beta2,kmul(K13L,JacPDstandardNth1beta3))),ToReal(2),JacPDstandardNth11alpha)))))))); + kmadd(G111,JacPDstandardNth1alpha,kmadd(beta1L,JacPDstandardNth1K11,kmadd(G211,JacPDstandardNth2alpha,kmadd(beta2L,JacPDstandardNth2K11,kmadd(G311,JacPDstandardNth3alpha,kmadd(beta3L,JacPDstandardNth3K11,kmadd(alphaL,kadd(R11,kmadd(K11L,trK,kmul(kmadd(K11L,Km11,kmadd(K12L,Km21,kmul(K13L,Km31))),ToReal(-2)))),kmsub(kmadd(K11L,JacPDstandardNth1beta1,kmadd(K12L,JacPDstandardNth1beta2,kmul(K13L,JacPDstandardNth1beta3))),ToReal(2),JacPDstandardNth11alpha)))))))); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED K12rhsL = - kmadd(G112,JacPDstandardNth1alpha,kmadd(K22L,JacPDstandardNth1beta2,kmadd(K23L,JacPDstandardNth1beta3,kmadd(beta1L,JacPDstandardNth1K12,kmadd(G212,JacPDstandardNth2alpha,kmadd(K11L,JacPDstandardNth2beta1,kmadd(K12L,kadd(JacPDstandardNth1beta1,JacPDstandardNth2beta2),kmadd(K13L,JacPDstandardNth2beta3,kmadd(beta2L,JacPDstandardNth2K12,kmadd(G312,JacPDstandardNth3alpha,kmadd(beta3L,JacPDstandardNth3K12,kmsub(alphaL,kadd(R12,kmadd(kmadd(K11L,Km12,kmul(K13L,Km32)),ToReal(-2),kmul(K12L,kmadd(Km22,ToReal(-2),trK)))),JacPDstandardNth12alpha)))))))))))); + kmadd(G112,JacPDstandardNth1alpha,kmadd(K22L,JacPDstandardNth1beta2,kmadd(K23L,JacPDstandardNth1beta3,kmadd(beta1L,JacPDstandardNth1K12,kmadd(G212,JacPDstandardNth2alpha,kmadd(K11L,JacPDstandardNth2beta1,kmadd(K12L,kadd(JacPDstandardNth1beta1,JacPDstandardNth2beta2),kmadd(K13L,JacPDstandardNth2beta3,kmadd(beta2L,JacPDstandardNth2K12,kmadd(G312,JacPDstandardNth3alpha,kmadd(beta3L,JacPDstandardNth3K12,kmsub(alphaL,kadd(R12,kmadd(K12L,trK,kmul(kmadd(K11L,Km12,kmadd(K12L,Km22,kmul(K13L,Km32))),ToReal(-2)))),JacPDstandardNth12alpha)))))))))))); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED K13rhsL = kmadd(G113,JacPDstandardNth1alpha,kmadd(K23L,JacPDstandardNth1beta2,kmadd(K33L,JacPDstandardNth1beta3,kmadd(beta1L,JacPDstandardNth1K13,kmadd(G213,JacPDstandardNth2alpha,kmadd(beta2L,JacPDstandardNth2K13,kmadd(G313,JacPDstandardNth3alpha,kmadd(K11L,JacPDstandardNth3beta1,kmadd(K12L,JacPDstandardNth3beta2,kmadd(K13L,kadd(JacPDstandardNth1beta1,JacPDstandardNth3beta3),kmadd(beta3L,JacPDstandardNth3K13,kmsub(alphaL,kadd(R13,kmadd(K13L,trK,kmul(kmadd(K11L,Km13,kmadd(K12L,Km23,kmul(K13L,Km33))),ToReal(-2)))),JacPDstandardNth13alpha)))))))))))); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED K22rhsL = - kmadd(G122,JacPDstandardNth1alpha,kmadd(beta1L,JacPDstandardNth1K22,kmadd(G222,JacPDstandardNth2alpha,kmadd(beta2L,JacPDstandardNth2K22,kmadd(G322,JacPDstandardNth3alpha,kmadd(beta3L,JacPDstandardNth3K22,kmadd(alphaL,kadd(R22,kmadd(kmadd(K12L,Km12,kmul(K23L,Km32)),ToReal(-2),kmul(K22L,kmadd(Km22,ToReal(-2),trK)))),kmsub(kmadd(K12L,JacPDstandardNth2beta1,kmadd(K22L,JacPDstandardNth2beta2,kmul(K23L,JacPDstandardNth2beta3))),ToReal(2),JacPDstandardNth22alpha)))))))); + kmadd(G122,JacPDstandardNth1alpha,kmadd(beta1L,JacPDstandardNth1K22,kmadd(G222,JacPDstandardNth2alpha,kmadd(beta2L,JacPDstandardNth2K22,kmadd(G322,JacPDstandardNth3alpha,kmadd(beta3L,JacPDstandardNth3K22,kmadd(alphaL,kadd(R22,kmadd(K22L,trK,kmul(kmadd(K12L,Km12,kmadd(K22L,Km22,kmul(K23L,Km32))),ToReal(-2)))),kmsub(kmadd(K12L,JacPDstandardNth2beta1,kmadd(K22L,JacPDstandardNth2beta2,kmul(K23L,JacPDstandardNth2beta3))),ToReal(2),JacPDstandardNth22alpha)))))))); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED K23rhsL = kmadd(G123,JacPDstandardNth1alpha,kmadd(beta1L,JacPDstandardNth1K23,kmadd(G223,JacPDstandardNth2alpha,kmadd(K13L,JacPDstandardNth2beta1,kmadd(K33L,JacPDstandardNth2beta3,kmadd(beta2L,JacPDstandardNth2K23,kmadd(G323,JacPDstandardNth3alpha,kmadd(K12L,JacPDstandardNth3beta1,kmadd(K22L,JacPDstandardNth3beta2,kmadd(K23L,kadd(JacPDstandardNth2beta2,JacPDstandardNth3beta3),kmadd(beta3L,JacPDstandardNth3K23,kmsub(alphaL,kadd(R23,kmadd(K23L,trK,kmul(kmadd(K12L,Km13,kmadd(K22L,Km23,kmul(K23L,Km33))),ToReal(-2)))),JacPDstandardNth23alpha)))))))))))); diff --git a/ML_ADM/src/ML_ADM_constraints.cc b/ML_ADM/src/ML_ADM_constraints.cc index f209776..1ff7b5c 100644 --- a/ML_ADM/src/ML_ADM_constraints.cc +++ b/ML_ADM/src/ML_ADM_constraints.cc @@ -1162,22 +1162,22 @@ static void ML_ADM_constraints_Body(cGH const * restrict const cctkGH, int const kmul(ToReal(0.5),kmadd(gu33,JacPDstandardNth3g33,kmadd(gu31,kmsub(JacPDstandardNth3g13,ToReal(2),JacPDstandardNth1g33),kmul(gu32,kmsub(JacPDstandardNth3g23,ToReal(2),JacPDstandardNth2g33))))); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED R11 = - kmul(ToReal(0.5),knmsub(gu22,kadd(JacPDstandardNth11g22,kmadd(JacPDstandardNth21g12,ToReal(-2),JacPDstandardNth22g11)),kmadd(gu12,ksub(JacPDstandardNth21g11,JacPDstandardNth12g11),kmadd(gu13,ksub(JacPDstandardNth31g11,JacPDstandardNth13g11),kmadd(gu23,kadd(JacPDstandardNth21g13,ksub(JacPDstandardNth31g12,kadd(JacPDstandardNth23g11,JacPDstandardNth11g23))),kmadd(gu32,kadd(JacPDstandardNth21g13,ksub(JacPDstandardNth31g12,kadd(JacPDstandardNth32g11,JacPDstandardNth11g23))),kmadd(ToReal(2),kmadd(gu23,kmadd(g13L,kmul(G112,G313),kmadd(g23L,kmul(G212,G313),kmadd(g33L,kmul(G312,G313),kmadd(g11L,kmsub(G112,G113,kmul(G111,G123)),kmadd(g12L,kmadd(G113,G212,kmsub(G112,G213,kmadd(G111,G223,kmul(G123,G211)))),kmadd(g22L,kmsub(G212,G213,kmul(G211,G223)),knmsub(G311,kmadd(g13L,G123,kmadd(g33L,G323,kmul(g23L,G223))),kmadd(g13L,kmsub(G113,G312,kmul(G111,G323)),kmul(g23L,kmsub(G213,G312,kmul(G211,G323))))))))))),kmadd(gu32,kmadd(g13L,kmul(G112,G313),kmadd(g23L,kmul(G212,G313),kmadd(g33L,kmul(G312,G313),kmadd(g11L,kmsub(G112,G113,kmul(G111,G123)),kmadd(g12L,kmadd(G113,G212,kmsub(G112,G213,kmadd(G111,G223,kmul(G123,G211)))),kmadd(g22L,kmsub(G212,G213,kmul(G211,G223)),knmsub(G311,kmadd(g13L,G123,kmadd(g33L,G323,kmul(g23L,G223))),kmadd(g13L,kmsub(G113,G312,kmul(G111,G323)),kmul(g23L,kmsub(G213,G312,kmul(G211,G323))))))))))),kmul(gu22,kmadd(g22L,kmul(G212,G212),knmsub(g12L,kmadd(G122,G211,kmadd(G111,G222,kmul(G112,kmul(G212,ToReal(-2))))),kmadd(g11L,kmsub(G112,G112,kmul(G111,G122)),knmsub(G222,kmadd(g23L,G311,kmul(g22L,G211)),kmadd(g33L,kmsub(G312,G312,kmul(G311,G322)),kmadd(g13L,knmsub(G122,G311,kmsub(G112,kmul(G312,ToReal(2)),kmul(G111,G322))),kmul(g23L,kmsub(G212,kmul(G312,ToReal(2)),kmul(G211,G322)))))))))))),kmul(gu33,ksub(ksub(kmadd(JacPDstandardNth31g13,ToReal(2),kmul(ToReal(2),kmadd(g22L,kmul(G213,G213),knmsub(g12L,kmadd(G133,G211,kmadd(G111,G233,kmul(G113,kmul(G213,ToReal(-2))))),kmadd(g11L,kmsub(G113,G113,kmul(G111,G133)),knmsub(G233,kmadd(g23L,G311,kmul(g22L,G211)),kmadd(g33L,kmsub(G313,G313,kmul(G311,G333)),kmadd(g13L,knmsub(G133,G311,kmsub(G113,kmul(G313,ToReal(2)),kmul(G111,G333))),kmul(g23L,kmsub(G213,kmul(G313,ToReal(2)),kmul(G211,G333))))))))))),JacPDstandardNth33g11),JacPDstandardNth11g33))))))))); + kmul(ToReal(0.5),kmadd(gu23,JacPDstandardNth31g12,kmadd(gu32,JacPDstandardNth31g12,kmadd(gu12,ksub(JacPDstandardNth21g11,JacPDstandardNth12g11),kmadd(gu13,ksub(JacPDstandardNth31g11,JacPDstandardNth13g11),kmadd(kmadd(g23L,kmul(G212,kmul(G313,gu23)),kmadd(g33L,kmul(G312,kmul(G313,gu23)),kmadd(g22L,kmul(G212,kmul(G213,gu32)),kmadd(g13L,kmul(G113,kmul(G312,gu32)),kmadd(g23L,kmul(G213,kmul(G312,gu32)),kmadd(g13L,kmul(G112,kmul(G313,gu32)),kmadd(g23L,kmul(G212,kmul(G313,gu32)),kmul(gu22,JacPDstandardNth21g12)))))))),ToReal(2),kmadd(g22L,kmul(gu22,kmul(kmul(G212,G212),ToReal(2))),kmadd(g22L,kmul(gu33,kmul(kmul(G213,G213),ToReal(2))),kmadd(g33L,kmul(gu22,kmul(kmul(G312,G312),ToReal(2))),kmadd(g33L,kmul(gu33,kmul(kmul(G313,G313),ToReal(2))),kmadd(g11L,kmul(kmadd(G112,kmul(G113,kadd(gu23,gu32)),kmadd(gu22,kmul(G112,G112),kmsub(gu33,kmul(G113,G113),kmul(G111,kmadd(G122,gu22,kmadd(G123,kadd(gu23,gu32),kmul(G133,gu33))))))),ToReal(2)),kmadd(gu23,kadd(JacPDstandardNth21g13,kmadd(g13L,kmul(G123,kmul(G311,ToReal(-2))),kmadd(g23L,kmul(G223,kmul(G311,ToReal(-2))),kmadd(g13L,kmul(G111,kmul(G323,ToReal(-2))),ksub(ksub(kmadd(g22L,kmul(G212,kmul(G213,ToReal(2))),kmadd(g13L,kmul(G113,kmul(G312,ToReal(2))),kmadd(g23L,kmul(G213,kmul(G312,ToReal(2))),kmul(g13L,kmul(G112,kmul(G313,ToReal(2))))))),JacPDstandardNth23g11),JacPDstandardNth11g23))))),kmadd(gu32,kadd(JacPDstandardNth21g13,kmadd(g13L,kmul(G111,kmul(G323,ToReal(-2))),ksub(kmsub(g33L,kmul(G312,kmul(G313,ToReal(2))),JacPDstandardNth32g11),JacPDstandardNth11g23))),kmadd(gu33,ksub(kmsub(JacPDstandardNth31g13,ToReal(2),JacPDstandardNth33g11),JacPDstandardNth11g33),kmadd(ToReal(-2),kmadd(g22L,kmul(G211,kmul(G222,gu22)),kmadd(g13L,kmul(G122,kmul(G311,gu22)),kmadd(g23L,kmul(G211,kmul(G322,gu22)),kmadd(g33L,kmul(G311,kmul(G322,gu22)),kmadd(g22L,kmul(G211,kmul(G223,gu23)),kmadd(g23L,kmul(G211,kmul(G323,gu23)),kmadd(g33L,kmul(G311,kmul(G323,gu23)),kmadd(g22L,kmul(G211,kmul(G223,gu32)),kmadd(g13L,kmul(G123,kmul(G311,gu32)),kmadd(g23L,kmul(G223,kmul(G311,gu32)),kmadd(g23L,kmul(G211,kmul(G323,gu32)),kmadd(g33L,kmul(G311,kmul(G323,gu32)),kmadd(g22L,kmul(G211,kmul(G233,gu33)),kmadd(g13L,kmul(G133,kmul(G311,gu33)),kmadd(g23L,kmul(G233,kmul(G311,gu33)),kmadd(g13L,kmul(G111,kmul(G333,gu33)),kmadd(g23L,kmul(G211,kmul(G333,gu33)),kmadd(g33L,kmul(G311,kmul(G333,gu33)),kmul(g12L,kmadd(G122,kmul(G211,gu22),kmadd(G111,kmul(G222,gu22),kmadd(G111,kmul(G223,gu23),kmadd(G111,kmul(G223,gu32),kmadd(G133,kmul(G211,gu33),kmadd(G111,kmul(G233,gu33),kmadd(G113,kmul(G213,kmul(gu33,ToReal(-2))),kmadd(gu23,kmsub(G123,G211,kmul(G113,G212)),kmsub(gu32,kmsub(G123,G211,kmul(G113,G212)),kmul(G112,kmadd(G213,kadd(gu23,gu32),kmul(G212,kmul(gu22,ToReal(2))))))))))))))))))))))))))))))))),kmadd(kmadd(g13L,kmul(G113,kmul(G313,gu33)),kmul(g23L,kmul(G213,kmul(G313,gu33)))),ToReal(4),kmul(gu22,kmadd(g23L,kmul(G222,kmul(G311,ToReal(-2))),kmadd(g13L,kmul(G111,kmul(G322,ToReal(-2))),ksub(ksub(kmadd(g13L,kmul(G112,kmul(G312,ToReal(4))),kmul(g23L,kmul(G212,kmul(G312,ToReal(4))))),JacPDstandardNth22g11),JacPDstandardNth11g22)))))))))))))))))))); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED R12 = - kmul(ToReal(0.5),kmadd(gu12,kadd(JacPDstandardNth11g22,kmadd(JacPDstandardNth12g12,ToReal(-2),JacPDstandardNth22g11)),kmadd(gu22,ksub(JacPDstandardNth21g22,JacPDstandardNth12g22),kmadd(gu13,kadd(JacPDstandardNth11g23,ksub(JacPDstandardNth32g11,kadd(JacPDstandardNth13g12,JacPDstandardNth12g13))),kmadd(gu23,kadd(JacPDstandardNth21g23,ksub(JacPDstandardNth32g12,kadd(JacPDstandardNth23g12,JacPDstandardNth12g23))),kmadd(gu32,kadd(JacPDstandardNth22g13,ksub(JacPDstandardNth31g22,kadd(JacPDstandardNth32g12,JacPDstandardNth12g23))),kmadd(gu33,kadd(JacPDstandardNth31g23,ksub(JacPDstandardNth32g13,kadd(JacPDstandardNth33g12,JacPDstandardNth12g33))),kmul(kmadd(gu21,kmadd(G122,kmadd(g11L,G111,kmadd(g12L,G211,kmul(g13L,G311))),kmadd(G222,kmadd(g22L,G211,kmul(g23L,G311)),kmadd(kmadd(g13L,G111,kmadd(g23L,G211,kmul(g33L,G311))),G322,kmadd(kmadd(g13L,kmul(G112,G312),kmul(g23L,kmul(G212,G312))),ToReal(-2),kmsub(g12L,kmadd(G111,G222,kmul(G112,kmul(G212,ToReal(-2)))),kmadd(g11L,kmul(G112,G112),kmadd(g33L,kmul(G312,G312),kmul(g22L,kmul(G212,G212))))))))),kmadd(gu31,kmadd(g13L,kmul(G111,G323),kmadd(g23L,kmul(G211,G323),kmadd(G311,kmadd(g13L,G123,kmadd(g23L,G223,kmul(g33L,G323))),kmadd(g11L,kmsub(G111,G123,kmul(G112,G113)),knmsub(G313,kmadd(g23L,G212,kmul(g13L,G112)),kmadd(g12L,kmadd(G123,G211,kmsub(G111,G223,kmadd(G112,G213,kmul(G113,G212)))),kmsub(g22L,kmsub(G211,G223,kmul(G212,G213)),kmul(G312,kmadd(g13L,G113,kmadd(g33L,G313,kmul(g23L,G213))))))))))),kmadd(gu23,kmadd(g13L,kmul(G113,G322),kmadd(g23L,kmul(G213,G322),kmadd(g33L,kmul(G313,G322),kmadd(g11L,kmsub(G113,G122,kmul(G112,G123)),kmadd(g12L,kmadd(G122,G213,kmsub(G113,G222,kmadd(G112,G223,kmul(G123,G212)))),kmadd(g22L,kmsub(G213,G222,kmul(G212,G223)),knmsub(G312,kmadd(g13L,G123,kmadd(g33L,G323,kmul(g23L,G223))),kmadd(g13L,kmsub(G122,G313,kmul(G112,G323)),kmul(g23L,kmsub(G222,G313,kmul(G212,G323))))))))))),kmul(gu33,kmadd(g13L,kmul(G113,G323),kmadd(g23L,kmul(G213,G323),kmadd(g33L,kmul(G313,G323),kmadd(g11L,kmsub(G113,G123,kmul(G112,G133)),kmadd(g12L,kmadd(G123,G213,kmsub(G113,G223,kmadd(G112,G233,kmul(G133,G212)))),kmadd(g22L,kmsub(G213,G223,kmul(G212,G233)),knmsub(G312,kmadd(g13L,G133,kmadd(g33L,G333,kmul(g23L,G233))),kmadd(g13L,kmsub(G123,G313,kmul(G112,G333)),kmul(g23L,kmsub(G223,G313,kmul(G212,G333))))))))))))))),ToReal(2))))))))); + kmul(ToReal(0.5),kmadd(gu23,JacPDstandardNth21g23,kmadd(gu12,kadd(JacPDstandardNth11g22,JacPDstandardNth22g11),kmadd(gu32,JacPDstandardNth31g22,kmadd(gu33,JacPDstandardNth31g23,kmadd(gu13,JacPDstandardNth32g11,kmadd(gu23,JacPDstandardNth32g12,kmadd(gu33,JacPDstandardNth32g13,kmadd(g33L,kmul(gu21,kmul(kmul(G312,G312),ToReal(-2))),kmadd(gu22,ksub(JacPDstandardNth21g22,JacPDstandardNth12g22),kmadd(gu13,ksub(JacPDstandardNth11g23,kadd(JacPDstandardNth13g12,JacPDstandardNth12g13)),kmadd(gu32,ksub(JacPDstandardNth22g13,kadd(JacPDstandardNth32g12,JacPDstandardNth12g23)),kmadd(ToReal(-2),kmadd(g22L,kmul(G212,kmul(G223,gu23)),kmadd(g13L,kmul(G123,kmul(G312,gu23)),kmadd(g23L,kmul(G212,kmul(G323,gu23)),kmadd(g33L,kmul(G312,kmul(G323,gu23)),kmadd(g22L,kmul(G212,kmul(G213,gu31)),kmadd(g13L,kmul(G113,kmul(G312,gu31)),kmadd(g23L,kmul(G213,kmul(G312,gu31)),kmadd(g13L,kmul(G112,kmul(G313,gu31)),kmadd(g23L,kmul(G212,kmul(G313,gu31)),kmadd(g22L,kmul(G212,kmul(G233,gu33)),kmadd(g13L,kmul(G133,kmul(G312,gu33)),kmadd(g23L,kmul(G233,kmul(G312,gu33)),kmadd(g23L,kmul(G212,kmul(G333,gu33)),kmadd(g33L,kmul(G312,kmul(G333,gu33)),kmadd(gu12,JacPDstandardNth12g12,kmul(g11L,kmadd(G112,kmadd(G123,gu23,kmadd(G113,gu31,kmul(G133,gu33))),kmsub(gu21,kmul(G112,G112),kmadd(G113,kmadd(G122,gu23,kmul(G123,gu33)),kmul(G111,kmadd(G122,gu21,kmul(G123,gu31)))))))))))))))))))))),kmadd(gu21,kmadd(g13L,kmul(G112,kmul(G312,ToReal(-4))),kmadd(g23L,kmul(G212,kmul(G312,ToReal(-4))),kmadd(g22L,kmul(kmul(G212,G212),ToReal(-2)),kmadd(g23L,kmul(G222,kmul(G311,ToReal(2))),kmul(g13L,kmul(G111,kmul(G322,ToReal(2)))))))),kmadd(gu23,kmadd(g23L,kmul(G223,kmul(G312,ToReal(-2))),kmadd(g13L,kmul(G112,kmul(G323,ToReal(-2))),ksub(ksub(kmadd(g13L,kmul(G122,kmul(G313,ToReal(2))),kmadd(g23L,kmul(G222,kmul(G313,ToReal(2))),kmul(g13L,kmul(G113,kmul(G322,ToReal(2)))))),JacPDstandardNth23g12),JacPDstandardNth12g23))),kmadd(gu31,kmadd(g33L,kmul(G312,kmul(G313,ToReal(-2))),kmul(g13L,kmul(G111,kmul(G323,ToReal(2))))),kmadd(gu33,kmadd(g13L,kmul(G112,kmul(G333,ToReal(-2))),ksub(kmsub(g33L,kmul(G313,kmul(G323,ToReal(2))),JacPDstandardNth33g12),JacPDstandardNth12g33)),kmul(ToReal(2),kmadd(g22L,kmul(G211,kmul(G222,gu21)),kmadd(g13L,kmul(G122,kmul(G311,gu21)),kmadd(g23L,kmul(G211,kmul(G322,gu21)),kmadd(g33L,kmul(G311,kmul(G322,gu21)),kmadd(g22L,kmul(G213,kmul(G222,gu23)),kmadd(g23L,kmul(G213,kmul(G322,gu23)),kmadd(g33L,kmul(G313,kmul(G322,gu23)),kmadd(g22L,kmul(G211,kmul(G223,gu31)),kmadd(g13L,kmul(G123,kmul(G311,gu31)),kmadd(g23L,kmul(G223,kmul(G311,gu31)),kmadd(g23L,kmul(G211,kmul(G323,gu31)),kmadd(g33L,kmul(G311,kmul(G323,gu31)),kmadd(g22L,kmul(G213,kmul(G223,gu33)),kmadd(g13L,kmul(G123,kmul(G313,gu33)),kmadd(g23L,kmul(G223,kmul(G313,gu33)),kmadd(g13L,kmul(G113,kmul(G323,gu33)),kmadd(g23L,kmul(G213,kmul(G323,gu33)),kmul(g12L,kmadd(G111,kmul(G222,gu21),kmadd(G122,kmadd(G211,gu21,kmul(G213,gu23)),kmadd(G111,kmul(G223,gu31),kmadd(G113,kmul(G223,gu33),kmadd(gu31,kmsub(G123,G211,kmul(G113,G212)),kmadd(gu23,kmsub(G113,G222,kmul(G123,G212)),kmsub(gu33,kmsub(G123,G213,kmul(G133,G212)),kmul(G112,kmadd(G223,gu23,kmadd(G213,gu31,kmadd(G233,gu33,kmul(G212,kmul(gu21,ToReal(2)))))))))))))))))))))))))))))))))))))))))))))))))); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED R13 = - kmul(ToReal(0.5),kmadd(gu13,kadd(JacPDstandardNth11g33,kmadd(JacPDstandardNth13g13,ToReal(-2),JacPDstandardNth33g11)),kmadd(gu12,kadd(JacPDstandardNth11g23,ksub(JacPDstandardNth23g11,kadd(JacPDstandardNth13g12,JacPDstandardNth12g13))),kmadd(gu33,ksub(JacPDstandardNth31g33,JacPDstandardNth13g33),kmadd(gu22,kadd(JacPDstandardNth21g23,ksub(JacPDstandardNth23g12,kadd(JacPDstandardNth22g13,JacPDstandardNth13g22))),kmadd(gu23,kadd(JacPDstandardNth21g33,ksub(JacPDstandardNth33g12,kadd(JacPDstandardNth23g13,JacPDstandardNth13g23))),kmadd(gu32,kadd(JacPDstandardNth23g13,ksub(JacPDstandardNth31g23,kadd(JacPDstandardNth32g13,JacPDstandardNth13g23))),kmul(kmadd(gu31,kmadd(G133,kmadd(g11L,G111,kmadd(g12L,G211,kmul(g13L,G311))),kmadd(G233,kmadd(g22L,G211,kmul(g23L,G311)),kmadd(kmadd(g13L,G111,kmadd(g23L,G211,kmul(g33L,G311))),G333,kmadd(kmadd(g13L,kmul(G113,G313),kmul(g23L,kmul(G213,G313))),ToReal(-2),kmsub(g12L,kmadd(G111,G233,kmul(G113,kmul(G213,ToReal(-2)))),kmadd(g11L,kmul(G113,G113),kmadd(g33L,kmul(G313,G313),kmul(g22L,kmul(G213,G213))))))))),kmadd(gu21,kmadd(g13L,kmul(G111,G323),kmadd(g23L,kmul(G211,G323),kmadd(G311,kmadd(g13L,G123,kmadd(g23L,G223,kmul(g33L,G323))),kmadd(g11L,kmsub(G111,G123,kmul(G112,G113)),knmsub(G313,kmadd(g23L,G212,kmul(g13L,G112)),kmadd(g12L,kmadd(G123,G211,kmsub(G111,G223,kmadd(G112,G213,kmul(G113,G212)))),kmsub(g22L,kmsub(G211,G223,kmul(G212,G213)),kmul(G312,kmadd(g13L,G113,kmadd(g33L,G313,kmul(g23L,G213))))))))))),kmadd(gu22,kmadd(g13L,kmul(G112,G323),kmadd(g23L,kmul(G212,G323),kmadd(G312,kmadd(g13L,G123,kmadd(g23L,G223,kmul(g33L,G323))),kmadd(g11L,kmsub(G112,G123,kmul(G113,G122)),knmsub(G322,kmadd(g23L,G213,kmul(g13L,G113)),kmadd(g12L,kmadd(G123,G212,kmsub(G112,G223,kmadd(G113,G222,kmul(G122,G213)))),kmsub(g22L,kmsub(G212,G223,kmul(G213,G222)),kmul(G313,kmadd(g13L,G122,kmadd(g33L,G322,kmul(g23L,G222))))))))))),kmul(gu32,kmadd(g13L,kmul(G112,G333),kmadd(g23L,kmul(G212,G333),kmadd(G312,kmadd(g13L,G133,kmadd(g23L,G233,kmul(g33L,G333))),kmadd(g11L,kmsub(G112,G133,kmul(G113,G123)),knmsub(G323,kmadd(g23L,G213,kmul(g13L,G113)),kmadd(g12L,kmadd(G133,G212,kmsub(G112,G233,kmadd(G113,G223,kmul(G123,G213)))),kmsub(g22L,kmsub(G212,G233,kmul(G213,G223)),kmul(G313,kmadd(g13L,G123,kmadd(g33L,G323,kmul(g23L,G223))))))))))))))),ToReal(2))))))))); + kmul(ToReal(0.5),kmadd(gu22,JacPDstandardNth23g12,kmadd(gu32,JacPDstandardNth23g13,kmadd(gu32,JacPDstandardNth31g23,kmadd(gu13,kadd(JacPDstandardNth11g33,JacPDstandardNth33g11),kmadd(gu23,JacPDstandardNth33g12,kmadd(kmadd(g13L,kmul(G113,kmul(G312,gu21)),kmadd(g23L,kmul(G213,kmul(G312,gu21)),kmadd(g13L,kmul(G112,kmul(G313,gu21)),kmadd(g23L,kmul(G212,kmul(G313,gu21)),kmadd(g13L,kmul(G122,kmul(G313,gu22)),kmadd(g23L,kmul(G222,kmul(G313,gu22)),kmadd(g13L,kmul(G113,kmul(G322,gu22)),kmadd(g23L,kmul(G213,kmul(G322,gu22)),kmadd(g33L,kmul(G313,kmul(G322,gu22)),kmadd(g22L,kmul(G213,kmul(G223,gu32)),kmadd(g23L,kmul(G213,kmul(G323,gu32)),kmadd(g33L,kmul(G313,kmul(G323,gu32)),kmadd(gu13,JacPDstandardNth13g13,kmul(g33L,kmul(gu31,kmul(G313,G313)))))))))))))))),ToReal(-2),kmadd(gu12,kadd(JacPDstandardNth11g23,ksub(JacPDstandardNth23g11,kadd(JacPDstandardNth13g12,JacPDstandardNth12g13))),kmadd(gu33,ksub(JacPDstandardNth31g33,JacPDstandardNth13g33),kmadd(gu22,ksub(JacPDstandardNth21g23,kadd(JacPDstandardNth22g13,JacPDstandardNth13g22)),kmadd(gu23,ksub(JacPDstandardNth21g33,kadd(JacPDstandardNth23g13,JacPDstandardNth13g23)),kmadd(gu21,kmadd(g33L,kmul(G312,kmul(G313,ToReal(-2))),kmul(g13L,kmul(G111,kmul(G323,ToReal(2))))),kmadd(gu31,kmadd(g13L,kmul(G113,kmul(G313,ToReal(-4))),kmadd(g23L,kmul(G213,kmul(G313,ToReal(-4))),kmadd(g22L,kmul(kmul(G213,G213),ToReal(-2)),kmadd(g23L,kmul(G233,kmul(G311,ToReal(2))),kmul(g13L,kmul(G111,kmul(G333,ToReal(2)))))))),kmadd(gu32,kmadd(g13L,kmul(G123,kmul(G313,ToReal(-2))),kmadd(g23L,kmul(G223,kmul(G313,ToReal(-2))),kmadd(g13L,kmul(G113,kmul(G323,ToReal(-2))),ksub(ksub(kmadd(g23L,kmul(G233,kmul(G312,ToReal(2))),kmul(g13L,kmul(G112,kmul(G333,ToReal(2))))),JacPDstandardNth32g13),JacPDstandardNth13g23)))),kmadd(g22L,kmadd(G212,kmul(G213,kmul(gu21,ToReal(-2))),kmadd(G213,kmul(G222,kmul(gu22,ToReal(-2))),kmadd(G211,kmul(G223,kmul(gu21,ToReal(2))),kmadd(G211,kmul(G233,kmul(gu31,ToReal(2))),kmul(G212,kmul(G233,kmul(gu32,ToReal(2)))))))),kmul(ToReal(2),kmadd(g13L,kmul(G123,kmul(G311,gu21)),kmadd(g23L,kmul(G223,kmul(G311,gu21)),kmadd(g23L,kmul(G211,kmul(G323,gu21)),kmadd(g33L,kmul(G311,kmul(G323,gu21)),kmadd(g22L,kmul(G212,kmul(G223,gu22)),kmadd(g13L,kmul(G123,kmul(G312,gu22)),kmadd(g23L,kmul(G223,kmul(G312,gu22)),kmadd(g13L,kmul(G112,kmul(G323,gu22)),kmadd(g23L,kmul(G212,kmul(G323,gu22)),kmadd(g33L,kmul(G312,kmul(G323,gu22)),kmadd(g13L,kmul(G133,kmul(G311,gu31)),kmadd(g23L,kmul(G211,kmul(G333,gu31)),kmadd(g33L,kmul(G311,kmul(G333,gu31)),kmadd(g13L,kmul(G133,kmul(G312,gu32)),kmadd(g23L,kmul(G212,kmul(G333,gu32)),kmadd(g33L,kmul(G312,kmul(G333,gu32)),kmadd(g11L,kmadd(G111,kmadd(G123,gu21,kmul(G133,gu31)),kmsub(G112,kmadd(G123,gu22,kmsub(G133,gu32,kmul(G113,gu21))),kmul(G113,kmadd(G122,gu22,kmadd(G113,gu31,kmul(G123,gu32)))))),kmul(g12L,kmadd(G133,kmul(G211,gu31),kmadd(G111,kmul(G233,gu31),kmadd(G133,kmul(G212,gu32),kmadd(G112,kmul(G233,gu32),kmadd(gu21,kmsub(G111,G223,kmul(G112,G213)),kmadd(gu22,kmsub(G112,G223,kmul(G122,G213)),kmsub(G123,kmadd(G211,gu21,kmsub(G212,gu22,kmul(G213,gu32))),kmul(G113,kmadd(G212,gu21,kmadd(G222,gu22,kmadd(G223,gu32,kmul(G213,kmul(gu31,ToReal(2)))))))))))))))))))))))))))))))))))))))))))))))); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED R22 = - kmul(ToReal(0.5),knmsub(gu11,kadd(JacPDstandardNth11g22,kmadd(JacPDstandardNth12g12,ToReal(-2),JacPDstandardNth22g11)),kmadd(gu21,ksub(JacPDstandardNth12g22,JacPDstandardNth21g22),kmadd(gu13,kadd(JacPDstandardNth12g23,ksub(JacPDstandardNth32g12,kadd(JacPDstandardNth22g13,JacPDstandardNth13g22))),kmadd(gu23,ksub(JacPDstandardNth32g22,JacPDstandardNth23g22),kmadd(gu31,kadd(JacPDstandardNth12g23,ksub(JacPDstandardNth32g12,kadd(JacPDstandardNth31g22,JacPDstandardNth22g13))),kmadd(ToReal(2),kmadd(gu13,kmadd(g13L,kmul(G112,G323),kmadd(g23L,kmul(G212,G323),kmadd(g11L,kmsub(G112,G123,kmul(G113,G122)),kmadd(g12L,kmadd(G123,G212,kmsub(G112,G223,kmadd(G113,G222,kmul(G122,G213)))),kmadd(g22L,kmsub(G212,G223,kmul(G213,G222)),kmadd(g13L,kmsub(G123,G312,kmadd(G113,G322,kmul(G122,G313))),kmadd(g23L,kmsub(G223,G312,kmadd(G213,G322,kmul(G222,G313))),kmul(g33L,kmsub(G312,G323,kmul(G313,G322)))))))))),kmadd(gu31,kmadd(g13L,kmul(G112,G323),kmadd(g23L,kmul(G212,G323),kmadd(g11L,kmsub(G112,G123,kmul(G113,G122)),kmadd(g12L,kmadd(G123,G212,kmsub(G112,G223,kmadd(G113,G222,kmul(G122,G213)))),kmadd(g22L,kmsub(G212,G223,kmul(G213,G222)),kmadd(g13L,kmsub(G123,G312,kmadd(G113,G322,kmul(G122,G313))),kmadd(g23L,kmsub(G223,G312,kmadd(G213,G322,kmul(G222,G313))),kmul(g33L,kmsub(G312,G323,kmul(G313,G322)))))))))),kmul(gu11,kmadd(g22L,kmul(G212,G212),knmsub(g12L,kmadd(G122,G211,kmadd(G111,G222,kmul(G112,kmul(G212,ToReal(-2))))),kmadd(g11L,kmsub(G112,G112,kmul(G111,G122)),knmsub(G222,kmadd(g23L,G311,kmul(g22L,G211)),kmadd(g33L,kmsub(G312,G312,kmul(G311,G322)),kmadd(g13L,knmsub(G122,G311,kmsub(G112,kmul(G312,ToReal(2)),kmul(G111,G322))),kmul(g23L,kmsub(G212,kmul(G312,ToReal(2)),kmul(G211,G322)))))))))))),kmul(gu33,ksub(ksub(kmadd(JacPDstandardNth32g23,ToReal(2),kmul(ToReal(2),kmadd(g22L,kmul(G223,G223),knmsub(g12L,kmadd(G133,G222,kmadd(G122,G233,kmul(G123,kmul(G223,ToReal(-2))))),kmadd(g11L,kmsub(G123,G123,kmul(G122,G133)),knmsub(G233,kmadd(g23L,G322,kmul(g22L,G222)),kmadd(g33L,kmsub(G323,G323,kmul(G322,G333)),kmadd(g13L,knmsub(G133,G322,kmsub(G123,kmul(G323,ToReal(2)),kmul(G122,G333))),kmul(g23L,kmsub(G223,kmul(G323,ToReal(2)),kmul(G222,G333))))))))))),JacPDstandardNth33g22),JacPDstandardNth22g33))))))))); + kmul(ToReal(0.5),kmadd(gu31,JacPDstandardNth32g12,kmadd(gu21,ksub(JacPDstandardNth12g22,JacPDstandardNth21g22),kmadd(gu23,ksub(JacPDstandardNth32g22,JacPDstandardNth23g22),knmsub(gu33,kadd(JacPDstandardNth33g22,JacPDstandardNth22g33),kmadd(kmadd(g22L,kmul(G212,kmul(G223,gu13)),kmadd(g13L,kmul(G123,kmul(G312,gu13)),kmadd(g23L,kmul(G212,kmul(G323,gu13)),kmadd(g33L,kmul(G312,kmul(G323,gu13)),kmadd(g22L,kmul(G212,kmul(G223,gu31)),kmadd(g13L,kmul(G123,kmul(G312,gu31)),kmadd(g23L,kmul(G223,kmul(G312,gu31)),kmadd(g23L,kmul(G212,kmul(G323,gu31)),kmadd(g33L,kmul(G312,kmul(G323,gu31)),kmadd(gu11,JacPDstandardNth12g12,kmul(gu33,JacPDstandardNth32g23))))))))))),ToReal(2),kmadd(g22L,kmul(gu33,kmul(kmul(G223,G223),ToReal(2))),kmadd(g33L,kmul(gu11,kmul(kmul(G312,G312),ToReal(2))),kmadd(g33L,kmul(gu33,kmul(kmul(G323,G323),ToReal(2))),kmadd(g11L,kmul(kmadd(G112,kmul(G123,kadd(gu13,gu31)),kmadd(gu11,kmul(G112,G112),kmsub(gu33,kmul(G123,G123),kmul(G122,kmadd(G111,gu11,kmadd(G113,kadd(gu31,gu13),kmul(G133,gu33))))))),ToReal(2)),kmadd(gu31,kadd(JacPDstandardNth12g23,kmadd(g33L,kmul(G313,kmul(G322,ToReal(-2))),ksub(kmsub(g13L,kmul(G112,kmul(G323,ToReal(2))),JacPDstandardNth31g22),JacPDstandardNth22g13))),kmadd(gu13,kadd(JacPDstandardNth12g23,kadd(JacPDstandardNth32g12,kmadd(g13L,kmul(G122,kmul(G313,ToReal(-2))),kmadd(g23L,kmul(G222,kmul(G313,ToReal(-2))),kmadd(g13L,kmul(G113,kmul(G322,ToReal(-2))),ksub(ksub(kmadd(g23L,kmul(G223,kmul(G312,ToReal(2))),kmul(g13L,kmul(G112,kmul(G323,ToReal(2))))),JacPDstandardNth22g13),JacPDstandardNth13g22)))))),kmadd(ToReal(-2),kmadd(g22L,kmul(G211,kmul(G222,gu11)),kmadd(g13L,kmul(G122,kmul(G311,gu11)),kmadd(g23L,kmul(G211,kmul(G322,gu11)),kmadd(g33L,kmul(G311,kmul(G322,gu11)),kmadd(g22L,kmul(G213,kmul(G222,gu13)),kmadd(g23L,kmul(G213,kmul(G322,gu13)),kmadd(g33L,kmul(G313,kmul(G322,gu13)),kmadd(g22L,kmul(G213,kmul(G222,gu31)),kmadd(g13L,kmul(G122,kmul(G313,gu31)),kmadd(g23L,kmul(G222,kmul(G313,gu31)),kmadd(g13L,kmul(G113,kmul(G322,gu31)),kmadd(g23L,kmul(G213,kmul(G322,gu31)),kmadd(g22L,kmul(G222,kmul(G233,gu33)),kmadd(g13L,kmul(G133,kmul(G322,gu33)),kmadd(g23L,kmul(G233,kmul(G322,gu33)),kmadd(g13L,kmul(G122,kmul(G333,gu33)),kmadd(g23L,kmul(G222,kmul(G333,gu33)),kmadd(g33L,kmul(G322,kmul(G333,gu33)),kmul(g12L,kmadd(G111,kmul(G222,gu11),kmadd(G133,kmul(G222,gu33),kmadd(G122,kmadd(G211,gu11,kmadd(G213,kadd(gu13,gu31),kmul(G233,gu33))),kmadd(G123,kmul(G223,kmul(gu33,ToReal(-2))),kmadd(gu13,kmsub(G113,G222,kmul(G123,G212)),kmsub(gu31,kmsub(G113,G222,kmul(G123,G212)),kmul(G112,kmadd(G223,kadd(gu13,gu31),kmul(G212,kmul(gu11,ToReal(2)))))))))))))))))))))))))))))),kmadd(kmadd(g13L,kmul(G123,kmul(G323,gu33)),kmul(g23L,kmul(G223,kmul(G323,gu33)))),ToReal(4),kmul(gu11,kmadd(g23L,kmul(G222,kmul(G311,ToReal(-2))),kmadd(g13L,kmul(G111,kmul(G322,ToReal(-2))),ksub(ksub(kmadd(g22L,kmul(kmul(G212,G212),ToReal(2)),kmadd(g13L,kmul(G112,kmul(G312,ToReal(4))),kmul(g23L,kmul(G212,kmul(G312,ToReal(4)))))),JacPDstandardNth22g11),JacPDstandardNth11g22)))))))))))))))))); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED R23 = - kmul(ToReal(0.5),kmadd(gu23,kadd(JacPDstandardNth22g33,kmadd(JacPDstandardNth23g23,ToReal(-2),JacPDstandardNth33g22)),kmadd(gu11,kadd(JacPDstandardNth12g13,ksub(JacPDstandardNth13g12,kadd(JacPDstandardNth23g11,JacPDstandardNth11g23))),kmadd(gu21,kadd(JacPDstandardNth13g22,ksub(JacPDstandardNth22g13,kadd(JacPDstandardNth23g12,JacPDstandardNth21g23))),kmadd(gu13,kadd(JacPDstandardNth12g33,ksub(JacPDstandardNth33g12,kadd(JacPDstandardNth23g13,JacPDstandardNth13g23))),kmadd(gu33,ksub(JacPDstandardNth32g33,JacPDstandardNth23g33),kmadd(gu31,kadd(JacPDstandardNth13g23,ksub(JacPDstandardNth32g13,kadd(JacPDstandardNth31g23,JacPDstandardNth23g13))),kmul(kmadd(gu32,kmadd(G133,kmadd(g11L,G122,kmadd(g12L,G222,kmul(g13L,G322))),kmadd(G233,kmadd(g22L,G222,kmul(g23L,G322)),kmadd(kmadd(g13L,G122,kmadd(g23L,G222,kmul(g33L,G322))),G333,kmadd(kmadd(g13L,kmul(G123,G323),kmul(g23L,kmul(G223,G323))),ToReal(-2),kmsub(g12L,kmadd(G122,G233,kmul(G123,kmul(G223,ToReal(-2)))),kmadd(g11L,kmul(G123,G123),kmadd(g33L,kmul(G323,G323),kmul(g22L,kmul(G223,G223))))))))),kmadd(gu31,kmadd(g13L,kmul(G112,G333),kmadd(g23L,kmul(G212,G333),kmadd(G312,kmadd(g13L,G133,kmadd(g23L,G233,kmul(g33L,G333))),kmadd(g11L,kmsub(G112,G133,kmul(G113,G123)),knmsub(G323,kmadd(g23L,G213,kmul(g13L,G113)),kmadd(g12L,kmadd(G133,G212,kmsub(G112,G233,kmadd(G113,G223,kmul(G123,G213)))),kmsub(g22L,kmsub(G212,G233,kmul(G213,G223)),kmul(G313,kmadd(g13L,G123,kmadd(g33L,G323,kmul(g23L,G223))))))))))),kmadd(gu11,kmadd(g13L,kmul(G112,G313),kmadd(g23L,kmul(G212,G313),kmadd(g33L,kmul(G312,G313),kmadd(g11L,kmsub(G112,G113,kmul(G111,G123)),kmadd(g12L,kmadd(G113,G212,kmsub(G112,G213,kmadd(G111,G223,kmul(G123,G211)))),kmadd(g22L,kmsub(G212,G213,kmul(G211,G223)),knmsub(G311,kmadd(g13L,G123,kmadd(g33L,G323,kmul(g23L,G223))),kmadd(g13L,kmsub(G113,G312,kmul(G111,G323)),kmul(g23L,kmsub(G213,G312,kmul(G211,G323))))))))))),kmul(gu12,kmadd(g13L,kmul(G113,G322),kmadd(g23L,kmul(G213,G322),kmadd(g33L,kmul(G313,G322),kmadd(g11L,kmsub(G113,G122,kmul(G112,G123)),kmadd(g12L,kmadd(G122,G213,kmsub(G113,G222,kmadd(G112,G223,kmul(G123,G212)))),kmadd(g22L,kmsub(G213,G222,kmul(G212,G223)),knmsub(G312,kmadd(g13L,G123,kmadd(g33L,G323,kmul(g23L,G223))),kmadd(g13L,kmsub(G122,G313,kmul(G112,G323)),kmul(g23L,kmsub(G222,G313,kmul(G212,G323))))))))))))))),ToReal(2))))))))); + kmul(ToReal(0.5),kmadd(gu23,kadd(JacPDstandardNth22g33,JacPDstandardNth33g22),kmadd(kmadd(kmadd(G223,kmadd(g22L,G212,kmul(g23L,G312)),kmadd(g23L,kmul(G212,G323),kmul(G312,kmadd(g13L,G123,kmul(g33L,G323))))),gu12,kmadd(kmadd(g13L,G123,kmul(g23L,G223)),kmadd(G311,gu11,kmul(G313,gu31)),kmadd(G323,kmadd(kmadd(g23L,G211,kmul(g33L,G311)),gu11,kmadd(kmadd(g23L,G213,kmul(g33L,G313)),gu31,kmul(g13L,kmadd(G112,gu12,kmul(G113,gu31))))),kmadd(gu23,JacPDstandardNth23g23,kmul(g33L,kmul(gu32,kmul(G323,G323))))))),ToReal(-2),kmadd(gu21,kadd(JacPDstandardNth13g22,ksub(JacPDstandardNth22g13,kadd(JacPDstandardNth23g12,JacPDstandardNth21g23))),kmadd(gu13,kadd(JacPDstandardNth12g33,ksub(JacPDstandardNth33g12,kadd(JacPDstandardNth23g13,JacPDstandardNth13g23))),kmadd(gu33,ksub(JacPDstandardNth32g33,JacPDstandardNth23g33),kmadd(gu31,kadd(JacPDstandardNth13g23,ksub(JacPDstandardNth32g13,kadd(JacPDstandardNth31g23,JacPDstandardNth23g13))),kmadd(gu11,kadd(JacPDstandardNth12g13,kadd(JacPDstandardNth13g12,kmadd(g13L,kmul(G111,kmul(G323,ToReal(-2))),ksub(kmsub(g33L,kmul(G312,kmul(G313,ToReal(2))),JacPDstandardNth23g11),JacPDstandardNth11g23)))),kmadd(gu32,kmadd(g22L,kmul(kmul(G223,G223),ToReal(-2)),kmadd(g23L,kmadd(G223,kmul(G323,ToReal(-4)),kmul(G233,kmul(G322,ToReal(2)))),kmul(g13L,kmadd(G123,kmul(G323,ToReal(-4)),kmul(G122,kmul(G333,ToReal(2))))))),kmadd(g22L,kmadd(G222,kmul(G233,kmul(gu32,ToReal(2))),kmadd(gu11,kmadd(G211,kmul(G223,ToReal(-2)),kmul(G212,kmul(G213,ToReal(2)))),kmul(G213,kmadd(G223,kmul(gu31,ToReal(-2)),kmul(G222,kmul(gu12,ToReal(2))))))),kmul(ToReal(2),kmadd(kmadd(g23L,kmul(G222,G313),kmul(g13L,kmul(G113,G322))),gu12,kmadd(G313,kmadd(kmadd(g13L,G112,kmul(g23L,G212)),gu11,kmul(g13L,kmul(G122,gu12))),kmadd(G233,kmul(kmadd(g22L,G212,kmul(g23L,G312)),gu31),kmadd(G312,kmadd(kmadd(g13L,G113,kmul(g23L,G213)),gu11,kmul(kmadd(g13L,G133,kmul(g33L,G333)),gu31)),kmadd(G333,kmadd(kmadd(g13L,G112,kmul(g23L,G212)),gu31,kmul(g23L,kmul(G222,gu32))),kmadd(G322,kmadd(kmadd(g23L,G213,kmul(g33L,G313)),gu12,kmul(kmadd(g13L,G133,kmul(g33L,G333)),gu32)),kmadd(g11L,kmadd(G113,kmul(G122,gu12),kmadd(G112,kmadd(G113,gu11,kmsub(G133,gu31,kmul(G123,gu12))),kmsub(gu32,kmsub(G122,G133,kmul(G123,G123)),kmul(G123,kmadd(G113,gu31,kmul(G111,gu11)))))),kmul(g12L,kmadd(kmadd(G133,G222,kmul(G122,G233)),gu32,kmadd(gu11,kmadd(G113,G212,kmsub(G112,G213,kmul(G111,G223))),kmadd(gu12,kmadd(G122,G213,kmsub(G113,G222,kmul(G112,G223))),kmsub(gu31,kmadd(G133,G212,kmsub(G112,G233,kmul(G113,G223))),kmul(G123,kmadd(G211,gu11,kmadd(G212,gu12,kmadd(G213,gu31,kmul(G223,kmul(gu32,ToReal(2)))))))))))))))))))))))))))))); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED R33 = - kmul(ToReal(0.5),knmsub(gu11,kadd(JacPDstandardNth11g33,kmadd(JacPDstandardNth13g13,ToReal(-2),JacPDstandardNth33g11)),kmadd(gu31,ksub(JacPDstandardNth13g33,JacPDstandardNth31g33),kmadd(gu32,ksub(JacPDstandardNth23g33,JacPDstandardNth32g33),kmadd(gu12,kadd(JacPDstandardNth13g23,ksub(JacPDstandardNth23g13,kadd(JacPDstandardNth33g12,JacPDstandardNth12g33))),kmadd(gu21,kadd(JacPDstandardNth13g23,ksub(JacPDstandardNth23g13,kadd(JacPDstandardNth33g12,JacPDstandardNth21g33))),kmadd(ToReal(2),kmadd(gu12,kmadd(g13L,kmul(G113,G323),kmadd(g23L,kmul(G213,G323),kmadd(g33L,kmul(G313,G323),kmadd(g11L,kmsub(G113,G123,kmul(G112,G133)),kmadd(g12L,kmadd(G123,G213,kmsub(G113,G223,kmadd(G112,G233,kmul(G133,G212)))),kmadd(g22L,kmsub(G213,G223,kmul(G212,G233)),knmsub(G312,kmadd(g13L,G133,kmadd(g33L,G333,kmul(g23L,G233))),kmadd(g13L,kmsub(G123,G313,kmul(G112,G333)),kmul(g23L,kmsub(G223,G313,kmul(G212,G333))))))))))),kmadd(gu21,kmadd(g13L,kmul(G113,G323),kmadd(g23L,kmul(G213,G323),kmadd(g33L,kmul(G313,G323),kmadd(g11L,kmsub(G113,G123,kmul(G112,G133)),kmadd(g12L,kmadd(G123,G213,kmsub(G113,G223,kmadd(G112,G233,kmul(G133,G212)))),kmadd(g22L,kmsub(G213,G223,kmul(G212,G233)),knmsub(G312,kmadd(g13L,G133,kmadd(g33L,G333,kmul(g23L,G233))),kmadd(g13L,kmsub(G123,G313,kmul(G112,G333)),kmul(g23L,kmsub(G223,G313,kmul(G212,G333))))))))))),kmul(gu11,kmadd(g22L,kmul(G213,G213),knmsub(g12L,kmadd(G133,G211,kmadd(G111,G233,kmul(G113,kmul(G213,ToReal(-2))))),kmadd(g11L,kmsub(G113,G113,kmul(G111,G133)),knmsub(G233,kmadd(g23L,G311,kmul(g22L,G211)),kmadd(g33L,kmsub(G313,G313,kmul(G311,G333)),kmadd(g13L,knmsub(G133,G311,kmsub(G113,kmul(G313,ToReal(2)),kmul(G111,G333))),kmul(g23L,kmsub(G213,kmul(G313,ToReal(2)),kmul(G211,G333)))))))))))),kmul(gu22,ksub(ksub(kmadd(JacPDstandardNth23g23,ToReal(2),kmul(ToReal(2),kmadd(g22L,kmul(G223,G223),knmsub(g12L,kmadd(G133,G222,kmadd(G122,G233,kmul(G123,kmul(G223,ToReal(-2))))),kmadd(g11L,kmsub(G123,G123,kmul(G122,G133)),knmsub(G233,kmadd(g23L,G322,kmul(g22L,G222)),kmadd(g33L,kmsub(G323,G323,kmul(G322,G333)),kmadd(g13L,knmsub(G133,G322,kmsub(G123,kmul(G323,ToReal(2)),kmul(G122,G333))),kmul(g23L,kmsub(G223,kmul(G323,ToReal(2)),kmul(G222,G333))))))))))),JacPDstandardNth33g22),JacPDstandardNth22g33))))))))); + kmul(ToReal(0.5),kmadd(gu31,ksub(JacPDstandardNth13g33,JacPDstandardNth31g33),kmadd(gu32,ksub(JacPDstandardNth23g33,JacPDstandardNth32g33),knmsub(gu22,kadd(JacPDstandardNth33g22,JacPDstandardNth22g33),kmadd(kmadd(g23L,kmul(G213,kmul(G323,gu12)),kmadd(g33L,kmul(G313,kmul(G323,gu12)),kmadd(g22L,kmul(G213,kmul(G223,gu21)),kmadd(g13L,kmul(G123,kmul(G313,gu21)),kmadd(g23L,kmul(G223,kmul(G313,gu21)),kmadd(g13L,kmul(G113,kmul(G323,gu21)),kmadd(g23L,kmul(G213,kmul(G323,gu21)),kmadd(gu11,JacPDstandardNth13g13,kmul(gu22,JacPDstandardNth23g23))))))))),ToReal(2),kmadd(g22L,kmul(gu11,kmul(kmul(G213,G213),ToReal(2))),kmadd(g22L,kmul(gu22,kmul(kmul(G223,G223),ToReal(2))),kmadd(g33L,kmul(gu11,kmul(kmul(G313,G313),ToReal(2))),kmadd(g33L,kmul(gu22,kmul(kmul(G323,G323),ToReal(2))),kmadd(g11L,kmul(kmadd(G113,kmul(G123,kadd(gu12,gu21)),kmadd(gu11,kmul(G113,G113),kmsub(gu22,kmul(G123,G123),kmul(G133,kmadd(G111,gu11,kmadd(G112,kadd(gu21,gu12),kmul(G122,gu22))))))),ToReal(2)),kmadd(gu12,kadd(JacPDstandardNth13g23,kadd(JacPDstandardNth23g13,kmadd(g13L,kmul(G133,kmul(G312,ToReal(-2))),kmadd(g23L,kmul(G233,kmul(G312,ToReal(-2))),kmadd(g13L,kmul(G112,kmul(G333,ToReal(-2))),ksub(ksub(kmadd(g22L,kmul(G213,kmul(G223,ToReal(2))),kmadd(g13L,kmul(G123,kmul(G313,ToReal(2))),kmadd(g23L,kmul(G223,kmul(G313,ToReal(2))),kmul(g13L,kmul(G113,kmul(G323,ToReal(2))))))),JacPDstandardNth33g12),JacPDstandardNth12g33)))))),kmadd(gu21,kadd(JacPDstandardNth13g23,kadd(JacPDstandardNth23g13,kmadd(g13L,kmul(G112,kmul(G333,ToReal(-2))),ksub(kmsub(g33L,kmul(G313,kmul(G323,ToReal(2))),JacPDstandardNth33g12),JacPDstandardNth21g33)))),kmadd(ToReal(-2),kmadd(g22L,kmul(G211,kmul(G233,gu11)),kmadd(g13L,kmul(G133,kmul(G311,gu11)),kmadd(g23L,kmul(G211,kmul(G333,gu11)),kmadd(g33L,kmul(G311,kmul(G333,gu11)),kmadd(g22L,kmul(G212,kmul(G233,gu12)),kmadd(g23L,kmul(G212,kmul(G333,gu12)),kmadd(g33L,kmul(G312,kmul(G333,gu12)),kmadd(g22L,kmul(G212,kmul(G233,gu21)),kmadd(g13L,kmul(G133,kmul(G312,gu21)),kmadd(g23L,kmul(G233,kmul(G312,gu21)),kmadd(g23L,kmul(G212,kmul(G333,gu21)),kmadd(g33L,kmul(G312,kmul(G333,gu21)),kmadd(g22L,kmul(G222,kmul(G233,gu22)),kmadd(g13L,kmul(G133,kmul(G322,gu22)),kmadd(g23L,kmul(G233,kmul(G322,gu22)),kmadd(g13L,kmul(G122,kmul(G333,gu22)),kmadd(g23L,kmul(G222,kmul(G333,gu22)),kmadd(g33L,kmul(G322,kmul(G333,gu22)),kmul(g12L,kmadd(G111,kmul(G233,gu11),kmadd(G122,kmul(G233,gu22),kmadd(G133,kmadd(G211,gu11,kmadd(G212,kadd(gu12,gu21),kmul(G222,gu22))),kmadd(G123,kmul(G223,kmul(gu22,ToReal(-2))),kmadd(gu12,kmsub(G112,G233,kmul(G123,G213)),kmsub(gu21,kmsub(G112,G233,kmul(G123,G213)),kmul(G113,kmadd(G223,kadd(gu12,gu21),kmul(G213,kmul(gu11,ToReal(2)))))))))))))))))))))))))))))),kmadd(kmadd(g13L,kmul(G123,kmul(G323,gu22)),kmul(g23L,kmul(G223,kmul(G323,gu22)))),ToReal(4),kmul(gu11,kmadd(g23L,kmul(G233,kmul(G311,ToReal(-2))),kmadd(g13L,kmul(G111,kmul(G333,ToReal(-2))),ksub(ksub(kmadd(g13L,kmul(G113,kmul(G313,ToReal(4))),kmul(g23L,kmul(G213,kmul(G313,ToReal(4))))),JacPDstandardNth33g11),JacPDstandardNth11g33)))))))))))))))))); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED trR = kmadd(gu11,R11,kmadd(kadd(gu12,gu21),R12,kmadd(kadd(gu13,gu31),R13,kmadd(gu22,R22,kmadd(kadd(gu23,gu32),R23,kmul(gu33,R33)))))); diff --git a/ML_ADMConstraints/schedule.ccl b/ML_ADMConstraints/schedule.ccl index d0c358e..cfbd37e 100644 --- a/ML_ADMConstraints/schedule.ccl +++ b/ML_ADMConstraints/schedule.ccl @@ -47,32 +47,32 @@ schedule group ML_ADMConstraints_group in MoL_PseudoEvolution after MoL_PostStep schedule ML_ADMConstraints in ML_ADMConstraints_group { LANG: C + READS: ADMBase::kxx + READS: ADMBase::kxy + READS: ADMBase::kxz + READS: ADMBase::kyy + READS: ADMBase::kyz + READS: ADMBase::kzz READS: ADMBase::alp + READS: ADMBase::gxx + READS: ADMBase::gxy + READS: ADMBase::gxz + READS: ADMBase::gyy + READS: ADMBase::gyz + READS: ADMBase::gzz READS: ADMBase::betax READS: ADMBase::betay READS: ADMBase::betaz READS: TmunuBase::eTtt - READS: TmunuBase::eTtx - READS: TmunuBase::eTty - READS: TmunuBase::eTtz READS: TmunuBase::eTxx READS: TmunuBase::eTxy READS: TmunuBase::eTxz READS: TmunuBase::eTyy READS: TmunuBase::eTyz READS: TmunuBase::eTzz - READS: ADMBase::gxx - READS: ADMBase::gxy - READS: ADMBase::gxz - READS: ADMBase::gyy - READS: ADMBase::gyz - READS: ADMBase::gzz - READS: ADMBase::kxx - READS: ADMBase::kxy - READS: ADMBase::kxz - READS: ADMBase::kyy - READS: ADMBase::kyz - READS: ADMBase::kzz + READS: TmunuBase::eTtx + READS: TmunuBase::eTty + READS: TmunuBase::eTtz WRITES: ML_ADMConstraints::H WRITES: ML_ADMConstraints::M1 WRITES: ML_ADMConstraints::M2 diff --git a/ML_ADMConstraints/src/ML_ADMConstraints.cc b/ML_ADMConstraints/src/ML_ADMConstraints.cc index e0ce995..a9f6db1 100644 --- a/ML_ADMConstraints/src/ML_ADMConstraints.cc +++ b/ML_ADMConstraints/src/ML_ADMConstraints.cc @@ -1080,158 +1080,176 @@ static void ML_ADMConstraints_Body(cGH const * restrict const cctkGH, int const gu32*(-JacPDstandardNth2gzz + 2*JacPDstandardNth3gyz) + gu33*JacPDstandardNth3gzz); - CCTK_REAL CCTK_ATTRIBUTE_UNUSED R11 = - 0.5*(gu21*(-JacPDstandardNth12gxx + JacPDstandardNth21gxx) + + CCTK_REAL CCTK_ATTRIBUTE_UNUSED R11 = 0.5*(gxyL*(-4*G123*G211*gu32 + + 4*G113*G212*gu32 + 4*G113*G213*gu33) + gxzL*(-4*G123*G311*gu32 + + 4*G113*G312*gu32 - 2*G133*G311*gu33) + 4*(gyzL*G213*G312*gu32 + + gyzL*G212*G313*gu32 + G112*(gxyL*G212*gu22 + gxzL*G312*gu22 + + gxxL*G113*gu32 + gxyL*G213*gu32 + gxzL*G313*gu32) + + gxzL*G113*G313*gu33) + gyzL*(4*G212*G312*gu22 - 2*G211*G322*gu22 - + 4*G223*G311*gu32 - 2*G233*G311*gu33 + 4*G213*G313*gu33 - + 2*G211*G333*gu33) - 2*(gxyL*G122*G211*gu22 + gyyL*G211*G222*gu22 + + gxzL*G122*G311*gu22 + gyzL*G222*G311*gu22 + gzzL*G311*G322*gu22 + + gxyL*G133*G211*gu33 + gzzL*G311*G333*gu33 + G111*(gxxL*G122*gu22 + + gxyL*G222*gu22 + gxzL*G322*gu22 + 2*gxxL*G123*gu32 + + 2*gxyL*G223*gu32 + 2*gxzL*G323*gu32 + gxxL*G133*gu33 + + gxyL*G233*gu33 + gxzL*G333*gu33)) + gu21*(-JacPDstandardNth12gxx + + JacPDstandardNth21gxx) + gu22*(-JacPDstandardNth11gyy + + 2*JacPDstandardNth21gxy - JacPDstandardNth22gxx) + + gu32*(4*gzzL*G312*G313 - 4*gyzL*G211*G323 - 4*gzzL*G311*G323 - + 2*JacPDstandardNth11gyz - JacPDstandardNth23gxx) + gu31*(-JacPDstandardNth13gxx + JacPDstandardNth31gxx) + - gu32*(4*(gxyL*(-(G123*G211) + G113*G212) + gyyL*(G212*G213 - - G211*G223) + gxzL*G113*G312 + gyzL*G212*G313 + gzzL*G312*G313 + - G112*(gxxL*G113 + gxyL*G213 + gxzL*G313) - G111*(gxxL*G123 + - gxyL*G223 + gxzL*G323) + G311*(-(gxzL*G123) - gyzL*G223 - - gzzL*G323) + gyzL*(G213*G312 - G211*G323)) - JacPDstandardNth11gyz - + JacPDstandardNth21gxz - JacPDstandardNth23gxx + - JacPDstandardNth31gxy) + gu32*(-JacPDstandardNth11gyz + - JacPDstandardNth21gxz + JacPDstandardNth31gxy - JacPDstandardNth32gxx) - + gu22*(-JacPDstandardNth11gyy + 2*JacPDstandardNth21gxy - - JacPDstandardNth22gxx + 2*(2*gyzL*G212*G312 + 2*G112*(gxyL*G212 + - gxzL*G312) - G111*(gxxL*G122 + gxyL*G222 + gxzL*G322) + - G211*(-(gxyL*G122) - gyyL*G222 - gyzL*G322) + - G311*(-(gxzL*G122) - gyzL*G222 - gzzL*G322) + gxxL*SQR(G112) + - gyyL*SQR(G212) + gzzL*SQR(G312))) + gu33*(-JacPDstandardNth11gzz + - 2*JacPDstandardNth31gxz - JacPDstandardNth33gxx + 2*(2*gyzL*G213*G313 - + 2*G113*(gxyL*G213 + gxzL*G313) - G111*(gxxL*G133 + gxyL*G233 - + gxzL*G333) + G211*(-(gxyL*G133) - gyyL*G233 - gyzL*G333) + - G311*(-(gxzL*G133) - gyzL*G233 - gzzL*G333) + gxxL*SQR(G113) + - gyyL*SQR(G213) + gzzL*SQR(G313)))); + 2*gu32*JacPDstandardNth31gxy + gu32*(2*JacPDstandardNth21gxz - + JacPDstandardNth32gxx) + gu33*(-JacPDstandardNth11gzz + + 2*JacPDstandardNth31gxz - JacPDstandardNth33gxx) + + 2*gxxL*gu22*SQR(G112) + 2*gxxL*gu33*SQR(G113) + + gyyL*(4*G212*G213*gu32 - 4*G211*G223*gu32 - 2*G211*G233*gu33 + + 2*gu22*SQR(G212)) + 2*gyyL*gu33*SQR(G213) + 2*gzzL*gu22*SQR(G312) + + 2*gzzL*gu33*SQR(G313)); - CCTK_REAL CCTK_ATTRIBUTE_UNUSED R12 = - 0.5*(gu22*(-JacPDstandardNth12gyy + JacPDstandardNth21gyy) + - gu21*(JacPDstandardNth11gyy - 2*JacPDstandardNth12gxy + - JacPDstandardNth22gxx) + gu32*(-2*JacPDstandardNth12gyz + - JacPDstandardNth21gyz + JacPDstandardNth22gxz - JacPDstandardNth23gxy + - JacPDstandardNth31gyy) + gu31*(JacPDstandardNth11gyz - - JacPDstandardNth12gxz - JacPDstandardNth13gxy + JacPDstandardNth32gxx) - + gu33*(-JacPDstandardNth12gzz + JacPDstandardNth31gyz + - JacPDstandardNth32gxz - JacPDstandardNth33gxy) + 2*((gxyL*(G123*G211 - - G113*G212) + gyyL*(-(G212*G213) + G211*G223) - G112*(gxxL*G113 + - gxyL*G213 + gxzL*G313) + G312*(-(gxzL*G113) - gyzL*G213 - - gzzL*G313) + G111*(gxxL*G123 + gxyL*G223 + gxzL*G323) + - G311*(gxzL*G123 + gyzL*G223 + gzzL*G323) + gyzL*(-(G212*G313) + - G211*G323))*gu31 + (gxyL*(-(G123*G212) + G122*G213) + - gyyL*(G213*G222 - G212*G223) + gxzL*G122*G313 + gyzL*G213*G322 + - gzzL*G313*G322 + G113*(gxxL*G122 + gxyL*G222 + gxzL*G322) - - G112*(gxxL*G123 + gxyL*G223 + gxzL*G323) + G312*(-(gxzL*G123) - - gyzL*G223 - gzzL*G323) + gyzL*(G222*G313 - G212*G323))*gu32 + - (gxyL*(-(G133*G212) + G123*G213) + gyyL*(G213*G223 - G212*G233) + - gxzL*G123*G313 + gyzL*G213*G323 + gzzL*G313*G323 + - G113*(gxxL*G123 + gxyL*G223 + gxzL*G323) - G112*(gxxL*G133 + - gxyL*G233 + gxzL*G333) + G312*(-(gxzL*G133) - gyzL*G233 - - gzzL*G333) + gyzL*(G223*G313 - G212*G333))*gu33 + - gu21*(-2*(gyzL*G212*G312 + G112*(gxyL*G212 + gxzL*G312)) + - G111*(gxxL*G122 + gxyL*G222 + gxzL*G322) + G211*(gxyL*G122 + - gyyL*G222 + gyzL*G322) + G311*(gxzL*G122 + gyzL*G222 + - gzzL*G322) - gxxL*SQR(G112) - gyyL*SQR(G212) - - gzzL*SQR(G312)))); + CCTK_REAL CCTK_ATTRIBUTE_UNUSED R12 = 0.5*(gxzL*((-2*G123*G312 + + 2*G122*G313)*gu32 - 2*G133*G312*gu33) + gyzL*((-4*G212*G312 + + 2*G211*G322)*gu21 + G223*(2*G311*gu31 - 2*G312*gu32) - + 2*G233*G312*gu33) + 2*(G111*((gxxL*G122 + gxyL*G222 + + gxzL*G322)*gu21 + (gxxL*G123 + gxyL*G223 + gxzL*G323)*gu31) + + G211*((gxyL*G122 + gyyL*G222)*gu21 + (gxyL*G123 + gyyL*G223 + + gyzL*G323)*gu31) + G311*((gxzL*G122 + gyzL*G222 + + gzzL*G322)*gu21 + (gxzL*G123 + gzzL*G323)*gu31) + + (G122*(gxxL*G113 + gxyL*G213) + G222*(gxyL*G113 + gyzL*G313) + + (gxzL*G113 + gyzL*G213 + gzzL*G313)*G322)*gu32 + + (G123*(gxxL*G113 + gxyL*G213) + gxyL*G113*G223 + (gyzL*G213 + + gzzL*G313)*G323)*gu33) + gu22*(-JacPDstandardNth12gyy + + JacPDstandardNth21gyy) + gu32*(JacPDstandardNth21gyz + + JacPDstandardNth22gxz - JacPDstandardNth23gxy + JacPDstandardNth31gyy) + + gu31*(JacPDstandardNth11gyz - JacPDstandardNth12gxz - + JacPDstandardNth13gxy + JacPDstandardNth32gxx) + + gu33*(2*(gxzL*G123*G313 + gyzL*G223*G313 + gxzL*G113*G323) - + 2*(gyyL*G212*G233 + gyzL*G212*G333 + gzzL*G312*G333) - + JacPDstandardNth12gzz + JacPDstandardNth31gyz + JacPDstandardNth32gxz - + JacPDstandardNth33gxy) - 2*((gyzL*G212 + gzzL*G312)*G313*gu31 + + G312*((gxzL*G113 + gyzL*G213)*gu31 + gzzL*G323*gu32) + + G212*((gxyL*G113 + gyyL*G213)*gu31 + (gxyL*G123 + + gyzL*G323)*gu32 + gxyL*G133*gu33) + G112*(2*(gxyL*G212 + + gxzL*G312)*gu21 + (gxxL*G113 + gxyL*G213 + gxzL*G313)*gu31 + + (gxxL*G123 + gxyL*G223 + gxzL*G323)*gu32 + (gxxL*G133 + + gxyL*G233 + gxzL*G333)*gu33) + gu32*JacPDstandardNth12gyz + + gu21*(JacPDstandardNth12gxy + gxxL*SQR(G112))) + gyyL*((2*G213*G222 + - 2*G212*G223)*gu32 + 2*G213*G223*gu33 - 2*gu21*SQR(G212)) + + gu21*(JacPDstandardNth11gyy + JacPDstandardNth22gxx - + 2*gzzL*SQR(G312))); - CCTK_REAL CCTK_ATTRIBUTE_UNUSED R13 = 0.5*(gu21*(JacPDstandardNth11gyz - - JacPDstandardNth12gxz - JacPDstandardNth13gxy + - JacPDstandardNth23gxx) + gu22*(-JacPDstandardNth13gyy + - JacPDstandardNth21gyz - JacPDstandardNth22gxz + JacPDstandardNth23gxy) - + gu33*(-JacPDstandardNth13gzz + JacPDstandardNth31gzz) + - gu31*(JacPDstandardNth11gzz - 2*JacPDstandardNth13gxz + - JacPDstandardNth33gxx) + gu32*(-2*JacPDstandardNth13gyz + - JacPDstandardNth21gzz + JacPDstandardNth31gyz - JacPDstandardNth32gxz + - JacPDstandardNth33gxy) + 2*((gxyL*(G123*G211 - G113*G212) + - gyyL*(-(G212*G213) + G211*G223) - G112*(gxxL*G113 + gxyL*G213 + - gxzL*G313) + G312*(-(gxzL*G113) - gyzL*G213 - gzzL*G313) + - G111*(gxxL*G123 + gxyL*G223 + gxzL*G323) + G311*(gxzL*G123 + - gyzL*G223 + gzzL*G323) + gyzL*(-(G212*G313) + G211*G323))*gu21 + - (gxyL*(G123*G212 - G122*G213) + gyyL*(-(G213*G222) + G212*G223) - - G113*(gxxL*G122 + gxyL*G222 + gxzL*G322) + G313*(-(gxzL*G122) - - gyzL*G222 - gzzL*G322) + G112*(gxxL*G123 + gxyL*G223 + - gxzL*G323) + G312*(gxzL*G123 + gyzL*G223 + gzzL*G323) + - gyzL*(-(G213*G322) + G212*G323))*gu22 + (gxyL*(G133*G212 - - G123*G213) + gyyL*(-(G213*G223) + G212*G233) - G113*(gxxL*G123 + - gxyL*G223 + gxzL*G323) + G313*(-(gxzL*G123) - gyzL*G223 - - gzzL*G323) + G112*(gxxL*G133 + gxyL*G233 + gxzL*G333) + - G312*(gxzL*G133 + gyzL*G233 + gzzL*G333) + gyzL*(-(G213*G323) + - G212*G333))*gu32 + gu31*(-2*(gyzL*G213*G313 + G113*(gxyL*G213 + - gxzL*G313)) + G111*(gxxL*G133 + gxyL*G233 + gxzL*G333) + - G211*(gxyL*G133 + gyyL*G233 + gyzL*G333) + G311*(gxzL*G133 + - gyzL*G233 + gzzL*G333) - gxxL*SQR(G113) - gyyL*SQR(G213) - - gzzL*SQR(G313)))); + CCTK_REAL CCTK_ATTRIBUTE_UNUSED R13 = 0.5*(gxyL*(2*G123*G211*gu21 - + 2*G113*G212*gu21 + 2*G123*G212*gu22 + 2*G133*G211*gu31 + + 2*G133*G212*gu32) + 2*(gxzL*G123*G311*gu21 + gyzL*G223*G311*gu21 + + gyzL*G211*G323*gu21 + gzzL*G311*G323*gu21 + gyyL*G212*G223*gu22 + + gxzL*G123*G312*gu22 + gyzL*G223*G312*gu22 + gxzL*G133*G311*gu31 + + G111*(gxxL*G123*gu21 + gxyL*G223*gu21 + gxzL*G323*gu21 + + gxxL*G133*gu31 + gxyL*G233*gu31 + gxzL*G333*gu31) + + gyyL*G212*G233*gu32 + gxzL*G133*G312*gu32 + gyzL*G233*G312*gu32 + + G112*((-(gxxL*G113) - gxyL*G213 - gxzL*G313)*gu21 + + gxxL*G123*gu22 + gxyL*G223*gu22 + gxzL*G323*gu22 + + gxxL*G133*gu32 + gxyL*G233*gu32 + gxzL*G333*gu32)) + + gu22*(-2*gzzL*G313*G322 + 2*gyzL*G212*G323 + 2*gzzL*G312*G323 - + JacPDstandardNth13gyy + JacPDstandardNth21gyz - JacPDstandardNth22gxz) + + gu21*(JacPDstandardNth11gyz - JacPDstandardNth12gxz - + JacPDstandardNth13gxy + JacPDstandardNth23gxx) + + gu22*JacPDstandardNth23gxy + gu32*JacPDstandardNth31gyz + + gu33*(-JacPDstandardNth13gzz + JacPDstandardNth31gzz) + + gu32*(-2*gzzL*G313*G323 + 2*gyzL*G212*G333 + 2*gzzL*G312*G333 - + 2*JacPDstandardNth13gyz + JacPDstandardNth21gzz - + JacPDstandardNth32gxz) + gu31*(-4*gxyL*G113*G213 + 2*gyzL*G233*G311 + - 4*gxzL*G113*G313 - 4*gyzL*G213*G313 + 2*gyzL*G211*G333 + + 2*gzzL*G311*G333 + JacPDstandardNth11gzz + JacPDstandardNth33gxx) + + gu32*JacPDstandardNth33gxy - 2*(gxzL*G113*G312*gu21 + + gyzL*G213*G312*gu21 + gyzL*G212*G313*gu21 + gzzL*G312*G313*gu21 + + gxxL*G113*G122*gu22 + gxyL*G122*G213*gu22 + gxyL*G113*G222*gu22 + + gxzL*G122*G313*gu22 + gyzL*G222*G313*gu22 + gxzL*G113*G322*gu22 + + gyzL*G213*G322*gu22 + gxxL*G113*G123*gu32 + gxyL*G123*G213*gu32 + + gxyL*G113*G223*gu32 + gxzL*G123*G313*gu32 + gyzL*G223*G313*gu32 + + gxzL*G113*G323*gu32 + gyzL*G213*G323*gu32 + + gu31*JacPDstandardNth13gxz + gxxL*gu31*SQR(G113)) + + gyyL*(-2*G212*G213*gu21 + 2*G211*G223*gu21 - 2*G213*G222*gu22 + + 2*G211*G233*gu31 - 2*G213*G223*gu32 - 2*gu31*SQR(G213)) - + 2*gzzL*gu31*SQR(G313)); - CCTK_REAL CCTK_ATTRIBUTE_UNUSED R22 = 0.5*(gu21*(JacPDstandardNth12gyy - - JacPDstandardNth21gyy) + gu31*(4*(gxyL*(G123*G212 - G122*G213) + - gyyL*(-(G213*G222) + G212*G223) + gxzL*(G123*G312 - G122*G313) - - G113*(gxxL*G122 + gxyL*G222 + gxzL*G322) + gyzL*(G223*G312 - - G222*G313 - G213*G322) + gyzL*G212*G323 + G112*(gxxL*G123 + - gxyL*G223 + gxzL*G323) + gzzL*(-(G313*G322) + G312*G323)) + - JacPDstandardNth12gyz - JacPDstandardNth13gyy - JacPDstandardNth22gxz + - JacPDstandardNth32gxy) + gu31*(JacPDstandardNth12gyz - - JacPDstandardNth22gxz - JacPDstandardNth31gyy + JacPDstandardNth32gxy) - + gu32*(-JacPDstandardNth23gyy + JacPDstandardNth32gyy) + + CCTK_REAL CCTK_ATTRIBUTE_UNUSED R22 = 0.5*(-4*(gyzL*G213*G322*gu31 + + gzzL*G313*G322*gu31) + gyzL*(4*G212*G312*gu11 - 2*G211*G322*gu11 + + 4*G223*G312*gu31 - 2*G233*G322*gu33) + 4*(G112*(gxyL*G212*gu11 + + gxzL*G312*gu11 + gxxL*G123*gu31 + gxyL*G223*gu31 + + gxzL*G323*gu31) + gxyL*G123*G223*gu33 + gxzL*G123*G323*gu33 + + gyzL*G223*G323*gu33) + gu21*(JacPDstandardNth12gyy - + JacPDstandardNth21gyy) - 2*(gxyL*G122*G211*gu11 + + gyyL*G211*G222*gu11 + gxzL*G122*G311*gu11 + gyzL*G222*G311*gu11 + + gzzL*G311*G322*gu11 + G111*(gxxL*G122 + gxyL*G222 + + gxzL*G322)*gu11 + gxxL*G122*G133*gu33 + gxyL*G133*G222*gu33 + + gxyL*G122*G233*gu33 + gyyL*G222*G233*gu33 + gxzL*G133*G322*gu33 + + gxzL*G122*G333*gu33 + gyzL*G222*G333*gu33 + gzzL*G322*G333*gu33 + + gu31*JacPDstandardNth22gxz) + gu31*(-4*gxxL*G113*G122 + + 4*gxyL*G123*G212 - 4*gxyL*G122*G213 - 4*gxyL*G113*G222 - + 4*gyyL*G213*G222 + 4*gyyL*G212*G223 + 4*gxzL*G123*G312 - + 4*gxzL*G122*G313 - 4*gyzL*G222*G313 - 4*gxzL*G113*G322 + + 4*gyzL*G212*G323 + 4*gzzL*G312*G323 + 2*JacPDstandardNth12gyz - + JacPDstandardNth13gyy - JacPDstandardNth31gyy) + + 2*gu31*JacPDstandardNth32gxy + gu32*(-JacPDstandardNth23gyy + + JacPDstandardNth32gyy) + gu33*(-JacPDstandardNth22gzz + + 2*JacPDstandardNth32gyz - JacPDstandardNth33gyy) + gu11*(-JacPDstandardNth11gyy + 2*JacPDstandardNth12gxy - - JacPDstandardNth22gxx + 2*(2*gyzL*G212*G312 + 2*G112*(gxyL*G212 + - gxzL*G312) - G111*(gxxL*G122 + gxyL*G222 + gxzL*G322) + - G211*(-(gxyL*G122) - gyyL*G222 - gyzL*G322) + - G311*(-(gxzL*G122) - gyzL*G222 - gzzL*G322) + gxxL*SQR(G112) + - gyyL*SQR(G212) + gzzL*SQR(G312))) + gu33*(-JacPDstandardNth22gzz + - 2*JacPDstandardNth32gyz - JacPDstandardNth33gyy + 2*(2*gyzL*G223*G323 - + 2*G123*(gxyL*G223 + gxzL*G323) - G122*(gxxL*G133 + gxyL*G233 - + gxzL*G333) + G222*(-(gxyL*G133) - gyyL*G233 - gyzL*G333) + - G322*(-(gxzL*G133) - gyzL*G233 - gzzL*G333) + gxxL*SQR(G123) + - gyyL*SQR(G223) + gzzL*SQR(G323)))); + JacPDstandardNth22gxx + 2*gxxL*SQR(G112)) + 2*gxxL*gu33*SQR(G123) + + 2*gyyL*gu11*SQR(G212) + 2*gyyL*gu33*SQR(G223) + + 2*gzzL*gu11*SQR(G312) + 2*gzzL*gu33*SQR(G323)); - CCTK_REAL CCTK_ATTRIBUTE_UNUSED R23 = - 0.5*(gu11*(-JacPDstandardNth11gyz + JacPDstandardNth12gxz + + CCTK_REAL CCTK_ATTRIBUTE_UNUSED R23 = 0.5*(gxyL*((-2*G123*G211 + + 2*G113*G212)*gu11 + G212*(-2*G123*gu21 + 2*G133*gu31) + (-4*G123*G223 + + 2*(G133*G222 + G122*G233))*gu32) + 2*(G313*((gyzL*G212 + + gzzL*G312)*gu11 + gxzL*G122*gu21) + G312*((gxzL*G113 + + gyzL*G213)*gu11 + gxzL*G133*gu31) + G112*((gxxL*G113 + + gxyL*G213 + gxzL*G313)*gu11 + (-(gxxL*G123) - gxyL*G223 - + gxzL*G323)*gu21 + gxxL*G133*gu31 + gxyL*G233*gu31 + + gxzL*G333*gu31) + (gxzL*G122 + gzzL*G322)*G333*gu32 + + G122*((gxxL*G113 + gxyL*G213)*gu21 + gxxL*G133*gu32) + + G322*((gxzL*G113 + gyzL*G213)*gu21 + gxzL*G133*gu32) + + G233*((gyyL*G212 + gyzL*G312)*gu31 + gyzL*G322*gu32) + + G222*((gxyL*G113 + gyzL*G313)*gu21 + gyzL*G333*gu32)) + + gu11*(-JacPDstandardNth11gyz + JacPDstandardNth12gxz + JacPDstandardNth13gxy - JacPDstandardNth23gxx) + - gu21*(JacPDstandardNth13gyy - JacPDstandardNth21gyz + - JacPDstandardNth22gxz - JacPDstandardNth23gxy) + - gu33*(-JacPDstandardNth23gzz + JacPDstandardNth32gzz) + - gu31*(JacPDstandardNth12gzz - 2*JacPDstandardNth23gxz - + gu21*(-2*gyzL*G212*G323 + gzzL*(2*G313*G322 - 2*G312*G323) + + JacPDstandardNth13gyy - JacPDstandardNth21gyz + JacPDstandardNth22gxz - + JacPDstandardNth23gxy) + gu33*(-JacPDstandardNth23gzz + + JacPDstandardNth32gzz) + gu31*(2*gyzL*G212*G333 + + gzzL*(-2*G313*G323 + 2*G312*G333) + JacPDstandardNth12gzz - JacPDstandardNth31gyz + JacPDstandardNth32gxz + JacPDstandardNth33gxy) - + gu32*(JacPDstandardNth22gzz - 2*JacPDstandardNth23gyz + - JacPDstandardNth33gyy) + 2*((gxyL*(-(G123*G211) + G113*G212) + - gyyL*(G212*G213 - G211*G223) + gxzL*G113*G312 + gyzL*G212*G313 + - gzzL*G312*G313 + G112*(gxxL*G113 + gxyL*G213 + gxzL*G313) - - G111*(gxxL*G123 + gxyL*G223 + gxzL*G323) + G311*(-(gxzL*G123) - - gyzL*G223 - gzzL*G323) + gyzL*(G213*G312 - G211*G323))*gu11 + - (gxyL*(-(G123*G212) + G122*G213) + gyyL*(G213*G222 - G212*G223) + - gxzL*G122*G313 + gyzL*G213*G322 + gzzL*G313*G322 + - G113*(gxxL*G122 + gxyL*G222 + gxzL*G322) - G112*(gxxL*G123 + - gxyL*G223 + gxzL*G323) + G312*(-(gxzL*G123) - gyzL*G223 - - gzzL*G323) + gyzL*(G222*G313 - G212*G323))*gu21 + - (gxyL*(G133*G212 - G123*G213) + gyyL*(-(G213*G223) + G212*G233) - - G113*(gxxL*G123 + gxyL*G223 + gxzL*G323) + G313*(-(gxzL*G123) - - gyzL*G223 - gzzL*G323) + G112*(gxxL*G133 + gxyL*G233 + - gxzL*G333) + G312*(gxzL*G133 + gyzL*G233 + gzzL*G333) + - gyzL*(-(G213*G323) + G212*G333))*gu31 + gu32*(-2*(gyzL*G223*G323 + - G123*(gxyL*G223 + gxzL*G323)) + G122*(gxxL*G133 + gxyL*G233 + - gxzL*G333) + G222*(gxyL*G133 + gyyL*G233 + gyzL*G333) + - G322*(gxzL*G133 + gyzL*G233 + gzzL*G333) - gxxL*SQR(G123) - - gyyL*SQR(G223) - gzzL*SQR(G323)))); + - 2*(((gxzL*G123 + gyzL*G223)*G311 + G111*(gxxL*G123 + + gxyL*G223 + gxzL*G323))*gu11 + (gxzL*G123*G312 + + G223*(gyyL*G212 + gyzL*G312))*gu21 + G323*((gyzL*G211 + + gzzL*G311)*gu11 + gxzL*G113*gu31) + gu31*(G123*(gxxL*G113 + + gxyL*G213 + gxzL*G313) + G223*(gxyL*G113 + gyzL*G313) + + gyzL*G213*G323 + JacPDstandardNth23gxz) + gxxL*gu32*SQR(G123)) + + gyyL*((2*G212*G213 - 2*G211*G223)*gu11 + G213*(2*G222*gu21 - + 2*G223*gu31) + gu32*(2*G222*G233 - 2*SQR(G223))) + + gu32*(-4*(gxzL*G123 + gyzL*G223)*G323 + JacPDstandardNth22gzz + + JacPDstandardNth33gyy - 2*(JacPDstandardNth23gyz + gzzL*SQR(G323)))); - CCTK_REAL CCTK_ATTRIBUTE_UNUSED R33 = 0.5*(gu31*(JacPDstandardNth13gzz - - JacPDstandardNth31gzz) + gu32*(JacPDstandardNth23gzz - - JacPDstandardNth32gzz) + gu21*(4*(gxyL*(-(G133*G212) + G123*G213) + - gyyL*(G213*G223 - G212*G233) + gxzL*G123*G313 + gyzL*G213*G323 + - gzzL*G313*G323 + G113*(gxxL*G123 + gxyL*G223 + gxzL*G323) - - G112*(gxxL*G133 + gxyL*G233 + gxzL*G333) + G312*(-(gxzL*G133) - - gyzL*G233 - gzzL*G333) + gyzL*(G223*G313 - G212*G333)) - - JacPDstandardNth12gzz + JacPDstandardNth13gyz + JacPDstandardNth23gxz - - JacPDstandardNth33gxy) + gu21*(JacPDstandardNth13gyz - - JacPDstandardNth21gzz + JacPDstandardNth23gxz - JacPDstandardNth33gxy) - + gu11*(-JacPDstandardNth11gzz + 2*JacPDstandardNth13gxz - - JacPDstandardNth33gxx + 2*(2*gyzL*G213*G313 + 2*G113*(gxyL*G213 + - gxzL*G313) - G111*(gxxL*G133 + gxyL*G233 + gxzL*G333) + - G211*(-(gxyL*G133) - gyyL*G233 - gyzL*G333) + - G311*(-(gxzL*G133) - gyzL*G233 - gzzL*G333) + gxxL*SQR(G113) + - gyyL*SQR(G213) + gzzL*SQR(G313))) + gu22*(-JacPDstandardNth22gzz + - 2*JacPDstandardNth23gyz - JacPDstandardNth33gyy + 2*(2*gyzL*G223*G323 - + 2*G123*(gxyL*G223 + gxzL*G323) - G122*(gxxL*G133 + gxyL*G233 - + gxzL*G333) + G222*(-(gxyL*G133) - gyyL*G233 - gyzL*G333) + - G322*(-(gxzL*G133) - gyzL*G233 - gzzL*G333) + gxxL*SQR(G123) + - gyyL*SQR(G223) + gzzL*SQR(G323)))); + CCTK_REAL CCTK_ATTRIBUTE_UNUSED R33 = 0.5*(4*((G213*(gxyL*G123 + + gyyL*G223) + gyzL*(G223*G313 + G213*G323))*gu21 + + G113*((gxyL*G213 + gxzL*G313)*gu11 + (gxxL*G123 + gxyL*G223 + + gxzL*G323)*gu21) + gxyL*G123*G223*gu22) + gxzL*((-4*G133*G312 + + 4*G123*G313)*gu21 - 2*G133*G322*gu22) + gyzL*((4*G213*G313 - + 2*G211*G333)*gu11 + G233*(-4*G312*gu21 - 2*G322*gu22)) + + gu21*(4*gzzL*G313*G323 - 4*(G133*(gxxL*G112 + gxyL*G212) + + (gxyL*G112 + gyyL*G212)*G233 + gxzL*G112*G333 + (gyzL*G212 + + gzzL*G312)*G333) - JacPDstandardNth12gzz + 2*JacPDstandardNth13gyz - + JacPDstandardNth21gzz + 2*JacPDstandardNth23gxz) + + gu31*(JacPDstandardNth13gzz - JacPDstandardNth31gzz) + + gu32*(JacPDstandardNth23gzz - JacPDstandardNth32gzz) - + 2*((G211*(gxyL*G133 + gyyL*G233) + G111*(gxxL*G133 + gxyL*G233 + + gxzL*G333) + G311*(gxzL*G133 + gyzL*G233 + gzzL*G333))*gu11 + + (G133*(gxxL*G122 + gxyL*G222) + gxyL*G122*G233 + (gyzL*G222 + + gzzL*G322)*G333)*gu22 + gu21*JacPDstandardNth33gxy) + + gu11*(-JacPDstandardNth11gzz + 2*JacPDstandardNth13gxz - + JacPDstandardNth33gxx + 2*gxxL*SQR(G113) + 2*gyyL*SQR(G213) + + 2*gzzL*SQR(G313)) + gu22*(4*(gxzL*G123*G323 + gyzL*G223*G323) - + 2*(gyyL*G222*G233 + gxzL*G122*G333) - JacPDstandardNth22gzz + + 2*JacPDstandardNth23gyz - JacPDstandardNth33gyy + 2*gxxL*SQR(G123) + + 2*gyyL*SQR(G223) + 2*gzzL*SQR(G323))); CCTK_REAL CCTK_ATTRIBUTE_UNUSED trR = gu11*R11 + gu22*R22 + 2*(gu21*R12 + gu31*R13 + gu32*R23) + gu33*R33; diff --git a/ML_ADMQuantities/schedule.ccl b/ML_ADMQuantities/schedule.ccl index 5253582..8adebf2 100644 --- a/ML_ADMQuantities/schedule.ccl +++ b/ML_ADMQuantities/schedule.ccl @@ -47,40 +47,40 @@ schedule group ML_ADMQuantities_group in MoL_PseudoEvolution after MoL_PostStep schedule ML_ADMQuantities in ML_ADMQuantities_group { LANG: C - READS: ML_BSSN::alpha + READS: grid::x + READS: grid::y + READS: grid::z READS: ML_BSSN::At11 READS: ML_BSSN::At12 READS: ML_BSSN::At13 READS: ML_BSSN::At22 READS: ML_BSSN::At23 READS: ML_BSSN::At33 + READS: ML_BSSN::Xt1 + READS: ML_BSSN::Xt2 + READS: ML_BSSN::Xt3 + READS: ML_BSSN::alpha + READS: ML_BSSN::phi + READS: ML_BSSN::gt11 + READS: ML_BSSN::gt12 + READS: ML_BSSN::gt13 + READS: ML_BSSN::gt22 + READS: ML_BSSN::gt23 + READS: ML_BSSN::gt33 READS: ML_BSSN::beta1 READS: ML_BSSN::beta2 READS: ML_BSSN::beta3 + READS: ML_BSSN::trK READS: TmunuBase::eTtt - READS: TmunuBase::eTtx - READS: TmunuBase::eTty - READS: TmunuBase::eTtz READS: TmunuBase::eTxx READS: TmunuBase::eTxy READS: TmunuBase::eTxz READS: TmunuBase::eTyy READS: TmunuBase::eTyz READS: TmunuBase::eTzz - READS: ML_BSSN::gt11 - READS: ML_BSSN::gt12 - READS: ML_BSSN::gt13 - READS: ML_BSSN::gt22 - READS: ML_BSSN::gt23 - READS: ML_BSSN::gt33 - READS: ML_BSSN::phi - READS: ML_BSSN::trK - READS: grid::x - READS: ML_BSSN::Xt1 - READS: ML_BSSN::Xt2 - READS: ML_BSSN::Xt3 - READS: grid::y - READS: grid::z + READS: TmunuBase::eTtx + READS: TmunuBase::eTty + READS: TmunuBase::eTtz WRITES: ML_ADMQuantities::Jadm1 WRITES: ML_ADMQuantities::Jadm2 WRITES: ML_ADMQuantities::Jadm3 diff --git a/ML_ADMQuantities/src/ML_ADMQuantities.cc b/ML_ADMQuantities/src/ML_ADMQuantities.cc index d7e40ab..4879023 100644 --- a/ML_ADMQuantities/src/ML_ADMQuantities.cc +++ b/ML_ADMQuantities/src/ML_ADMQuantities.cc @@ -1074,11 +1074,11 @@ static void ML_ADMQuantities_Body(cGH const * restrict const cctkGH, int const d JacPDstandardNth1gt22*SQR(gtu21) - JacPDstandardNth1gt33*SQR(gtu31); CCTK_REAL CCTK_ATTRIBUTE_UNUSED dgtu211 = - -(gtu11*(gtu21*JacPDstandardNth1gt11 + gtu22*JacPDstandardNth1gt12 + - gtu32*JacPDstandardNth1gt13)) - gtu21*(gtu31*JacPDstandardNth1gt13 + - gtu22*JacPDstandardNth1gt22 + gtu32*JacPDstandardNth1gt23) - - gtu31*(gtu22*JacPDstandardNth1gt23 + gtu32*JacPDstandardNth1gt33) - - JacPDstandardNth1gt12*SQR(gtu21); + -(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 = -(gtu11*(gtu31*JacPDstandardNth1gt11 + gtu32*JacPDstandardNth1gt12 + @@ -1110,11 +1110,11 @@ static void ML_ADMQuantities_Body(cGH const * restrict const cctkGH, int const d JacPDstandardNth2gt22*SQR(gtu21) - JacPDstandardNth2gt33*SQR(gtu31); CCTK_REAL CCTK_ATTRIBUTE_UNUSED dgtu212 = - -(gtu11*(gtu21*JacPDstandardNth2gt11 + gtu22*JacPDstandardNth2gt12 + - gtu32*JacPDstandardNth2gt13)) - gtu21*(gtu31*JacPDstandardNth2gt13 + - gtu22*JacPDstandardNth2gt22 + gtu32*JacPDstandardNth2gt23) - - gtu31*(gtu22*JacPDstandardNth2gt23 + gtu32*JacPDstandardNth2gt33) - - JacPDstandardNth2gt12*SQR(gtu21); + -(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 = -(gtu11*(gtu31*JacPDstandardNth2gt11 + gtu32*JacPDstandardNth2gt12 + @@ -1146,11 +1146,11 @@ static void ML_ADMQuantities_Body(cGH const * restrict const cctkGH, int const d JacPDstandardNth3gt22*SQR(gtu21) - JacPDstandardNth3gt33*SQR(gtu31); CCTK_REAL CCTK_ATTRIBUTE_UNUSED dgtu213 = - -(gtu11*(gtu21*JacPDstandardNth3gt11 + gtu22*JacPDstandardNth3gt12 + - gtu32*JacPDstandardNth3gt13)) - gtu21*(gtu31*JacPDstandardNth3gt13 + - gtu22*JacPDstandardNth3gt22 + gtu32*JacPDstandardNth3gt23) - - gtu31*(gtu22*JacPDstandardNth3gt23 + gtu32*JacPDstandardNth3gt33) - - JacPDstandardNth3gt12*SQR(gtu21); + -(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 = -(gtu11*(gtu31*JacPDstandardNth3gt11 + gtu32*JacPDstandardNth3gt12 + @@ -1365,18 +1365,18 @@ static void ML_ADMQuantities_Body(cGH const * restrict const cctkGH, int const d CCTK_REAL CCTK_ATTRIBUTE_UNUSED Xtn3 = Gt311*gtu11 + Gt322*gtu22 + 2*(Gt312*gtu21 + Gt313*gtu31 + Gt323*gtu32) + Gt333*gtu33; - CCTK_REAL CCTK_ATTRIBUTE_UNUSED Rt11 = 0.5*(6*(Gt111*Gtlu111 + - Gt112*Gtlu112 + Gt113*Gtlu113) + 4*(Gt211*Gtlu121 + Gt212*Gtlu122 + - Gt213*Gtlu123 + Gt311*Gtlu131 + Gt312*Gtlu132 + Gt313*Gtlu133) - - gtu11*JacPDstandardNth11gt11 + 2*(Gt211*Gtlu211 + Gt212*Gtlu212 + - Gt213*Gtlu213 + Gt311*Gtlu311 + Gt312*Gtlu312 + Gt313*Gtlu313 + - gt11L*JacPDstandardNth1Xt1) + 2*gt12L*JacPDstandardNth1Xt2 + - 2*gt13L*JacPDstandardNth1Xt3 + gtu21*(-JacPDstandardNth12gt11 - + 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) - - gtu33*JacPDstandardNth33gt11 + 2*Gtl111*Xtn1 + 2*Gtl112*Xtn2 + - 2*Gtl113*Xtn3); + 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 + @@ -1411,18 +1411,18 @@ static void ML_ADMQuantities_Body(cGH const * restrict const cctkGH, int const d Gtl113*Xtn1 + Gtl311*Xtn1 + Gtl123*Xtn2 + Gtl312*Xtn2 + Gtl133*Xtn3 + Gtl313*Xtn3); - CCTK_REAL CCTK_ATTRIBUTE_UNUSED Rt22 = 0.5*(6*(Gt212*Gtlu221 + - Gt222*Gtlu222 + Gt223*Gtlu223) + 4*(Gt123*Gtlu213 + Gt312*Gtlu231 + - Gt322*Gtlu232 + Gt323*Gtlu233) - gtu11*JacPDstandardNth11gt22 + - gtu21*(-JacPDstandardNth12gt22 - JacPDstandardNth21gt22) - - gtu22*JacPDstandardNth22gt22 + 2*(Gt123*Gtlu123 + Gt112*(Gtlu121 + - 2*Gtlu211) + Gt122*(Gtlu122 + 2*Gtlu212) + Gt312*Gtlu321 + - Gt322*Gtlu322 + Gt323*Gtlu323 + gt12L*JacPDstandardNth2Xt1) + - 2*gt22L*JacPDstandardNth2Xt2 + 2*gt23L*JacPDstandardNth2Xt3 + + 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) - - gtu33*JacPDstandardNth33gt22 + 2*Gtl212*Xtn1 + 2*Gtl222*Xtn2 + - 2*Gtl223*Xtn3); + 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 + @@ -1441,18 +1441,17 @@ static void ML_ADMQuantities_Body(cGH const * restrict const cctkGH, int const d Gtl213*Xtn1 + Gtl312*Xtn1 + Gtl223*Xtn2 + Gtl322*Xtn2 + Gtl233*Xtn3 + Gtl323*Xtn3); - CCTK_REAL CCTK_ATTRIBUTE_UNUSED Rt33 = 0.5*(4*(Gt133*Gtlu313 + - Gt213*Gtlu321 + Gt223*Gtlu322 + Gt233*Gtlu323) + 6*(Gt313*Gtlu331 + - Gt323*Gtlu332 + Gt333*Gtlu333) - gtu11*JacPDstandardNth11gt33 + - gtu21*(-JacPDstandardNth12gt33 - JacPDstandardNth21gt33) - + 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 + - 2*(Gt133*Gtlu133 + Gt213*Gtlu231 + Gt223*Gtlu232 + Gt233*Gtlu233 + - Gt113*(Gtlu131 + 2*Gtlu311) + Gt123*(Gtlu132 + 2*Gtlu312) + - gt13L*JacPDstandardNth3Xt1) + 2*gt23L*JacPDstandardNth3Xt2 + - 2*gt33L*JacPDstandardNth3Xt3 + 2*Gtl313*Xtn1 + 2*Gtl323*Xtn2 + - 2*Gtl333*Xtn3); + 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; @@ -1502,57 +1501,63 @@ static void ML_ADMQuantities_Body(cGH const * restrict const cctkGH, int const d beta2L*eTyzL + beta3L*eTzzL)*INV(alphaL); CCTK_REAL CCTK_ATTRIBUTE_UNUSED MadmL = - -0.0625*INV(Pi)*((Gt111*Gtlu111 + Gt112*Gtlu112 + Gt113*Gtlu113 + - Gt211*Gtlu121 + Gt212*Gtlu122 + Gt213*Gtlu123 + Gt311*Gtlu131 + - Gt312*Gtlu132 + Gt313*Gtlu133)*gtu11 + (Gt122*Gtlu112 + Gt123*Gtlu113 + - Gt222*Gtlu122 + Gt223*Gtlu123 + Gt322*Gtlu132 + Gt323*Gtlu133 + - Gt111*Gtlu211 + Gt112*(Gtlu111 + Gtlu212) + Gt113*Gtlu213 + - Gt211*Gtlu221 + Gt212*(Gtlu121 + Gtlu222) + Gt213*Gtlu223 + - Gt311*Gtlu231 + Gt312*(Gtlu131 + Gtlu232) + Gt313*Gtlu233)*gtu21 + - (Gt112*Gtlu211 + Gt122*Gtlu212 + Gt123*Gtlu213 + Gt212*Gtlu221 + - Gt222*Gtlu222 + Gt223*Gtlu223 + Gt312*Gtlu231 + Gt322*Gtlu232 + - Gt323*Gtlu233)*gtu22 + (Gt123*Gtlu112 + Gt133*Gtlu113 + Gt223*Gtlu122 + - Gt233*Gtlu123 + Gt323*Gtlu132 + Gt333*Gtlu133 + Gt111*Gtlu311 + - Gt112*Gtlu312 + Gt113*(Gtlu111 + Gtlu313) + Gt211*Gtlu321 + - Gt212*Gtlu322 + Gt213*(Gtlu121 + Gtlu323) + Gt311*Gtlu331 + - Gt312*Gtlu332 + Gt313*(Gtlu131 + Gtlu333))*gtu31 + (Gt113*Gtlu211 + - Gt133*Gtlu213 + Gt213*Gtlu221 + Gt233*Gtlu223 + Gt313*Gtlu231 + - Gt333*Gtlu233 + Gt112*Gtlu311 + Gt122*Gtlu312 + Gt123*(Gtlu212 + - Gtlu313) + Gt212*Gtlu321 + Gt222*Gtlu322 + Gt223*(Gtlu222 + Gtlu323) + - Gt312*Gtlu331 + Gt322*Gtlu332 + Gt323*(Gtlu232 + Gtlu333))*gtu32 + - (Gt113*Gtlu311 + Gt123*Gtlu312 + Gt133*Gtlu313 + Gt213*Gtlu321 + - Gt223*Gtlu322 + Gt233*Gtlu323 + Gt313*Gtlu331 + Gt323*Gtlu332 + - Gt333*Gtlu333)*gtu33 + (-1 + ephi)*trRt - pow(ephi,5)*(2*(Atm12*Atm21 + - Atm13*Atm31 + Atm23*Atm32) + 16*Pi*rho - + -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) + + Gtlu121*(Gt211*gtu11 + Gt212*gtu21 + Gt213*gtu31) + + Gtlu122*(Gt212*gtu11 + Gt222*gtu21 + Gt223*gtu31) + + Gtlu123*(Gt213*gtu11 + Gt223*gtu21 + Gt233*gtu31) + + Gtlu131*(Gt311*gtu11 + Gt312*gtu21 + Gt313*gtu31) + + Gtlu132*(Gt312*gtu11 + Gt322*gtu21 + Gt323*gtu31) + + Gtlu133*(Gt313*gtu11 + Gt323*gtu21 + Gt333*gtu31) + + Gtlu211*(Gt111*gtu21 + Gt112*gtu22 + Gt113*gtu32) + + Gtlu212*(Gt112*gtu21 + Gt122*gtu22 + Gt123*gtu32) + + Gtlu213*(Gt113*gtu21 + Gt123*gtu22 + Gt133*gtu32) + + Gtlu221*(Gt211*gtu21 + Gt212*gtu22 + Gt213*gtu32) + + Gtlu222*(Gt212*gtu21 + Gt222*gtu22 + Gt223*gtu32) + + Gtlu223*(Gt213*gtu21 + Gt223*gtu22 + Gt233*gtu32) + + Gtlu231*(Gt311*gtu21 + Gt312*gtu22 + Gt313*gtu32) + + Gtlu232*(Gt312*gtu21 + Gt322*gtu22 + Gt323*gtu32) + + Gtlu233*(Gt313*gtu21 + Gt323*gtu22 + Gt333*gtu32) + + Gtlu311*(Gt111*gtu31 + Gt112*gtu32 + Gt113*gtu33) + + Gtlu312*(Gt112*gtu31 + Gt122*gtu32 + Gt123*gtu33) + + Gtlu313*(Gt113*gtu31 + Gt123*gtu32 + Gt133*gtu33) + + Gtlu321*(Gt211*gtu31 + Gt212*gtu32 + Gt213*gtu33) + + Gtlu322*(Gt212*gtu31 + Gt222*gtu32 + Gt223*gtu33) + + 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))); - CCTK_REAL CCTK_ATTRIBUTE_UNUSED Jadm1L = - 0.0208333333333333333333333333333*(At11L*(3*zL*dgtu112 - - 3*yL*dgtu113) + At22L*(3*zL*dgtu222 - 3*yL*dgtu223) + 6*(Atm23 - + zL*(At12L*dgtu212 + At13L*dgtu312 + At23L*dgtu322)) - - 6*(Atm32 + yL*(At12L*dgtu213 + At13L*dgtu313 + At23L*dgtu323)) - + At33L*(3*zL*dgtu332 - 3*yL*dgtu333) + - zL*(-4*JacPDstandardNth2trK - 48*Pi*S2) + - yL*(4*JacPDstandardNth3trK + 48*Pi*S3))*INV(Pi)*pow(ephi,6); - - CCTK_REAL CCTK_ATTRIBUTE_UNUSED Jadm2L = - 0.0208333333333333333333333333333*(At11L*(-3*zL*dgtu111 + - 3*xL*dgtu113) + At22L*(-3*zL*dgtu221 + 3*xL*dgtu223) - 6*(Atm13 - + zL*(At12L*dgtu211 + At13L*dgtu311 + At23L*dgtu321)) + - 6*(Atm31 + xL*(At12L*dgtu213 + At13L*dgtu313 + At23L*dgtu323)) - + At33L*(-3*zL*dgtu331 + 3*xL*dgtu333) + - zL*(4*JacPDstandardNth1trK + 48*Pi*S1) + - xL*(-4*JacPDstandardNth3trK - 48*Pi*S3))*INV(Pi)*pow(ephi,6); - - CCTK_REAL CCTK_ATTRIBUTE_UNUSED Jadm3L = - 0.0208333333333333333333333333333*(At11L*(3*yL*dgtu111 - - 3*xL*dgtu112) + At22L*(3*yL*dgtu221 - 3*xL*dgtu222) + 6*(Atm12 - + yL*(At12L*dgtu211 + At13L*dgtu311 + At23L*dgtu321)) - - 6*(Atm21 + xL*(At12L*dgtu212 + At13L*dgtu312 + At23L*dgtu322)) - + At33L*(3*yL*dgtu331 - 3*xL*dgtu332) + - yL*(-4*JacPDstandardNth1trK - 48*Pi*S1) + - xL*(4*JacPDstandardNth2trK + 48*Pi*S2))*INV(Pi)*pow(ephi,6); + 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 + + 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); /* Copy local copies back to grid functions */ Jadm1[index] = Jadm1L; @@ -1581,7 +1586,6 @@ extern "C" void ML_ADMQuantities(CCTK_ARGUMENTS) const char *const groups[] = { "grid::coordinates", - "Grid::coordinates", "ML_BSSN::ML_curv", "ML_BSSN::ML_Gamma", "ML_BSSN::ML_lapse", @@ -1591,7 +1595,7 @@ extern "C" void ML_ADMQuantities(CCTK_ARGUMENTS) "ML_BSSN::ML_trace_curv", "ML_ADMQuantities::ML_Jadm", "ML_ADMQuantities::ML_Madm"}; - GenericFD_AssertGroupStorage(cctkGH, "ML_ADMQuantities", 11, groups); + GenericFD_AssertGroupStorage(cctkGH, "ML_ADMQuantities", 10, groups); GenericFD_EnsureStencilFits(cctkGH, "ML_ADMQuantities", 2, 2, 2); diff --git a/ML_WaveToy/src/WT_Standing.cc b/ML_WaveToy/src/WT_Standing.cc index b1a640a..7af03f1 100644 --- a/ML_WaveToy/src/WT_Standing.cc +++ b/ML_WaveToy/src/WT_Standing.cc @@ -100,7 +100,7 @@ static void WT_Standing_Body(cGH const * restrict const cctkGH, int const dir, i /* Calculate temporaries and grid functions */ CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED kvec = - ToReal(Pi*ScalarINV(width)); + ToReal(3.14159265358979323846264338328*ScalarINV(width)); CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED omega = ksqrt(kmul(kmul(kvec,kvec),ToReal(3))); diff --git a/ML_WaveToy_CL/src/WT_CL_Standing.cc b/ML_WaveToy_CL/src/WT_CL_Standing.cc index 1d84d1c..7b35e3f 100644 --- a/ML_WaveToy_CL/src/WT_CL_Standing.cc +++ b/ML_WaveToy_CL/src/WT_CL_Standing.cc @@ -102,7 +102,7 @@ static void WT_CL_Standing_Body(cGH const * restrict const cctkGH, int const dir " \n" " /* Calculate temporaries and grid functions */\n" " CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED kvec = \n" - " ToReal(Pi*ScalarINV(width));\n" + " ToReal(3.14159265358979323846264338328*ScalarINV(width));\n" " \n" " CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED omega = \n" " ksqrt(kmul(kmul(kvec,kvec),ToReal(3)));\n" |