diff options
Diffstat (limited to 'ML_ADM/src/ML_ADM_convertToADMBase.cc')
-rw-r--r-- | ML_ADM/src/ML_ADM_convertToADMBase.cc | 59 |
1 files changed, 32 insertions, 27 deletions
diff --git a/ML_ADM/src/ML_ADM_convertToADMBase.cc b/ML_ADM/src/ML_ADM_convertToADMBase.cc index 4bed36b..b1e128b 100644 --- a/ML_ADM/src/ML_ADM_convertToADMBase.cc +++ b/ML_ADM/src/ML_ADM_convertToADMBase.cc @@ -51,13 +51,18 @@ static void ML_ADM_convertToADMBase_Body(const cGH* restrict const cctkGH, const 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,45 +244,45 @@ static void ML_ADM_convertToADMBase_Body(const cGH* restrict const cctkGH, const } /* Calculate temporaries and grid functions */ - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED gxxL = g11L; + CCTK_REAL_VEC gxxL CCTK_ATTRIBUTE_UNUSED = g11L; - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED gxyL = g12L; + CCTK_REAL_VEC gxyL CCTK_ATTRIBUTE_UNUSED = g12L; - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED gxzL = g13L; + CCTK_REAL_VEC gxzL CCTK_ATTRIBUTE_UNUSED = g13L; - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED gyyL = g22L; + CCTK_REAL_VEC gyyL CCTK_ATTRIBUTE_UNUSED = g22L; - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED gyzL = g23L; + CCTK_REAL_VEC gyzL CCTK_ATTRIBUTE_UNUSED = g23L; - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED gzzL = g33L; + CCTK_REAL_VEC gzzL CCTK_ATTRIBUTE_UNUSED = g33L; - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED kxxL = K11L; + CCTK_REAL_VEC kxxL CCTK_ATTRIBUTE_UNUSED = K11L; - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED kxyL = K12L; + CCTK_REAL_VEC kxyL CCTK_ATTRIBUTE_UNUSED = K12L; - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED kxzL = K13L; + CCTK_REAL_VEC kxzL CCTK_ATTRIBUTE_UNUSED = K13L; - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED kyyL = K22L; + CCTK_REAL_VEC kyyL CCTK_ATTRIBUTE_UNUSED = K22L; - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED kyzL = K23L; + CCTK_REAL_VEC kyzL CCTK_ATTRIBUTE_UNUSED = K23L; - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED kzzL = K33L; + CCTK_REAL_VEC kzzL CCTK_ATTRIBUTE_UNUSED = K33L; - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED alpL = alphaL; + CCTK_REAL_VEC alpL CCTK_ATTRIBUTE_UNUSED = alphaL; - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED dtalpL = ToReal(0); + CCTK_REAL_VEC dtalpL CCTK_ATTRIBUTE_UNUSED = ToReal(0); - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED betaxL = beta1L; + CCTK_REAL_VEC betaxL CCTK_ATTRIBUTE_UNUSED = beta1L; - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED betayL = beta2L; + CCTK_REAL_VEC betayL CCTK_ATTRIBUTE_UNUSED = beta2L; - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED betazL = beta3L; + CCTK_REAL_VEC betazL CCTK_ATTRIBUTE_UNUSED = beta3L; - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED dtbetaxL = ToReal(0); + CCTK_REAL_VEC dtbetaxL CCTK_ATTRIBUTE_UNUSED = ToReal(0); - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED dtbetayL = ToReal(0); + CCTK_REAL_VEC dtbetayL CCTK_ATTRIBUTE_UNUSED = ToReal(0); - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED dtbetazL = ToReal(0); + CCTK_REAL_VEC dtbetazL CCTK_ATTRIBUTE_UNUSED = ToReal(0); /* Copy local copies back to grid functions */ vec_store_partial_prepare(i,vecimin,vecimax); |