aboutsummaryrefslogtreecommitdiff
path: root/ML_CCZ4/src/ML_CCZ4_InitGamma.cc
diff options
context:
space:
mode:
authorIan Hinder <ian.hinder@aei.mpg.de>2013-09-10 16:03:25 +0200
committerIan Hinder <ian.hinder@aei.mpg.de>2013-09-10 16:03:25 +0200
commit8df3aae833eb6bf334b0decb57ebf9e1baf33ddf (patch)
treebd9193699a88b7a5743bdee1d58ebf203ef7b94b /ML_CCZ4/src/ML_CCZ4_InitGamma.cc
parent51f5f80e196f950622f63977ae344580faf1cb2e (diff)
Regenerate code (Kranc commit 10daf4a83a783a812c7b5a3aba0df115f0eac98f)
Diffstat (limited to 'ML_CCZ4/src/ML_CCZ4_InitGamma.cc')
-rw-r--r--ML_CCZ4/src/ML_CCZ4_InitGamma.cc19
1 files changed, 12 insertions, 7 deletions
diff --git a/ML_CCZ4/src/ML_CCZ4_InitGamma.cc b/ML_CCZ4/src/ML_CCZ4_InitGamma.cc
index cf4904b..e16fb52 100644
--- a/ML_CCZ4/src/ML_CCZ4_InitGamma.cc
+++ b/ML_CCZ4/src/ML_CCZ4_InitGamma.cc
@@ -51,13 +51,18 @@ static void ML_CCZ4_InitGamma_Body(const cGH* restrict const cctkGH, const int d
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 p1o1024dx CCTK_ATTRIBUTE_UNUSED = kdiv(ToReal(0.0009765625),dx);