diff options
Diffstat (limited to 'ML_BSSN/src/ML_BSSN_convertFromADMBase.cc')
-rw-r--r-- | ML_BSSN/src/ML_BSSN_convertFromADMBase.cc | 58 |
1 files changed, 29 insertions, 29 deletions
diff --git a/ML_BSSN/src/ML_BSSN_convertFromADMBase.cc b/ML_BSSN/src/ML_BSSN_convertFromADMBase.cc index 5bafb53..b843763 100644 --- a/ML_BSSN/src/ML_BSSN_convertFromADMBase.cc +++ b/ML_BSSN/src/ML_BSSN_convertFromADMBase.cc @@ -283,37 +283,37 @@ static void ML_BSSN_convertFromADMBase_Body(const cGH* restrict const cctkGH, co } /* Calculate temporaries and grid functions */ - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED g11 = gxxL; + CCTK_REAL_VEC g11 CCTK_ATTRIBUTE_UNUSED = gxxL; - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED g12 = gxyL; + CCTK_REAL_VEC g12 CCTK_ATTRIBUTE_UNUSED = gxyL; - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED g13 = gxzL; + CCTK_REAL_VEC g13 CCTK_ATTRIBUTE_UNUSED = gxzL; - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED g22 = gyyL; + CCTK_REAL_VEC g22 CCTK_ATTRIBUTE_UNUSED = gyyL; - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED g23 = gyzL; + CCTK_REAL_VEC g23 CCTK_ATTRIBUTE_UNUSED = gyzL; - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED g33 = gzzL; + CCTK_REAL_VEC g33 CCTK_ATTRIBUTE_UNUSED = gzzL; - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED detg = + CCTK_REAL_VEC detg CCTK_ATTRIBUTE_UNUSED = knmsub(g22,kmul(g13,g13),knmsub(g11,kmul(g23,g23),kmadd(g33,kmsub(g11,g22,kmul(g12,g12)),kmul(g12,kmul(g13,kmul(g23,ToReal(2))))))); - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED gu11 = + CCTK_REAL_VEC gu11 CCTK_ATTRIBUTE_UNUSED = kdiv(kmsub(g22,g33,kmul(g23,g23)),detg); - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED gu12 = + CCTK_REAL_VEC gu12 CCTK_ATTRIBUTE_UNUSED = kdiv(kmsub(g13,g23,kmul(g12,g33)),detg); - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED gu13 = + CCTK_REAL_VEC gu13 CCTK_ATTRIBUTE_UNUSED = kdiv(kmsub(g12,g23,kmul(g13,g22)),detg); - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED gu22 = + CCTK_REAL_VEC gu22 CCTK_ATTRIBUTE_UNUSED = kdiv(kmsub(g11,g33,kmul(g13,g13)),detg); - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED gu23 = + CCTK_REAL_VEC gu23 CCTK_ATTRIBUTE_UNUSED = kdiv(kmsub(g12,g13,kmul(g11,g23)),detg); - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED gu33 = + CCTK_REAL_VEC gu33 CCTK_ATTRIBUTE_UNUSED = kdiv(kmsub(g11,g22,kmul(g12,g12)),detg); CCTK_REAL_VEC em4phi CCTK_ATTRIBUTE_UNUSED; @@ -331,46 +331,46 @@ static void ML_BSSN_convertFromADMBase_Body(const cGH* restrict const cctkGH, co em4phi = kexp(kmul(phiL,ToReal(-4))); } - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED gt11L = kmul(em4phi,g11); + CCTK_REAL_VEC gt11L CCTK_ATTRIBUTE_UNUSED = kmul(em4phi,g11); - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED gt12L = kmul(em4phi,g12); + CCTK_REAL_VEC gt12L CCTK_ATTRIBUTE_UNUSED = kmul(em4phi,g12); - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED gt13L = kmul(em4phi,g13); + CCTK_REAL_VEC gt13L CCTK_ATTRIBUTE_UNUSED = kmul(em4phi,g13); - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED gt22L = kmul(em4phi,g22); + CCTK_REAL_VEC gt22L CCTK_ATTRIBUTE_UNUSED = kmul(em4phi,g22); - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED gt23L = kmul(em4phi,g23); + CCTK_REAL_VEC gt23L CCTK_ATTRIBUTE_UNUSED = kmul(em4phi,g23); - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED gt33L = kmul(em4phi,g33); + CCTK_REAL_VEC gt33L CCTK_ATTRIBUTE_UNUSED = kmul(em4phi,g33); trKL = kmadd(kxxL,gu11,kmadd(kyyL,gu22,kmadd(kzzL,gu33,kmul(kmadd(kxyL,gu12,kmadd(kxzL,gu13,kmul(kyzL,gu23))),ToReal(2))))); - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED At11L = + CCTK_REAL_VEC At11L CCTK_ATTRIBUTE_UNUSED = kmul(em4phi,kmadd(trKL,kmul(g11,ToReal(-0.333333333333333333333333333333)),kxxL)); - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED At12L = + CCTK_REAL_VEC At12L CCTK_ATTRIBUTE_UNUSED = kmul(em4phi,kmadd(trKL,kmul(g12,ToReal(-0.333333333333333333333333333333)),kxyL)); - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED At13L = + CCTK_REAL_VEC At13L CCTK_ATTRIBUTE_UNUSED = kmul(em4phi,kmadd(trKL,kmul(g13,ToReal(-0.333333333333333333333333333333)),kxzL)); - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED At22L = + CCTK_REAL_VEC At22L CCTK_ATTRIBUTE_UNUSED = kmul(em4phi,kmadd(trKL,kmul(g22,ToReal(-0.333333333333333333333333333333)),kyyL)); - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED At23L = + CCTK_REAL_VEC At23L CCTK_ATTRIBUTE_UNUSED = kmul(em4phi,kmadd(trKL,kmul(g23,ToReal(-0.333333333333333333333333333333)),kyzL)); - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED At33L = + CCTK_REAL_VEC At33L CCTK_ATTRIBUTE_UNUSED = kmul(em4phi,kmadd(trKL,kmul(g33,ToReal(-0.333333333333333333333333333333)),kzzL)); - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED alphaL = alpL; + CCTK_REAL_VEC alphaL CCTK_ATTRIBUTE_UNUSED = alpL; - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED beta1L = betaxL; + CCTK_REAL_VEC beta1L CCTK_ATTRIBUTE_UNUSED = betaxL; - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED beta2L = betayL; + CCTK_REAL_VEC beta2L CCTK_ATTRIBUTE_UNUSED = betayL; - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED beta3L = betazL; + CCTK_REAL_VEC beta3L CCTK_ATTRIBUTE_UNUSED = betazL; /* Copy local copies back to grid functions */ vec_store_partial_prepare(i,vecimin,vecimax); |