Regenerate code (Kranc commit 10daf4a83a783a812c7b5a3aba0df115f0eac98f)

6 files changed, 84 insertions, 54 deletions

diff --git a/ML_WaveToy_CL/src/WT_CL_Dirichlet.cc b/ML_WaveToy_CL/src/WT_CL_Dirichlet.cc
index d83016a..49cc6ee 100644
--- a/ML_WaveToy_CL/src/WT_CL_Dirichlet.cc
+++ b/ML_WaveToy_CL/src/WT_CL_Dirichlet.cc
@@ -70,13 +70,18 @@ static void WT_CL_Dirichlet_Body(const cGH* restrict const cctkGH, const int dir
   "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"
@@ -122,9 +127,9 @@ static void WT_CL_Dirichlet_Body(const cGH* restrict const cctkGH, const int dir
   "  /* Precompute derivatives */\n"
   "  \n"
   "  /* Calculate temporaries and grid functions */\n"
-  "  CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED urhsL = ToReal(0);\n"
+  "  CCTK_REAL_VEC urhsL CCTK_ATTRIBUTE_UNUSED = ToReal(0);\n"
   "  \n"
-  "  CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED rhorhsL = ToReal(0);\n"
+  "  CCTK_REAL_VEC rhorhsL 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"
diff --git a/ML_WaveToy_CL/src/WT_CL_Energy.cc b/ML_WaveToy_CL/src/WT_CL_Energy.cc
index 4ed7d60..c8beec9 100644
--- a/ML_WaveToy_CL/src/WT_CL_Energy.cc
+++ b/ML_WaveToy_CL/src/WT_CL_Energy.cc
@@ -67,13 +67,18 @@ static void WT_CL_Energy_Body(const cGH* restrict const cctkGH, const int dir, c
   "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"
@@ -124,7 +129,7 @@ static void WT_CL_Energy_Body(const cGH* restrict const cctkGH, const int dir, c
   "  const CCTK_REAL_VEC PDstandardNth3u CCTK_ATTRIBUTE_UNUSED = PDstandardNth3(&u[index]);\n"
   "  \n"
   "  /* Calculate temporaries and grid functions */\n"
-  "  CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED epsL = \n"
+  "  CCTK_REAL_VEC epsL CCTK_ATTRIBUTE_UNUSED = \n"
   "    kmul(kmadd(rhoL,rhoL,kmadd(PDstandardNth1u,PDstandardNth1u,kmadd(PDstandardNth2u,PDstandardNth2u,kmul(PDstandardNth3u,PDstandardNth3u)))),ToReal(0.5));\n"
   "  \n"
   "  /* Copy local copies back to grid functions */\n"
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"
diff --git a/ML_WaveToy_CL/src/WT_CL_Gaussian.cc b/ML_WaveToy_CL/src/WT_CL_Gaussian.cc
index 9cfe036..4d29eff 100644
--- a/ML_WaveToy_CL/src/WT_CL_Gaussian.cc
+++ b/ML_WaveToy_CL/src/WT_CL_Gaussian.cc
@@ -53,13 +53,18 @@ static void WT_CL_Gaussian_Body(const cGH* restrict const cctkGH, const int dir,
   "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"
@@ -106,10 +111,10 @@ static void WT_CL_Gaussian_Body(const cGH* restrict const cctkGH, const int dir,
   "  /* Precompute derivatives */\n"
   "  \n"
   "  /* Calculate temporaries and grid functions */\n"
-  "  CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED uL = \n"
+  "  CCTK_REAL_VEC uL CCTK_ATTRIBUTE_UNUSED = \n"
   "    kmul(kexp(kmul(kmul(kmul(rL,rL),ToReal(-0.5)),ToReal(ScalarINV(ScalarSQR(width))))),ToReal(amplitude));\n"
   "  \n"
-  "  CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED rhoL = ToReal(0);\n"
+  "  CCTK_REAL_VEC rhoL 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"
diff --git a/ML_WaveToy_CL/src/WT_CL_RHS.cc b/ML_WaveToy_CL/src/WT_CL_RHS.cc
index 128496a..454ca3a 100644
--- a/ML_WaveToy_CL/src/WT_CL_RHS.cc
+++ b/ML_WaveToy_CL/src/WT_CL_RHS.cc
@@ -70,13 +70,18 @@ static void WT_CL_RHS_Body(const cGH* restrict const cctkGH, const int dir, cons
   "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"
@@ -127,9 +132,9 @@ static void WT_CL_RHS_Body(const cGH* restrict const cctkGH, const int dir, cons
   "  const CCTK_REAL_VEC PDstandardNth33u CCTK_ATTRIBUTE_UNUSED = PDstandardNth33(&u[index]);\n"
   "  \n"
   "  /* Calculate temporaries and grid functions */\n"
-  "  CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED urhsL = rhoL;\n"
+  "  CCTK_REAL_VEC urhsL CCTK_ATTRIBUTE_UNUSED = rhoL;\n"
   "  \n"
-  "  CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED rhorhsL = \n"
+  "  CCTK_REAL_VEC rhorhsL CCTK_ATTRIBUTE_UNUSED = \n"
   "    kadd(PDstandardNth11u,kadd(PDstandardNth22u,PDstandardNth33u));\n"
   "  \n"
   "  /* Copy local copies back to grid functions */\n"
diff --git a/ML_WaveToy_CL/src/WT_CL_Standing.cc b/ML_WaveToy_CL/src/WT_CL_Standing.cc
index 6d590b8..9a154f7 100644
--- a/ML_WaveToy_CL/src/WT_CL_Standing.cc
+++ b/ML_WaveToy_CL/src/WT_CL_Standing.cc
@@ -53,13 +53,18 @@ static void WT_CL_Standing_Body(const cGH* restrict const cctkGH, const int dir,
   "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"
@@ -108,16 +113,16 @@ static void WT_CL_Standing_Body(const cGH* restrict const cctkGH, const int dir,
   "  /* Precompute derivatives */\n"
   "  \n"
   "  /* Calculate temporaries and grid functions */\n"
-  "  CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED kvec = \n"
+  "  CCTK_REAL_VEC kvec CCTK_ATTRIBUTE_UNUSED = \n"
   "    ToReal(3.14159265358979323846264338328*ScalarINV(width));\n"
   "  \n"
-  "  CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED omega = \n"
+  "  CCTK_REAL_VEC omega CCTK_ATTRIBUTE_UNUSED = \n"
   "    ksqrt(kmul(kmul(kvec,kvec),ToReal(3)));\n"
   "  \n"
-  "  CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED uL = \n"
+  "  CCTK_REAL_VEC uL CCTK_ATTRIBUTE_UNUSED = \n"
   "    kmul(kcos(kmul(xL,kvec)),kmul(kcos(kmul(yL,kvec)),kmul(kcos(kmul(zL,kvec)),kmul(kcos(kmul(omega,t)),ToReal(amplitude)))));\n"
   "  \n"
-  "  CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED rhoL = \n"
+  "  CCTK_REAL_VEC rhoL CCTK_ATTRIBUTE_UNUSED = \n"
   "    kneg(kmul(omega,kmul(kcos(kmul(xL,kvec)),kmul(kcos(kmul(yL,kvec)),kmul(kcos(kmul(zL,kvec)),kmul(ksin(kmul(omega,t)),ToReal(amplitude)))))));\n"
   "  \n"
   "  /* Copy local copies back to grid functions */\n"