diff options
Diffstat (limited to 'ML_ADM/src/ML_ADM_constraints_boundary.cc')
| -rw-r--r-- | ML_ADM/src/ML_ADM_constraints_boundary.cc | 27 |
1 files changed, 16 insertions, 11 deletions
diff --git a/ML_ADM/src/ML_ADM_constraints_boundary.cc b/ML_ADM/src/ML_ADM_constraints_boundary.cc index 7f93d73..cfe9564 100644 --- a/ML_ADM/src/ML_ADM_constraints_boundary.cc +++ b/ML_ADM/src/ML_ADM_constraints_boundary.cc @@ -68,13 +68,18 @@ static void ML_ADM_constraints_boundary_Body(const cGH* restrict const cctkGH, c 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); @@ -240,13 +245,13 @@ static void ML_ADM_constraints_boundary_Body(const cGH* restrict const cctkGH, c } /* Calculate temporaries and grid functions */ - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED HL = ToReal(0); + CCTK_REAL_VEC HL CCTK_ATTRIBUTE_UNUSED = ToReal(0); - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED M1L = ToReal(0); + CCTK_REAL_VEC M1L CCTK_ATTRIBUTE_UNUSED = ToReal(0); - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED M2L = ToReal(0); + CCTK_REAL_VEC M2L CCTK_ATTRIBUTE_UNUSED = ToReal(0); - CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED M3L = ToReal(0); + CCTK_REAL_VEC M3L CCTK_ATTRIBUTE_UNUSED = ToReal(0); /* Copy local copies back to grid functions */ vec_store_partial_prepare(i,vecimin,vecimax); |