diff options
Diffstat (limited to 'ML_BSSN_UPW/src/ML_BSSN_UPW_RHS2.cc')
| -rw-r--r-- | ML_BSSN_UPW/src/ML_BSSN_UPW_RHS2.cc | 66 |
1 files changed, 33 insertions, 33 deletions
diff --git a/ML_BSSN_UPW/src/ML_BSSN_UPW_RHS2.cc b/ML_BSSN_UPW/src/ML_BSSN_UPW_RHS2.cc index 4178637..50107ad 100644 --- a/ML_BSSN_UPW/src/ML_BSSN_UPW_RHS2.cc +++ b/ML_BSSN_UPW/src/ML_BSSN_UPW_RHS2.cc @@ -1626,6 +1626,35 @@ static void ML_BSSN_UPW_RHS2_Body(cGH const * restrict const cctkGH, int const d CCTK_REAL_VEC Xtn3 = kmadd(Gt311,gtu11,kmadd(Gt322,gtu22,kmadd(Gt333,gtu33,kmul(kmadd(Gt312,gtu12,kmadd(Gt313,gtu13,kmul(Gt323,gtu23))),ToReal(2))))); + CCTK_REAL_VEC e4phi = IfThen(conformalMethod == + 1,INV(SQR(phiL)),kexp(kmul(phiL,ToReal(4)))); + + CCTK_REAL_VEC em4phi = INV(e4phi); + + CCTK_REAL_VEC g11 = kmul(gt11L,e4phi); + + CCTK_REAL_VEC g12 = kmul(gt12L,e4phi); + + CCTK_REAL_VEC g13 = kmul(gt13L,e4phi); + + CCTK_REAL_VEC g22 = kmul(gt22L,e4phi); + + CCTK_REAL_VEC g23 = kmul(gt23L,e4phi); + + CCTK_REAL_VEC g33 = kmul(gt33L,e4phi); + + CCTK_REAL_VEC gu11 = kmul(em4phi,gtu11); + + CCTK_REAL_VEC gu12 = kmul(em4phi,gtu12); + + CCTK_REAL_VEC gu13 = kmul(em4phi,gtu13); + + CCTK_REAL_VEC gu22 = kmul(em4phi,gtu22); + + CCTK_REAL_VEC gu23 = kmul(em4phi,gtu23); + + CCTK_REAL_VEC gu33 = kmul(em4phi,gtu33); + CCTK_REAL_VEC Rt11 = kmul(ToReal(0.5),knmsub(gtu11,JacPDstandardNth11gt11,knmsub(gtu22,JacPDstandardNth22gt11,knmsub(gtu33,JacPDstandardNth33gt11,knmsub(gtu12,kadd(JacPDstandardNth21gt11,JacPDstandardNth12gt11),knmsub(gtu13,kadd(JacPDstandardNth31gt11,JacPDstandardNth13gt11),knmsub(gtu23,kadd(JacPDstandardNth32gt11,JacPDstandardNth23gt11),kmadd(kmadd(Gt211,Gtlu211,kmadd(Gt212,Gtlu212,kmadd(Gt213,Gtlu213,kmadd(Gt311,Gtlu311,kmadd(Gt312,Gtlu312,kmadd(Gt313,Gtlu313,kmul(gt11L,JacPDstandardNth1Xt1))))))),ToReal(2),kmadd(gt12L,kmul(JacPDstandardNth1Xt2,ToReal(2)),kmadd(gt13L,kmul(JacPDstandardNth1Xt3,ToReal(2)),kmadd(Gtl111,kmul(Xtn1,ToReal(2)),kmadd(Gtl112,kmul(Xtn2,ToReal(2)),kmadd(Gtl113,kmul(Xtn3,ToReal(2)),kmadd(kmadd(Gt211,Gtlu121,kmadd(Gt212,Gtlu122,kmadd(Gt213,Gtlu123,kmadd(Gt311,Gtlu131,kmadd(Gt312,Gtlu132,kmul(Gt313,Gtlu133)))))),ToReal(4),kmul(kmadd(Gt111,Gtlu111,kmadd(Gt112,Gtlu112,kmul(Gt113,Gtlu113))),ToReal(6)))))))))))))))); @@ -1644,8 +1673,8 @@ static void ML_BSSN_UPW_RHS2_Body(cGH const * restrict const cctkGH, int const d CCTK_REAL_VEC Rt33 = kmul(ToReal(0.5),knmsub(gtu11,JacPDstandardNth11gt33,knmsub(gtu22,JacPDstandardNth22gt33,knmsub(gtu33,JacPDstandardNth33gt33,knmsub(gtu12,kadd(JacPDstandardNth21gt33,JacPDstandardNth12gt33),knmsub(gtu13,kadd(JacPDstandardNth31gt33,JacPDstandardNth13gt33),knmsub(gtu23,kadd(JacPDstandardNth32gt33,JacPDstandardNth23gt33),kmadd(gt23L,kmul(JacPDstandardNth3Xt2,ToReal(2)),kmadd(gt33L,kmul(JacPDstandardNth3Xt3,ToReal(2)),kmadd(Gtl313,kmul(Xtn1,ToReal(2)),kmadd(Gtl323,kmul(Xtn2,ToReal(2)),kmadd(Gtl333,kmul(Xtn3,ToReal(2)),kmadd(ToReal(2),kmadd(Gt133,Gtlu133,kmadd(Gt213,Gtlu231,kmadd(Gt223,Gtlu232,kmadd(Gt233,Gtlu233,kmadd(gt13L,JacPDstandardNth3Xt1,kmadd(Gt113,kmadd(Gtlu311,ToReal(2),Gtlu131),kmul(Gt123,kmadd(Gtlu312,ToReal(2),Gtlu132)))))))),kmadd(kmadd(Gt133,Gtlu313,kmadd(Gt213,Gtlu321,kmadd(Gt223,Gtlu322,kmul(Gt233,Gtlu323)))),ToReal(4),kmul(kmadd(Gt313,Gtlu331,kmadd(Gt323,Gtlu332,kmul(Gt333,Gtlu333))),ToReal(6)))))))))))))))); - CCTK_REAL_VEC fac1 = - IfThen(conformalMethod,kmul(INV(phiL),ToReal(-0.5)),ToReal(1)); + CCTK_REAL_VEC fac1 = IfThen(conformalMethod == + 1,kmul(INV(phiL),ToReal(-0.5)),ToReal(1)); CCTK_REAL_VEC cdphi1 = kmul(fac1,JacPDstandardNth1phi); @@ -1653,8 +1682,8 @@ static void ML_BSSN_UPW_RHS2_Body(cGH const * restrict const cctkGH, int const d CCTK_REAL_VEC cdphi3 = kmul(fac1,JacPDstandardNth3phi); - CCTK_REAL_VEC fac2 = - IfThen(conformalMethod,kmul(INV(SQR(phiL)),ToReal(0.5)),ToReal(0)); + CCTK_REAL_VEC fac2 = IfThen(conformalMethod == + 1,kmul(INV(SQR(phiL)),ToReal(0.5)),ToReal(0)); CCTK_REAL_VEC cdphi211 = kmadd(fac2,SQR(JacPDstandardNth1phi),kmul(fac1,ksub(JacPDstandardNth11phi,kmadd(Gt111,JacPDstandardNth1phi,kmadd(Gt311,JacPDstandardNth3phi,kmul(Gt211,JacPDstandardNth2phi)))))); @@ -1719,35 +1748,6 @@ static void ML_BSSN_UPW_RHS2_Body(cGH const * restrict const cctkGH, int const d CCTK_REAL_VEC Atm33 = kmadd(At13L,gtu13,kmadd(At23L,gtu23,kmul(At33L,gtu33))); - CCTK_REAL_VEC e4phi = - IfThen(conformalMethod,INV(SQR(phiL)),kexp(kmul(phiL,ToReal(4)))); - - CCTK_REAL_VEC em4phi = INV(e4phi); - - CCTK_REAL_VEC g11 = kmul(gt11L,e4phi); - - CCTK_REAL_VEC g12 = kmul(gt12L,e4phi); - - CCTK_REAL_VEC g13 = kmul(gt13L,e4phi); - - CCTK_REAL_VEC g22 = kmul(gt22L,e4phi); - - CCTK_REAL_VEC g23 = kmul(gt23L,e4phi); - - CCTK_REAL_VEC g33 = kmul(gt33L,e4phi); - - CCTK_REAL_VEC gu11 = kmul(em4phi,gtu11); - - CCTK_REAL_VEC gu12 = kmul(em4phi,gtu12); - - CCTK_REAL_VEC gu13 = kmul(em4phi,gtu13); - - CCTK_REAL_VEC gu22 = kmul(em4phi,gtu22); - - CCTK_REAL_VEC gu23 = kmul(em4phi,gtu23); - - CCTK_REAL_VEC gu33 = kmul(em4phi,gtu33); - CCTK_REAL_VEC R11 = kadd(Rphi11,Rt11); CCTK_REAL_VEC R12 = kadd(Rphi12,Rt12); |