aboutsummaryrefslogtreecommitdiff
path: root/Examples/Burgers
diff options
context:
space:
mode:
Diffstat (limited to 'Examples/Burgers')
-rw-r--r--Examples/Burgers/interface.ccl2
-rw-r--r--Examples/Burgers/schedule.ccl37
-rw-r--r--Examples/Burgers/src/Boundaries.cc10
-rw-r--r--Examples/Burgers/src/Differencing.h60
-rw-r--r--Examples/Burgers/src/RegisterMoL.cc2
-rw-r--r--Examples/Burgers/src/Startup.cc2
-rw-r--r--Examples/Burgers/src/burgers_flux_1.cc67
-rw-r--r--Examples/Burgers/src/burgers_initial_shock.cc60
-rw-r--r--Examples/Burgers/src/burgers_initial_sine.cc59
-rw-r--r--Examples/Burgers/src/burgers_reconstruct_1.cc74
-rw-r--r--Examples/Burgers/src/burgers_rhs_1.cc64
-rw-r--r--Examples/Burgers/src/burgers_zero_rhs.cc56
12 files changed, 243 insertions, 250 deletions
diff --git a/Examples/Burgers/interface.ccl b/Examples/Burgers/interface.ccl
index fadd7cc..e1dbc48 100644
--- a/Examples/Burgers/interface.ccl
+++ b/Examples/Burgers/interface.ccl
@@ -52,7 +52,7 @@ CCTK_REAL u_group type=GF timelevels=3 tags='tensortypealias="Scalar" tensorweig
} "u_group"
public:
-CCTK_REAL u_grouprhs type=GF timelevels=3 tags='tensortypealias="Scalar" tensorweight=0'
+CCTK_REAL u_grouprhs type=GF timelevels=3 tags='tensortypealias="Scalar" tensorweight=0 Prolongation="None"'
{
urhs
} "u_grouprhs"
diff --git a/Examples/Burgers/schedule.ccl b/Examples/Burgers/schedule.ccl
index a52f5d3..ad2286a 100644
--- a/Examples/Burgers/schedule.ccl
+++ b/Examples/Burgers/schedule.ccl
@@ -60,8 +60,8 @@ if (CCTK_EQUALS(initial_data, "sine"))
schedule burgers_initial_sine at CCTK_INITIAL as burgers_initial
{
LANG: C
- READS: grid::x
- WRITES: Burgers::u
+ READS: grid::x(Everywhere)
+ WRITES: Burgers::u(Everywhere)
} "burgers_initial_sine"
}
@@ -71,28 +71,27 @@ if (CCTK_EQUALS(initial_data, "shock"))
schedule burgers_initial_shock at CCTK_INITIAL as burgers_initial
{
LANG: C
- READS: grid::x
- WRITES: Burgers::u
+ READS: grid::x(Everywhere)
+ WRITES: Burgers::u(Everywhere)
} "burgers_initial_shock"
}
schedule burgers_zero_rhs in MoL_CalcRHS
{
LANG: C
- WRITES: Burgers::urhs
+ WRITES: Burgers::urhs(Everywhere)
} "burgers_zero_rhs"
schedule group burgers_reconstruct_1_group in MoL_CalcRHS after burgers_zero_rhs
{
- # no language specified
} "burgers_reconstruct_1"
schedule burgers_reconstruct_1 in burgers_reconstruct_1_group
{
LANG: C
- READS: Burgers::u
- WRITES: Burgers::uLeft
- WRITES: Burgers::uR
+ READS: Burgers::u(Everywhere)
+ WRITES: Burgers::uLeft(Interior)
+ WRITES: Burgers::uR(Interior)
} "burgers_reconstruct_1"
schedule burgers_reconstruct_1_SelectBCs in burgers_reconstruct_1_bc_group
@@ -105,30 +104,26 @@ schedule burgers_reconstruct_1_SelectBCs in burgers_reconstruct_1_bc_group
schedule group ApplyBCs as burgers_reconstruct_1_ApplyBCs in burgers_reconstruct_1_bc_group after burgers_reconstruct_1_SelectBCs
{
- # no language specified
} "Apply BCs for groups set in burgers_reconstruct_1"
schedule group burgers_reconstruct_1_bc_group in burgers_reconstruct_1_group after burgers_reconstruct_1
{
- # no language specified
} "burgers_reconstruct_1"
schedule group burgers_reconstruct_1_bc_group in MoL_PseudoEvolutionBoundaries after MoL_PostStep
{
- # no language specified
} "burgers_reconstruct_1"
schedule group burgers_flux_1_group in MoL_CalcRHS after burgers_reconstruct_1
{
- # no language specified
} "burgers_flux_1"
schedule burgers_flux_1 in burgers_flux_1_group
{
LANG: C
- READS: Burgers::uLeft
- READS: Burgers::uR
- WRITES: Burgers::uF
+ READS: Burgers::uLeft(Everywhere)
+ READS: Burgers::uR(Everywhere)
+ WRITES: Burgers::uF(Interior)
} "burgers_flux_1"
schedule burgers_flux_1_SelectBCs in burgers_flux_1_bc_group
@@ -140,25 +135,22 @@ schedule burgers_flux_1_SelectBCs in burgers_flux_1_bc_group
schedule group ApplyBCs as burgers_flux_1_ApplyBCs in burgers_flux_1_bc_group after burgers_flux_1_SelectBCs
{
- # no language specified
} "Apply BCs for groups set in burgers_flux_1"
schedule group burgers_flux_1_bc_group in burgers_flux_1_group after burgers_flux_1
{
- # no language specified
} "burgers_flux_1"
schedule group burgers_flux_1_bc_group in MoL_PseudoEvolutionBoundaries after MoL_PostStep
{
- # no language specified
} "burgers_flux_1"
schedule burgers_rhs_1 in MoL_CalcRHS after burgers_flux_1
{
LANG: C
- READS: Burgers::uF
- READS: Burgers::urhs
- WRITES: Burgers::urhs
+ READS: Burgers::uF(Everywhere)
+ READS: Burgers::urhs(Everywhere)
+ WRITES: Burgers::urhs(Interior)
} "burgers_rhs_1"
schedule Burgers_SelectBoundConds in MoL_PostStep
@@ -182,5 +174,4 @@ schedule Burgers_RegisterVars in MoL_Register
schedule group ApplyBCs as Burgers_ApplyBCs in MoL_PostStep after Burgers_SelectBoundConds
{
- # no language specified
} "Apply boundary conditions controlled by thorn Boundary"
diff --git a/Examples/Burgers/src/Boundaries.cc b/Examples/Burgers/src/Boundaries.cc
index b0ba04d..876a25b 100644
--- a/Examples/Burgers/src/Boundaries.cc
+++ b/Examples/Burgers/src/Boundaries.cc
@@ -30,7 +30,7 @@ extern "C" void Burgers_SelectBoundConds(CCTK_ARGUMENTS)
DECLARE_CCTK_ARGUMENTS;
DECLARE_CCTK_PARAMETERS;
- CCTK_INT ierr = 0;
+ CCTK_INT ierr CCTK_ATTRIBUTE_UNUSED = 0;
if (CCTK_EQUALS(u_group_bound, "none" ) ||
CCTK_EQUALS(u_group_bound, "static") ||
@@ -57,7 +57,7 @@ extern "C" void Burgers_SelectBoundConds(CCTK_ARGUMENTS)
if (CCTK_EQUALS(u_group_bound, "radiative"))
{
/* select radiation boundary condition */
- static CCTK_INT handle_u_group_bound = -1;
+ static CCTK_INT handle_u_group_bound CCTK_ATTRIBUTE_UNUSED = -1;
if (handle_u_group_bound < 0) handle_u_group_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE);
if (handle_u_group_bound < 0) CCTK_WARN(0, "could not create table!");
if (Util_TableSetReal(handle_u_group_bound , u_group_bound_limit, "LIMIT") < 0)
@@ -76,7 +76,7 @@ extern "C" void Burgers_SelectBoundConds(CCTK_ARGUMENTS)
if (CCTK_EQUALS(u_bound, "radiative"))
{
/* select radiation boundary condition */
- static CCTK_INT handle_u_bound = -1;
+ static CCTK_INT handle_u_bound CCTK_ATTRIBUTE_UNUSED = -1;
if (handle_u_bound < 0) handle_u_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE);
if (handle_u_bound < 0) CCTK_WARN(0, "could not create table!");
if (Util_TableSetReal(handle_u_bound , u_bound_limit, "LIMIT") < 0)
@@ -95,7 +95,7 @@ extern "C" void Burgers_SelectBoundConds(CCTK_ARGUMENTS)
if (CCTK_EQUALS(u_group_bound, "scalar"))
{
/* select scalar boundary condition */
- static CCTK_INT handle_u_group_bound = -1;
+ static CCTK_INT handle_u_group_bound CCTK_ATTRIBUTE_UNUSED = -1;
if (handle_u_group_bound < 0) handle_u_group_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE);
if (handle_u_group_bound < 0) CCTK_WARN(0, "could not create table!");
if (Util_TableSetReal(handle_u_group_bound ,u_group_bound_scalar, "SCALAR") < 0)
@@ -112,7 +112,7 @@ extern "C" void Burgers_SelectBoundConds(CCTK_ARGUMENTS)
if (CCTK_EQUALS(u_bound, "scalar"))
{
/* select scalar boundary condition */
- static CCTK_INT handle_u_bound = -1;
+ static CCTK_INT handle_u_bound CCTK_ATTRIBUTE_UNUSED = -1;
if (handle_u_bound < 0) handle_u_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE);
if (handle_u_bound < 0) CCTK_WARN(0, "could not create table!");
if (Util_TableSetReal(handle_u_bound ,u_bound_scalar, "SCALAR") < 0)
diff --git a/Examples/Burgers/src/Differencing.h b/Examples/Burgers/src/Differencing.h
index dec7ce7..ea2f017 100644
--- a/Examples/Burgers/src/Differencing.h
+++ b/Examples/Burgers/src/Differencing.h
@@ -1,108 +1,108 @@
#ifndef KRANC_DIFF_FUNCTIONS
-# define PDplus1(u) ((-KRANC_GFOFFSET3D(u,0,0,0) + KRANC_GFOFFSET3D(u,1,0,0))*p1odx)
+# define PDplus1(u) ((-1.*KRANC_GFOFFSET3D(u,0,0,0) + KRANC_GFOFFSET3D(u,1,0,0))*p1odx)
#else
# define PDplus1(u) (PDplus1_impl(u,p1odx,cdj,cdk))
static CCTK_REAL PDplus1_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1odx, ptrdiff_t const cdj, ptrdiff_t const cdk) CCTK_ATTRIBUTE_NOINLINE CCTK_ATTRIBUTE_UNUSED;
static CCTK_REAL PDplus1_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1odx, ptrdiff_t const cdj, ptrdiff_t const cdk)
{
- ptrdiff_t const cdi=sizeof(CCTK_REAL);
- return (-KRANC_GFOFFSET3D(u,0,0,0) + KRANC_GFOFFSET3D(u,1,0,0))*p1odx;
+ ptrdiff_t const cdi CCTK_ATTRIBUTE_UNUSED = sizeof(CCTK_REAL);
+ return (-1.*KRANC_GFOFFSET3D(u,0,0,0) + KRANC_GFOFFSET3D(u,1,0,0))*p1odx;
}
#endif
#ifndef KRANC_DIFF_FUNCTIONS
-# define PDplus2(u) ((-KRANC_GFOFFSET3D(u,0,0,0) + KRANC_GFOFFSET3D(u,0,1,0))*p1ody)
+# define PDplus2(u) ((-1.*KRANC_GFOFFSET3D(u,0,0,0) + KRANC_GFOFFSET3D(u,0,1,0))*p1ody)
#else
# define PDplus2(u) (PDplus2_impl(u,p1ody,cdj,cdk))
static CCTK_REAL PDplus2_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1ody, ptrdiff_t const cdj, ptrdiff_t const cdk) CCTK_ATTRIBUTE_NOINLINE CCTK_ATTRIBUTE_UNUSED;
static CCTK_REAL PDplus2_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1ody, ptrdiff_t const cdj, ptrdiff_t const cdk)
{
- ptrdiff_t const cdi=sizeof(CCTK_REAL);
- return (-KRANC_GFOFFSET3D(u,0,0,0) + KRANC_GFOFFSET3D(u,0,1,0))*p1ody;
+ ptrdiff_t const cdi CCTK_ATTRIBUTE_UNUSED = sizeof(CCTK_REAL);
+ return (-1.*KRANC_GFOFFSET3D(u,0,0,0) + KRANC_GFOFFSET3D(u,0,1,0))*p1ody;
}
#endif
#ifndef KRANC_DIFF_FUNCTIONS
-# define PDplus3(u) ((-KRANC_GFOFFSET3D(u,0,0,0) + KRANC_GFOFFSET3D(u,0,0,1))*p1odz)
+# define PDplus3(u) ((-1.*KRANC_GFOFFSET3D(u,0,0,0) + KRANC_GFOFFSET3D(u,0,0,1))*p1odz)
#else
# define PDplus3(u) (PDplus3_impl(u,p1odz,cdj,cdk))
static CCTK_REAL PDplus3_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1odz, ptrdiff_t const cdj, ptrdiff_t const cdk) CCTK_ATTRIBUTE_NOINLINE CCTK_ATTRIBUTE_UNUSED;
static CCTK_REAL PDplus3_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1odz, ptrdiff_t const cdj, ptrdiff_t const cdk)
{
- ptrdiff_t const cdi=sizeof(CCTK_REAL);
- return (-KRANC_GFOFFSET3D(u,0,0,0) + KRANC_GFOFFSET3D(u,0,0,1))*p1odz;
+ ptrdiff_t const cdi CCTK_ATTRIBUTE_UNUSED = sizeof(CCTK_REAL);
+ return (-1.*KRANC_GFOFFSET3D(u,0,0,0) + KRANC_GFOFFSET3D(u,0,0,1))*p1odz;
}
#endif
#ifndef KRANC_DIFF_FUNCTIONS
-# define DiffPlus1(u) ((-KRANC_GFOFFSET3D(u,0,0,0) + KRANC_GFOFFSET3D(u,1,0,0))*p1o1)
+# define DiffPlus1(u) ((-1.*KRANC_GFOFFSET3D(u,0,0,0) + KRANC_GFOFFSET3D(u,1,0,0))*p1o1)
#else
# define DiffPlus1(u) (DiffPlus1_impl(u,p1o1,cdj,cdk))
static CCTK_REAL DiffPlus1_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1o1, ptrdiff_t const cdj, ptrdiff_t const cdk) CCTK_ATTRIBUTE_NOINLINE CCTK_ATTRIBUTE_UNUSED;
static CCTK_REAL DiffPlus1_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1o1, ptrdiff_t const cdj, ptrdiff_t const cdk)
{
- ptrdiff_t const cdi=sizeof(CCTK_REAL);
- return (-KRANC_GFOFFSET3D(u,0,0,0) + KRANC_GFOFFSET3D(u,1,0,0))*p1o1;
+ ptrdiff_t const cdi CCTK_ATTRIBUTE_UNUSED = sizeof(CCTK_REAL);
+ return (-1.*KRANC_GFOFFSET3D(u,0,0,0) + KRANC_GFOFFSET3D(u,1,0,0))*p1o1;
}
#endif
#ifndef KRANC_DIFF_FUNCTIONS
-# define DiffPlus2(u) ((-KRANC_GFOFFSET3D(u,0,0,0) + KRANC_GFOFFSET3D(u,0,1,0))*p1o1)
+# define DiffPlus2(u) ((-1.*KRANC_GFOFFSET3D(u,0,0,0) + KRANC_GFOFFSET3D(u,0,1,0))*p1o1)
#else
# define DiffPlus2(u) (DiffPlus2_impl(u,p1o1,cdj,cdk))
static CCTK_REAL DiffPlus2_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1o1, ptrdiff_t const cdj, ptrdiff_t const cdk) CCTK_ATTRIBUTE_NOINLINE CCTK_ATTRIBUTE_UNUSED;
static CCTK_REAL DiffPlus2_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1o1, ptrdiff_t const cdj, ptrdiff_t const cdk)
{
- ptrdiff_t const cdi=sizeof(CCTK_REAL);
- return (-KRANC_GFOFFSET3D(u,0,0,0) + KRANC_GFOFFSET3D(u,0,1,0))*p1o1;
+ ptrdiff_t const cdi CCTK_ATTRIBUTE_UNUSED = sizeof(CCTK_REAL);
+ return (-1.*KRANC_GFOFFSET3D(u,0,0,0) + KRANC_GFOFFSET3D(u,0,1,0))*p1o1;
}
#endif
#ifndef KRANC_DIFF_FUNCTIONS
-# define DiffPlus3(u) ((-KRANC_GFOFFSET3D(u,0,0,0) + KRANC_GFOFFSET3D(u,0,0,1))*p1o1)
+# define DiffPlus3(u) ((-1.*KRANC_GFOFFSET3D(u,0,0,0) + KRANC_GFOFFSET3D(u,0,0,1))*p1o1)
#else
# define DiffPlus3(u) (DiffPlus3_impl(u,p1o1,cdj,cdk))
static CCTK_REAL DiffPlus3_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1o1, ptrdiff_t const cdj, ptrdiff_t const cdk) CCTK_ATTRIBUTE_NOINLINE CCTK_ATTRIBUTE_UNUSED;
static CCTK_REAL DiffPlus3_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1o1, ptrdiff_t const cdj, ptrdiff_t const cdk)
{
- ptrdiff_t const cdi=sizeof(CCTK_REAL);
- return (-KRANC_GFOFFSET3D(u,0,0,0) + KRANC_GFOFFSET3D(u,0,0,1))*p1o1;
+ ptrdiff_t const cdi CCTK_ATTRIBUTE_UNUSED = sizeof(CCTK_REAL);
+ return (-1.*KRANC_GFOFFSET3D(u,0,0,0) + KRANC_GFOFFSET3D(u,0,0,1))*p1o1;
}
#endif
#ifndef KRANC_DIFF_FUNCTIONS
-# define DiffMinus1(u) ((KRANC_GFOFFSET3D(u,0,0,0) - KRANC_GFOFFSET3D(u,-1,0,0))*p1o1)
+# define DiffMinus1(u) ((KRANC_GFOFFSET3D(u,0,0,0) - 1.*KRANC_GFOFFSET3D(u,-1,0,0))*p1o1)
#else
# define DiffMinus1(u) (DiffMinus1_impl(u,p1o1,cdj,cdk))
static CCTK_REAL DiffMinus1_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1o1, ptrdiff_t const cdj, ptrdiff_t const cdk) CCTK_ATTRIBUTE_NOINLINE CCTK_ATTRIBUTE_UNUSED;
static CCTK_REAL DiffMinus1_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1o1, ptrdiff_t const cdj, ptrdiff_t const cdk)
{
- ptrdiff_t const cdi=sizeof(CCTK_REAL);
- return (KRANC_GFOFFSET3D(u,0,0,0) - KRANC_GFOFFSET3D(u,-1,0,0))*p1o1;
+ ptrdiff_t const cdi CCTK_ATTRIBUTE_UNUSED = sizeof(CCTK_REAL);
+ return (KRANC_GFOFFSET3D(u,0,0,0) - 1.*KRANC_GFOFFSET3D(u,-1,0,0))*p1o1;
}
#endif
#ifndef KRANC_DIFF_FUNCTIONS
-# define DiffMinus2(u) ((KRANC_GFOFFSET3D(u,0,0,0) - KRANC_GFOFFSET3D(u,0,-1,0))*p1o1)
+# define DiffMinus2(u) ((KRANC_GFOFFSET3D(u,0,0,0) - 1.*KRANC_GFOFFSET3D(u,0,-1,0))*p1o1)
#else
# define DiffMinus2(u) (DiffMinus2_impl(u,p1o1,cdj,cdk))
static CCTK_REAL DiffMinus2_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1o1, ptrdiff_t const cdj, ptrdiff_t const cdk) CCTK_ATTRIBUTE_NOINLINE CCTK_ATTRIBUTE_UNUSED;
static CCTK_REAL DiffMinus2_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1o1, ptrdiff_t const cdj, ptrdiff_t const cdk)
{
- ptrdiff_t const cdi=sizeof(CCTK_REAL);
- return (KRANC_GFOFFSET3D(u,0,0,0) - KRANC_GFOFFSET3D(u,0,-1,0))*p1o1;
+ ptrdiff_t const cdi CCTK_ATTRIBUTE_UNUSED = sizeof(CCTK_REAL);
+ return (KRANC_GFOFFSET3D(u,0,0,0) - 1.*KRANC_GFOFFSET3D(u,0,-1,0))*p1o1;
}
#endif
#ifndef KRANC_DIFF_FUNCTIONS
-# define DiffMinus3(u) ((KRANC_GFOFFSET3D(u,0,0,0) - KRANC_GFOFFSET3D(u,0,0,-1))*p1o1)
+# define DiffMinus3(u) ((KRANC_GFOFFSET3D(u,0,0,0) - 1.*KRANC_GFOFFSET3D(u,0,0,-1))*p1o1)
#else
# define DiffMinus3(u) (DiffMinus3_impl(u,p1o1,cdj,cdk))
static CCTK_REAL DiffMinus3_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1o1, ptrdiff_t const cdj, ptrdiff_t const cdk) CCTK_ATTRIBUTE_NOINLINE CCTK_ATTRIBUTE_UNUSED;
static CCTK_REAL DiffMinus3_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1o1, ptrdiff_t const cdj, ptrdiff_t const cdk)
{
- ptrdiff_t const cdi=sizeof(CCTK_REAL);
- return (KRANC_GFOFFSET3D(u,0,0,0) - KRANC_GFOFFSET3D(u,0,0,-1))*p1o1;
+ ptrdiff_t const cdi CCTK_ATTRIBUTE_UNUSED = sizeof(CCTK_REAL);
+ return (KRANC_GFOFFSET3D(u,0,0,0) - 1.*KRANC_GFOFFSET3D(u,0,0,-1))*p1o1;
}
#endif
@@ -113,7 +113,7 @@ static CCTK_REAL DiffMinus3_impl(CCTK_REAL const* restrict const u, CCTK_REAL co
static CCTK_REAL ShiftMinus1_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1o1, ptrdiff_t const cdj, ptrdiff_t const cdk) CCTK_ATTRIBUTE_NOINLINE CCTK_ATTRIBUTE_UNUSED;
static CCTK_REAL ShiftMinus1_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1o1, ptrdiff_t const cdj, ptrdiff_t const cdk)
{
- ptrdiff_t const cdi=sizeof(CCTK_REAL);
+ ptrdiff_t const cdi CCTK_ATTRIBUTE_UNUSED = sizeof(CCTK_REAL);
return KRANC_GFOFFSET3D(u,-1,0,0)*p1o1;
}
#endif
@@ -125,7 +125,7 @@ static CCTK_REAL ShiftMinus1_impl(CCTK_REAL const* restrict const u, CCTK_REAL c
static CCTK_REAL ShiftMinus2_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1o1, ptrdiff_t const cdj, ptrdiff_t const cdk) CCTK_ATTRIBUTE_NOINLINE CCTK_ATTRIBUTE_UNUSED;
static CCTK_REAL ShiftMinus2_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1o1, ptrdiff_t const cdj, ptrdiff_t const cdk)
{
- ptrdiff_t const cdi=sizeof(CCTK_REAL);
+ ptrdiff_t const cdi CCTK_ATTRIBUTE_UNUSED = sizeof(CCTK_REAL);
return KRANC_GFOFFSET3D(u,0,-1,0)*p1o1;
}
#endif
@@ -137,7 +137,7 @@ static CCTK_REAL ShiftMinus2_impl(CCTK_REAL const* restrict const u, CCTK_REAL c
static CCTK_REAL ShiftMinus3_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1o1, ptrdiff_t const cdj, ptrdiff_t const cdk) CCTK_ATTRIBUTE_NOINLINE CCTK_ATTRIBUTE_UNUSED;
static CCTK_REAL ShiftMinus3_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1o1, ptrdiff_t const cdj, ptrdiff_t const cdk)
{
- ptrdiff_t const cdi=sizeof(CCTK_REAL);
+ ptrdiff_t const cdi CCTK_ATTRIBUTE_UNUSED = sizeof(CCTK_REAL);
return KRANC_GFOFFSET3D(u,0,0,-1)*p1o1;
}
#endif
diff --git a/Examples/Burgers/src/RegisterMoL.cc b/Examples/Burgers/src/RegisterMoL.cc
index 3863790..cb22a57 100644
--- a/Examples/Burgers/src/RegisterMoL.cc
+++ b/Examples/Burgers/src/RegisterMoL.cc
@@ -9,7 +9,7 @@ extern "C" void Burgers_RegisterVars(CCTK_ARGUMENTS)
DECLARE_CCTK_ARGUMENTS;
DECLARE_CCTK_PARAMETERS;
- CCTK_INT ierr = 0;
+ CCTK_INT ierr CCTK_ATTRIBUTE_UNUSED = 0;
/* Register all the evolved grid functions with MoL */
ierr += MoLRegisterEvolved(CCTK_VarIndex("Burgers::u"), CCTK_VarIndex("Burgers::urhs"));
diff --git a/Examples/Burgers/src/Startup.cc b/Examples/Burgers/src/Startup.cc
index f092a5e..0102836 100644
--- a/Examples/Burgers/src/Startup.cc
+++ b/Examples/Burgers/src/Startup.cc
@@ -4,7 +4,7 @@
extern "C" int Burgers_Startup(void)
{
- const char * banner = "Burgers";
+ const char * banner CCTK_ATTRIBUTE_UNUSED = "Burgers";
CCTK_RegisterBanner(banner);
return 0;
}
diff --git a/Examples/Burgers/src/burgers_flux_1.cc b/Examples/Burgers/src/burgers_flux_1.cc
index a0a2300..acd9e91 100644
--- a/Examples/Burgers/src/burgers_flux_1.cc
+++ b/Examples/Burgers/src/burgers_flux_1.cc
@@ -27,7 +27,7 @@ extern "C" void burgers_flux_1_SelectBCs(CCTK_ARGUMENTS)
DECLARE_CCTK_ARGUMENTS;
DECLARE_CCTK_PARAMETERS;
- CCTK_INT ierr = 0;
+ CCTK_INT ierr CCTK_ATTRIBUTE_UNUSED = 0;
ierr = Boundary_SelectGroupForBC(cctkGH, CCTK_ALL_FACES, GenericFD_GetBoundaryWidth(cctkGH), -1 /* no table */, "Burgers::uF_group","flat");
if (ierr < 0)
CCTK_WARN(1, "Failed to register flat BC for Burgers::uF_group.");
@@ -43,34 +43,34 @@ static void burgers_flux_1_Body(cGH const * restrict const cctkGH, int const dir
/* Include user-supplied include files */
/* Initialise finite differencing variables */
- ptrdiff_t const di = 1;
- ptrdiff_t const dj = CCTK_GFINDEX3D(cctkGH,0,1,0) - CCTK_GFINDEX3D(cctkGH,0,0,0);
- ptrdiff_t const dk = CCTK_GFINDEX3D(cctkGH,0,0,1) - CCTK_GFINDEX3D(cctkGH,0,0,0);
- ptrdiff_t const cdi = sizeof(CCTK_REAL) * di;
- ptrdiff_t const cdj = sizeof(CCTK_REAL) * dj;
- ptrdiff_t const cdk = sizeof(CCTK_REAL) * dk;
- CCTK_REAL const dx = ToReal(CCTK_DELTA_SPACE(0));
- CCTK_REAL const dy = ToReal(CCTK_DELTA_SPACE(1));
- CCTK_REAL const dz = ToReal(CCTK_DELTA_SPACE(2));
- CCTK_REAL const dt = ToReal(CCTK_DELTA_TIME);
- CCTK_REAL const t = ToReal(cctk_time);
- CCTK_REAL const dxi = INV(dx);
- CCTK_REAL const dyi = INV(dy);
- CCTK_REAL const dzi = INV(dz);
- CCTK_REAL const khalf = 0.5;
- CCTK_REAL const kthird = 1/3.0;
- CCTK_REAL const ktwothird = 2.0/3.0;
- CCTK_REAL const kfourthird = 4.0/3.0;
- CCTK_REAL const keightthird = 8.0/3.0;
- CCTK_REAL const hdxi = 0.5 * dxi;
- CCTK_REAL const hdyi = 0.5 * dyi;
- CCTK_REAL const hdzi = 0.5 * dzi;
+ ptrdiff_t /*const*/ di CCTK_ATTRIBUTE_UNUSED = 1;
+ ptrdiff_t /*const*/ dj CCTK_ATTRIBUTE_UNUSED = CCTK_GFINDEX3D(cctkGH,0,1,0) - CCTK_GFINDEX3D(cctkGH,0,0,0);
+ ptrdiff_t /*const*/ dk CCTK_ATTRIBUTE_UNUSED = CCTK_GFINDEX3D(cctkGH,0,0,1) - CCTK_GFINDEX3D(cctkGH,0,0,0);
+ ptrdiff_t /*const*/ cdi CCTK_ATTRIBUTE_UNUSED = sizeof(CCTK_REAL) * di;
+ ptrdiff_t /*const*/ cdj CCTK_ATTRIBUTE_UNUSED = sizeof(CCTK_REAL) * dj;
+ ptrdiff_t /*const*/ cdk CCTK_ATTRIBUTE_UNUSED = sizeof(CCTK_REAL) * dk;
+ CCTK_REAL /*const*/ dx CCTK_ATTRIBUTE_UNUSED = ToReal(CCTK_DELTA_SPACE(0));
+ CCTK_REAL /*const*/ dy CCTK_ATTRIBUTE_UNUSED = ToReal(CCTK_DELTA_SPACE(1));
+ CCTK_REAL /*const*/ dz CCTK_ATTRIBUTE_UNUSED = ToReal(CCTK_DELTA_SPACE(2));
+ CCTK_REAL /*const*/ dt CCTK_ATTRIBUTE_UNUSED = ToReal(CCTK_DELTA_TIME);
+ CCTK_REAL /*const*/ t CCTK_ATTRIBUTE_UNUSED = ToReal(cctk_time);
+ CCTK_REAL /*const*/ dxi CCTK_ATTRIBUTE_UNUSED = INV(dx);
+ CCTK_REAL /*const*/ dyi CCTK_ATTRIBUTE_UNUSED = INV(dy);
+ CCTK_REAL /*const*/ dzi CCTK_ATTRIBUTE_UNUSED = INV(dz);
+ CCTK_REAL /*const*/ khalf CCTK_ATTRIBUTE_UNUSED = 0.5;
+ CCTK_REAL /*const*/ kthird CCTK_ATTRIBUTE_UNUSED = 1/3.0;
+ CCTK_REAL /*const*/ ktwothird CCTK_ATTRIBUTE_UNUSED = 2.0/3.0;
+ CCTK_REAL /*const*/ kfourthird CCTK_ATTRIBUTE_UNUSED = 4.0/3.0;
+ CCTK_REAL /*const*/ keightthird CCTK_ATTRIBUTE_UNUSED = 8.0/3.0;
+ CCTK_REAL /*const*/ hdxi CCTK_ATTRIBUTE_UNUSED = 0.5 * dxi;
+ CCTK_REAL /*const*/ hdyi CCTK_ATTRIBUTE_UNUSED = 0.5 * dyi;
+ CCTK_REAL /*const*/ hdzi CCTK_ATTRIBUTE_UNUSED = 0.5 * dzi;
/* Initialize predefined quantities */
- CCTK_REAL const p1o1 = 1;
- CCTK_REAL const p1odx = INV(dx);
- CCTK_REAL const p1ody = INV(dy);
- CCTK_REAL const p1odz = INV(dz);
+ CCTK_REAL /*const*/ p1o1 CCTK_ATTRIBUTE_UNUSED = 1.;
+ CCTK_REAL /*const*/ p1odx CCTK_ATTRIBUTE_UNUSED = INV(dx);
+ CCTK_REAL /*const*/ p1ody CCTK_ATTRIBUTE_UNUSED = INV(dy);
+ CCTK_REAL /*const*/ p1odz CCTK_ATTRIBUTE_UNUSED = INV(dz);
/* Assign local copies of arrays functions */
@@ -86,22 +86,23 @@ static void burgers_flux_1_Body(cGH const * restrict const cctkGH, int const dir
i,j,k, imin[0],imin[1],imin[2], imax[0],imax[1],imax[2],
cctk_ash[0],cctk_ash[1],cctk_ash[2])
{
- ptrdiff_t const index = di*i + dj*j + dk*k;
+ ptrdiff_t /*const*/ index CCTK_ATTRIBUTE_UNUSED = di*i + dj*j + dk*k;
/* Assign local copies of grid functions */
- CCTK_REAL uLeftL = uLeft[index];
- CCTK_REAL uRL = uR[index];
+ CCTK_REAL uLeftL CCTK_ATTRIBUTE_UNUSED = uLeft[index];
+ CCTK_REAL uRL CCTK_ATTRIBUTE_UNUSED = uR[index];
/* Include user supplied include files */
/* Precompute derivatives */
- CCTK_REAL const ShiftMinus1uR = ShiftMinus1(&uR[index]);
+ CCTK_REAL /*const*/ ShiftMinus1uR CCTK_ATTRIBUTE_UNUSED = ShiftMinus1(&uR[index]);
/* Calculate temporaries and grid functions */
- CCTK_REAL uFL = 0.25*(SQR(ShiftMinus1uR) + uLeftL*(uLeftL -
- 2*ToReal(alpha)) + 2*ShiftMinus1uR*ToReal(alpha));
+ CCTK_REAL CCTK_ATTRIBUTE_UNUSED uFL = 0.25*(SQR(ShiftMinus1uR) +
+ uLeftL*(uLeftL - 2.*ToReal(alpha)) +
+ 2.*ShiftMinus1uR*ToReal(alpha));
/* Copy local copies back to grid functions */
uF[index] = uFL;
diff --git a/Examples/Burgers/src/burgers_initial_shock.cc b/Examples/Burgers/src/burgers_initial_shock.cc
index 22dc55a..1551662 100644
--- a/Examples/Burgers/src/burgers_initial_shock.cc
+++ b/Examples/Burgers/src/burgers_initial_shock.cc
@@ -31,34 +31,34 @@ static void burgers_initial_shock_Body(cGH const * restrict const cctkGH, int co
/* Include user-supplied include files */
/* Initialise finite differencing variables */
- ptrdiff_t const di = 1;
- ptrdiff_t const dj = CCTK_GFINDEX3D(cctkGH,0,1,0) - CCTK_GFINDEX3D(cctkGH,0,0,0);
- ptrdiff_t const dk = CCTK_GFINDEX3D(cctkGH,0,0,1) - CCTK_GFINDEX3D(cctkGH,0,0,0);
- ptrdiff_t const cdi = sizeof(CCTK_REAL) * di;
- ptrdiff_t const cdj = sizeof(CCTK_REAL) * dj;
- ptrdiff_t const cdk = sizeof(CCTK_REAL) * dk;
- CCTK_REAL const dx = ToReal(CCTK_DELTA_SPACE(0));
- CCTK_REAL const dy = ToReal(CCTK_DELTA_SPACE(1));
- CCTK_REAL const dz = ToReal(CCTK_DELTA_SPACE(2));
- CCTK_REAL const dt = ToReal(CCTK_DELTA_TIME);
- CCTK_REAL const t = ToReal(cctk_time);
- CCTK_REAL const dxi = INV(dx);
- CCTK_REAL const dyi = INV(dy);
- CCTK_REAL const dzi = INV(dz);
- CCTK_REAL const khalf = 0.5;
- CCTK_REAL const kthird = 1/3.0;
- CCTK_REAL const ktwothird = 2.0/3.0;
- CCTK_REAL const kfourthird = 4.0/3.0;
- CCTK_REAL const keightthird = 8.0/3.0;
- CCTK_REAL const hdxi = 0.5 * dxi;
- CCTK_REAL const hdyi = 0.5 * dyi;
- CCTK_REAL const hdzi = 0.5 * dzi;
+ ptrdiff_t /*const*/ di CCTK_ATTRIBUTE_UNUSED = 1;
+ ptrdiff_t /*const*/ dj CCTK_ATTRIBUTE_UNUSED = CCTK_GFINDEX3D(cctkGH,0,1,0) - CCTK_GFINDEX3D(cctkGH,0,0,0);
+ ptrdiff_t /*const*/ dk CCTK_ATTRIBUTE_UNUSED = CCTK_GFINDEX3D(cctkGH,0,0,1) - CCTK_GFINDEX3D(cctkGH,0,0,0);
+ ptrdiff_t /*const*/ cdi CCTK_ATTRIBUTE_UNUSED = sizeof(CCTK_REAL) * di;
+ ptrdiff_t /*const*/ cdj CCTK_ATTRIBUTE_UNUSED = sizeof(CCTK_REAL) * dj;
+ ptrdiff_t /*const*/ cdk CCTK_ATTRIBUTE_UNUSED = sizeof(CCTK_REAL) * dk;
+ CCTK_REAL /*const*/ dx CCTK_ATTRIBUTE_UNUSED = ToReal(CCTK_DELTA_SPACE(0));
+ CCTK_REAL /*const*/ dy CCTK_ATTRIBUTE_UNUSED = ToReal(CCTK_DELTA_SPACE(1));
+ CCTK_REAL /*const*/ dz CCTK_ATTRIBUTE_UNUSED = ToReal(CCTK_DELTA_SPACE(2));
+ CCTK_REAL /*const*/ dt CCTK_ATTRIBUTE_UNUSED = ToReal(CCTK_DELTA_TIME);
+ CCTK_REAL /*const*/ t CCTK_ATTRIBUTE_UNUSED = ToReal(cctk_time);
+ CCTK_REAL /*const*/ dxi CCTK_ATTRIBUTE_UNUSED = INV(dx);
+ CCTK_REAL /*const*/ dyi CCTK_ATTRIBUTE_UNUSED = INV(dy);
+ CCTK_REAL /*const*/ dzi CCTK_ATTRIBUTE_UNUSED = INV(dz);
+ CCTK_REAL /*const*/ khalf CCTK_ATTRIBUTE_UNUSED = 0.5;
+ CCTK_REAL /*const*/ kthird CCTK_ATTRIBUTE_UNUSED = 1/3.0;
+ CCTK_REAL /*const*/ ktwothird CCTK_ATTRIBUTE_UNUSED = 2.0/3.0;
+ CCTK_REAL /*const*/ kfourthird CCTK_ATTRIBUTE_UNUSED = 4.0/3.0;
+ CCTK_REAL /*const*/ keightthird CCTK_ATTRIBUTE_UNUSED = 8.0/3.0;
+ CCTK_REAL /*const*/ hdxi CCTK_ATTRIBUTE_UNUSED = 0.5 * dxi;
+ CCTK_REAL /*const*/ hdyi CCTK_ATTRIBUTE_UNUSED = 0.5 * dyi;
+ CCTK_REAL /*const*/ hdzi CCTK_ATTRIBUTE_UNUSED = 0.5 * dzi;
/* Initialize predefined quantities */
- CCTK_REAL const p1o1 = 1;
- CCTK_REAL const p1odx = INV(dx);
- CCTK_REAL const p1ody = INV(dy);
- CCTK_REAL const p1odz = INV(dz);
+ CCTK_REAL /*const*/ p1o1 CCTK_ATTRIBUTE_UNUSED = 1.;
+ CCTK_REAL /*const*/ p1odx CCTK_ATTRIBUTE_UNUSED = INV(dx);
+ CCTK_REAL /*const*/ p1ody CCTK_ATTRIBUTE_UNUSED = INV(dy);
+ CCTK_REAL /*const*/ p1odz CCTK_ATTRIBUTE_UNUSED = INV(dz);
/* Assign local copies of arrays functions */
@@ -74,11 +74,11 @@ static void burgers_initial_shock_Body(cGH const * restrict const cctkGH, int co
i,j,k, imin[0],imin[1],imin[2], imax[0],imax[1],imax[2],
cctk_ash[0],cctk_ash[1],cctk_ash[2])
{
- ptrdiff_t const index = di*i + dj*j + dk*k;
+ ptrdiff_t /*const*/ index CCTK_ATTRIBUTE_UNUSED = di*i + dj*j + dk*k;
/* Assign local copies of grid functions */
- CCTK_REAL xL = x[index];
+ CCTK_REAL xL CCTK_ATTRIBUTE_UNUSED = x[index];
/* Include user supplied include files */
@@ -86,8 +86,8 @@ static void burgers_initial_shock_Body(cGH const * restrict const cctkGH, int co
/* Precompute derivatives */
/* Calculate temporaries and grid functions */
- CCTK_REAL uL = ToReal(uL0) + StepFunction(-0.5 + xL)*(-ToReal(uL0) +
- ToReal(uR0));
+ CCTK_REAL CCTK_ATTRIBUTE_UNUSED uL = ToReal(uL0) + StepFunction(-0.5 +
+ xL)*(-1.*ToReal(uL0) + ToReal(uR0));
/* Copy local copies back to grid functions */
u[index] = uL;
diff --git a/Examples/Burgers/src/burgers_initial_sine.cc b/Examples/Burgers/src/burgers_initial_sine.cc
index 0e5daaa..e5e23b9 100644
--- a/Examples/Burgers/src/burgers_initial_sine.cc
+++ b/Examples/Burgers/src/burgers_initial_sine.cc
@@ -31,34 +31,34 @@ static void burgers_initial_sine_Body(cGH const * restrict const cctkGH, int con
/* Include user-supplied include files */
/* Initialise finite differencing variables */
- ptrdiff_t const di = 1;
- ptrdiff_t const dj = CCTK_GFINDEX3D(cctkGH,0,1,0) - CCTK_GFINDEX3D(cctkGH,0,0,0);
- ptrdiff_t const dk = CCTK_GFINDEX3D(cctkGH,0,0,1) - CCTK_GFINDEX3D(cctkGH,0,0,0);
- ptrdiff_t const cdi = sizeof(CCTK_REAL) * di;
- ptrdiff_t const cdj = sizeof(CCTK_REAL) * dj;
- ptrdiff_t const cdk = sizeof(CCTK_REAL) * dk;
- CCTK_REAL const dx = ToReal(CCTK_DELTA_SPACE(0));
- CCTK_REAL const dy = ToReal(CCTK_DELTA_SPACE(1));
- CCTK_REAL const dz = ToReal(CCTK_DELTA_SPACE(2));
- CCTK_REAL const dt = ToReal(CCTK_DELTA_TIME);
- CCTK_REAL const t = ToReal(cctk_time);
- CCTK_REAL const dxi = INV(dx);
- CCTK_REAL const dyi = INV(dy);
- CCTK_REAL const dzi = INV(dz);
- CCTK_REAL const khalf = 0.5;
- CCTK_REAL const kthird = 1/3.0;
- CCTK_REAL const ktwothird = 2.0/3.0;
- CCTK_REAL const kfourthird = 4.0/3.0;
- CCTK_REAL const keightthird = 8.0/3.0;
- CCTK_REAL const hdxi = 0.5 * dxi;
- CCTK_REAL const hdyi = 0.5 * dyi;
- CCTK_REAL const hdzi = 0.5 * dzi;
+ ptrdiff_t /*const*/ di CCTK_ATTRIBUTE_UNUSED = 1;
+ ptrdiff_t /*const*/ dj CCTK_ATTRIBUTE_UNUSED = CCTK_GFINDEX3D(cctkGH,0,1,0) - CCTK_GFINDEX3D(cctkGH,0,0,0);
+ ptrdiff_t /*const*/ dk CCTK_ATTRIBUTE_UNUSED = CCTK_GFINDEX3D(cctkGH,0,0,1) - CCTK_GFINDEX3D(cctkGH,0,0,0);
+ ptrdiff_t /*const*/ cdi CCTK_ATTRIBUTE_UNUSED = sizeof(CCTK_REAL) * di;
+ ptrdiff_t /*const*/ cdj CCTK_ATTRIBUTE_UNUSED = sizeof(CCTK_REAL) * dj;
+ ptrdiff_t /*const*/ cdk CCTK_ATTRIBUTE_UNUSED = sizeof(CCTK_REAL) * dk;
+ CCTK_REAL /*const*/ dx CCTK_ATTRIBUTE_UNUSED = ToReal(CCTK_DELTA_SPACE(0));
+ CCTK_REAL /*const*/ dy CCTK_ATTRIBUTE_UNUSED = ToReal(CCTK_DELTA_SPACE(1));
+ CCTK_REAL /*const*/ dz CCTK_ATTRIBUTE_UNUSED = ToReal(CCTK_DELTA_SPACE(2));
+ CCTK_REAL /*const*/ dt CCTK_ATTRIBUTE_UNUSED = ToReal(CCTK_DELTA_TIME);
+ CCTK_REAL /*const*/ t CCTK_ATTRIBUTE_UNUSED = ToReal(cctk_time);
+ CCTK_REAL /*const*/ dxi CCTK_ATTRIBUTE_UNUSED = INV(dx);
+ CCTK_REAL /*const*/ dyi CCTK_ATTRIBUTE_UNUSED = INV(dy);
+ CCTK_REAL /*const*/ dzi CCTK_ATTRIBUTE_UNUSED = INV(dz);
+ CCTK_REAL /*const*/ khalf CCTK_ATTRIBUTE_UNUSED = 0.5;
+ CCTK_REAL /*const*/ kthird CCTK_ATTRIBUTE_UNUSED = 1/3.0;
+ CCTK_REAL /*const*/ ktwothird CCTK_ATTRIBUTE_UNUSED = 2.0/3.0;
+ CCTK_REAL /*const*/ kfourthird CCTK_ATTRIBUTE_UNUSED = 4.0/3.0;
+ CCTK_REAL /*const*/ keightthird CCTK_ATTRIBUTE_UNUSED = 8.0/3.0;
+ CCTK_REAL /*const*/ hdxi CCTK_ATTRIBUTE_UNUSED = 0.5 * dxi;
+ CCTK_REAL /*const*/ hdyi CCTK_ATTRIBUTE_UNUSED = 0.5 * dyi;
+ CCTK_REAL /*const*/ hdzi CCTK_ATTRIBUTE_UNUSED = 0.5 * dzi;
/* Initialize predefined quantities */
- CCTK_REAL const p1o1 = 1;
- CCTK_REAL const p1odx = INV(dx);
- CCTK_REAL const p1ody = INV(dy);
- CCTK_REAL const p1odz = INV(dz);
+ CCTK_REAL /*const*/ p1o1 CCTK_ATTRIBUTE_UNUSED = 1.;
+ CCTK_REAL /*const*/ p1odx CCTK_ATTRIBUTE_UNUSED = INV(dx);
+ CCTK_REAL /*const*/ p1ody CCTK_ATTRIBUTE_UNUSED = INV(dy);
+ CCTK_REAL /*const*/ p1odz CCTK_ATTRIBUTE_UNUSED = INV(dz);
/* Assign local copies of arrays functions */
@@ -74,11 +74,11 @@ static void burgers_initial_sine_Body(cGH const * restrict const cctkGH, int con
i,j,k, imin[0],imin[1],imin[2], imax[0],imax[1],imax[2],
cctk_ash[0],cctk_ash[1],cctk_ash[2])
{
- ptrdiff_t const index = di*i + dj*j + dk*k;
+ ptrdiff_t /*const*/ index CCTK_ATTRIBUTE_UNUSED = di*i + dj*j + dk*k;
/* Assign local copies of grid functions */
- CCTK_REAL xL = x[index];
+ CCTK_REAL xL CCTK_ATTRIBUTE_UNUSED = x[index];
/* Include user supplied include files */
@@ -86,7 +86,8 @@ static void burgers_initial_sine_Body(cGH const * restrict const cctkGH, int con
/* Precompute derivatives */
/* Calculate temporaries and grid functions */
- CCTK_REAL uL = 1 + sin(2*xL*Pi)*ToReal(amp);
+ CCTK_REAL CCTK_ATTRIBUTE_UNUSED uL = 1. +
+ sin(6.283185307179586*xL)*ToReal(amp);
/* Copy local copies back to grid functions */
u[index] = uL;
diff --git a/Examples/Burgers/src/burgers_reconstruct_1.cc b/Examples/Burgers/src/burgers_reconstruct_1.cc
index 2a14321..3ce2fc9 100644
--- a/Examples/Burgers/src/burgers_reconstruct_1.cc
+++ b/Examples/Burgers/src/burgers_reconstruct_1.cc
@@ -27,7 +27,7 @@ extern "C" void burgers_reconstruct_1_SelectBCs(CCTK_ARGUMENTS)
DECLARE_CCTK_ARGUMENTS;
DECLARE_CCTK_PARAMETERS;
- CCTK_INT ierr = 0;
+ CCTK_INT ierr CCTK_ATTRIBUTE_UNUSED = 0;
ierr = Boundary_SelectGroupForBC(cctkGH, CCTK_ALL_FACES, GenericFD_GetBoundaryWidth(cctkGH), -1 /* no table */, "Burgers::uLeft_group","flat");
if (ierr < 0)
CCTK_WARN(1, "Failed to register flat BC for Burgers::uLeft_group.");
@@ -46,34 +46,34 @@ static void burgers_reconstruct_1_Body(cGH const * restrict const cctkGH, int co
/* Include user-supplied include files */
/* Initialise finite differencing variables */
- ptrdiff_t const di = 1;
- ptrdiff_t const dj = CCTK_GFINDEX3D(cctkGH,0,1,0) - CCTK_GFINDEX3D(cctkGH,0,0,0);
- ptrdiff_t const dk = CCTK_GFINDEX3D(cctkGH,0,0,1) - CCTK_GFINDEX3D(cctkGH,0,0,0);
- ptrdiff_t const cdi = sizeof(CCTK_REAL) * di;
- ptrdiff_t const cdj = sizeof(CCTK_REAL) * dj;
- ptrdiff_t const cdk = sizeof(CCTK_REAL) * dk;
- CCTK_REAL const dx = ToReal(CCTK_DELTA_SPACE(0));
- CCTK_REAL const dy = ToReal(CCTK_DELTA_SPACE(1));
- CCTK_REAL const dz = ToReal(CCTK_DELTA_SPACE(2));
- CCTK_REAL const dt = ToReal(CCTK_DELTA_TIME);
- CCTK_REAL const t = ToReal(cctk_time);
- CCTK_REAL const dxi = INV(dx);
- CCTK_REAL const dyi = INV(dy);
- CCTK_REAL const dzi = INV(dz);
- CCTK_REAL const khalf = 0.5;
- CCTK_REAL const kthird = 1/3.0;
- CCTK_REAL const ktwothird = 2.0/3.0;
- CCTK_REAL const kfourthird = 4.0/3.0;
- CCTK_REAL const keightthird = 8.0/3.0;
- CCTK_REAL const hdxi = 0.5 * dxi;
- CCTK_REAL const hdyi = 0.5 * dyi;
- CCTK_REAL const hdzi = 0.5 * dzi;
+ ptrdiff_t /*const*/ di CCTK_ATTRIBUTE_UNUSED = 1;
+ ptrdiff_t /*const*/ dj CCTK_ATTRIBUTE_UNUSED = CCTK_GFINDEX3D(cctkGH,0,1,0) - CCTK_GFINDEX3D(cctkGH,0,0,0);
+ ptrdiff_t /*const*/ dk CCTK_ATTRIBUTE_UNUSED = CCTK_GFINDEX3D(cctkGH,0,0,1) - CCTK_GFINDEX3D(cctkGH,0,0,0);
+ ptrdiff_t /*const*/ cdi CCTK_ATTRIBUTE_UNUSED = sizeof(CCTK_REAL) * di;
+ ptrdiff_t /*const*/ cdj CCTK_ATTRIBUTE_UNUSED = sizeof(CCTK_REAL) * dj;
+ ptrdiff_t /*const*/ cdk CCTK_ATTRIBUTE_UNUSED = sizeof(CCTK_REAL) * dk;
+ CCTK_REAL /*const*/ dx CCTK_ATTRIBUTE_UNUSED = ToReal(CCTK_DELTA_SPACE(0));
+ CCTK_REAL /*const*/ dy CCTK_ATTRIBUTE_UNUSED = ToReal(CCTK_DELTA_SPACE(1));
+ CCTK_REAL /*const*/ dz CCTK_ATTRIBUTE_UNUSED = ToReal(CCTK_DELTA_SPACE(2));
+ CCTK_REAL /*const*/ dt CCTK_ATTRIBUTE_UNUSED = ToReal(CCTK_DELTA_TIME);
+ CCTK_REAL /*const*/ t CCTK_ATTRIBUTE_UNUSED = ToReal(cctk_time);
+ CCTK_REAL /*const*/ dxi CCTK_ATTRIBUTE_UNUSED = INV(dx);
+ CCTK_REAL /*const*/ dyi CCTK_ATTRIBUTE_UNUSED = INV(dy);
+ CCTK_REAL /*const*/ dzi CCTK_ATTRIBUTE_UNUSED = INV(dz);
+ CCTK_REAL /*const*/ khalf CCTK_ATTRIBUTE_UNUSED = 0.5;
+ CCTK_REAL /*const*/ kthird CCTK_ATTRIBUTE_UNUSED = 1/3.0;
+ CCTK_REAL /*const*/ ktwothird CCTK_ATTRIBUTE_UNUSED = 2.0/3.0;
+ CCTK_REAL /*const*/ kfourthird CCTK_ATTRIBUTE_UNUSED = 4.0/3.0;
+ CCTK_REAL /*const*/ keightthird CCTK_ATTRIBUTE_UNUSED = 8.0/3.0;
+ CCTK_REAL /*const*/ hdxi CCTK_ATTRIBUTE_UNUSED = 0.5 * dxi;
+ CCTK_REAL /*const*/ hdyi CCTK_ATTRIBUTE_UNUSED = 0.5 * dyi;
+ CCTK_REAL /*const*/ hdzi CCTK_ATTRIBUTE_UNUSED = 0.5 * dzi;
/* Initialize predefined quantities */
- CCTK_REAL const p1o1 = 1;
- CCTK_REAL const p1odx = INV(dx);
- CCTK_REAL const p1ody = INV(dy);
- CCTK_REAL const p1odz = INV(dz);
+ CCTK_REAL /*const*/ p1o1 CCTK_ATTRIBUTE_UNUSED = 1.;
+ CCTK_REAL /*const*/ p1odx CCTK_ATTRIBUTE_UNUSED = INV(dx);
+ CCTK_REAL /*const*/ p1ody CCTK_ATTRIBUTE_UNUSED = INV(dy);
+ CCTK_REAL /*const*/ p1odz CCTK_ATTRIBUTE_UNUSED = INV(dz);
/* Assign local copies of arrays functions */
@@ -89,30 +89,30 @@ static void burgers_reconstruct_1_Body(cGH const * restrict const cctkGH, int co
i,j,k, imin[0],imin[1],imin[2], imax[0],imax[1],imax[2],
cctk_ash[0],cctk_ash[1],cctk_ash[2])
{
- ptrdiff_t const index = di*i + dj*j + dk*k;
+ ptrdiff_t /*const*/ index CCTK_ATTRIBUTE_UNUSED = di*i + dj*j + dk*k;
/* Assign local copies of grid functions */
- CCTK_REAL uL = u[index];
+ CCTK_REAL uL CCTK_ATTRIBUTE_UNUSED = u[index];
/* Include user supplied include files */
/* Precompute derivatives */
- CCTK_REAL const DiffPlus1u = DiffPlus1(&u[index]);
- CCTK_REAL const DiffMinus1u = DiffMinus1(&u[index]);
+ CCTK_REAL /*const*/ DiffPlus1u CCTK_ATTRIBUTE_UNUSED = DiffPlus1(&u[index]);
+ CCTK_REAL /*const*/ DiffMinus1u CCTK_ATTRIBUTE_UNUSED = DiffMinus1(&u[index]);
/* Calculate temporaries and grid functions */
- CCTK_REAL slopeL = DiffMinus1u;
+ CCTK_REAL CCTK_ATTRIBUTE_UNUSED slopeL = DiffMinus1u;
- CCTK_REAL slopeR = DiffPlus1u;
+ CCTK_REAL CCTK_ATTRIBUTE_UNUSED slopeR = DiffPlus1u;
- CCTK_REAL slope = IfThen(slopeL*slopeR < 0,0,IfThen(fabs(slopeL) <
- fabs(slopeR),slopeL,slopeR));
+ CCTK_REAL CCTK_ATTRIBUTE_UNUSED slope = IfThen(slopeL*slopeR <
+ 0.,0.,IfThen(fabs(slopeL) < fabs(slopeR),slopeL,slopeR));
- CCTK_REAL uLeftL = uL - 0.5*slope;
+ CCTK_REAL CCTK_ATTRIBUTE_UNUSED uLeftL = uL - 0.5*slope;
- CCTK_REAL uRL = uL + 0.5*slope;
+ CCTK_REAL CCTK_ATTRIBUTE_UNUSED uRL = uL + 0.5*slope;
/* Copy local copies back to grid functions */
uLeft[index] = uLeftL;
diff --git a/Examples/Burgers/src/burgers_rhs_1.cc b/Examples/Burgers/src/burgers_rhs_1.cc
index 644aca7..409831a 100644
--- a/Examples/Burgers/src/burgers_rhs_1.cc
+++ b/Examples/Burgers/src/burgers_rhs_1.cc
@@ -27,7 +27,7 @@ extern "C" void burgers_rhs_1_SelectBCs(CCTK_ARGUMENTS)
DECLARE_CCTK_ARGUMENTS;
DECLARE_CCTK_PARAMETERS;
- CCTK_INT ierr = 0;
+ CCTK_INT ierr CCTK_ATTRIBUTE_UNUSED = 0;
ierr = Boundary_SelectGroupForBC(cctkGH, CCTK_ALL_FACES, GenericFD_GetBoundaryWidth(cctkGH), -1 /* no table */, "Burgers::u_grouprhs","flat");
if (ierr < 0)
CCTK_WARN(1, "Failed to register flat BC for Burgers::u_grouprhs.");
@@ -43,34 +43,34 @@ static void burgers_rhs_1_Body(cGH const * restrict const cctkGH, int const dir,
/* Include user-supplied include files */
/* Initialise finite differencing variables */
- ptrdiff_t const di = 1;
- ptrdiff_t const dj = CCTK_GFINDEX3D(cctkGH,0,1,0) - CCTK_GFINDEX3D(cctkGH,0,0,0);
- ptrdiff_t const dk = CCTK_GFINDEX3D(cctkGH,0,0,1) - CCTK_GFINDEX3D(cctkGH,0,0,0);
- ptrdiff_t const cdi = sizeof(CCTK_REAL) * di;
- ptrdiff_t const cdj = sizeof(CCTK_REAL) * dj;
- ptrdiff_t const cdk = sizeof(CCTK_REAL) * dk;
- CCTK_REAL const dx = ToReal(CCTK_DELTA_SPACE(0));
- CCTK_REAL const dy = ToReal(CCTK_DELTA_SPACE(1));
- CCTK_REAL const dz = ToReal(CCTK_DELTA_SPACE(2));
- CCTK_REAL const dt = ToReal(CCTK_DELTA_TIME);
- CCTK_REAL const t = ToReal(cctk_time);
- CCTK_REAL const dxi = INV(dx);
- CCTK_REAL const dyi = INV(dy);
- CCTK_REAL const dzi = INV(dz);
- CCTK_REAL const khalf = 0.5;
- CCTK_REAL const kthird = 1/3.0;
- CCTK_REAL const ktwothird = 2.0/3.0;
- CCTK_REAL const kfourthird = 4.0/3.0;
- CCTK_REAL const keightthird = 8.0/3.0;
- CCTK_REAL const hdxi = 0.5 * dxi;
- CCTK_REAL const hdyi = 0.5 * dyi;
- CCTK_REAL const hdzi = 0.5 * dzi;
+ ptrdiff_t /*const*/ di CCTK_ATTRIBUTE_UNUSED = 1;
+ ptrdiff_t /*const*/ dj CCTK_ATTRIBUTE_UNUSED = CCTK_GFINDEX3D(cctkGH,0,1,0) - CCTK_GFINDEX3D(cctkGH,0,0,0);
+ ptrdiff_t /*const*/ dk CCTK_ATTRIBUTE_UNUSED = CCTK_GFINDEX3D(cctkGH,0,0,1) - CCTK_GFINDEX3D(cctkGH,0,0,0);
+ ptrdiff_t /*const*/ cdi CCTK_ATTRIBUTE_UNUSED = sizeof(CCTK_REAL) * di;
+ ptrdiff_t /*const*/ cdj CCTK_ATTRIBUTE_UNUSED = sizeof(CCTK_REAL) * dj;
+ ptrdiff_t /*const*/ cdk CCTK_ATTRIBUTE_UNUSED = sizeof(CCTK_REAL) * dk;
+ CCTK_REAL /*const*/ dx CCTK_ATTRIBUTE_UNUSED = ToReal(CCTK_DELTA_SPACE(0));
+ CCTK_REAL /*const*/ dy CCTK_ATTRIBUTE_UNUSED = ToReal(CCTK_DELTA_SPACE(1));
+ CCTK_REAL /*const*/ dz CCTK_ATTRIBUTE_UNUSED = ToReal(CCTK_DELTA_SPACE(2));
+ CCTK_REAL /*const*/ dt CCTK_ATTRIBUTE_UNUSED = ToReal(CCTK_DELTA_TIME);
+ CCTK_REAL /*const*/ t CCTK_ATTRIBUTE_UNUSED = ToReal(cctk_time);
+ CCTK_REAL /*const*/ dxi CCTK_ATTRIBUTE_UNUSED = INV(dx);
+ CCTK_REAL /*const*/ dyi CCTK_ATTRIBUTE_UNUSED = INV(dy);
+ CCTK_REAL /*const*/ dzi CCTK_ATTRIBUTE_UNUSED = INV(dz);
+ CCTK_REAL /*const*/ khalf CCTK_ATTRIBUTE_UNUSED = 0.5;
+ CCTK_REAL /*const*/ kthird CCTK_ATTRIBUTE_UNUSED = 1/3.0;
+ CCTK_REAL /*const*/ ktwothird CCTK_ATTRIBUTE_UNUSED = 2.0/3.0;
+ CCTK_REAL /*const*/ kfourthird CCTK_ATTRIBUTE_UNUSED = 4.0/3.0;
+ CCTK_REAL /*const*/ keightthird CCTK_ATTRIBUTE_UNUSED = 8.0/3.0;
+ CCTK_REAL /*const*/ hdxi CCTK_ATTRIBUTE_UNUSED = 0.5 * dxi;
+ CCTK_REAL /*const*/ hdyi CCTK_ATTRIBUTE_UNUSED = 0.5 * dyi;
+ CCTK_REAL /*const*/ hdzi CCTK_ATTRIBUTE_UNUSED = 0.5 * dzi;
/* Initialize predefined quantities */
- CCTK_REAL const p1o1 = 1;
- CCTK_REAL const p1odx = INV(dx);
- CCTK_REAL const p1ody = INV(dy);
- CCTK_REAL const p1odz = INV(dz);
+ CCTK_REAL /*const*/ p1o1 CCTK_ATTRIBUTE_UNUSED = 1.;
+ CCTK_REAL /*const*/ p1odx CCTK_ATTRIBUTE_UNUSED = INV(dx);
+ CCTK_REAL /*const*/ p1ody CCTK_ATTRIBUTE_UNUSED = INV(dy);
+ CCTK_REAL /*const*/ p1odz CCTK_ATTRIBUTE_UNUSED = INV(dz);
/* Assign local copies of arrays functions */
@@ -86,21 +86,21 @@ static void burgers_rhs_1_Body(cGH const * restrict const cctkGH, int const dir,
i,j,k, imin[0],imin[1],imin[2], imax[0],imax[1],imax[2],
cctk_ash[0],cctk_ash[1],cctk_ash[2])
{
- ptrdiff_t const index = di*i + dj*j + dk*k;
+ ptrdiff_t /*const*/ index CCTK_ATTRIBUTE_UNUSED = di*i + dj*j + dk*k;
/* Assign local copies of grid functions */
- CCTK_REAL uFL = uF[index];
- CCTK_REAL urhsL = urhs[index];
+ CCTK_REAL uFL CCTK_ATTRIBUTE_UNUSED = uF[index];
+ CCTK_REAL urhsL CCTK_ATTRIBUTE_UNUSED = urhs[index];
/* Include user supplied include files */
/* Precompute derivatives */
- CCTK_REAL const PDplus1uF = PDplus1(&uF[index]);
+ CCTK_REAL /*const*/ PDplus1uF CCTK_ATTRIBUTE_UNUSED = PDplus1(&uF[index]);
/* Calculate temporaries and grid functions */
- urhsL = urhsL - PDplus1uF;
+ urhsL = urhsL - 1.*PDplus1uF;
/* Copy local copies back to grid functions */
urhs[index] = urhsL;
diff --git a/Examples/Burgers/src/burgers_zero_rhs.cc b/Examples/Burgers/src/burgers_zero_rhs.cc
index 689339b..0d5b35e 100644
--- a/Examples/Burgers/src/burgers_zero_rhs.cc
+++ b/Examples/Burgers/src/burgers_zero_rhs.cc
@@ -31,34 +31,34 @@ static void burgers_zero_rhs_Body(cGH const * restrict const cctkGH, int const d
/* Include user-supplied include files */
/* Initialise finite differencing variables */
- ptrdiff_t const di = 1;
- ptrdiff_t const dj = CCTK_GFINDEX3D(cctkGH,0,1,0) - CCTK_GFINDEX3D(cctkGH,0,0,0);
- ptrdiff_t const dk = CCTK_GFINDEX3D(cctkGH,0,0,1) - CCTK_GFINDEX3D(cctkGH,0,0,0);
- ptrdiff_t const cdi = sizeof(CCTK_REAL) * di;
- ptrdiff_t const cdj = sizeof(CCTK_REAL) * dj;
- ptrdiff_t const cdk = sizeof(CCTK_REAL) * dk;
- CCTK_REAL const dx = ToReal(CCTK_DELTA_SPACE(0));
- CCTK_REAL const dy = ToReal(CCTK_DELTA_SPACE(1));
- CCTK_REAL const dz = ToReal(CCTK_DELTA_SPACE(2));
- CCTK_REAL const dt = ToReal(CCTK_DELTA_TIME);
- CCTK_REAL const t = ToReal(cctk_time);
- CCTK_REAL const dxi = INV(dx);
- CCTK_REAL const dyi = INV(dy);
- CCTK_REAL const dzi = INV(dz);
- CCTK_REAL const khalf = 0.5;
- CCTK_REAL const kthird = 1/3.0;
- CCTK_REAL const ktwothird = 2.0/3.0;
- CCTK_REAL const kfourthird = 4.0/3.0;
- CCTK_REAL const keightthird = 8.0/3.0;
- CCTK_REAL const hdxi = 0.5 * dxi;
- CCTK_REAL const hdyi = 0.5 * dyi;
- CCTK_REAL const hdzi = 0.5 * dzi;
+ ptrdiff_t /*const*/ di CCTK_ATTRIBUTE_UNUSED = 1;
+ ptrdiff_t /*const*/ dj CCTK_ATTRIBUTE_UNUSED = CCTK_GFINDEX3D(cctkGH,0,1,0) - CCTK_GFINDEX3D(cctkGH,0,0,0);
+ ptrdiff_t /*const*/ dk CCTK_ATTRIBUTE_UNUSED = CCTK_GFINDEX3D(cctkGH,0,0,1) - CCTK_GFINDEX3D(cctkGH,0,0,0);
+ ptrdiff_t /*const*/ cdi CCTK_ATTRIBUTE_UNUSED = sizeof(CCTK_REAL) * di;
+ ptrdiff_t /*const*/ cdj CCTK_ATTRIBUTE_UNUSED = sizeof(CCTK_REAL) * dj;
+ ptrdiff_t /*const*/ cdk CCTK_ATTRIBUTE_UNUSED = sizeof(CCTK_REAL) * dk;
+ CCTK_REAL /*const*/ dx CCTK_ATTRIBUTE_UNUSED = ToReal(CCTK_DELTA_SPACE(0));
+ CCTK_REAL /*const*/ dy CCTK_ATTRIBUTE_UNUSED = ToReal(CCTK_DELTA_SPACE(1));
+ CCTK_REAL /*const*/ dz CCTK_ATTRIBUTE_UNUSED = ToReal(CCTK_DELTA_SPACE(2));
+ CCTK_REAL /*const*/ dt CCTK_ATTRIBUTE_UNUSED = ToReal(CCTK_DELTA_TIME);
+ CCTK_REAL /*const*/ t CCTK_ATTRIBUTE_UNUSED = ToReal(cctk_time);
+ CCTK_REAL /*const*/ dxi CCTK_ATTRIBUTE_UNUSED = INV(dx);
+ CCTK_REAL /*const*/ dyi CCTK_ATTRIBUTE_UNUSED = INV(dy);
+ CCTK_REAL /*const*/ dzi CCTK_ATTRIBUTE_UNUSED = INV(dz);
+ CCTK_REAL /*const*/ khalf CCTK_ATTRIBUTE_UNUSED = 0.5;
+ CCTK_REAL /*const*/ kthird CCTK_ATTRIBUTE_UNUSED = 1/3.0;
+ CCTK_REAL /*const*/ ktwothird CCTK_ATTRIBUTE_UNUSED = 2.0/3.0;
+ CCTK_REAL /*const*/ kfourthird CCTK_ATTRIBUTE_UNUSED = 4.0/3.0;
+ CCTK_REAL /*const*/ keightthird CCTK_ATTRIBUTE_UNUSED = 8.0/3.0;
+ CCTK_REAL /*const*/ hdxi CCTK_ATTRIBUTE_UNUSED = 0.5 * dxi;
+ CCTK_REAL /*const*/ hdyi CCTK_ATTRIBUTE_UNUSED = 0.5 * dyi;
+ CCTK_REAL /*const*/ hdzi CCTK_ATTRIBUTE_UNUSED = 0.5 * dzi;
/* Initialize predefined quantities */
- CCTK_REAL const p1o1 = 1;
- CCTK_REAL const p1odx = INV(dx);
- CCTK_REAL const p1ody = INV(dy);
- CCTK_REAL const p1odz = INV(dz);
+ CCTK_REAL /*const*/ p1o1 CCTK_ATTRIBUTE_UNUSED = 1.;
+ CCTK_REAL /*const*/ p1odx CCTK_ATTRIBUTE_UNUSED = INV(dx);
+ CCTK_REAL /*const*/ p1ody CCTK_ATTRIBUTE_UNUSED = INV(dy);
+ CCTK_REAL /*const*/ p1odz CCTK_ATTRIBUTE_UNUSED = INV(dz);
/* Assign local copies of arrays functions */
@@ -74,7 +74,7 @@ static void burgers_zero_rhs_Body(cGH const * restrict const cctkGH, int const d
i,j,k, imin[0],imin[1],imin[2], imax[0],imax[1],imax[2],
cctk_ash[0],cctk_ash[1],cctk_ash[2])
{
- ptrdiff_t const index = di*i + dj*j + dk*k;
+ ptrdiff_t /*const*/ index CCTK_ATTRIBUTE_UNUSED = di*i + dj*j + dk*k;
/* Assign local copies of grid functions */
@@ -85,7 +85,7 @@ static void burgers_zero_rhs_Body(cGH const * restrict const cctkGH, int const d
/* Precompute derivatives */
/* Calculate temporaries and grid functions */
- CCTK_REAL urhsL = 0;
+ CCTK_REAL CCTK_ATTRIBUTE_UNUSED urhsL = 0.;
/* Copy local copies back to grid functions */
urhs[index] = urhsL;