diff options
Diffstat (limited to 'ML_BSSN/src/ML_BSSN_constraints2.cc')
| -rw-r--r-- | ML_BSSN/src/ML_BSSN_constraints2.cc | 116 |
1 files changed, 58 insertions, 58 deletions
diff --git a/ML_BSSN/src/ML_BSSN_constraints2.cc b/ML_BSSN/src/ML_BSSN_constraints2.cc index f33cab2..29c7116 100644 --- a/ML_BSSN/src/ML_BSSN_constraints2.cc +++ b/ML_BSSN/src/ML_BSSN_constraints2.cc @@ -817,176 +817,176 @@ static void ML_BSSN_constraints2_Body(const cGH* restrict const cctkGH, const in JacPDstandardNth3trK = PDstandardNth3trK; } - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED detgt = ToReal(1); + CCTK_REAL_VEC detgt CCTK_ATTRIBUTE_UNUSED = ToReal(1); - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED gtu11 = + CCTK_REAL_VEC gtu11 CCTK_ATTRIBUTE_UNUSED = kdiv(kmsub(gt22L,gt33L,kmul(gt23L,gt23L)),detgt); - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED gtu12 = + CCTK_REAL_VEC gtu12 CCTK_ATTRIBUTE_UNUSED = kdiv(kmsub(gt13L,gt23L,kmul(gt12L,gt33L)),detgt); - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED gtu13 = + CCTK_REAL_VEC gtu13 CCTK_ATTRIBUTE_UNUSED = kdiv(kmsub(gt12L,gt23L,kmul(gt13L,gt22L)),detgt); - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED gtu22 = + CCTK_REAL_VEC gtu22 CCTK_ATTRIBUTE_UNUSED = kdiv(kmsub(gt11L,gt33L,kmul(gt13L,gt13L)),detgt); - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED gtu23 = + CCTK_REAL_VEC gtu23 CCTK_ATTRIBUTE_UNUSED = kdiv(kmsub(gt12L,gt13L,kmul(gt11L,gt23L)),detgt); - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED gtu33 = + CCTK_REAL_VEC gtu33 CCTK_ATTRIBUTE_UNUSED = kdiv(kmsub(gt11L,gt22L,kmul(gt12L,gt12L)),detgt); - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED Gtl111 = + CCTK_REAL_VEC Gtl111 CCTK_ATTRIBUTE_UNUSED = kmul(JacPDstandardNth1gt11,ToReal(0.5)); - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED Gtl112 = + CCTK_REAL_VEC Gtl112 CCTK_ATTRIBUTE_UNUSED = kmul(JacPDstandardNth2gt11,ToReal(0.5)); - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED Gtl113 = + CCTK_REAL_VEC Gtl113 CCTK_ATTRIBUTE_UNUSED = kmul(JacPDstandardNth3gt11,ToReal(0.5)); - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED Gtl122 = + CCTK_REAL_VEC Gtl122 CCTK_ATTRIBUTE_UNUSED = kmadd(JacPDstandardNth1gt22,ToReal(-0.5),JacPDstandardNth2gt12); - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED Gtl123 = + CCTK_REAL_VEC Gtl123 CCTK_ATTRIBUTE_UNUSED = kmul(kadd(JacPDstandardNth2gt13,ksub(JacPDstandardNth3gt12,JacPDstandardNth1gt23)),ToReal(0.5)); - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED Gtl133 = + CCTK_REAL_VEC Gtl133 CCTK_ATTRIBUTE_UNUSED = kmadd(JacPDstandardNth1gt33,ToReal(-0.5),JacPDstandardNth3gt13); - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED Gtl211 = + CCTK_REAL_VEC Gtl211 CCTK_ATTRIBUTE_UNUSED = kmadd(JacPDstandardNth2gt11,ToReal(-0.5),JacPDstandardNth1gt12); - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED Gtl212 = + CCTK_REAL_VEC Gtl212 CCTK_ATTRIBUTE_UNUSED = kmul(JacPDstandardNth1gt22,ToReal(0.5)); - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED Gtl213 = + CCTK_REAL_VEC Gtl213 CCTK_ATTRIBUTE_UNUSED = kmul(kadd(JacPDstandardNth1gt23,ksub(JacPDstandardNth3gt12,JacPDstandardNth2gt13)),ToReal(0.5)); - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED Gtl222 = + CCTK_REAL_VEC Gtl222 CCTK_ATTRIBUTE_UNUSED = kmul(JacPDstandardNth2gt22,ToReal(0.5)); - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED Gtl223 = + CCTK_REAL_VEC Gtl223 CCTK_ATTRIBUTE_UNUSED = kmul(JacPDstandardNth3gt22,ToReal(0.5)); - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED Gtl233 = + CCTK_REAL_VEC Gtl233 CCTK_ATTRIBUTE_UNUSED = kmadd(JacPDstandardNth2gt33,ToReal(-0.5),JacPDstandardNth3gt23); - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED Gtl311 = + CCTK_REAL_VEC Gtl311 CCTK_ATTRIBUTE_UNUSED = kmadd(JacPDstandardNth3gt11,ToReal(-0.5),JacPDstandardNth1gt13); - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED Gtl312 = + CCTK_REAL_VEC Gtl312 CCTK_ATTRIBUTE_UNUSED = kmul(kadd(JacPDstandardNth1gt23,ksub(JacPDstandardNth2gt13,JacPDstandardNth3gt12)),ToReal(0.5)); - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED Gtl313 = + CCTK_REAL_VEC Gtl313 CCTK_ATTRIBUTE_UNUSED = kmul(JacPDstandardNth1gt33,ToReal(0.5)); - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED Gtl322 = + CCTK_REAL_VEC Gtl322 CCTK_ATTRIBUTE_UNUSED = kmadd(JacPDstandardNth3gt22,ToReal(-0.5),JacPDstandardNth2gt23); - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED Gtl323 = + CCTK_REAL_VEC Gtl323 CCTK_ATTRIBUTE_UNUSED = kmul(JacPDstandardNth2gt33,ToReal(0.5)); - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED Gtl333 = + CCTK_REAL_VEC Gtl333 CCTK_ATTRIBUTE_UNUSED = kmul(JacPDstandardNth3gt33,ToReal(0.5)); - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED Gt111 = + CCTK_REAL_VEC Gt111 CCTK_ATTRIBUTE_UNUSED = kmadd(Gtl111,gtu11,kmadd(Gtl211,gtu12,kmul(Gtl311,gtu13))); - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED Gt211 = + CCTK_REAL_VEC Gt211 CCTK_ATTRIBUTE_UNUSED = kmadd(Gtl111,gtu12,kmadd(Gtl211,gtu22,kmul(Gtl311,gtu23))); - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED Gt311 = + CCTK_REAL_VEC Gt311 CCTK_ATTRIBUTE_UNUSED = kmadd(Gtl111,gtu13,kmadd(Gtl211,gtu23,kmul(Gtl311,gtu33))); - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED Gt112 = + CCTK_REAL_VEC Gt112 CCTK_ATTRIBUTE_UNUSED = kmadd(Gtl112,gtu11,kmadd(Gtl212,gtu12,kmul(Gtl312,gtu13))); - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED Gt212 = + CCTK_REAL_VEC Gt212 CCTK_ATTRIBUTE_UNUSED = kmadd(Gtl112,gtu12,kmadd(Gtl212,gtu22,kmul(Gtl312,gtu23))); - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED Gt312 = + CCTK_REAL_VEC Gt312 CCTK_ATTRIBUTE_UNUSED = kmadd(Gtl112,gtu13,kmadd(Gtl212,gtu23,kmul(Gtl312,gtu33))); - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED Gt113 = + CCTK_REAL_VEC Gt113 CCTK_ATTRIBUTE_UNUSED = kmadd(Gtl113,gtu11,kmadd(Gtl213,gtu12,kmul(Gtl313,gtu13))); - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED Gt213 = + CCTK_REAL_VEC Gt213 CCTK_ATTRIBUTE_UNUSED = kmadd(Gtl113,gtu12,kmadd(Gtl213,gtu22,kmul(Gtl313,gtu23))); - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED Gt313 = + CCTK_REAL_VEC Gt313 CCTK_ATTRIBUTE_UNUSED = kmadd(Gtl113,gtu13,kmadd(Gtl213,gtu23,kmul(Gtl313,gtu33))); - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED Gt122 = + CCTK_REAL_VEC Gt122 CCTK_ATTRIBUTE_UNUSED = kmadd(Gtl122,gtu11,kmadd(Gtl222,gtu12,kmul(Gtl322,gtu13))); - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED Gt222 = + CCTK_REAL_VEC Gt222 CCTK_ATTRIBUTE_UNUSED = kmadd(Gtl122,gtu12,kmadd(Gtl222,gtu22,kmul(Gtl322,gtu23))); - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED Gt322 = + CCTK_REAL_VEC Gt322 CCTK_ATTRIBUTE_UNUSED = kmadd(Gtl122,gtu13,kmadd(Gtl222,gtu23,kmul(Gtl322,gtu33))); - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED Gt123 = + CCTK_REAL_VEC Gt123 CCTK_ATTRIBUTE_UNUSED = kmadd(Gtl123,gtu11,kmadd(Gtl223,gtu12,kmul(Gtl323,gtu13))); - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED Gt223 = + CCTK_REAL_VEC Gt223 CCTK_ATTRIBUTE_UNUSED = kmadd(Gtl123,gtu12,kmadd(Gtl223,gtu22,kmul(Gtl323,gtu23))); - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED Gt323 = + CCTK_REAL_VEC Gt323 CCTK_ATTRIBUTE_UNUSED = kmadd(Gtl123,gtu13,kmadd(Gtl223,gtu23,kmul(Gtl323,gtu33))); - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED Gt133 = + CCTK_REAL_VEC Gt133 CCTK_ATTRIBUTE_UNUSED = kmadd(Gtl133,gtu11,kmadd(Gtl233,gtu12,kmul(Gtl333,gtu13))); - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED Gt233 = + CCTK_REAL_VEC Gt233 CCTK_ATTRIBUTE_UNUSED = kmadd(Gtl133,gtu12,kmadd(Gtl233,gtu22,kmul(Gtl333,gtu23))); - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED Gt333 = + CCTK_REAL_VEC Gt333 CCTK_ATTRIBUTE_UNUSED = kmadd(Gtl133,gtu13,kmadd(Gtl233,gtu23,kmul(Gtl333,gtu33))); - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED fac1 = IfThen(conformalMethod == + CCTK_REAL_VEC fac1 CCTK_ATTRIBUTE_UNUSED = IfThen(conformalMethod == 1,kdiv(ToReal(-0.5),phiL),ToReal(1)); - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED cdphi1 = + CCTK_REAL_VEC cdphi1 CCTK_ATTRIBUTE_UNUSED = kmul(fac1,JacPDstandardNth1phi); - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED cdphi2 = + CCTK_REAL_VEC cdphi2 CCTK_ATTRIBUTE_UNUSED = kmul(fac1,JacPDstandardNth2phi); - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED cdphi3 = + CCTK_REAL_VEC cdphi3 CCTK_ATTRIBUTE_UNUSED = kmul(fac1,JacPDstandardNth3phi); - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED S1 = + CCTK_REAL_VEC S1 CCTK_ATTRIBUTE_UNUSED = kdiv(kmadd(beta1L,eTxxL,kmadd(beta2L,eTxyL,kmsub(beta3L,eTxzL,eTtxL))),alphaL); - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED S2 = + CCTK_REAL_VEC S2 CCTK_ATTRIBUTE_UNUSED = kdiv(kmadd(beta1L,eTxyL,kmadd(beta2L,eTyyL,kmsub(beta3L,eTyzL,eTtyL))),alphaL); - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED S3 = + CCTK_REAL_VEC S3 CCTK_ATTRIBUTE_UNUSED = kdiv(kmadd(beta1L,eTxzL,kmadd(beta2L,eTyzL,kmsub(beta3L,eTzzL,eTtzL))),alphaL); - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED M1L = + CCTK_REAL_VEC M1L CCTK_ATTRIBUTE_UNUSED = kmadd(gtu11,JacPDstandardNth1At11,kmadd(gtu12,JacPDstandardNth2At11,kmadd(gtu22,JacPDstandardNth2At12,kmadd(gtu23,JacPDstandardNth2At13,kmadd(gtu13,JacPDstandardNth3At11,kmadd(gtu23,JacPDstandardNth3At12,kmadd(gtu33,JacPDstandardNth3At13,kmadd(gtu12,kadd(JacPDstandardNth1At12,kmsub(At13L,kmul(Gt312,ToReal(-3)),kmul(At22L,Gt211))),kmadd(gtu13,kadd(JacPDstandardNth1At13,kmsub(At13L,kmul(Gt313,ToReal(-3)),kmul(At23L,Gt211))),knmsub(gtu22,kmadd(At23L,Gt312,kmul(At22L,Gt212)),kmadd(Gt311,kmsub(At13L,kmul(gtu11,ToReal(-2)),kmadd(At33L,gtu13,kmul(At23L,gtu12))),kmadd(At13L,kmsub(Gt323,kmul(gtu23,ToReal(-2)),kmul(Gt113,gtu33)),kmadd(JacPDstandardNth1trK,ToReal(-0.666666666666666666666666666667),knmsub(At12L,kmadd(Gt111,gtu12,kmadd(Gt112,gtu22,kmadd(Gt222,gtu22,kmadd(Gt113,gtu23,kmadd(Gt233,gtu33,kmadd(cdphi1,kmul(gtu12,ToReal(-6)),kmadd(cdphi2,kmul(gtu22,ToReal(-6)),kmadd(cdphi3,kmul(gtu23,ToReal(-6)),kmadd(Gt211,kmul(gtu11,ToReal(2)),kmadd(Gt223,kmul(gtu23,ToReal(2)),kmadd(Gt212,kmul(gtu12,ToReal(3)),kmul(Gt213,kmul(gtu13,ToReal(3)))))))))))))),kmadd(gtu23,knmsub(At22L,Gt213,knmsub(At33L,Gt312,kmsub(At13L,kmul(cdphi2,ToReal(6)),kmul(At23L,kadd(Gt313,Gt212))))),kmadd(gtu33,knmsub(At23L,Gt213,kmsub(At13L,kmul(cdphi3,ToReal(6)),kmul(At33L,Gt313))),kmadd(At11L,kmadd(Gt112,kmul(gtu12,ToReal(-3)),kmadd(Gt113,kmul(gtu13,ToReal(-3)),kmadd(Gt111,kmul(gtu11,ToReal(-2)),kmadd(Gt123,kmul(gtu23,ToReal(-2)),knmsub(Gt122,gtu22,knmsub(Gt133,gtu33,kmadd(cdphi1,kmul(gtu11,ToReal(6)),kmadd(cdphi2,kmul(gtu12,ToReal(6)),kmul(cdphi3,kmul(gtu13,ToReal(6))))))))))),kmadd(At13L,knmsub(Gt322,gtu22,knmsub(Gt112,gtu23,kmsub(gtu13,kmsub(cdphi1,ToReal(6),Gt111),kmul(Gt333,gtu33)))),kmul(S1,ToReal(-25.1327412287183459077011470662)))))))))))))))))))); - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED M2L = + CCTK_REAL_VEC M2L CCTK_ATTRIBUTE_UNUSED = kmadd(gtu11,JacPDstandardNth1At12,kmadd(gtu12,JacPDstandardNth1At22,kmadd(gtu13,JacPDstandardNth1At23,kmadd(gtu12,JacPDstandardNth2At12,kmadd(gtu22,JacPDstandardNth2At22,kmadd(gtu23,JacPDstandardNth2At23,kmadd(gtu13,JacPDstandardNth3At12,kmadd(gtu23,JacPDstandardNth3At22,kmadd(gtu33,JacPDstandardNth3At23,kmadd(At22L,kmul(Gt222,kmul(gtu22,ToReal(-2))),knmsub(At11L,kmadd(Gt112,gtu11,kmadd(Gt122,gtu12,kmul(Gt123,gtu13))),knmsub(gtu11,kmadd(At22L,Gt211,kmadd(At13L,Gt312,kmul(At23L,Gt311))),kmadd(Gt212,kmsub(At22L,kmul(gtu12,ToReal(-3)),kmul(At23L,gtu13)),kmadd(Gt312,kmsub(At23L,kmul(gtu12,ToReal(-3)),kmul(At33L,gtu13)),kmadd(Gt322,kmsub(At23L,kmul(gtu22,ToReal(-2)),kmul(At33L,gtu23)),kmadd(Gt223,kmsub(At22L,kmul(gtu23,ToReal(-3)),kmul(At23L,gtu33)),kmadd(Gt323,kmsub(At23L,kmul(gtu23,ToReal(-3)),kmul(At33L,gtu33)),kmadd(At22L,kmsub(Gt213,kmul(gtu13,ToReal(-2)),kmul(Gt233,gtu33)),kmadd(At23L,kmsub(Gt313,kmul(gtu13,ToReal(-2)),kmadd(Gt333,gtu33,kmul(Gt222,gtu23))),kmadd(JacPDstandardNth2trK,ToReal(-0.666666666666666666666666666667),kmadd(At22L,kmul(cdphi2,kmul(gtu22,ToReal(6))),kmadd(At22L,kmul(cdphi3,kmul(gtu23,ToReal(6))),kmadd(gtu12,kmsub(At22L,kmul(cdphi1,ToReal(6)),kmul(At13L,Gt322)),kmadd(gtu13,kmsub(At23L,kmul(cdphi1,ToReal(6)),kmul(At13L,kadd(Gt323,Gt112))),kmadd(gtu23,kmsub(At23L,kmul(cdphi2,ToReal(6)),kmul(At13L,Gt122)),kmadd(gtu33,kmsub(At23L,kmul(cdphi3,ToReal(6)),kmul(At13L,Gt123)),kmadd(At12L,kmadd(Gt112,kmul(gtu12,ToReal(-3)),kmadd(Gt123,kmul(gtu23,ToReal(-3)),kmadd(Gt113,kmul(gtu13,ToReal(-2)),kmadd(Gt122,kmul(gtu22,ToReal(-2)),knmsub(Gt133,gtu33,kmadd(gtu11,ksub(kmsub(cdphi1,ToReal(6),Gt212),Gt111),kmadd(gtu12,kmsub(cdphi2,ToReal(6),Gt222),kmul(gtu13,kmsub(cdphi3,ToReal(6),Gt223))))))))),kmul(S2,ToReal(-25.1327412287183459077011470662))))))))))))))))))))))))))))); - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED M3L = + CCTK_REAL_VEC M3L CCTK_ATTRIBUTE_UNUSED = kmadd(gtu11,JacPDstandardNth1At13,kmadd(gtu12,JacPDstandardNth1At23,kmadd(gtu13,JacPDstandardNth1At33,kmadd(gtu12,JacPDstandardNth2At13,kmadd(gtu22,JacPDstandardNth2At23,kmadd(gtu23,JacPDstandardNth2At33,kmadd(gtu13,JacPDstandardNth3At13,kmadd(gtu23,JacPDstandardNth3At23,kmadd(gtu33,JacPDstandardNth3At33,kmadd(At33L,kmul(Gt313,kmul(gtu13,ToReal(-3))),kmadd(At23L,kmul(Gt223,kmul(gtu23,ToReal(-3))),kmadd(At33L,kmul(Gt323,kmul(gtu23,ToReal(-3))),kmadd(At23L,kmul(Gt233,kmul(gtu33,ToReal(-2))),kmadd(At33L,kmul(Gt333,kmul(gtu33,ToReal(-2))),knmsub(At11L,kmadd(Gt113,gtu11,kmadd(Gt123,gtu12,kmul(Gt133,gtu13))),kmadd(gtu12,kmsub(At23L,kmul(Gt212,ToReal(-2)),kmul(At12L,Gt223)),knmsub(gtu11,kmadd(At23L,Gt211,kmadd(At33L,Gt311,kmul(At12L,Gt213))),kmadd(At33L,kmsub(Gt312,kmul(gtu12,ToReal(-2)),kmul(Gt322,gtu22)),kmadd(At23L,kmsub(Gt213,kmul(gtu13,ToReal(-3)),kmadd(Gt333,gtu23,kmul(Gt222,gtu22))),kmadd(JacPDstandardNth3trK,ToReal(-0.666666666666666666666666666667),kmadd(At33L,kmul(cdphi3,kmul(gtu33,ToReal(6))),kmadd(gtu13,kmsub(At33L,kmul(cdphi1,ToReal(6)),kmul(At12L,Gt233)),kmadd(gtu23,kmsub(At33L,kmul(cdphi2,ToReal(6)),kmul(At12L,Gt133)),kmadd(gtu23,kmsub(At23L,kmul(cdphi3,ToReal(6)),kmul(At22L,Gt233)),kmadd(gtu12,knmsub(At12L,Gt113,kmsub(At23L,kmsub(cdphi1,ToReal(6),Gt313),kmul(At22L,Gt213))),kmadd(gtu22,knmsub(At12L,Gt123,kmsub(At23L,kmsub(cdphi2,ToReal(6),Gt323),kmul(At22L,Gt223))),kmadd(At13L,kmadd(Gt113,kmul(gtu13,ToReal(-3)),kmadd(Gt123,kmul(gtu23,ToReal(-3)),kmadd(Gt112,kmul(gtu12,ToReal(-2)),kmadd(Gt133,kmul(gtu33,ToReal(-2)),knmsub(Gt122,gtu22,kmadd(gtu11,ksub(kmsub(cdphi1,ToReal(6),Gt313),Gt111),kmadd(gtu12,kmsub(cdphi2,ToReal(6),Gt323),kmul(gtu13,kmsub(cdphi3,ToReal(6),Gt333))))))))),kmul(S3,ToReal(-25.1327412287183459077011470662))))))))))))))))))))))))))))); - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED cSL = klog(detgt); + CCTK_REAL_VEC cSL CCTK_ATTRIBUTE_UNUSED = klog(detgt); - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED cXt1L = + CCTK_REAL_VEC cXt1L CCTK_ATTRIBUTE_UNUSED = kmadd(Gt111,gtu11,kmadd(Gt122,gtu22,kmadd(Gt133,gtu33,kmsub(kmadd(Gt112,gtu12,kmadd(Gt113,gtu13,kmul(Gt123,gtu23))),ToReal(2),Xt1L)))); - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED cXt2L = + CCTK_REAL_VEC cXt2L CCTK_ATTRIBUTE_UNUSED = kmadd(Gt211,gtu11,kmadd(Gt222,gtu22,kmadd(Gt233,gtu33,kmsub(kmadd(Gt212,gtu12,kmadd(Gt213,gtu13,kmul(Gt223,gtu23))),ToReal(2),Xt2L)))); - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED cXt3L = + CCTK_REAL_VEC cXt3L CCTK_ATTRIBUTE_UNUSED = kmadd(Gt311,gtu11,kmadd(Gt322,gtu22,kmadd(Gt333,gtu33,kmsub(kmadd(Gt312,gtu12,kmadd(Gt313,gtu13,kmul(Gt323,gtu23))),ToReal(2),Xt3L)))); - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED cAL = + CCTK_REAL_VEC cAL CCTK_ATTRIBUTE_UNUSED = kmadd(At11L,gtu11,kmadd(At22L,gtu22,kmadd(At33L,gtu33,kmul(kmadd(At12L,gtu12,kmadd(At13L,gtu13,kmul(At23L,gtu23))),ToReal(2))))); /* Copy local copies back to grid functions */ |