diff options
Diffstat (limited to 'Examples/Burgers')
-rw-r--r-- | Examples/Burgers/interface.ccl | 2 | ||||
-rw-r--r-- | Examples/Burgers/schedule.ccl | 37 | ||||
-rw-r--r-- | Examples/Burgers/src/Boundaries.cc | 10 | ||||
-rw-r--r-- | Examples/Burgers/src/Differencing.h | 60 | ||||
-rw-r--r-- | Examples/Burgers/src/RegisterMoL.cc | 2 | ||||
-rw-r--r-- | Examples/Burgers/src/Startup.cc | 2 | ||||
-rw-r--r-- | Examples/Burgers/src/burgers_flux_1.cc | 67 | ||||
-rw-r--r-- | Examples/Burgers/src/burgers_initial_shock.cc | 60 | ||||
-rw-r--r-- | Examples/Burgers/src/burgers_initial_sine.cc | 59 | ||||
-rw-r--r-- | Examples/Burgers/src/burgers_reconstruct_1.cc | 74 | ||||
-rw-r--r-- | Examples/Burgers/src/burgers_rhs_1.cc | 64 | ||||
-rw-r--r-- | Examples/Burgers/src/burgers_zero_rhs.cc | 56 |
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; |