diff options
Diffstat (limited to 'ML_WaveToy_CL/src/WT_CL_EnergyBoundary.cc')
| -rw-r--r-- | ML_WaveToy_CL/src/WT_CL_EnergyBoundary.cc | 21 |
1 files changed, 13 insertions, 8 deletions
diff --git a/ML_WaveToy_CL/src/WT_CL_EnergyBoundary.cc b/ML_WaveToy_CL/src/WT_CL_EnergyBoundary.cc index 8161612..d9fb6f8 100644 --- a/ML_WaveToy_CL/src/WT_CL_EnergyBoundary.cc +++ b/ML_WaveToy_CL/src/WT_CL_EnergyBoundary.cc @@ -67,13 +67,18 @@ static void WT_CL_EnergyBoundary_Body(const cGH* restrict const cctkGH, const in "const CCTK_REAL_VEC dyi CCTK_ATTRIBUTE_UNUSED = INV(dy);\n" "const CCTK_REAL_VEC dzi CCTK_ATTRIBUTE_UNUSED = INV(dz);\n" "const CCTK_REAL_VEC khalf CCTK_ATTRIBUTE_UNUSED = ToReal(0.5);\n" - "const CCTK_REAL_VEC kthird CCTK_ATTRIBUTE_UNUSED = ToReal(1.0/3.0);\n" - "const CCTK_REAL_VEC ktwothird CCTK_ATTRIBUTE_UNUSED = ToReal(2.0/3.0);\n" - "const CCTK_REAL_VEC kfourthird CCTK_ATTRIBUTE_UNUSED = ToReal(4.0/3.0);\n" - "const CCTK_REAL_VEC keightthird CCTK_ATTRIBUTE_UNUSED = ToReal(8.0/3.0);\n" - "const CCTK_REAL_VEC hdxi CCTK_ATTRIBUTE_UNUSED = kmul(ToReal(0.5), dxi);\n" - "const CCTK_REAL_VEC hdyi CCTK_ATTRIBUTE_UNUSED = kmul(ToReal(0.5), dyi);\n" - "const CCTK_REAL_VEC hdzi CCTK_ATTRIBUTE_UNUSED = kmul(ToReal(0.5), dzi);\n" + "const CCTK_REAL_VEC kthird CCTK_ATTRIBUTE_UNUSED = \n" + " ToReal(0.333333333333333333333333333333);\n" + "const CCTK_REAL_VEC ktwothird CCTK_ATTRIBUTE_UNUSED = \n" + " ToReal(0.666666666666666666666666666667);\n" + "const CCTK_REAL_VEC kfourthird CCTK_ATTRIBUTE_UNUSED = \n" + " ToReal(1.33333333333333333333333333333);\n" + "const CCTK_REAL_VEC hdxi CCTK_ATTRIBUTE_UNUSED = \n" + " kmul(dxi,ToReal(0.5));\n" + "const CCTK_REAL_VEC hdyi CCTK_ATTRIBUTE_UNUSED = \n" + " kmul(dyi,ToReal(0.5));\n" + "const CCTK_REAL_VEC hdzi CCTK_ATTRIBUTE_UNUSED = \n" + " kmul(dzi,ToReal(0.5));\n" "\n" "/* Initialize predefined quantities */\n" "const CCTK_REAL_VEC p1o12dx CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(0.0833333333333333333333333333333),dx);\n" @@ -119,7 +124,7 @@ static void WT_CL_EnergyBoundary_Body(const cGH* restrict const cctkGH, const in " /* Precompute derivatives */\n" " \n" " /* Calculate temporaries and grid functions */\n" - " CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED epsL = ToReal(0);\n" + " CCTK_REAL_VEC epsL CCTK_ATTRIBUTE_UNUSED = ToReal(0);\n" " \n" " /* Copy local copies back to grid functions */\n" " vec_store_partial_prepare(i,lc_imin,lc_imax);\n" |