diff options
Diffstat (limited to 'ML_ADM/src/ML_ADM_convertFromADMBase.cc')
-rw-r--r-- | ML_ADM/src/ML_ADM_convertFromADMBase.cc | 51 |
1 files changed, 28 insertions, 23 deletions
diff --git a/ML_ADM/src/ML_ADM_convertFromADMBase.cc b/ML_ADM/src/ML_ADM_convertFromADMBase.cc index 5043bec..f8006f6 100644 --- a/ML_ADM/src/ML_ADM_convertFromADMBase.cc +++ b/ML_ADM/src/ML_ADM_convertFromADMBase.cc @@ -51,13 +51,18 @@ static void ML_ADM_convertFromADMBase_Body(const cGH* restrict const cctkGH, con const CCTK_REAL_VEC dyi CCTK_ATTRIBUTE_UNUSED = INV(dy); const CCTK_REAL_VEC dzi CCTK_ATTRIBUTE_UNUSED = INV(dz); const CCTK_REAL_VEC khalf CCTK_ATTRIBUTE_UNUSED = ToReal(0.5); - const CCTK_REAL_VEC kthird CCTK_ATTRIBUTE_UNUSED = ToReal(1.0/3.0); - const CCTK_REAL_VEC ktwothird CCTK_ATTRIBUTE_UNUSED = ToReal(2.0/3.0); - const CCTK_REAL_VEC kfourthird CCTK_ATTRIBUTE_UNUSED = ToReal(4.0/3.0); - const CCTK_REAL_VEC keightthird CCTK_ATTRIBUTE_UNUSED = ToReal(8.0/3.0); - const CCTK_REAL_VEC hdxi CCTK_ATTRIBUTE_UNUSED = kmul(ToReal(0.5), dxi); - const CCTK_REAL_VEC hdyi CCTK_ATTRIBUTE_UNUSED = kmul(ToReal(0.5), dyi); - const CCTK_REAL_VEC hdzi CCTK_ATTRIBUTE_UNUSED = kmul(ToReal(0.5), dzi); + const CCTK_REAL_VEC kthird CCTK_ATTRIBUTE_UNUSED = + ToReal(0.333333333333333333333333333333); + const CCTK_REAL_VEC ktwothird CCTK_ATTRIBUTE_UNUSED = + ToReal(0.666666666666666666666666666667); + const CCTK_REAL_VEC kfourthird CCTK_ATTRIBUTE_UNUSED = + ToReal(1.33333333333333333333333333333); + const CCTK_REAL_VEC hdxi CCTK_ATTRIBUTE_UNUSED = + kmul(dxi,ToReal(0.5)); + const CCTK_REAL_VEC hdyi CCTK_ATTRIBUTE_UNUSED = + kmul(dyi,ToReal(0.5)); + const CCTK_REAL_VEC hdzi CCTK_ATTRIBUTE_UNUSED = + kmul(dzi,ToReal(0.5)); /* Initialize predefined quantities */ const CCTK_REAL_VEC p1o12dx CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(0.0833333333333333333333333333333),dx); @@ -239,37 +244,37 @@ static void ML_ADM_convertFromADMBase_Body(const cGH* restrict const cctkGH, con } /* Calculate temporaries and grid functions */ - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED g11L = gxxL; + CCTK_REAL_VEC g11L CCTK_ATTRIBUTE_UNUSED = gxxL; - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED g12L = gxyL; + CCTK_REAL_VEC g12L CCTK_ATTRIBUTE_UNUSED = gxyL; - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED g13L = gxzL; + CCTK_REAL_VEC g13L CCTK_ATTRIBUTE_UNUSED = gxzL; - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED g22L = gyyL; + CCTK_REAL_VEC g22L CCTK_ATTRIBUTE_UNUSED = gyyL; - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED g23L = gyzL; + CCTK_REAL_VEC g23L CCTK_ATTRIBUTE_UNUSED = gyzL; - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED g33L = gzzL; + CCTK_REAL_VEC g33L CCTK_ATTRIBUTE_UNUSED = gzzL; - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED K11L = kxxL; + CCTK_REAL_VEC K11L CCTK_ATTRIBUTE_UNUSED = kxxL; - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED K12L = kxyL; + CCTK_REAL_VEC K12L CCTK_ATTRIBUTE_UNUSED = kxyL; - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED K13L = kxzL; + CCTK_REAL_VEC K13L CCTK_ATTRIBUTE_UNUSED = kxzL; - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED K22L = kyyL; + CCTK_REAL_VEC K22L CCTK_ATTRIBUTE_UNUSED = kyyL; - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED K23L = kyzL; + CCTK_REAL_VEC K23L CCTK_ATTRIBUTE_UNUSED = kyzL; - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED K33L = kzzL; + CCTK_REAL_VEC K33L CCTK_ATTRIBUTE_UNUSED = 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); |