diff options
Diffstat (limited to 'Examples/Wave')
-rw-r--r-- | Examples/Wave/interface.ccl | 2 | ||||
-rw-r--r-- | Examples/Wave/schedule.ccl | 91 | ||||
-rw-r--r-- | Examples/Wave/src/Boundaries.cc | 14 | ||||
-rw-r--r-- | Examples/Wave/src/Differencing.h | 264 | ||||
-rw-r--r-- | Examples/Wave/src/RegisterMoL.cc | 2 | ||||
-rw-r--r-- | Examples/Wave/src/Startup.cc | 2 | ||||
-rw-r--r-- | Examples/Wave/src/wave_boundary.cc | 139 | ||||
-rw-r--r-- | Examples/Wave/src/wave_calc_errors.cc | 114 | ||||
-rw-r--r-- | Examples/Wave/src/wave_calc_norm.cc | 133 | ||||
-rw-r--r-- | Examples/Wave/src/wave_evolve.cc | 120 | ||||
-rw-r--r-- | Examples/Wave/src/wave_exact_gaussian.cc | 133 | ||||
-rw-r--r-- | Examples/Wave/src/wave_exact_sine.cc | 122 | ||||
-rw-r--r-- | Examples/Wave/src/wave_import_exact.cc | 110 |
13 files changed, 624 insertions, 622 deletions
diff --git a/Examples/Wave/interface.ccl b/Examples/Wave/interface.ccl index 89fef10..a07844c 100644 --- a/Examples/Wave/interface.ccl +++ b/Examples/Wave/interface.ccl @@ -57,7 +57,7 @@ CCTK_REAL evolved type=GF timelevels=3 tags='' } "evolved" public: -CCTK_REAL evolvedrhs type=GF timelevels=3 tags='' +CCTK_REAL evolvedrhs type=GF timelevels=3 tags='Prolongation="None"' { phirhs, pirhs diff --git a/Examples/Wave/schedule.ccl b/Examples/Wave/schedule.ccl index 719d86f..6d2a02c 100644 --- a/Examples/Wave/schedule.ccl +++ b/Examples/Wave/schedule.ccl @@ -84,11 +84,11 @@ if (CCTK_EQUALS(initial_data, "sine")) schedule wave_exact_sine AT INITIAL before import_exact before import_exact { LANG: C - READS: grid::x - READS: grid::y - READS: grid::z - WRITES: Wave::phiExact - WRITES: Wave::piExact + READS: grid::x(Everywhere) + READS: grid::y(Everywhere) + READS: grid::z(Everywhere) + WRITES: Wave::phiExact(Everywhere) + WRITES: Wave::piExact(Everywhere) } "wave_exact_sine" } @@ -98,11 +98,11 @@ if (CCTK_EQUALS(initial_data, "sine")) schedule wave_exact_sine AT POSTSTEP before calc_errors before import_exact { LANG: C - READS: grid::x - READS: grid::y - READS: grid::z - WRITES: Wave::phiExact - WRITES: Wave::piExact + READS: grid::x(Everywhere) + READS: grid::y(Everywhere) + READS: grid::z(Everywhere) + WRITES: Wave::phiExact(Everywhere) + WRITES: Wave::piExact(Everywhere) } "wave_exact_sine" } @@ -112,9 +112,9 @@ if (CCTK_EQUALS(initial_data, "gaussian")) schedule wave_exact_gaussian AT INITIAL before import_exact before import_exact { LANG: C - READS: grid::r - WRITES: Wave::phiExact - WRITES: Wave::piExact + READS: grid::r(Everywhere) + WRITES: Wave::phiExact(Everywhere) + WRITES: Wave::piExact(Everywhere) } "wave_exact_gaussian" } @@ -124,52 +124,52 @@ if (CCTK_EQUALS(initial_data, "gaussian")) schedule wave_exact_gaussian AT POSTSTEP before calc_errors before import_exact { LANG: C - READS: grid::r - WRITES: Wave::phiExact - WRITES: Wave::piExact + READS: grid::r(Everywhere) + WRITES: Wave::phiExact(Everywhere) + WRITES: Wave::piExact(Everywhere) } "wave_exact_gaussian" } schedule wave_import_exact at INITIAL as import_exact { LANG: C - READS: Wave::phiExact - READS: Wave::piExact - WRITES: Wave::phi - WRITES: Wave::pi + READS: Wave::phiExact(Everywhere) + READS: Wave::piExact(Everywhere) + WRITES: Wave::phi(Everywhere) + WRITES: Wave::pi(Everywhere) } "wave_import_exact" schedule wave_evolve in MoL_CalcRHS as evolve { LANG: C - READS: Wave::phi - READS: Wave::pi - WRITES: Wave::phirhs - WRITES: Wave::pirhs + READS: Wave::phi(Everywhere) + READS: Wave::pi(Everywhere) + WRITES: Wave::phirhs(Interior) + WRITES: Wave::pirhs(Interior) } "wave_evolve" schedule wave_calc_errors at ANALYSIS as calc_errors { LANG: C - READS: Wave::phi - READS: Wave::phiExact - READS: Wave::pi - READS: Wave::piExact - WRITES: Wave::phiError - WRITES: Wave::piError + READS: Wave::phi(Everywhere) + READS: Wave::pi(Everywhere) + READS: Wave::phiExact(Everywhere) + READS: Wave::piExact(Everywhere) + WRITES: Wave::phiError(Everywhere) + WRITES: Wave::piError(Everywhere) } "wave_calc_errors" schedule wave_calc_norm at ANALYSIS as calc_norm { LANG: C SYNC: norms - READS: Wave::phi - READS: Wave::phiError - READS: Wave::pi - READS: Wave::piError - WRITES: Wave::EL2 - WRITES: Wave::VDP - WRITES: Wave::VL2 + READS: Wave::phiError(Everywhere) + READS: Wave::piError(Everywhere) + READS: Wave::phi(Everywhere) + READS: Wave::pi(Everywhere) + WRITES: Wave::EL2(Interior) + WRITES: Wave::VDP(Interior) + WRITES: Wave::VL2(Interior) } "wave_calc_norm" @@ -178,14 +178,14 @@ if (CCTK_EQUALS(boundary_condition, "radiative")) schedule wave_boundary in MoL_RHSBoundaries { LANG: C - READS: Wave::phi - READS: Wave::pi - READS: grid::r - READS: grid::x - READS: grid::y - READS: grid::z - WRITES: Wave::phirhs - WRITES: Wave::pirhs + READS: grid::x(Everywhere) + READS: grid::y(Everywhere) + READS: grid::z(Everywhere) + READS: grid::r(Everywhere) + READS: Wave::phi(Everywhere) + READS: Wave::pi(Everywhere) + WRITES: Wave::phirhs(Boundary) + WRITES: Wave::pirhs(Boundary) } "wave_boundary" } @@ -210,5 +210,4 @@ schedule Wave_RegisterVars in MoL_Register schedule group ApplyBCs as Wave_ApplyBCs in MoL_PostStep after Wave_SelectBoundConds { - # no language specified } "Apply boundary conditions controlled by thorn Boundary" diff --git a/Examples/Wave/src/Boundaries.cc b/Examples/Wave/src/Boundaries.cc index 672d7f3..2ea1b08 100644 --- a/Examples/Wave/src/Boundaries.cc +++ b/Examples/Wave/src/Boundaries.cc @@ -30,7 +30,7 @@ extern "C" void Wave_SelectBoundConds(CCTK_ARGUMENTS) DECLARE_CCTK_ARGUMENTS; DECLARE_CCTK_PARAMETERS; - CCTK_INT ierr = 0; + CCTK_INT ierr CCTK_ATTRIBUTE_UNUSED = 0; if (CCTK_EQUALS(evolved_bound, "none" ) || CCTK_EQUALS(evolved_bound, "static") || @@ -68,7 +68,7 @@ extern "C" void Wave_SelectBoundConds(CCTK_ARGUMENTS) if (CCTK_EQUALS(evolved_bound, "radiative")) { /* select radiation boundary condition */ - static CCTK_INT handle_evolved_bound = -1; + static CCTK_INT handle_evolved_bound CCTK_ATTRIBUTE_UNUSED = -1; if (handle_evolved_bound < 0) handle_evolved_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE); if (handle_evolved_bound < 0) CCTK_WARN(0, "could not create table!"); if (Util_TableSetReal(handle_evolved_bound , evolved_bound_limit, "LIMIT") < 0) @@ -87,7 +87,7 @@ extern "C" void Wave_SelectBoundConds(CCTK_ARGUMENTS) if (CCTK_EQUALS(phi_bound, "radiative")) { /* select radiation boundary condition */ - static CCTK_INT handle_phi_bound = -1; + static CCTK_INT handle_phi_bound CCTK_ATTRIBUTE_UNUSED = -1; if (handle_phi_bound < 0) handle_phi_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE); if (handle_phi_bound < 0) CCTK_WARN(0, "could not create table!"); if (Util_TableSetReal(handle_phi_bound , phi_bound_limit, "LIMIT") < 0) @@ -106,7 +106,7 @@ extern "C" void Wave_SelectBoundConds(CCTK_ARGUMENTS) if (CCTK_EQUALS(pi_bound, "radiative")) { /* select radiation boundary condition */ - static CCTK_INT handle_pi_bound = -1; + static CCTK_INT handle_pi_bound CCTK_ATTRIBUTE_UNUSED = -1; if (handle_pi_bound < 0) handle_pi_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE); if (handle_pi_bound < 0) CCTK_WARN(0, "could not create table!"); if (Util_TableSetReal(handle_pi_bound , pi_bound_limit, "LIMIT") < 0) @@ -125,7 +125,7 @@ extern "C" void Wave_SelectBoundConds(CCTK_ARGUMENTS) if (CCTK_EQUALS(evolved_bound, "scalar")) { /* select scalar boundary condition */ - static CCTK_INT handle_evolved_bound = -1; + static CCTK_INT handle_evolved_bound CCTK_ATTRIBUTE_UNUSED = -1; if (handle_evolved_bound < 0) handle_evolved_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE); if (handle_evolved_bound < 0) CCTK_WARN(0, "could not create table!"); if (Util_TableSetReal(handle_evolved_bound ,evolved_bound_scalar, "SCALAR") < 0) @@ -142,7 +142,7 @@ extern "C" void Wave_SelectBoundConds(CCTK_ARGUMENTS) if (CCTK_EQUALS(phi_bound, "scalar")) { /* select scalar boundary condition */ - static CCTK_INT handle_phi_bound = -1; + static CCTK_INT handle_phi_bound CCTK_ATTRIBUTE_UNUSED = -1; if (handle_phi_bound < 0) handle_phi_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE); if (handle_phi_bound < 0) CCTK_WARN(0, "could not create table!"); if (Util_TableSetReal(handle_phi_bound ,phi_bound_scalar, "SCALAR") < 0) @@ -159,7 +159,7 @@ extern "C" void Wave_SelectBoundConds(CCTK_ARGUMENTS) if (CCTK_EQUALS(pi_bound, "scalar")) { /* select scalar boundary condition */ - static CCTK_INT handle_pi_bound = -1; + static CCTK_INT handle_pi_bound CCTK_ATTRIBUTE_UNUSED = -1; if (handle_pi_bound < 0) handle_pi_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE); if (handle_pi_bound < 0) CCTK_WARN(0, "could not create table!"); if (Util_TableSetReal(handle_pi_bound ,pi_bound_scalar, "SCALAR") < 0) diff --git a/Examples/Wave/src/Differencing.h b/Examples/Wave/src/Differencing.h index 28eb5d5..4032807 100644 --- a/Examples/Wave/src/Differencing.h +++ b/Examples/Wave/src/Differencing.h @@ -1,528 +1,528 @@ #ifndef KRANC_DIFF_FUNCTIONS -# define PDstandardfdOrder21(u) ((-KRANC_GFOFFSET3D(u,-1,0,0) + KRANC_GFOFFSET3D(u,1,0,0))*p1o2dx) +# define PDstandardfdOrder21(u) ((-1.*KRANC_GFOFFSET3D(u,-1,0,0) + KRANC_GFOFFSET3D(u,1,0,0))*p1o2dx) #else # define PDstandardfdOrder21(u) (PDstandardfdOrder21_impl(u,p1o2dx,cdj,cdk)) static CCTK_REAL PDstandardfdOrder21_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1o2dx, ptrdiff_t const cdj, ptrdiff_t const cdk) CCTK_ATTRIBUTE_NOINLINE CCTK_ATTRIBUTE_UNUSED; static CCTK_REAL PDstandardfdOrder21_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1o2dx, ptrdiff_t const cdj, ptrdiff_t const cdk) { - ptrdiff_t const cdi=sizeof(CCTK_REAL); - return (-KRANC_GFOFFSET3D(u,-1,0,0) + KRANC_GFOFFSET3D(u,1,0,0))*p1o2dx; + ptrdiff_t const cdi CCTK_ATTRIBUTE_UNUSED = sizeof(CCTK_REAL); + return (-1.*KRANC_GFOFFSET3D(u,-1,0,0) + KRANC_GFOFFSET3D(u,1,0,0))*p1o2dx; } #endif #ifndef KRANC_DIFF_FUNCTIONS -# define PDstandardfdOrder22(u) ((-KRANC_GFOFFSET3D(u,0,-1,0) + KRANC_GFOFFSET3D(u,0,1,0))*p1o2dy) +# define PDstandardfdOrder22(u) ((-1.*KRANC_GFOFFSET3D(u,0,-1,0) + KRANC_GFOFFSET3D(u,0,1,0))*p1o2dy) #else # define PDstandardfdOrder22(u) (PDstandardfdOrder22_impl(u,p1o2dy,cdj,cdk)) static CCTK_REAL PDstandardfdOrder22_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1o2dy, ptrdiff_t const cdj, ptrdiff_t const cdk) CCTK_ATTRIBUTE_NOINLINE CCTK_ATTRIBUTE_UNUSED; static CCTK_REAL PDstandardfdOrder22_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1o2dy, ptrdiff_t const cdj, ptrdiff_t const cdk) { - ptrdiff_t const cdi=sizeof(CCTK_REAL); - return (-KRANC_GFOFFSET3D(u,0,-1,0) + KRANC_GFOFFSET3D(u,0,1,0))*p1o2dy; + ptrdiff_t const cdi CCTK_ATTRIBUTE_UNUSED = sizeof(CCTK_REAL); + return (-1.*KRANC_GFOFFSET3D(u,0,-1,0) + KRANC_GFOFFSET3D(u,0,1,0))*p1o2dy; } #endif #ifndef KRANC_DIFF_FUNCTIONS -# define PDstandardfdOrder23(u) ((-KRANC_GFOFFSET3D(u,0,0,-1) + KRANC_GFOFFSET3D(u,0,0,1))*p1o2dz) +# define PDstandardfdOrder23(u) ((-1.*KRANC_GFOFFSET3D(u,0,0,-1) + KRANC_GFOFFSET3D(u,0,0,1))*p1o2dz) #else # define PDstandardfdOrder23(u) (PDstandardfdOrder23_impl(u,p1o2dz,cdj,cdk)) static CCTK_REAL PDstandardfdOrder23_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1o2dz, ptrdiff_t const cdj, ptrdiff_t const cdk) CCTK_ATTRIBUTE_NOINLINE CCTK_ATTRIBUTE_UNUSED; static CCTK_REAL PDstandardfdOrder23_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1o2dz, ptrdiff_t const cdj, ptrdiff_t const cdk) { - ptrdiff_t const cdi=sizeof(CCTK_REAL); - return (-KRANC_GFOFFSET3D(u,0,0,-1) + KRANC_GFOFFSET3D(u,0,0,1))*p1o2dz; + ptrdiff_t const cdi CCTK_ATTRIBUTE_UNUSED = sizeof(CCTK_REAL); + return (-1.*KRANC_GFOFFSET3D(u,0,0,-1) + KRANC_GFOFFSET3D(u,0,0,1))*p1o2dz; } #endif #ifndef KRANC_DIFF_FUNCTIONS -# define PDstandardfdOrder41(u) ((-8*KRANC_GFOFFSET3D(u,-1,0,0) + 8*KRANC_GFOFFSET3D(u,1,0,0) + KRANC_GFOFFSET3D(u,-2,0,0) - KRANC_GFOFFSET3D(u,2,0,0))*p1o12dx) +# define PDstandardfdOrder41(u) ((-8.*KRANC_GFOFFSET3D(u,-1,0,0) + 8.*KRANC_GFOFFSET3D(u,1,0,0) + KRANC_GFOFFSET3D(u,-2,0,0) - 1.*KRANC_GFOFFSET3D(u,2,0,0))*p1o12dx) #else # define PDstandardfdOrder41(u) (PDstandardfdOrder41_impl(u,p1o12dx,cdj,cdk)) static CCTK_REAL PDstandardfdOrder41_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1o12dx, ptrdiff_t const cdj, ptrdiff_t const cdk) CCTK_ATTRIBUTE_NOINLINE CCTK_ATTRIBUTE_UNUSED; static CCTK_REAL PDstandardfdOrder41_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1o12dx, ptrdiff_t const cdj, ptrdiff_t const cdk) { - ptrdiff_t const cdi=sizeof(CCTK_REAL); - return (-8*KRANC_GFOFFSET3D(u,-1,0,0) + 8*KRANC_GFOFFSET3D(u,1,0,0) + KRANC_GFOFFSET3D(u,-2,0,0) - KRANC_GFOFFSET3D(u,2,0,0))*p1o12dx; + ptrdiff_t const cdi CCTK_ATTRIBUTE_UNUSED = sizeof(CCTK_REAL); + return (-8.*KRANC_GFOFFSET3D(u,-1,0,0) + 8.*KRANC_GFOFFSET3D(u,1,0,0) + KRANC_GFOFFSET3D(u,-2,0,0) - 1.*KRANC_GFOFFSET3D(u,2,0,0))*p1o12dx; } #endif #ifndef KRANC_DIFF_FUNCTIONS -# define PDstandardfdOrder42(u) ((-8*KRANC_GFOFFSET3D(u,0,-1,0) + 8*KRANC_GFOFFSET3D(u,0,1,0) + KRANC_GFOFFSET3D(u,0,-2,0) - KRANC_GFOFFSET3D(u,0,2,0))*p1o12dy) +# define PDstandardfdOrder42(u) ((-8.*KRANC_GFOFFSET3D(u,0,-1,0) + 8.*KRANC_GFOFFSET3D(u,0,1,0) + KRANC_GFOFFSET3D(u,0,-2,0) - 1.*KRANC_GFOFFSET3D(u,0,2,0))*p1o12dy) #else # define PDstandardfdOrder42(u) (PDstandardfdOrder42_impl(u,p1o12dy,cdj,cdk)) static CCTK_REAL PDstandardfdOrder42_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1o12dy, ptrdiff_t const cdj, ptrdiff_t const cdk) CCTK_ATTRIBUTE_NOINLINE CCTK_ATTRIBUTE_UNUSED; static CCTK_REAL PDstandardfdOrder42_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1o12dy, ptrdiff_t const cdj, ptrdiff_t const cdk) { - ptrdiff_t const cdi=sizeof(CCTK_REAL); - return (-8*KRANC_GFOFFSET3D(u,0,-1,0) + 8*KRANC_GFOFFSET3D(u,0,1,0) + KRANC_GFOFFSET3D(u,0,-2,0) - KRANC_GFOFFSET3D(u,0,2,0))*p1o12dy; + ptrdiff_t const cdi CCTK_ATTRIBUTE_UNUSED = sizeof(CCTK_REAL); + return (-8.*KRANC_GFOFFSET3D(u,0,-1,0) + 8.*KRANC_GFOFFSET3D(u,0,1,0) + KRANC_GFOFFSET3D(u,0,-2,0) - 1.*KRANC_GFOFFSET3D(u,0,2,0))*p1o12dy; } #endif #ifndef KRANC_DIFF_FUNCTIONS -# define PDstandardfdOrder43(u) ((-8*KRANC_GFOFFSET3D(u,0,0,-1) + 8*KRANC_GFOFFSET3D(u,0,0,1) + KRANC_GFOFFSET3D(u,0,0,-2) - KRANC_GFOFFSET3D(u,0,0,2))*p1o12dz) +# define PDstandardfdOrder43(u) ((-8.*KRANC_GFOFFSET3D(u,0,0,-1) + 8.*KRANC_GFOFFSET3D(u,0,0,1) + KRANC_GFOFFSET3D(u,0,0,-2) - 1.*KRANC_GFOFFSET3D(u,0,0,2))*p1o12dz) #else # define PDstandardfdOrder43(u) (PDstandardfdOrder43_impl(u,p1o12dz,cdj,cdk)) static CCTK_REAL PDstandardfdOrder43_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1o12dz, ptrdiff_t const cdj, ptrdiff_t const cdk) CCTK_ATTRIBUTE_NOINLINE CCTK_ATTRIBUTE_UNUSED; static CCTK_REAL PDstandardfdOrder43_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1o12dz, ptrdiff_t const cdj, ptrdiff_t const cdk) { - ptrdiff_t const cdi=sizeof(CCTK_REAL); - return (-8*KRANC_GFOFFSET3D(u,0,0,-1) + 8*KRANC_GFOFFSET3D(u,0,0,1) + KRANC_GFOFFSET3D(u,0,0,-2) - KRANC_GFOFFSET3D(u,0,0,2))*p1o12dz; + ptrdiff_t const cdi CCTK_ATTRIBUTE_UNUSED = sizeof(CCTK_REAL); + return (-8.*KRANC_GFOFFSET3D(u,0,0,-1) + 8.*KRANC_GFOFFSET3D(u,0,0,1) + KRANC_GFOFFSET3D(u,0,0,-2) - 1.*KRANC_GFOFFSET3D(u,0,0,2))*p1o12dz; } #endif #ifndef KRANC_DIFF_FUNCTIONS -# define PDstandardfdOrder211(u) ((-2*KRANC_GFOFFSET3D(u,0,0,0) + KRANC_GFOFFSET3D(u,-1,0,0) + KRANC_GFOFFSET3D(u,1,0,0))*p1odx2) +# define PDstandardfdOrder211(u) ((-2.*KRANC_GFOFFSET3D(u,0,0,0) + KRANC_GFOFFSET3D(u,-1,0,0) + KRANC_GFOFFSET3D(u,1,0,0))*p1odx2) #else # define PDstandardfdOrder211(u) (PDstandardfdOrder211_impl(u,p1odx2,cdj,cdk)) static CCTK_REAL PDstandardfdOrder211_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1odx2, ptrdiff_t const cdj, ptrdiff_t const cdk) CCTK_ATTRIBUTE_NOINLINE CCTK_ATTRIBUTE_UNUSED; static CCTK_REAL PDstandardfdOrder211_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1odx2, ptrdiff_t const cdj, ptrdiff_t const cdk) { - ptrdiff_t const cdi=sizeof(CCTK_REAL); - return (-2*KRANC_GFOFFSET3D(u,0,0,0) + KRANC_GFOFFSET3D(u,-1,0,0) + KRANC_GFOFFSET3D(u,1,0,0))*p1odx2; + ptrdiff_t const cdi CCTK_ATTRIBUTE_UNUSED = sizeof(CCTK_REAL); + return (-2.*KRANC_GFOFFSET3D(u,0,0,0) + KRANC_GFOFFSET3D(u,-1,0,0) + KRANC_GFOFFSET3D(u,1,0,0))*p1odx2; } #endif #ifndef KRANC_DIFF_FUNCTIONS -# define PDstandardfdOrder222(u) ((-2*KRANC_GFOFFSET3D(u,0,0,0) + KRANC_GFOFFSET3D(u,0,-1,0) + KRANC_GFOFFSET3D(u,0,1,0))*p1ody2) +# define PDstandardfdOrder222(u) ((-2.*KRANC_GFOFFSET3D(u,0,0,0) + KRANC_GFOFFSET3D(u,0,-1,0) + KRANC_GFOFFSET3D(u,0,1,0))*p1ody2) #else # define PDstandardfdOrder222(u) (PDstandardfdOrder222_impl(u,p1ody2,cdj,cdk)) static CCTK_REAL PDstandardfdOrder222_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1ody2, ptrdiff_t const cdj, ptrdiff_t const cdk) CCTK_ATTRIBUTE_NOINLINE CCTK_ATTRIBUTE_UNUSED; static CCTK_REAL PDstandardfdOrder222_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1ody2, ptrdiff_t const cdj, ptrdiff_t const cdk) { - ptrdiff_t const cdi=sizeof(CCTK_REAL); - return (-2*KRANC_GFOFFSET3D(u,0,0,0) + KRANC_GFOFFSET3D(u,0,-1,0) + KRANC_GFOFFSET3D(u,0,1,0))*p1ody2; + ptrdiff_t const cdi CCTK_ATTRIBUTE_UNUSED = sizeof(CCTK_REAL); + return (-2.*KRANC_GFOFFSET3D(u,0,0,0) + KRANC_GFOFFSET3D(u,0,-1,0) + KRANC_GFOFFSET3D(u,0,1,0))*p1ody2; } #endif #ifndef KRANC_DIFF_FUNCTIONS -# define PDstandardfdOrder233(u) ((-2*KRANC_GFOFFSET3D(u,0,0,0) + KRANC_GFOFFSET3D(u,0,0,-1) + KRANC_GFOFFSET3D(u,0,0,1))*p1odz2) +# define PDstandardfdOrder233(u) ((-2.*KRANC_GFOFFSET3D(u,0,0,0) + KRANC_GFOFFSET3D(u,0,0,-1) + KRANC_GFOFFSET3D(u,0,0,1))*p1odz2) #else # define PDstandardfdOrder233(u) (PDstandardfdOrder233_impl(u,p1odz2,cdj,cdk)) static CCTK_REAL PDstandardfdOrder233_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1odz2, ptrdiff_t const cdj, ptrdiff_t const cdk) CCTK_ATTRIBUTE_NOINLINE CCTK_ATTRIBUTE_UNUSED; static CCTK_REAL PDstandardfdOrder233_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1odz2, ptrdiff_t const cdj, ptrdiff_t const cdk) { - ptrdiff_t const cdi=sizeof(CCTK_REAL); - return (-2*KRANC_GFOFFSET3D(u,0,0,0) + KRANC_GFOFFSET3D(u,0,0,-1) + KRANC_GFOFFSET3D(u,0,0,1))*p1odz2; + ptrdiff_t const cdi CCTK_ATTRIBUTE_UNUSED = sizeof(CCTK_REAL); + return (-2.*KRANC_GFOFFSET3D(u,0,0,0) + KRANC_GFOFFSET3D(u,0,0,-1) + KRANC_GFOFFSET3D(u,0,0,1))*p1odz2; } #endif #ifndef KRANC_DIFF_FUNCTIONS -# define PDstandardfdOrder411(u) ((30*KRANC_GFOFFSET3D(u,0,0,0) - 16*(KRANC_GFOFFSET3D(u,-1,0,0) + KRANC_GFOFFSET3D(u,1,0,0)) + KRANC_GFOFFSET3D(u,-2,0,0) + KRANC_GFOFFSET3D(u,2,0,0))*pm1o12dx2) +# define PDstandardfdOrder411(u) ((30.*KRANC_GFOFFSET3D(u,0,0,0) - 16.*(KRANC_GFOFFSET3D(u,-1,0,0) + KRANC_GFOFFSET3D(u,1,0,0)) + KRANC_GFOFFSET3D(u,-2,0,0) + KRANC_GFOFFSET3D(u,2,0,0))*pm1o12dx2) #else # define PDstandardfdOrder411(u) (PDstandardfdOrder411_impl(u,pm1o12dx2,cdj,cdk)) static CCTK_REAL PDstandardfdOrder411_impl(CCTK_REAL const* restrict const u, CCTK_REAL const pm1o12dx2, ptrdiff_t const cdj, ptrdiff_t const cdk) CCTK_ATTRIBUTE_NOINLINE CCTK_ATTRIBUTE_UNUSED; static CCTK_REAL PDstandardfdOrder411_impl(CCTK_REAL const* restrict const u, CCTK_REAL const pm1o12dx2, ptrdiff_t const cdj, ptrdiff_t const cdk) { - ptrdiff_t const cdi=sizeof(CCTK_REAL); - return (30*KRANC_GFOFFSET3D(u,0,0,0) - 16*(KRANC_GFOFFSET3D(u,-1,0,0) + KRANC_GFOFFSET3D(u,1,0,0)) + KRANC_GFOFFSET3D(u,-2,0,0) + KRANC_GFOFFSET3D(u,2,0,0))*pm1o12dx2; + ptrdiff_t const cdi CCTK_ATTRIBUTE_UNUSED = sizeof(CCTK_REAL); + return (30.*KRANC_GFOFFSET3D(u,0,0,0) - 16.*(KRANC_GFOFFSET3D(u,-1,0,0) + KRANC_GFOFFSET3D(u,1,0,0)) + KRANC_GFOFFSET3D(u,-2,0,0) + KRANC_GFOFFSET3D(u,2,0,0))*pm1o12dx2; } #endif #ifndef KRANC_DIFF_FUNCTIONS -# define PDstandardfdOrder422(u) ((30*KRANC_GFOFFSET3D(u,0,0,0) - 16*(KRANC_GFOFFSET3D(u,0,-1,0) + KRANC_GFOFFSET3D(u,0,1,0)) + KRANC_GFOFFSET3D(u,0,-2,0) + KRANC_GFOFFSET3D(u,0,2,0))*pm1o12dy2) +# define PDstandardfdOrder422(u) ((30.*KRANC_GFOFFSET3D(u,0,0,0) - 16.*(KRANC_GFOFFSET3D(u,0,-1,0) + KRANC_GFOFFSET3D(u,0,1,0)) + KRANC_GFOFFSET3D(u,0,-2,0) + KRANC_GFOFFSET3D(u,0,2,0))*pm1o12dy2) #else # define PDstandardfdOrder422(u) (PDstandardfdOrder422_impl(u,pm1o12dy2,cdj,cdk)) static CCTK_REAL PDstandardfdOrder422_impl(CCTK_REAL const* restrict const u, CCTK_REAL const pm1o12dy2, ptrdiff_t const cdj, ptrdiff_t const cdk) CCTK_ATTRIBUTE_NOINLINE CCTK_ATTRIBUTE_UNUSED; static CCTK_REAL PDstandardfdOrder422_impl(CCTK_REAL const* restrict const u, CCTK_REAL const pm1o12dy2, ptrdiff_t const cdj, ptrdiff_t const cdk) { - ptrdiff_t const cdi=sizeof(CCTK_REAL); - return (30*KRANC_GFOFFSET3D(u,0,0,0) - 16*(KRANC_GFOFFSET3D(u,0,-1,0) + KRANC_GFOFFSET3D(u,0,1,0)) + KRANC_GFOFFSET3D(u,0,-2,0) + KRANC_GFOFFSET3D(u,0,2,0))*pm1o12dy2; + ptrdiff_t const cdi CCTK_ATTRIBUTE_UNUSED = sizeof(CCTK_REAL); + return (30.*KRANC_GFOFFSET3D(u,0,0,0) - 16.*(KRANC_GFOFFSET3D(u,0,-1,0) + KRANC_GFOFFSET3D(u,0,1,0)) + KRANC_GFOFFSET3D(u,0,-2,0) + KRANC_GFOFFSET3D(u,0,2,0))*pm1o12dy2; } #endif #ifndef KRANC_DIFF_FUNCTIONS -# define PDstandardfdOrder433(u) ((30*KRANC_GFOFFSET3D(u,0,0,0) - 16*(KRANC_GFOFFSET3D(u,0,0,-1) + KRANC_GFOFFSET3D(u,0,0,1)) + KRANC_GFOFFSET3D(u,0,0,-2) + KRANC_GFOFFSET3D(u,0,0,2))*pm1o12dz2) +# define PDstandardfdOrder433(u) ((30.*KRANC_GFOFFSET3D(u,0,0,0) - 16.*(KRANC_GFOFFSET3D(u,0,0,-1) + KRANC_GFOFFSET3D(u,0,0,1)) + KRANC_GFOFFSET3D(u,0,0,-2) + KRANC_GFOFFSET3D(u,0,0,2))*pm1o12dz2) #else # define PDstandardfdOrder433(u) (PDstandardfdOrder433_impl(u,pm1o12dz2,cdj,cdk)) static CCTK_REAL PDstandardfdOrder433_impl(CCTK_REAL const* restrict const u, CCTK_REAL const pm1o12dz2, ptrdiff_t const cdj, ptrdiff_t const cdk) CCTK_ATTRIBUTE_NOINLINE CCTK_ATTRIBUTE_UNUSED; static CCTK_REAL PDstandardfdOrder433_impl(CCTK_REAL const* restrict const u, CCTK_REAL const pm1o12dz2, ptrdiff_t const cdj, ptrdiff_t const cdk) { - ptrdiff_t const cdi=sizeof(CCTK_REAL); - return (30*KRANC_GFOFFSET3D(u,0,0,0) - 16*(KRANC_GFOFFSET3D(u,0,0,-1) + KRANC_GFOFFSET3D(u,0,0,1)) + KRANC_GFOFFSET3D(u,0,0,-2) + KRANC_GFOFFSET3D(u,0,0,2))*pm1o12dz2; + ptrdiff_t const cdi CCTK_ATTRIBUTE_UNUSED = sizeof(CCTK_REAL); + return (30.*KRANC_GFOFFSET3D(u,0,0,0) - 16.*(KRANC_GFOFFSET3D(u,0,0,-1) + KRANC_GFOFFSET3D(u,0,0,1)) + KRANC_GFOFFSET3D(u,0,0,-2) + KRANC_GFOFFSET3D(u,0,0,2))*pm1o12dz2; } #endif #ifndef KRANC_DIFF_FUNCTIONS -# define PDstandardfdOrder212(u) ((KRANC_GFOFFSET3D(u,-1,-1,0) - KRANC_GFOFFSET3D(u,-1,1,0) - KRANC_GFOFFSET3D(u,1,-1,0) + KRANC_GFOFFSET3D(u,1,1,0))*p1o4dxdy) +# define PDstandardfdOrder212(u) ((KRANC_GFOFFSET3D(u,-1,-1,0) - 1.*(KRANC_GFOFFSET3D(u,-1,1,0) + KRANC_GFOFFSET3D(u,1,-1,0)) + KRANC_GFOFFSET3D(u,1,1,0))*p1o4dxdy) #else # define PDstandardfdOrder212(u) (PDstandardfdOrder212_impl(u,p1o4dxdy,cdj,cdk)) static CCTK_REAL PDstandardfdOrder212_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1o4dxdy, ptrdiff_t const cdj, ptrdiff_t const cdk) CCTK_ATTRIBUTE_NOINLINE CCTK_ATTRIBUTE_UNUSED; static CCTK_REAL PDstandardfdOrder212_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1o4dxdy, ptrdiff_t const cdj, ptrdiff_t const cdk) { - ptrdiff_t const cdi=sizeof(CCTK_REAL); - return (KRANC_GFOFFSET3D(u,-1,-1,0) - KRANC_GFOFFSET3D(u,-1,1,0) - KRANC_GFOFFSET3D(u,1,-1,0) + KRANC_GFOFFSET3D(u,1,1,0))*p1o4dxdy; + ptrdiff_t const cdi CCTK_ATTRIBUTE_UNUSED = sizeof(CCTK_REAL); + return (KRANC_GFOFFSET3D(u,-1,-1,0) - 1.*(KRANC_GFOFFSET3D(u,-1,1,0) + KRANC_GFOFFSET3D(u,1,-1,0)) + KRANC_GFOFFSET3D(u,1,1,0))*p1o4dxdy; } #endif #ifndef KRANC_DIFF_FUNCTIONS -# define PDstandardfdOrder213(u) ((KRANC_GFOFFSET3D(u,-1,0,-1) - KRANC_GFOFFSET3D(u,-1,0,1) - KRANC_GFOFFSET3D(u,1,0,-1) + KRANC_GFOFFSET3D(u,1,0,1))*p1o4dxdz) +# define PDstandardfdOrder213(u) ((KRANC_GFOFFSET3D(u,-1,0,-1) - 1.*(KRANC_GFOFFSET3D(u,-1,0,1) + KRANC_GFOFFSET3D(u,1,0,-1)) + KRANC_GFOFFSET3D(u,1,0,1))*p1o4dxdz) #else # define PDstandardfdOrder213(u) (PDstandardfdOrder213_impl(u,p1o4dxdz,cdj,cdk)) static CCTK_REAL PDstandardfdOrder213_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1o4dxdz, ptrdiff_t const cdj, ptrdiff_t const cdk) CCTK_ATTRIBUTE_NOINLINE CCTK_ATTRIBUTE_UNUSED; static CCTK_REAL PDstandardfdOrder213_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1o4dxdz, ptrdiff_t const cdj, ptrdiff_t const cdk) { - ptrdiff_t const cdi=sizeof(CCTK_REAL); - return (KRANC_GFOFFSET3D(u,-1,0,-1) - KRANC_GFOFFSET3D(u,-1,0,1) - KRANC_GFOFFSET3D(u,1,0,-1) + KRANC_GFOFFSET3D(u,1,0,1))*p1o4dxdz; + ptrdiff_t const cdi CCTK_ATTRIBUTE_UNUSED = sizeof(CCTK_REAL); + return (KRANC_GFOFFSET3D(u,-1,0,-1) - 1.*(KRANC_GFOFFSET3D(u,-1,0,1) + KRANC_GFOFFSET3D(u,1,0,-1)) + KRANC_GFOFFSET3D(u,1,0,1))*p1o4dxdz; } #endif #ifndef KRANC_DIFF_FUNCTIONS -# define PDstandardfdOrder221(u) ((KRANC_GFOFFSET3D(u,-1,-1,0) - KRANC_GFOFFSET3D(u,-1,1,0) - KRANC_GFOFFSET3D(u,1,-1,0) + KRANC_GFOFFSET3D(u,1,1,0))*p1o4dxdy) +# define PDstandardfdOrder221(u) ((KRANC_GFOFFSET3D(u,-1,-1,0) - 1.*(KRANC_GFOFFSET3D(u,-1,1,0) + KRANC_GFOFFSET3D(u,1,-1,0)) + KRANC_GFOFFSET3D(u,1,1,0))*p1o4dxdy) #else # define PDstandardfdOrder221(u) (PDstandardfdOrder221_impl(u,p1o4dxdy,cdj,cdk)) static CCTK_REAL PDstandardfdOrder221_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1o4dxdy, ptrdiff_t const cdj, ptrdiff_t const cdk) CCTK_ATTRIBUTE_NOINLINE CCTK_ATTRIBUTE_UNUSED; static CCTK_REAL PDstandardfdOrder221_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1o4dxdy, ptrdiff_t const cdj, ptrdiff_t const cdk) { - ptrdiff_t const cdi=sizeof(CCTK_REAL); - return (KRANC_GFOFFSET3D(u,-1,-1,0) - KRANC_GFOFFSET3D(u,-1,1,0) - KRANC_GFOFFSET3D(u,1,-1,0) + KRANC_GFOFFSET3D(u,1,1,0))*p1o4dxdy; + ptrdiff_t const cdi CCTK_ATTRIBUTE_UNUSED = sizeof(CCTK_REAL); + return (KRANC_GFOFFSET3D(u,-1,-1,0) - 1.*(KRANC_GFOFFSET3D(u,-1,1,0) + KRANC_GFOFFSET3D(u,1,-1,0)) + KRANC_GFOFFSET3D(u,1,1,0))*p1o4dxdy; } #endif #ifndef KRANC_DIFF_FUNCTIONS -# define PDstandardfdOrder223(u) ((KRANC_GFOFFSET3D(u,0,-1,-1) - KRANC_GFOFFSET3D(u,0,-1,1) - KRANC_GFOFFSET3D(u,0,1,-1) + KRANC_GFOFFSET3D(u,0,1,1))*p1o4dydz) +# define PDstandardfdOrder223(u) ((KRANC_GFOFFSET3D(u,0,-1,-1) - 1.*(KRANC_GFOFFSET3D(u,0,-1,1) + KRANC_GFOFFSET3D(u,0,1,-1)) + KRANC_GFOFFSET3D(u,0,1,1))*p1o4dydz) #else # define PDstandardfdOrder223(u) (PDstandardfdOrder223_impl(u,p1o4dydz,cdj,cdk)) static CCTK_REAL PDstandardfdOrder223_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1o4dydz, ptrdiff_t const cdj, ptrdiff_t const cdk) CCTK_ATTRIBUTE_NOINLINE CCTK_ATTRIBUTE_UNUSED; static CCTK_REAL PDstandardfdOrder223_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1o4dydz, ptrdiff_t const cdj, ptrdiff_t const cdk) { - ptrdiff_t const cdi=sizeof(CCTK_REAL); - return (KRANC_GFOFFSET3D(u,0,-1,-1) - KRANC_GFOFFSET3D(u,0,-1,1) - KRANC_GFOFFSET3D(u,0,1,-1) + KRANC_GFOFFSET3D(u,0,1,1))*p1o4dydz; + ptrdiff_t const cdi CCTK_ATTRIBUTE_UNUSED = sizeof(CCTK_REAL); + return (KRANC_GFOFFSET3D(u,0,-1,-1) - 1.*(KRANC_GFOFFSET3D(u,0,-1,1) + KRANC_GFOFFSET3D(u,0,1,-1)) + KRANC_GFOFFSET3D(u,0,1,1))*p1o4dydz; } #endif #ifndef KRANC_DIFF_FUNCTIONS -# define PDstandardfdOrder231(u) ((KRANC_GFOFFSET3D(u,-1,0,-1) - KRANC_GFOFFSET3D(u,-1,0,1) - KRANC_GFOFFSET3D(u,1,0,-1) + KRANC_GFOFFSET3D(u,1,0,1))*p1o4dxdz) +# define PDstandardfdOrder231(u) ((KRANC_GFOFFSET3D(u,-1,0,-1) - 1.*(KRANC_GFOFFSET3D(u,-1,0,1) + KRANC_GFOFFSET3D(u,1,0,-1)) + KRANC_GFOFFSET3D(u,1,0,1))*p1o4dxdz) #else # define PDstandardfdOrder231(u) (PDstandardfdOrder231_impl(u,p1o4dxdz,cdj,cdk)) static CCTK_REAL PDstandardfdOrder231_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1o4dxdz, ptrdiff_t const cdj, ptrdiff_t const cdk) CCTK_ATTRIBUTE_NOINLINE CCTK_ATTRIBUTE_UNUSED; static CCTK_REAL PDstandardfdOrder231_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1o4dxdz, ptrdiff_t const cdj, ptrdiff_t const cdk) { - ptrdiff_t const cdi=sizeof(CCTK_REAL); - return (KRANC_GFOFFSET3D(u,-1,0,-1) - KRANC_GFOFFSET3D(u,-1,0,1) - KRANC_GFOFFSET3D(u,1,0,-1) + KRANC_GFOFFSET3D(u,1,0,1))*p1o4dxdz; + ptrdiff_t const cdi CCTK_ATTRIBUTE_UNUSED = sizeof(CCTK_REAL); + return (KRANC_GFOFFSET3D(u,-1,0,-1) - 1.*(KRANC_GFOFFSET3D(u,-1,0,1) + KRANC_GFOFFSET3D(u,1,0,-1)) + KRANC_GFOFFSET3D(u,1,0,1))*p1o4dxdz; } #endif #ifndef KRANC_DIFF_FUNCTIONS -# define PDstandardfdOrder232(u) ((KRANC_GFOFFSET3D(u,0,-1,-1) - KRANC_GFOFFSET3D(u,0,-1,1) - KRANC_GFOFFSET3D(u,0,1,-1) + KRANC_GFOFFSET3D(u,0,1,1))*p1o4dydz) +# define PDstandardfdOrder232(u) ((KRANC_GFOFFSET3D(u,0,-1,-1) - 1.*(KRANC_GFOFFSET3D(u,0,-1,1) + KRANC_GFOFFSET3D(u,0,1,-1)) + KRANC_GFOFFSET3D(u,0,1,1))*p1o4dydz) #else # define PDstandardfdOrder232(u) (PDstandardfdOrder232_impl(u,p1o4dydz,cdj,cdk)) static CCTK_REAL PDstandardfdOrder232_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1o4dydz, ptrdiff_t const cdj, ptrdiff_t const cdk) CCTK_ATTRIBUTE_NOINLINE CCTK_ATTRIBUTE_UNUSED; static CCTK_REAL PDstandardfdOrder232_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1o4dydz, ptrdiff_t const cdj, ptrdiff_t const cdk) { - ptrdiff_t const cdi=sizeof(CCTK_REAL); - return (KRANC_GFOFFSET3D(u,0,-1,-1) - KRANC_GFOFFSET3D(u,0,-1,1) - KRANC_GFOFFSET3D(u,0,1,-1) + KRANC_GFOFFSET3D(u,0,1,1))*p1o4dydz; + ptrdiff_t const cdi CCTK_ATTRIBUTE_UNUSED = sizeof(CCTK_REAL); + return (KRANC_GFOFFSET3D(u,0,-1,-1) - 1.*(KRANC_GFOFFSET3D(u,0,-1,1) + KRANC_GFOFFSET3D(u,0,1,-1)) + KRANC_GFOFFSET3D(u,0,1,1))*p1o4dydz; } #endif #ifndef KRANC_DIFF_FUNCTIONS -# define PDstandardfdOrder412(u) ((-64*(KRANC_GFOFFSET3D(u,-1,1,0) + KRANC_GFOFFSET3D(u,1,-1,0)) + 64*(KRANC_GFOFFSET3D(u,-1,-1,0) + KRANC_GFOFFSET3D(u,1,1,0)) + 8*(KRANC_GFOFFSET3D(u,-1,2,0) + KRANC_GFOFFSET3D(u,1,-2,0) + KRANC_GFOFFSET3D(u,-2,1,0) + KRANC_GFOFFSET3D(u,2,-1,0)) - 8*(KRANC_GFOFFSET3D(u,-1,-2,0) + KRANC_GFOFFSET3D(u,1,2,0) + KRANC_GFOFFSET3D(u,-2,-1,0) + KRANC_GFOFFSET3D(u,2,1,0)) + KRANC_GFOFFSET3D(u,-2,-2,0) - KRANC_GFOFFSET3D(u,-2,2,0) - KRANC_GFOFFSET3D(u,2,-2,0) + KRANC_GFOFFSET3D(u,2,2,0))*p1o144dxdy) +# define PDstandardfdOrder412(u) ((-64.*(KRANC_GFOFFSET3D(u,-1,1,0) + KRANC_GFOFFSET3D(u,1,-1,0)) + 64.*(KRANC_GFOFFSET3D(u,-1,-1,0) + KRANC_GFOFFSET3D(u,1,1,0)) + 8.*(KRANC_GFOFFSET3D(u,-1,2,0) + KRANC_GFOFFSET3D(u,1,-2,0) + KRANC_GFOFFSET3D(u,-2,1,0) + KRANC_GFOFFSET3D(u,2,-1,0)) - 8.*(KRANC_GFOFFSET3D(u,-1,-2,0) + KRANC_GFOFFSET3D(u,1,2,0) + KRANC_GFOFFSET3D(u,-2,-1,0) + KRANC_GFOFFSET3D(u,2,1,0)) + KRANC_GFOFFSET3D(u,-2,-2,0) - 1.*(KRANC_GFOFFSET3D(u,-2,2,0) + KRANC_GFOFFSET3D(u,2,-2,0)) + KRANC_GFOFFSET3D(u,2,2,0))*p1o144dxdy) #else # define PDstandardfdOrder412(u) (PDstandardfdOrder412_impl(u,p1o144dxdy,cdj,cdk)) static CCTK_REAL PDstandardfdOrder412_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1o144dxdy, ptrdiff_t const cdj, ptrdiff_t const cdk) CCTK_ATTRIBUTE_NOINLINE CCTK_ATTRIBUTE_UNUSED; static CCTK_REAL PDstandardfdOrder412_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1o144dxdy, ptrdiff_t const cdj, ptrdiff_t const cdk) { - ptrdiff_t const cdi=sizeof(CCTK_REAL); - return (-64*(KRANC_GFOFFSET3D(u,-1,1,0) + KRANC_GFOFFSET3D(u,1,-1,0)) + 64*(KRANC_GFOFFSET3D(u,-1,-1,0) + KRANC_GFOFFSET3D(u,1,1,0)) + 8*(KRANC_GFOFFSET3D(u,-1,2,0) + KRANC_GFOFFSET3D(u,1,-2,0) + KRANC_GFOFFSET3D(u,-2,1,0) + KRANC_GFOFFSET3D(u,2,-1,0)) - 8*(KRANC_GFOFFSET3D(u,-1,-2,0) + KRANC_GFOFFSET3D(u,1,2,0) + KRANC_GFOFFSET3D(u,-2,-1,0) + KRANC_GFOFFSET3D(u,2,1,0)) + KRANC_GFOFFSET3D(u,-2,-2,0) - KRANC_GFOFFSET3D(u,-2,2,0) - KRANC_GFOFFSET3D(u,2,-2,0) + KRANC_GFOFFSET3D(u,2,2,0))*p1o144dxdy; + ptrdiff_t const cdi CCTK_ATTRIBUTE_UNUSED = sizeof(CCTK_REAL); + return (-64.*(KRANC_GFOFFSET3D(u,-1,1,0) + KRANC_GFOFFSET3D(u,1,-1,0)) + 64.*(KRANC_GFOFFSET3D(u,-1,-1,0) + KRANC_GFOFFSET3D(u,1,1,0)) + 8.*(KRANC_GFOFFSET3D(u,-1,2,0) + KRANC_GFOFFSET3D(u,1,-2,0) + KRANC_GFOFFSET3D(u,-2,1,0) + KRANC_GFOFFSET3D(u,2,-1,0)) - 8.*(KRANC_GFOFFSET3D(u,-1,-2,0) + KRANC_GFOFFSET3D(u,1,2,0) + KRANC_GFOFFSET3D(u,-2,-1,0) + KRANC_GFOFFSET3D(u,2,1,0)) + KRANC_GFOFFSET3D(u,-2,-2,0) - 1.*(KRANC_GFOFFSET3D(u,-2,2,0) + KRANC_GFOFFSET3D(u,2,-2,0)) + KRANC_GFOFFSET3D(u,2,2,0))*p1o144dxdy; } #endif #ifndef KRANC_DIFF_FUNCTIONS -# define PDstandardfdOrder413(u) ((-64*(KRANC_GFOFFSET3D(u,-1,0,1) + KRANC_GFOFFSET3D(u,1,0,-1)) + 64*(KRANC_GFOFFSET3D(u,-1,0,-1) + KRANC_GFOFFSET3D(u,1,0,1)) + 8*(KRANC_GFOFFSET3D(u,-1,0,2) + KRANC_GFOFFSET3D(u,1,0,-2) + KRANC_GFOFFSET3D(u,-2,0,1) + KRANC_GFOFFSET3D(u,2,0,-1)) - 8*(KRANC_GFOFFSET3D(u,-1,0,-2) + KRANC_GFOFFSET3D(u,1,0,2) + KRANC_GFOFFSET3D(u,-2,0,-1) + KRANC_GFOFFSET3D(u,2,0,1)) + KRANC_GFOFFSET3D(u,-2,0,-2) - KRANC_GFOFFSET3D(u,-2,0,2) - KRANC_GFOFFSET3D(u,2,0,-2) + KRANC_GFOFFSET3D(u,2,0,2))*p1o144dxdz) +# define PDstandardfdOrder413(u) ((-64.*(KRANC_GFOFFSET3D(u,-1,0,1) + KRANC_GFOFFSET3D(u,1,0,-1)) + 64.*(KRANC_GFOFFSET3D(u,-1,0,-1) + KRANC_GFOFFSET3D(u,1,0,1)) + 8.*(KRANC_GFOFFSET3D(u,-1,0,2) + KRANC_GFOFFSET3D(u,1,0,-2) + KRANC_GFOFFSET3D(u,-2,0,1) + KRANC_GFOFFSET3D(u,2,0,-1)) - 8.*(KRANC_GFOFFSET3D(u,-1,0,-2) + KRANC_GFOFFSET3D(u,1,0,2) + KRANC_GFOFFSET3D(u,-2,0,-1) + KRANC_GFOFFSET3D(u,2,0,1)) + KRANC_GFOFFSET3D(u,-2,0,-2) - 1.*(KRANC_GFOFFSET3D(u,-2,0,2) + KRANC_GFOFFSET3D(u,2,0,-2)) + KRANC_GFOFFSET3D(u,2,0,2))*p1o144dxdz) #else # define PDstandardfdOrder413(u) (PDstandardfdOrder413_impl(u,p1o144dxdz,cdj,cdk)) static CCTK_REAL PDstandardfdOrder413_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1o144dxdz, ptrdiff_t const cdj, ptrdiff_t const cdk) CCTK_ATTRIBUTE_NOINLINE CCTK_ATTRIBUTE_UNUSED; static CCTK_REAL PDstandardfdOrder413_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1o144dxdz, ptrdiff_t const cdj, ptrdiff_t const cdk) { - ptrdiff_t const cdi=sizeof(CCTK_REAL); - return (-64*(KRANC_GFOFFSET3D(u,-1,0,1) + KRANC_GFOFFSET3D(u,1,0,-1)) + 64*(KRANC_GFOFFSET3D(u,-1,0,-1) + KRANC_GFOFFSET3D(u,1,0,1)) + 8*(KRANC_GFOFFSET3D(u,-1,0,2) + KRANC_GFOFFSET3D(u,1,0,-2) + KRANC_GFOFFSET3D(u,-2,0,1) + KRANC_GFOFFSET3D(u,2,0,-1)) - 8*(KRANC_GFOFFSET3D(u,-1,0,-2) + KRANC_GFOFFSET3D(u,1,0,2) + KRANC_GFOFFSET3D(u,-2,0,-1) + KRANC_GFOFFSET3D(u,2,0,1)) + KRANC_GFOFFSET3D(u,-2,0,-2) - KRANC_GFOFFSET3D(u,-2,0,2) - KRANC_GFOFFSET3D(u,2,0,-2) + KRANC_GFOFFSET3D(u,2,0,2))*p1o144dxdz; + ptrdiff_t const cdi CCTK_ATTRIBUTE_UNUSED = sizeof(CCTK_REAL); + return (-64.*(KRANC_GFOFFSET3D(u,-1,0,1) + KRANC_GFOFFSET3D(u,1,0,-1)) + 64.*(KRANC_GFOFFSET3D(u,-1,0,-1) + KRANC_GFOFFSET3D(u,1,0,1)) + 8.*(KRANC_GFOFFSET3D(u,-1,0,2) + KRANC_GFOFFSET3D(u,1,0,-2) + KRANC_GFOFFSET3D(u,-2,0,1) + KRANC_GFOFFSET3D(u,2,0,-1)) - 8.*(KRANC_GFOFFSET3D(u,-1,0,-2) + KRANC_GFOFFSET3D(u,1,0,2) + KRANC_GFOFFSET3D(u,-2,0,-1) + KRANC_GFOFFSET3D(u,2,0,1)) + KRANC_GFOFFSET3D(u,-2,0,-2) - 1.*(KRANC_GFOFFSET3D(u,-2,0,2) + KRANC_GFOFFSET3D(u,2,0,-2)) + KRANC_GFOFFSET3D(u,2,0,2))*p1o144dxdz; } #endif #ifndef KRANC_DIFF_FUNCTIONS -# define PDstandardfdOrder421(u) ((-64*(KRANC_GFOFFSET3D(u,-1,1,0) + KRANC_GFOFFSET3D(u,1,-1,0)) + 64*(KRANC_GFOFFSET3D(u,-1,-1,0) + KRANC_GFOFFSET3D(u,1,1,0)) + 8*(KRANC_GFOFFSET3D(u,-1,2,0) + KRANC_GFOFFSET3D(u,1,-2,0) + KRANC_GFOFFSET3D(u,-2,1,0) + KRANC_GFOFFSET3D(u,2,-1,0)) - 8*(KRANC_GFOFFSET3D(u,-1,-2,0) + KRANC_GFOFFSET3D(u,1,2,0) + KRANC_GFOFFSET3D(u,-2,-1,0) + KRANC_GFOFFSET3D(u,2,1,0)) + KRANC_GFOFFSET3D(u,-2,-2,0) - KRANC_GFOFFSET3D(u,-2,2,0) - KRANC_GFOFFSET3D(u,2,-2,0) + KRANC_GFOFFSET3D(u,2,2,0))*p1o144dxdy) +# define PDstandardfdOrder421(u) ((-64.*(KRANC_GFOFFSET3D(u,-1,1,0) + KRANC_GFOFFSET3D(u,1,-1,0)) + 64.*(KRANC_GFOFFSET3D(u,-1,-1,0) + KRANC_GFOFFSET3D(u,1,1,0)) + 8.*(KRANC_GFOFFSET3D(u,-1,2,0) + KRANC_GFOFFSET3D(u,1,-2,0) + KRANC_GFOFFSET3D(u,-2,1,0) + KRANC_GFOFFSET3D(u,2,-1,0)) - 8.*(KRANC_GFOFFSET3D(u,-1,-2,0) + KRANC_GFOFFSET3D(u,1,2,0) + KRANC_GFOFFSET3D(u,-2,-1,0) + KRANC_GFOFFSET3D(u,2,1,0)) + KRANC_GFOFFSET3D(u,-2,-2,0) - 1.*(KRANC_GFOFFSET3D(u,-2,2,0) + KRANC_GFOFFSET3D(u,2,-2,0)) + KRANC_GFOFFSET3D(u,2,2,0))*p1o144dxdy) #else # define PDstandardfdOrder421(u) (PDstandardfdOrder421_impl(u,p1o144dxdy,cdj,cdk)) static CCTK_REAL PDstandardfdOrder421_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1o144dxdy, ptrdiff_t const cdj, ptrdiff_t const cdk) CCTK_ATTRIBUTE_NOINLINE CCTK_ATTRIBUTE_UNUSED; static CCTK_REAL PDstandardfdOrder421_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1o144dxdy, ptrdiff_t const cdj, ptrdiff_t const cdk) { - ptrdiff_t const cdi=sizeof(CCTK_REAL); - return (-64*(KRANC_GFOFFSET3D(u,-1,1,0) + KRANC_GFOFFSET3D(u,1,-1,0)) + 64*(KRANC_GFOFFSET3D(u,-1,-1,0) + KRANC_GFOFFSET3D(u,1,1,0)) + 8*(KRANC_GFOFFSET3D(u,-1,2,0) + KRANC_GFOFFSET3D(u,1,-2,0) + KRANC_GFOFFSET3D(u,-2,1,0) + KRANC_GFOFFSET3D(u,2,-1,0)) - 8*(KRANC_GFOFFSET3D(u,-1,-2,0) + KRANC_GFOFFSET3D(u,1,2,0) + KRANC_GFOFFSET3D(u,-2,-1,0) + KRANC_GFOFFSET3D(u,2,1,0)) + KRANC_GFOFFSET3D(u,-2,-2,0) - KRANC_GFOFFSET3D(u,-2,2,0) - KRANC_GFOFFSET3D(u,2,-2,0) + KRANC_GFOFFSET3D(u,2,2,0))*p1o144dxdy; + ptrdiff_t const cdi CCTK_ATTRIBUTE_UNUSED = sizeof(CCTK_REAL); + return (-64.*(KRANC_GFOFFSET3D(u,-1,1,0) + KRANC_GFOFFSET3D(u,1,-1,0)) + 64.*(KRANC_GFOFFSET3D(u,-1,-1,0) + KRANC_GFOFFSET3D(u,1,1,0)) + 8.*(KRANC_GFOFFSET3D(u,-1,2,0) + KRANC_GFOFFSET3D(u,1,-2,0) + KRANC_GFOFFSET3D(u,-2,1,0) + KRANC_GFOFFSET3D(u,2,-1,0)) - 8.*(KRANC_GFOFFSET3D(u,-1,-2,0) + KRANC_GFOFFSET3D(u,1,2,0) + KRANC_GFOFFSET3D(u,-2,-1,0) + KRANC_GFOFFSET3D(u,2,1,0)) + KRANC_GFOFFSET3D(u,-2,-2,0) - 1.*(KRANC_GFOFFSET3D(u,-2,2,0) + KRANC_GFOFFSET3D(u,2,-2,0)) + KRANC_GFOFFSET3D(u,2,2,0))*p1o144dxdy; } #endif #ifndef KRANC_DIFF_FUNCTIONS -# define PDstandardfdOrder423(u) ((-64*(KRANC_GFOFFSET3D(u,0,-1,1) + KRANC_GFOFFSET3D(u,0,1,-1)) + 64*(KRANC_GFOFFSET3D(u,0,-1,-1) + KRANC_GFOFFSET3D(u,0,1,1)) + 8*(KRANC_GFOFFSET3D(u,0,-1,2) + KRANC_GFOFFSET3D(u,0,1,-2) + KRANC_GFOFFSET3D(u,0,-2,1) + KRANC_GFOFFSET3D(u,0,2,-1)) - 8*(KRANC_GFOFFSET3D(u,0,-1,-2) + KRANC_GFOFFSET3D(u,0,1,2) + KRANC_GFOFFSET3D(u,0,-2,-1) + KRANC_GFOFFSET3D(u,0,2,1)) + KRANC_GFOFFSET3D(u,0,-2,-2) - KRANC_GFOFFSET3D(u,0,-2,2) - KRANC_GFOFFSET3D(u,0,2,-2) + KRANC_GFOFFSET3D(u,0,2,2))*p1o144dydz) +# define PDstandardfdOrder423(u) ((-64.*(KRANC_GFOFFSET3D(u,0,-1,1) + KRANC_GFOFFSET3D(u,0,1,-1)) + 64.*(KRANC_GFOFFSET3D(u,0,-1,-1) + KRANC_GFOFFSET3D(u,0,1,1)) + 8.*(KRANC_GFOFFSET3D(u,0,-1,2) + KRANC_GFOFFSET3D(u,0,1,-2) + KRANC_GFOFFSET3D(u,0,-2,1) + KRANC_GFOFFSET3D(u,0,2,-1)) - 8.*(KRANC_GFOFFSET3D(u,0,-1,-2) + KRANC_GFOFFSET3D(u,0,1,2) + KRANC_GFOFFSET3D(u,0,-2,-1) + KRANC_GFOFFSET3D(u,0,2,1)) + KRANC_GFOFFSET3D(u,0,-2,-2) - 1.*(KRANC_GFOFFSET3D(u,0,-2,2) + KRANC_GFOFFSET3D(u,0,2,-2)) + KRANC_GFOFFSET3D(u,0,2,2))*p1o144dydz) #else # define PDstandardfdOrder423(u) (PDstandardfdOrder423_impl(u,p1o144dydz,cdj,cdk)) static CCTK_REAL PDstandardfdOrder423_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1o144dydz, ptrdiff_t const cdj, ptrdiff_t const cdk) CCTK_ATTRIBUTE_NOINLINE CCTK_ATTRIBUTE_UNUSED; static CCTK_REAL PDstandardfdOrder423_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1o144dydz, ptrdiff_t const cdj, ptrdiff_t const cdk) { - ptrdiff_t const cdi=sizeof(CCTK_REAL); - return (-64*(KRANC_GFOFFSET3D(u,0,-1,1) + KRANC_GFOFFSET3D(u,0,1,-1)) + 64*(KRANC_GFOFFSET3D(u,0,-1,-1) + KRANC_GFOFFSET3D(u,0,1,1)) + 8*(KRANC_GFOFFSET3D(u,0,-1,2) + KRANC_GFOFFSET3D(u,0,1,-2) + KRANC_GFOFFSET3D(u,0,-2,1) + KRANC_GFOFFSET3D(u,0,2,-1)) - 8*(KRANC_GFOFFSET3D(u,0,-1,-2) + KRANC_GFOFFSET3D(u,0,1,2) + KRANC_GFOFFSET3D(u,0,-2,-1) + KRANC_GFOFFSET3D(u,0,2,1)) + KRANC_GFOFFSET3D(u,0,-2,-2) - KRANC_GFOFFSET3D(u,0,-2,2) - KRANC_GFOFFSET3D(u,0,2,-2) + KRANC_GFOFFSET3D(u,0,2,2))*p1o144dydz; + ptrdiff_t const cdi CCTK_ATTRIBUTE_UNUSED = sizeof(CCTK_REAL); + return (-64.*(KRANC_GFOFFSET3D(u,0,-1,1) + KRANC_GFOFFSET3D(u,0,1,-1)) + 64.*(KRANC_GFOFFSET3D(u,0,-1,-1) + KRANC_GFOFFSET3D(u,0,1,1)) + 8.*(KRANC_GFOFFSET3D(u,0,-1,2) + KRANC_GFOFFSET3D(u,0,1,-2) + KRANC_GFOFFSET3D(u,0,-2,1) + KRANC_GFOFFSET3D(u,0,2,-1)) - 8.*(KRANC_GFOFFSET3D(u,0,-1,-2) + KRANC_GFOFFSET3D(u,0,1,2) + KRANC_GFOFFSET3D(u,0,-2,-1) + KRANC_GFOFFSET3D(u,0,2,1)) + KRANC_GFOFFSET3D(u,0,-2,-2) - 1.*(KRANC_GFOFFSET3D(u,0,-2,2) + KRANC_GFOFFSET3D(u,0,2,-2)) + KRANC_GFOFFSET3D(u,0,2,2))*p1o144dydz; } #endif #ifndef KRANC_DIFF_FUNCTIONS -# define PDstandardfdOrder431(u) ((-64*(KRANC_GFOFFSET3D(u,-1,0,1) + KRANC_GFOFFSET3D(u,1,0,-1)) + 64*(KRANC_GFOFFSET3D(u,-1,0,-1) + KRANC_GFOFFSET3D(u,1,0,1)) + 8*(KRANC_GFOFFSET3D(u,-1,0,2) + KRANC_GFOFFSET3D(u,1,0,-2) + KRANC_GFOFFSET3D(u,-2,0,1) + KRANC_GFOFFSET3D(u,2,0,-1)) - 8*(KRANC_GFOFFSET3D(u,-1,0,-2) + KRANC_GFOFFSET3D(u,1,0,2) + KRANC_GFOFFSET3D(u,-2,0,-1) + KRANC_GFOFFSET3D(u,2,0,1)) + KRANC_GFOFFSET3D(u,-2,0,-2) - KRANC_GFOFFSET3D(u,-2,0,2) - KRANC_GFOFFSET3D(u,2,0,-2) + KRANC_GFOFFSET3D(u,2,0,2))*p1o144dxdz) +# define PDstandardfdOrder431(u) ((-64.*(KRANC_GFOFFSET3D(u,-1,0,1) + KRANC_GFOFFSET3D(u,1,0,-1)) + 64.*(KRANC_GFOFFSET3D(u,-1,0,-1) + KRANC_GFOFFSET3D(u,1,0,1)) + 8.*(KRANC_GFOFFSET3D(u,-1,0,2) + KRANC_GFOFFSET3D(u,1,0,-2) + KRANC_GFOFFSET3D(u,-2,0,1) + KRANC_GFOFFSET3D(u,2,0,-1)) - 8.*(KRANC_GFOFFSET3D(u,-1,0,-2) + KRANC_GFOFFSET3D(u,1,0,2) + KRANC_GFOFFSET3D(u,-2,0,-1) + KRANC_GFOFFSET3D(u,2,0,1)) + KRANC_GFOFFSET3D(u,-2,0,-2) - 1.*(KRANC_GFOFFSET3D(u,-2,0,2) + KRANC_GFOFFSET3D(u,2,0,-2)) + KRANC_GFOFFSET3D(u,2,0,2))*p1o144dxdz) #else # define PDstandardfdOrder431(u) (PDstandardfdOrder431_impl(u,p1o144dxdz,cdj,cdk)) static CCTK_REAL PDstandardfdOrder431_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1o144dxdz, ptrdiff_t const cdj, ptrdiff_t const cdk) CCTK_ATTRIBUTE_NOINLINE CCTK_ATTRIBUTE_UNUSED; static CCTK_REAL PDstandardfdOrder431_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1o144dxdz, ptrdiff_t const cdj, ptrdiff_t const cdk) { - ptrdiff_t const cdi=sizeof(CCTK_REAL); - return (-64*(KRANC_GFOFFSET3D(u,-1,0,1) + KRANC_GFOFFSET3D(u,1,0,-1)) + 64*(KRANC_GFOFFSET3D(u,-1,0,-1) + KRANC_GFOFFSET3D(u,1,0,1)) + 8*(KRANC_GFOFFSET3D(u,-1,0,2) + KRANC_GFOFFSET3D(u,1,0,-2) + KRANC_GFOFFSET3D(u,-2,0,1) + KRANC_GFOFFSET3D(u,2,0,-1)) - 8*(KRANC_GFOFFSET3D(u,-1,0,-2) + KRANC_GFOFFSET3D(u,1,0,2) + KRANC_GFOFFSET3D(u,-2,0,-1) + KRANC_GFOFFSET3D(u,2,0,1)) + KRANC_GFOFFSET3D(u,-2,0,-2) - KRANC_GFOFFSET3D(u,-2,0,2) - KRANC_GFOFFSET3D(u,2,0,-2) + KRANC_GFOFFSET3D(u,2,0,2))*p1o144dxdz; + ptrdiff_t const cdi CCTK_ATTRIBUTE_UNUSED = sizeof(CCTK_REAL); + return (-64.*(KRANC_GFOFFSET3D(u,-1,0,1) + KRANC_GFOFFSET3D(u,1,0,-1)) + 64.*(KRANC_GFOFFSET3D(u,-1,0,-1) + KRANC_GFOFFSET3D(u,1,0,1)) + 8.*(KRANC_GFOFFSET3D(u,-1,0,2) + KRANC_GFOFFSET3D(u,1,0,-2) + KRANC_GFOFFSET3D(u,-2,0,1) + KRANC_GFOFFSET3D(u,2,0,-1)) - 8.*(KRANC_GFOFFSET3D(u,-1,0,-2) + KRANC_GFOFFSET3D(u,1,0,2) + KRANC_GFOFFSET3D(u,-2,0,-1) + KRANC_GFOFFSET3D(u,2,0,1)) + KRANC_GFOFFSET3D(u,-2,0,-2) - 1.*(KRANC_GFOFFSET3D(u,-2,0,2) + KRANC_GFOFFSET3D(u,2,0,-2)) + KRANC_GFOFFSET3D(u,2,0,2))*p1o144dxdz; } #endif #ifndef KRANC_DIFF_FUNCTIONS -# define PDstandardfdOrder432(u) ((-64*(KRANC_GFOFFSET3D(u,0,-1,1) + KRANC_GFOFFSET3D(u,0,1,-1)) + 64*(KRANC_GFOFFSET3D(u,0,-1,-1) + KRANC_GFOFFSET3D(u,0,1,1)) + 8*(KRANC_GFOFFSET3D(u,0,-1,2) + KRANC_GFOFFSET3D(u,0,1,-2) + KRANC_GFOFFSET3D(u,0,-2,1) + KRANC_GFOFFSET3D(u,0,2,-1)) - 8*(KRANC_GFOFFSET3D(u,0,-1,-2) + KRANC_GFOFFSET3D(u,0,1,2) + KRANC_GFOFFSET3D(u,0,-2,-1) + KRANC_GFOFFSET3D(u,0,2,1)) + KRANC_GFOFFSET3D(u,0,-2,-2) - KRANC_GFOFFSET3D(u,0,-2,2) - KRANC_GFOFFSET3D(u,0,2,-2) + KRANC_GFOFFSET3D(u,0,2,2))*p1o144dydz) +# define PDstandardfdOrder432(u) ((-64.*(KRANC_GFOFFSET3D(u,0,-1,1) + KRANC_GFOFFSET3D(u,0,1,-1)) + 64.*(KRANC_GFOFFSET3D(u,0,-1,-1) + KRANC_GFOFFSET3D(u,0,1,1)) + 8.*(KRANC_GFOFFSET3D(u,0,-1,2) + KRANC_GFOFFSET3D(u,0,1,-2) + KRANC_GFOFFSET3D(u,0,-2,1) + KRANC_GFOFFSET3D(u,0,2,-1)) - 8.*(KRANC_GFOFFSET3D(u,0,-1,-2) + KRANC_GFOFFSET3D(u,0,1,2) + KRANC_GFOFFSET3D(u,0,-2,-1) + KRANC_GFOFFSET3D(u,0,2,1)) + KRANC_GFOFFSET3D(u,0,-2,-2) - 1.*(KRANC_GFOFFSET3D(u,0,-2,2) + KRANC_GFOFFSET3D(u,0,2,-2)) + KRANC_GFOFFSET3D(u,0,2,2))*p1o144dydz) #else # define PDstandardfdOrder432(u) (PDstandardfdOrder432_impl(u,p1o144dydz,cdj,cdk)) static CCTK_REAL PDstandardfdOrder432_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1o144dydz, ptrdiff_t const cdj, ptrdiff_t const cdk) CCTK_ATTRIBUTE_NOINLINE CCTK_ATTRIBUTE_UNUSED; static CCTK_REAL PDstandardfdOrder432_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1o144dydz, ptrdiff_t const cdj, ptrdiff_t const cdk) { - ptrdiff_t const cdi=sizeof(CCTK_REAL); - return (-64*(KRANC_GFOFFSET3D(u,0,-1,1) + KRANC_GFOFFSET3D(u,0,1,-1)) + 64*(KRANC_GFOFFSET3D(u,0,-1,-1) + KRANC_GFOFFSET3D(u,0,1,1)) + 8*(KRANC_GFOFFSET3D(u,0,-1,2) + KRANC_GFOFFSET3D(u,0,1,-2) + KRANC_GFOFFSET3D(u,0,-2,1) + KRANC_GFOFFSET3D(u,0,2,-1)) - 8*(KRANC_GFOFFSET3D(u,0,-1,-2) + KRANC_GFOFFSET3D(u,0,1,2) + KRANC_GFOFFSET3D(u,0,-2,-1) + KRANC_GFOFFSET3D(u,0,2,1)) + KRANC_GFOFFSET3D(u,0,-2,-2) - KRANC_GFOFFSET3D(u,0,-2,2) - KRANC_GFOFFSET3D(u,0,2,-2) + KRANC_GFOFFSET3D(u,0,2,2))*p1o144dydz; + ptrdiff_t const cdi CCTK_ATTRIBUTE_UNUSED = sizeof(CCTK_REAL); + return (-64.*(KRANC_GFOFFSET3D(u,0,-1,1) + KRANC_GFOFFSET3D(u,0,1,-1)) + 64.*(KRANC_GFOFFSET3D(u,0,-1,-1) + KRANC_GFOFFSET3D(u,0,1,1)) + 8.*(KRANC_GFOFFSET3D(u,0,-1,2) + KRANC_GFOFFSET3D(u,0,1,-2) + KRANC_GFOFFSET3D(u,0,-2,1) + KRANC_GFOFFSET3D(u,0,2,-1)) - 8.*(KRANC_GFOFFSET3D(u,0,-1,-2) + KRANC_GFOFFSET3D(u,0,1,2) + KRANC_GFOFFSET3D(u,0,-2,-1) + KRANC_GFOFFSET3D(u,0,2,1)) + KRANC_GFOFFSET3D(u,0,-2,-2) - 1.*(KRANC_GFOFFSET3D(u,0,-2,2) + KRANC_GFOFFSET3D(u,0,2,-2)) + KRANC_GFOFFSET3D(u,0,2,2))*p1o144dydz; } #endif #ifndef KRANC_DIFF_FUNCTIONS -# define PDonesided2nd1(u) ((3*KRANC_GFOFFSET3D(u,0,0,0) + KRANC_GFOFFSET3D(u,2*dir1,0,0) - 4*KRANC_GFOFFSET3D(u,dir1,0,0))*pm1o2dx*dir1) +# define PDonesided2nd1(u) ((3.*KRANC_GFOFFSET3D(u,0,0,0) + KRANC_GFOFFSET3D(u,2*dir1,0,0) - 4.*KRANC_GFOFFSET3D(u,dir1,0,0))*pm1o2dx*dir1) #else # define PDonesided2nd1(u) (PDonesided2nd1_impl(u,pm1o2dx,cdj,cdk,dir1,dir2,dir3)) static CCTK_REAL PDonesided2nd1_impl(CCTK_REAL const* restrict const u, CCTK_REAL const pm1o2dx, ptrdiff_t const cdj, ptrdiff_t const cdk, ptrdiff_t const dir1, ptrdiff_t const dir2, ptrdiff_t const dir3) CCTK_ATTRIBUTE_NOINLINE CCTK_ATTRIBUTE_UNUSED; static CCTK_REAL PDonesided2nd1_impl(CCTK_REAL const* restrict const u, CCTK_REAL const pm1o2dx, ptrdiff_t const cdj, ptrdiff_t const cdk, ptrdiff_t const dir1, ptrdiff_t const dir2, ptrdiff_t const dir3) { - ptrdiff_t const cdi=sizeof(CCTK_REAL); - return (3*KRANC_GFOFFSET3D(u,0,0,0) + KRANC_GFOFFSET3D(u,2*dir1,0,0) - 4*KRANC_GFOFFSET3D(u,dir1,0,0))*pm1o2dx*dir1; + ptrdiff_t const cdi CCTK_ATTRIBUTE_UNUSED = sizeof(CCTK_REAL); + return (3.*KRANC_GFOFFSET3D(u,0,0,0) + KRANC_GFOFFSET3D(u,2*dir1,0,0) - 4.*KRANC_GFOFFSET3D(u,dir1,0,0))*pm1o2dx*dir1; } #endif #ifndef KRANC_DIFF_FUNCTIONS -# define PDonesided2nd2(u) ((3*KRANC_GFOFFSET3D(u,0,0,0) + KRANC_GFOFFSET3D(u,0,2*dir2,0) - 4*KRANC_GFOFFSET3D(u,0,dir2,0))*pm1o2dy*dir2) +# define PDonesided2nd2(u) ((3.*KRANC_GFOFFSET3D(u,0,0,0) + KRANC_GFOFFSET3D(u,0,2*dir2,0) - 4.*KRANC_GFOFFSET3D(u,0,dir2,0))*pm1o2dy*dir2) #else # define PDonesided2nd2(u) (PDonesided2nd2_impl(u,pm1o2dy,cdj,cdk,dir1,dir2,dir3)) static CCTK_REAL PDonesided2nd2_impl(CCTK_REAL const* restrict const u, CCTK_REAL const pm1o2dy, ptrdiff_t const cdj, ptrdiff_t const cdk, ptrdiff_t const dir1, ptrdiff_t const dir2, ptrdiff_t const dir3) CCTK_ATTRIBUTE_NOINLINE CCTK_ATTRIBUTE_UNUSED; static CCTK_REAL PDonesided2nd2_impl(CCTK_REAL const* restrict const u, CCTK_REAL const pm1o2dy, ptrdiff_t const cdj, ptrdiff_t const cdk, ptrdiff_t const dir1, ptrdiff_t const dir2, ptrdiff_t const dir3) { - ptrdiff_t const cdi=sizeof(CCTK_REAL); - return (3*KRANC_GFOFFSET3D(u,0,0,0) + KRANC_GFOFFSET3D(u,0,2*dir2,0) - 4*KRANC_GFOFFSET3D(u,0,dir2,0))*pm1o2dy*dir2; + ptrdiff_t const cdi CCTK_ATTRIBUTE_UNUSED = sizeof(CCTK_REAL); + return (3.*KRANC_GFOFFSET3D(u,0,0,0) + KRANC_GFOFFSET3D(u,0,2*dir2,0) - 4.*KRANC_GFOFFSET3D(u,0,dir2,0))*pm1o2dy*dir2; } #endif #ifndef KRANC_DIFF_FUNCTIONS -# define PDonesided2nd3(u) ((3*KRANC_GFOFFSET3D(u,0,0,0) + KRANC_GFOFFSET3D(u,0,0,2*dir3) - 4*KRANC_GFOFFSET3D(u,0,0,dir3))*pm1o2dz*dir3) +# define PDonesided2nd3(u) ((3.*KRANC_GFOFFSET3D(u,0,0,0) + KRANC_GFOFFSET3D(u,0,0,2*dir3) - 4.*KRANC_GFOFFSET3D(u,0,0,dir3))*pm1o2dz*dir3) #else # define PDonesided2nd3(u) (PDonesided2nd3_impl(u,pm1o2dz,cdj,cdk,dir1,dir2,dir3)) static CCTK_REAL PDonesided2nd3_impl(CCTK_REAL const* restrict const u, CCTK_REAL const pm1o2dz, ptrdiff_t const cdj, ptrdiff_t const cdk, ptrdiff_t const dir1, ptrdiff_t const dir2, ptrdiff_t const dir3) CCTK_ATTRIBUTE_NOINLINE CCTK_ATTRIBUTE_UNUSED; static CCTK_REAL PDonesided2nd3_impl(CCTK_REAL const* restrict const u, CCTK_REAL const pm1o2dz, ptrdiff_t const cdj, ptrdiff_t const cdk, ptrdiff_t const dir1, ptrdiff_t const dir2, ptrdiff_t const dir3) { - ptrdiff_t const cdi=sizeof(CCTK_REAL); - return (3*KRANC_GFOFFSET3D(u,0,0,0) + KRANC_GFOFFSET3D(u,0,0,2*dir3) - 4*KRANC_GFOFFSET3D(u,0,0,dir3))*pm1o2dz*dir3; + ptrdiff_t const cdi CCTK_ATTRIBUTE_UNUSED = sizeof(CCTK_REAL); + return (3.*KRANC_GFOFFSET3D(u,0,0,0) + KRANC_GFOFFSET3D(u,0,0,2*dir3) - 4.*KRANC_GFOFFSET3D(u,0,0,dir3))*pm1o2dz*dir3; } #endif #ifndef KRANC_DIFF_FUNCTIONS -# define DissfdOrder2(u) (-(p1odxdydz*diss*((6*KRANC_GFOFFSET3D(u,0,0,0) - 4*(KRANC_GFOFFSET3D(u,0,0,-1) + KRANC_GFOFFSET3D(u,0,0,1)) + KRANC_GFOFFSET3D(u,0,0,-2) + KRANC_GFOFFSET3D(u,0,0,2))*INV(dxi*dyi) + ((6*KRANC_GFOFFSET3D(u,0,0,0) - 4*(KRANC_GFOFFSET3D(u,0,-1,0) + KRANC_GFOFFSET3D(u,0,1,0)) + KRANC_GFOFFSET3D(u,0,-2,0) + KRANC_GFOFFSET3D(u,0,2,0))*INV(dxi) + (6*KRANC_GFOFFSET3D(u,0,0,0) - 4*(KRANC_GFOFFSET3D(u,-1,0,0) + KRANC_GFOFFSET3D(u,1,0,0)) + KRANC_GFOFFSET3D(u,-2,0,0) + KRANC_GFOFFSET3D(u,2,0,0))*INV(dyi))*INV(dzi)))) +# define DissfdOrder2(u) (-1.*p1odxdydz*diss*((6.*KRANC_GFOFFSET3D(u,0,0,0) - 4.*(KRANC_GFOFFSET3D(u,0,0,-1) + KRANC_GFOFFSET3D(u,0,0,1)) + KRANC_GFOFFSET3D(u,0,0,-2) + KRANC_GFOFFSET3D(u,0,0,2))*INV(dxi*dyi) + ((6.*KRANC_GFOFFSET3D(u,0,0,0) - 4.*(KRANC_GFOFFSET3D(u,0,-1,0) + KRANC_GFOFFSET3D(u,0,1,0)) + KRANC_GFOFFSET3D(u,0,-2,0) + KRANC_GFOFFSET3D(u,0,2,0))*INV(dxi) + (6.*KRANC_GFOFFSET3D(u,0,0,0) - 4.*(KRANC_GFOFFSET3D(u,-1,0,0) + KRANC_GFOFFSET3D(u,1,0,0)) + KRANC_GFOFFSET3D(u,-2,0,0) + KRANC_GFOFFSET3D(u,2,0,0))*INV(dyi))*INV(dzi))) #else # define DissfdOrder2(u) (DissfdOrder2_impl(u,p1odxdydz,cdj,cdk)) static CCTK_REAL DissfdOrder2_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1odxdydz, ptrdiff_t const cdj, ptrdiff_t const cdk) CCTK_ATTRIBUTE_NOINLINE CCTK_ATTRIBUTE_UNUSED; static CCTK_REAL DissfdOrder2_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1odxdydz, ptrdiff_t const cdj, ptrdiff_t const cdk) { - ptrdiff_t const cdi=sizeof(CCTK_REAL); - return -(p1odxdydz*diss*((6*KRANC_GFOFFSET3D(u,0,0,0) - 4*(KRANC_GFOFFSET3D(u,0,0,-1) + KRANC_GFOFFSET3D(u,0,0,1)) + KRANC_GFOFFSET3D(u,0,0,-2) + KRANC_GFOFFSET3D(u,0,0,2))*INV(dxi*dyi) + ((6*KRANC_GFOFFSET3D(u,0,0,0) - 4*(KRANC_GFOFFSET3D(u,0,-1,0) + KRANC_GFOFFSET3D(u,0,1,0)) + KRANC_GFOFFSET3D(u,0,-2,0) + KRANC_GFOFFSET3D(u,0,2,0))*INV(dxi) + (6*KRANC_GFOFFSET3D(u,0,0,0) - 4*(KRANC_GFOFFSET3D(u,-1,0,0) + KRANC_GFOFFSET3D(u,1,0,0)) + KRANC_GFOFFSET3D(u,-2,0,0) + KRANC_GFOFFSET3D(u,2,0,0))*INV(dyi))*INV(dzi))); + ptrdiff_t const cdi CCTK_ATTRIBUTE_UNUSED = sizeof(CCTK_REAL); + return -1.*p1odxdydz*diss*((6.*KRANC_GFOFFSET3D(u,0,0,0) - 4.*(KRANC_GFOFFSET3D(u,0,0,-1) + KRANC_GFOFFSET3D(u,0,0,1)) + KRANC_GFOFFSET3D(u,0,0,-2) + KRANC_GFOFFSET3D(u,0,0,2))*INV(dxi*dyi) + ((6.*KRANC_GFOFFSET3D(u,0,0,0) - 4.*(KRANC_GFOFFSET3D(u,0,-1,0) + KRANC_GFOFFSET3D(u,0,1,0)) + KRANC_GFOFFSET3D(u,0,-2,0) + KRANC_GFOFFSET3D(u,0,2,0))*INV(dxi) + (6.*KRANC_GFOFFSET3D(u,0,0,0) - 4.*(KRANC_GFOFFSET3D(u,-1,0,0) + KRANC_GFOFFSET3D(u,1,0,0)) + KRANC_GFOFFSET3D(u,-2,0,0) + KRANC_GFOFFSET3D(u,2,0,0))*INV(dyi))*INV(dzi)); } #endif #ifndef KRANC_DIFF_FUNCTIONS -# define DissfdOrder4(u) (p1odxdydz*diss*((-20*KRANC_GFOFFSET3D(u,0,0,0) + 15*(KRANC_GFOFFSET3D(u,0,0,-1) + KRANC_GFOFFSET3D(u,0,0,1)) - 6*(KRANC_GFOFFSET3D(u,0,0,-2) + KRANC_GFOFFSET3D(u,0,0,2)) + KRANC_GFOFFSET3D(u,0,0,-3) + KRANC_GFOFFSET3D(u,0,0,3))*INV(dxi*dyi) + ((15*(KRANC_GFOFFSET3D(u,0,-1,0) + KRANC_GFOFFSET3D(u,0,1,0)) - 6*(KRANC_GFOFFSET3D(u,0,-2,0) + KRANC_GFOFFSET3D(u,0,2,0)) + KRANC_GFOFFSET3D(u,0,-3,0) + KRANC_GFOFFSET3D(u,0,3,0))*INV(dxi) + (15*(KRANC_GFOFFSET3D(u,-1,0,0) + KRANC_GFOFFSET3D(u,1,0,0)) - 6*(KRANC_GFOFFSET3D(u,-2,0,0) + KRANC_GFOFFSET3D(u,2,0,0)) + KRANC_GFOFFSET3D(u,-3,0,0) + KRANC_GFOFFSET3D(u,3,0,0))*INV(dyi) - 20*KRANC_GFOFFSET3D(u,0,0,0)*(INV(dxi) + INV(dyi)))*INV(dzi))) +# define DissfdOrder4(u) (p1odxdydz*diss*((-20.*KRANC_GFOFFSET3D(u,0,0,0) + 15.*(KRANC_GFOFFSET3D(u,0,0,-1) + KRANC_GFOFFSET3D(u,0,0,1)) - 6.*(KRANC_GFOFFSET3D(u,0,0,-2) + KRANC_GFOFFSET3D(u,0,0,2)) + KRANC_GFOFFSET3D(u,0,0,-3) + KRANC_GFOFFSET3D(u,0,0,3))*INV(dxi*dyi) + ((15.*(KRANC_GFOFFSET3D(u,0,-1,0) + KRANC_GFOFFSET3D(u,0,1,0)) - 6.*(KRANC_GFOFFSET3D(u,0,-2,0) + KRANC_GFOFFSET3D(u,0,2,0)) + KRANC_GFOFFSET3D(u,0,-3,0) + KRANC_GFOFFSET3D(u,0,3,0))*INV(dxi) + (15.*(KRANC_GFOFFSET3D(u,-1,0,0) + KRANC_GFOFFSET3D(u,1,0,0)) - 6.*(KRANC_GFOFFSET3D(u,-2,0,0) + KRANC_GFOFFSET3D(u,2,0,0)) + KRANC_GFOFFSET3D(u,-3,0,0) + KRANC_GFOFFSET3D(u,3,0,0))*INV(dyi) - 20.*KRANC_GFOFFSET3D(u,0,0,0)*(INV(dxi) + INV(dyi)))*INV(dzi))) #else # define DissfdOrder4(u) (DissfdOrder4_impl(u,p1odxdydz,cdj,cdk)) static CCTK_REAL DissfdOrder4_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1odxdydz, ptrdiff_t const cdj, ptrdiff_t const cdk) CCTK_ATTRIBUTE_NOINLINE CCTK_ATTRIBUTE_UNUSED; static CCTK_REAL DissfdOrder4_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1odxdydz, ptrdiff_t const cdj, ptrdiff_t const cdk) { - ptrdiff_t const cdi=sizeof(CCTK_REAL); - return p1odxdydz*diss*((-20*KRANC_GFOFFSET3D(u,0,0,0) + 15*(KRANC_GFOFFSET3D(u,0,0,-1) + KRANC_GFOFFSET3D(u,0,0,1)) - 6*(KRANC_GFOFFSET3D(u,0,0,-2) + KRANC_GFOFFSET3D(u,0,0,2)) + KRANC_GFOFFSET3D(u,0,0,-3) + KRANC_GFOFFSET3D(u,0,0,3))*INV(dxi*dyi) + ((15*(KRANC_GFOFFSET3D(u,0,-1,0) + KRANC_GFOFFSET3D(u,0,1,0)) - 6*(KRANC_GFOFFSET3D(u,0,-2,0) + KRANC_GFOFFSET3D(u,0,2,0)) + KRANC_GFOFFSET3D(u,0,-3,0) + KRANC_GFOFFSET3D(u,0,3,0))*INV(dxi) + (15*(KRANC_GFOFFSET3D(u,-1,0,0) + KRANC_GFOFFSET3D(u,1,0,0)) - 6*(KRANC_GFOFFSET3D(u,-2,0,0) + KRANC_GFOFFSET3D(u,2,0,0)) + KRANC_GFOFFSET3D(u,-3,0,0) + KRANC_GFOFFSET3D(u,3,0,0))*INV(dyi) - 20*KRANC_GFOFFSET3D(u,0,0,0)*(INV(dxi) + INV(dyi)))*INV(dzi)); + ptrdiff_t const cdi CCTK_ATTRIBUTE_UNUSED = sizeof(CCTK_REAL); + return p1odxdydz*diss*((-20.*KRANC_GFOFFSET3D(u,0,0,0) + 15.*(KRANC_GFOFFSET3D(u,0,0,-1) + KRANC_GFOFFSET3D(u,0,0,1)) - 6.*(KRANC_GFOFFSET3D(u,0,0,-2) + KRANC_GFOFFSET3D(u,0,0,2)) + KRANC_GFOFFSET3D(u,0,0,-3) + KRANC_GFOFFSET3D(u,0,0,3))*INV(dxi*dyi) + ((15.*(KRANC_GFOFFSET3D(u,0,-1,0) + KRANC_GFOFFSET3D(u,0,1,0)) - 6.*(KRANC_GFOFFSET3D(u,0,-2,0) + KRANC_GFOFFSET3D(u,0,2,0)) + KRANC_GFOFFSET3D(u,0,-3,0) + KRANC_GFOFFSET3D(u,0,3,0))*INV(dxi) + (15.*(KRANC_GFOFFSET3D(u,-1,0,0) + KRANC_GFOFFSET3D(u,1,0,0)) - 6.*(KRANC_GFOFFSET3D(u,-2,0,0) + KRANC_GFOFFSET3D(u,2,0,0)) + KRANC_GFOFFSET3D(u,-3,0,0) + KRANC_GFOFFSET3D(u,3,0,0))*INV(dyi) - 20.*KRANC_GFOFFSET3D(u,0,0,0)*(INV(dxi) + INV(dyi)))*INV(dzi)); } #endif #ifndef KRANC_DIFF_FUNCTIONS -# define PDzero1(u) ((-KRANC_GFOFFSET3D(u,-1,0,0) + KRANC_GFOFFSET3D(u,1,0,0))*p1o2dx) +# define PDzero1(u) ((-1.*KRANC_GFOFFSET3D(u,-1,0,0) + KRANC_GFOFFSET3D(u,1,0,0))*p1o2dx) #else # define PDzero1(u) (PDzero1_impl(u,p1o2dx,cdj,cdk)) static CCTK_REAL PDzero1_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1o2dx, ptrdiff_t const cdj, ptrdiff_t const cdk) CCTK_ATTRIBUTE_NOINLINE CCTK_ATTRIBUTE_UNUSED; static CCTK_REAL PDzero1_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1o2dx, ptrdiff_t const cdj, ptrdiff_t const cdk) { - ptrdiff_t const cdi=sizeof(CCTK_REAL); - return (-KRANC_GFOFFSET3D(u,-1,0,0) + KRANC_GFOFFSET3D(u,1,0,0))*p1o2dx; + ptrdiff_t const cdi CCTK_ATTRIBUTE_UNUSED = sizeof(CCTK_REAL); + return (-1.*KRANC_GFOFFSET3D(u,-1,0,0) + KRANC_GFOFFSET3D(u,1,0,0))*p1o2dx; } #endif #ifndef KRANC_DIFF_FUNCTIONS -# define PDzero2(u) ((-KRANC_GFOFFSET3D(u,0,-1,0) + KRANC_GFOFFSET3D(u,0,1,0))*p1o2dy) +# define PDzero2(u) ((-1.*KRANC_GFOFFSET3D(u,0,-1,0) + KRANC_GFOFFSET3D(u,0,1,0))*p1o2dy) #else # define PDzero2(u) (PDzero2_impl(u,p1o2dy,cdj,cdk)) static CCTK_REAL PDzero2_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1o2dy, ptrdiff_t const cdj, ptrdiff_t const cdk) CCTK_ATTRIBUTE_NOINLINE CCTK_ATTRIBUTE_UNUSED; static CCTK_REAL PDzero2_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1o2dy, ptrdiff_t const cdj, ptrdiff_t const cdk) { - ptrdiff_t const cdi=sizeof(CCTK_REAL); - return (-KRANC_GFOFFSET3D(u,0,-1,0) + KRANC_GFOFFSET3D(u,0,1,0))*p1o2dy; + ptrdiff_t const cdi CCTK_ATTRIBUTE_UNUSED = sizeof(CCTK_REAL); + return (-1.*KRANC_GFOFFSET3D(u,0,-1,0) + KRANC_GFOFFSET3D(u,0,1,0))*p1o2dy; } #endif #ifndef KRANC_DIFF_FUNCTIONS -# define PDzero3(u) ((-KRANC_GFOFFSET3D(u,0,0,-1) + KRANC_GFOFFSET3D(u,0,0,1))*p1o2dz) +# define PDzero3(u) ((-1.*KRANC_GFOFFSET3D(u,0,0,-1) + KRANC_GFOFFSET3D(u,0,0,1))*p1o2dz) #else # define PDzero3(u) (PDzero3_impl(u,p1o2dz,cdj,cdk)) static CCTK_REAL PDzero3_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1o2dz, ptrdiff_t const cdj, ptrdiff_t const cdk) CCTK_ATTRIBUTE_NOINLINE CCTK_ATTRIBUTE_UNUSED; static CCTK_REAL PDzero3_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1o2dz, ptrdiff_t const cdj, ptrdiff_t const cdk) { - ptrdiff_t const cdi=sizeof(CCTK_REAL); - return (-KRANC_GFOFFSET3D(u,0,0,-1) + KRANC_GFOFFSET3D(u,0,0,1))*p1o2dz; + ptrdiff_t const cdi CCTK_ATTRIBUTE_UNUSED = sizeof(CCTK_REAL); + return (-1.*KRANC_GFOFFSET3D(u,0,0,-1) + KRANC_GFOFFSET3D(u,0,0,1))*p1o2dz; } #endif #ifndef KRANC_DIFF_FUNCTIONS -# define PDzero11(u) ((-2*KRANC_GFOFFSET3D(u,0,0,0) + KRANC_GFOFFSET3D(u,-2,0,0) + KRANC_GFOFFSET3D(u,2,0,0))*p1o4dx2) +# define PDzero11(u) ((-2.*KRANC_GFOFFSET3D(u,0,0,0) + KRANC_GFOFFSET3D(u,-2,0,0) + KRANC_GFOFFSET3D(u,2,0,0))*p1o4dx2) #else # define PDzero11(u) (PDzero11_impl(u,p1o4dx2,cdj,cdk)) static CCTK_REAL PDzero11_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1o4dx2, ptrdiff_t const cdj, ptrdiff_t const cdk) CCTK_ATTRIBUTE_NOINLINE CCTK_ATTRIBUTE_UNUSED; static CCTK_REAL PDzero11_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1o4dx2, ptrdiff_t const cdj, ptrdiff_t const cdk) { - ptrdiff_t const cdi=sizeof(CCTK_REAL); - return (-2*KRANC_GFOFFSET3D(u,0,0,0) + KRANC_GFOFFSET3D(u,-2,0,0) + KRANC_GFOFFSET3D(u,2,0,0))*p1o4dx2; + ptrdiff_t const cdi CCTK_ATTRIBUTE_UNUSED = sizeof(CCTK_REAL); + return (-2.*KRANC_GFOFFSET3D(u,0,0,0) + KRANC_GFOFFSET3D(u,-2,0,0) + KRANC_GFOFFSET3D(u,2,0,0))*p1o4dx2; } #endif #ifndef KRANC_DIFF_FUNCTIONS -# define PDzero12(u) ((KRANC_GFOFFSET3D(u,-1,-1,0) - KRANC_GFOFFSET3D(u,-1,1,0) - KRANC_GFOFFSET3D(u,1,-1,0) + KRANC_GFOFFSET3D(u,1,1,0))*p1o4dxdy) +# define PDzero12(u) ((KRANC_GFOFFSET3D(u,-1,-1,0) - 1.*(KRANC_GFOFFSET3D(u,-1,1,0) + KRANC_GFOFFSET3D(u,1,-1,0)) + KRANC_GFOFFSET3D(u,1,1,0))*p1o4dxdy) #else # define PDzero12(u) (PDzero12_impl(u,p1o4dxdy,cdj,cdk)) static CCTK_REAL PDzero12_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1o4dxdy, ptrdiff_t const cdj, ptrdiff_t const cdk) CCTK_ATTRIBUTE_NOINLINE CCTK_ATTRIBUTE_UNUSED; static CCTK_REAL PDzero12_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1o4dxdy, ptrdiff_t const cdj, ptrdiff_t const cdk) { - ptrdiff_t const cdi=sizeof(CCTK_REAL); - return (KRANC_GFOFFSET3D(u,-1,-1,0) - KRANC_GFOFFSET3D(u,-1,1,0) - KRANC_GFOFFSET3D(u,1,-1,0) + KRANC_GFOFFSET3D(u,1,1,0))*p1o4dxdy; + ptrdiff_t const cdi CCTK_ATTRIBUTE_UNUSED = sizeof(CCTK_REAL); + return (KRANC_GFOFFSET3D(u,-1,-1,0) - 1.*(KRANC_GFOFFSET3D(u,-1,1,0) + KRANC_GFOFFSET3D(u,1,-1,0)) + KRANC_GFOFFSET3D(u,1,1,0))*p1o4dxdy; } #endif #ifndef KRANC_DIFF_FUNCTIONS -# define PDzero13(u) ((KRANC_GFOFFSET3D(u,-1,0,-1) - KRANC_GFOFFSET3D(u,-1,0,1) - KRANC_GFOFFSET3D(u,1,0,-1) + KRANC_GFOFFSET3D(u,1,0,1))*p1o4dxdz) +# define PDzero13(u) ((KRANC_GFOFFSET3D(u,-1,0,-1) - 1.*(KRANC_GFOFFSET3D(u,-1,0,1) + KRANC_GFOFFSET3D(u,1,0,-1)) + KRANC_GFOFFSET3D(u,1,0,1))*p1o4dxdz) #else # define PDzero13(u) (PDzero13_impl(u,p1o4dxdz,cdj,cdk)) static CCTK_REAL PDzero13_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1o4dxdz, ptrdiff_t const cdj, ptrdiff_t const cdk) CCTK_ATTRIBUTE_NOINLINE CCTK_ATTRIBUTE_UNUSED; static CCTK_REAL PDzero13_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1o4dxdz, ptrdiff_t const cdj, ptrdiff_t const cdk) { - ptrdiff_t const cdi=sizeof(CCTK_REAL); - return (KRANC_GFOFFSET3D(u,-1,0,-1) - KRANC_GFOFFSET3D(u,-1,0,1) - KRANC_GFOFFSET3D(u,1,0,-1) + KRANC_GFOFFSET3D(u,1,0,1))*p1o4dxdz; + ptrdiff_t const cdi CCTK_ATTRIBUTE_UNUSED = sizeof(CCTK_REAL); + return (KRANC_GFOFFSET3D(u,-1,0,-1) - 1.*(KRANC_GFOFFSET3D(u,-1,0,1) + KRANC_GFOFFSET3D(u,1,0,-1)) + KRANC_GFOFFSET3D(u,1,0,1))*p1o4dxdz; } #endif #ifndef KRANC_DIFF_FUNCTIONS -# define PDzero21(u) ((KRANC_GFOFFSET3D(u,-1,-1,0) - KRANC_GFOFFSET3D(u,-1,1,0) - KRANC_GFOFFSET3D(u,1,-1,0) + KRANC_GFOFFSET3D(u,1,1,0))*p1o4dxdy) +# define PDzero21(u) ((KRANC_GFOFFSET3D(u,-1,-1,0) - 1.*(KRANC_GFOFFSET3D(u,-1,1,0) + KRANC_GFOFFSET3D(u,1,-1,0)) + KRANC_GFOFFSET3D(u,1,1,0))*p1o4dxdy) #else # define PDzero21(u) (PDzero21_impl(u,p1o4dxdy,cdj,cdk)) static CCTK_REAL PDzero21_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1o4dxdy, ptrdiff_t const cdj, ptrdiff_t const cdk) CCTK_ATTRIBUTE_NOINLINE CCTK_ATTRIBUTE_UNUSED; static CCTK_REAL PDzero21_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1o4dxdy, ptrdiff_t const cdj, ptrdiff_t const cdk) { - ptrdiff_t const cdi=sizeof(CCTK_REAL); - return (KRANC_GFOFFSET3D(u,-1,-1,0) - KRANC_GFOFFSET3D(u,-1,1,0) - KRANC_GFOFFSET3D(u,1,-1,0) + KRANC_GFOFFSET3D(u,1,1,0))*p1o4dxdy; + ptrdiff_t const cdi CCTK_ATTRIBUTE_UNUSED = sizeof(CCTK_REAL); + return (KRANC_GFOFFSET3D(u,-1,-1,0) - 1.*(KRANC_GFOFFSET3D(u,-1,1,0) + KRANC_GFOFFSET3D(u,1,-1,0)) + KRANC_GFOFFSET3D(u,1,1,0))*p1o4dxdy; } #endif #ifndef KRANC_DIFF_FUNCTIONS -# define PDzero22(u) ((-2*KRANC_GFOFFSET3D(u,0,0,0) + KRANC_GFOFFSET3D(u,0,-2,0) + KRANC_GFOFFSET3D(u,0,2,0))*p1o4dy2) +# define PDzero22(u) ((-2.*KRANC_GFOFFSET3D(u,0,0,0) + KRANC_GFOFFSET3D(u,0,-2,0) + KRANC_GFOFFSET3D(u,0,2,0))*p1o4dy2) #else # define PDzero22(u) (PDzero22_impl(u,p1o4dy2,cdj,cdk)) static CCTK_REAL PDzero22_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1o4dy2, ptrdiff_t const cdj, ptrdiff_t const cdk) CCTK_ATTRIBUTE_NOINLINE CCTK_ATTRIBUTE_UNUSED; static CCTK_REAL PDzero22_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1o4dy2, ptrdiff_t const cdj, ptrdiff_t const cdk) { - ptrdiff_t const cdi=sizeof(CCTK_REAL); - return (-2*KRANC_GFOFFSET3D(u,0,0,0) + KRANC_GFOFFSET3D(u,0,-2,0) + KRANC_GFOFFSET3D(u,0,2,0))*p1o4dy2; + ptrdiff_t const cdi CCTK_ATTRIBUTE_UNUSED = sizeof(CCTK_REAL); + return (-2.*KRANC_GFOFFSET3D(u,0,0,0) + KRANC_GFOFFSET3D(u,0,-2,0) + KRANC_GFOFFSET3D(u,0,2,0))*p1o4dy2; } #endif #ifndef KRANC_DIFF_FUNCTIONS -# define PDzero23(u) ((KRANC_GFOFFSET3D(u,0,-1,-1) - KRANC_GFOFFSET3D(u,0,-1,1) - KRANC_GFOFFSET3D(u,0,1,-1) + KRANC_GFOFFSET3D(u,0,1,1))*p1o4dydz) +# define PDzero23(u) ((KRANC_GFOFFSET3D(u,0,-1,-1) - 1.*(KRANC_GFOFFSET3D(u,0,-1,1) + KRANC_GFOFFSET3D(u,0,1,-1)) + KRANC_GFOFFSET3D(u,0,1,1))*p1o4dydz) #else # define PDzero23(u) (PDzero23_impl(u,p1o4dydz,cdj,cdk)) static CCTK_REAL PDzero23_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1o4dydz, ptrdiff_t const cdj, ptrdiff_t const cdk) CCTK_ATTRIBUTE_NOINLINE CCTK_ATTRIBUTE_UNUSED; static CCTK_REAL PDzero23_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1o4dydz, ptrdiff_t const cdj, ptrdiff_t const cdk) { - ptrdiff_t const cdi=sizeof(CCTK_REAL); - return (KRANC_GFOFFSET3D(u,0,-1,-1) - KRANC_GFOFFSET3D(u,0,-1,1) - KRANC_GFOFFSET3D(u,0,1,-1) + KRANC_GFOFFSET3D(u,0,1,1))*p1o4dydz; + ptrdiff_t const cdi CCTK_ATTRIBUTE_UNUSED = sizeof(CCTK_REAL); + return (KRANC_GFOFFSET3D(u,0,-1,-1) - 1.*(KRANC_GFOFFSET3D(u,0,-1,1) + KRANC_GFOFFSET3D(u,0,1,-1)) + KRANC_GFOFFSET3D(u,0,1,1))*p1o4dydz; } #endif #ifndef KRANC_DIFF_FUNCTIONS -# define PDzero31(u) ((KRANC_GFOFFSET3D(u,-1,0,-1) - KRANC_GFOFFSET3D(u,-1,0,1) - KRANC_GFOFFSET3D(u,1,0,-1) + KRANC_GFOFFSET3D(u,1,0,1))*p1o4dxdz) +# define PDzero31(u) ((KRANC_GFOFFSET3D(u,-1,0,-1) - 1.*(KRANC_GFOFFSET3D(u,-1,0,1) + KRANC_GFOFFSET3D(u,1,0,-1)) + KRANC_GFOFFSET3D(u,1,0,1))*p1o4dxdz) #else # define PDzero31(u) (PDzero31_impl(u,p1o4dxdz,cdj,cdk)) static CCTK_REAL PDzero31_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1o4dxdz, ptrdiff_t const cdj, ptrdiff_t const cdk) CCTK_ATTRIBUTE_NOINLINE CCTK_ATTRIBUTE_UNUSED; static CCTK_REAL PDzero31_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1o4dxdz, ptrdiff_t const cdj, ptrdiff_t const cdk) { - ptrdiff_t const cdi=sizeof(CCTK_REAL); - return (KRANC_GFOFFSET3D(u,-1,0,-1) - KRANC_GFOFFSET3D(u,-1,0,1) - KRANC_GFOFFSET3D(u,1,0,-1) + KRANC_GFOFFSET3D(u,1,0,1))*p1o4dxdz; + ptrdiff_t const cdi CCTK_ATTRIBUTE_UNUSED = sizeof(CCTK_REAL); + return (KRANC_GFOFFSET3D(u,-1,0,-1) - 1.*(KRANC_GFOFFSET3D(u,-1,0,1) + KRANC_GFOFFSET3D(u,1,0,-1)) + KRANC_GFOFFSET3D(u,1,0,1))*p1o4dxdz; } #endif #ifndef KRANC_DIFF_FUNCTIONS -# define PDzero32(u) ((KRANC_GFOFFSET3D(u,0,-1,-1) - KRANC_GFOFFSET3D(u,0,-1,1) - KRANC_GFOFFSET3D(u,0,1,-1) + KRANC_GFOFFSET3D(u,0,1,1))*p1o4dydz) +# define PDzero32(u) ((KRANC_GFOFFSET3D(u,0,-1,-1) - 1.*(KRANC_GFOFFSET3D(u,0,-1,1) + KRANC_GFOFFSET3D(u,0,1,-1)) + KRANC_GFOFFSET3D(u,0,1,1))*p1o4dydz) #else # define PDzero32(u) (PDzero32_impl(u,p1o4dydz,cdj,cdk)) static CCTK_REAL PDzero32_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1o4dydz, ptrdiff_t const cdj, ptrdiff_t const cdk) CCTK_ATTRIBUTE_NOINLINE CCTK_ATTRIBUTE_UNUSED; static CCTK_REAL PDzero32_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1o4dydz, ptrdiff_t const cdj, ptrdiff_t const cdk) { - ptrdiff_t const cdi=sizeof(CCTK_REAL); - return (KRANC_GFOFFSET3D(u,0,-1,-1) - KRANC_GFOFFSET3D(u,0,-1,1) - KRANC_GFOFFSET3D(u,0,1,-1) + KRANC_GFOFFSET3D(u,0,1,1))*p1o4dydz; + ptrdiff_t const cdi CCTK_ATTRIBUTE_UNUSED = sizeof(CCTK_REAL); + return (KRANC_GFOFFSET3D(u,0,-1,-1) - 1.*(KRANC_GFOFFSET3D(u,0,-1,1) + KRANC_GFOFFSET3D(u,0,1,-1)) + KRANC_GFOFFSET3D(u,0,1,1))*p1o4dydz; } #endif #ifndef KRANC_DIFF_FUNCTIONS -# define PDzero33(u) ((-2*KRANC_GFOFFSET3D(u,0,0,0) + KRANC_GFOFFSET3D(u,0,0,-2) + KRANC_GFOFFSET3D(u,0,0,2))*p1o4dz2) +# define PDzero33(u) ((-2.*KRANC_GFOFFSET3D(u,0,0,0) + KRANC_GFOFFSET3D(u,0,0,-2) + KRANC_GFOFFSET3D(u,0,0,2))*p1o4dz2) #else # define PDzero33(u) (PDzero33_impl(u,p1o4dz2,cdj,cdk)) static CCTK_REAL PDzero33_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1o4dz2, ptrdiff_t const cdj, ptrdiff_t const cdk) CCTK_ATTRIBUTE_NOINLINE CCTK_ATTRIBUTE_UNUSED; static CCTK_REAL PDzero33_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1o4dz2, ptrdiff_t const cdj, ptrdiff_t const cdk) { - ptrdiff_t const cdi=sizeof(CCTK_REAL); - return (-2*KRANC_GFOFFSET3D(u,0,0,0) + KRANC_GFOFFSET3D(u,0,0,-2) + KRANC_GFOFFSET3D(u,0,0,2))*p1o4dz2; + ptrdiff_t const cdi CCTK_ATTRIBUTE_UNUSED = sizeof(CCTK_REAL); + return (-2.*KRANC_GFOFFSET3D(u,0,0,0) + KRANC_GFOFFSET3D(u,0,0,-2) + KRANC_GFOFFSET3D(u,0,0,2))*p1o4dz2; } #endif #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 diff --git a/Examples/Wave/src/RegisterMoL.cc b/Examples/Wave/src/RegisterMoL.cc index 623b8e3..4653774 100644 --- a/Examples/Wave/src/RegisterMoL.cc +++ b/Examples/Wave/src/RegisterMoL.cc @@ -9,7 +9,7 @@ extern "C" void Wave_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("Wave::phi"), CCTK_VarIndex("Wave::phirhs")); diff --git a/Examples/Wave/src/Startup.cc b/Examples/Wave/src/Startup.cc index 9281f39..70cc6b3 100644 --- a/Examples/Wave/src/Startup.cc +++ b/Examples/Wave/src/Startup.cc @@ -4,7 +4,7 @@ extern "C" int Wave_Startup(void) { - const char * banner = "Wave"; + const char * banner CCTK_ATTRIBUTE_UNUSED = "Wave"; CCTK_RegisterBanner(banner); return 0; } diff --git a/Examples/Wave/src/wave_boundary.cc b/Examples/Wave/src/wave_boundary.cc index 3ff86ec..794a9a1 100644 --- a/Examples/Wave/src/wave_boundary.cc +++ b/Examples/Wave/src/wave_boundary.cc @@ -27,7 +27,7 @@ extern "C" void wave_boundary_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 */, "Wave::evolvedrhs","flat"); if (ierr < 0) CCTK_WARN(1, "Failed to register flat BC for Wave::evolvedrhs."); @@ -43,58 +43,58 @@ static void wave_boundary_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 p1o12dx = 0.0833333333333333333333333333333*INV(dx); - CCTK_REAL const p1o12dy = 0.0833333333333333333333333333333*INV(dy); - CCTK_REAL const p1o12dz = 0.0833333333333333333333333333333*INV(dz); - CCTK_REAL const p1o144dxdy = 0.00694444444444444444444444444444*INV(dx*dy); - CCTK_REAL const p1o144dxdz = 0.00694444444444444444444444444444*INV(dx*dz); - CCTK_REAL const p1o144dydz = 0.00694444444444444444444444444444*INV(dy*dz); - CCTK_REAL const p1o2dx = 0.5*INV(dx); - CCTK_REAL const p1o2dy = 0.5*INV(dy); - CCTK_REAL const p1o2dz = 0.5*INV(dz); - CCTK_REAL const p1o4dx2 = 0.25*INV(SQR(dx)); - CCTK_REAL const p1o4dxdy = 0.25*INV(dx*dy); - CCTK_REAL const p1o4dxdz = 0.25*INV(dx*dz); - CCTK_REAL const p1o4dy2 = 0.25*INV(SQR(dy)); - CCTK_REAL const p1o4dydz = 0.25*INV(dy*dz); - CCTK_REAL const p1o4dz2 = 0.25*INV(SQR(dz)); - CCTK_REAL const p1odx = INV(dx); - CCTK_REAL const p1odx2 = INV(SQR(dx)); - CCTK_REAL const p1odxdydz = INV(dx*dy*dz); - CCTK_REAL const p1ody = INV(dy); - CCTK_REAL const p1ody2 = INV(SQR(dy)); - CCTK_REAL const p1odz = INV(dz); - CCTK_REAL const p1odz2 = INV(SQR(dz)); - CCTK_REAL const pm1o12dx2 = -0.0833333333333333333333333333333*INV(SQR(dx)); - CCTK_REAL const pm1o12dy2 = -0.0833333333333333333333333333333*INV(SQR(dy)); - CCTK_REAL const pm1o12dz2 = -0.0833333333333333333333333333333*INV(SQR(dz)); - CCTK_REAL const pm1o2dx = -0.5*INV(dx); - CCTK_REAL const pm1o2dy = -0.5*INV(dy); - CCTK_REAL const pm1o2dz = -0.5*INV(dz); + CCTK_REAL /*const*/ p1o12dx CCTK_ATTRIBUTE_UNUSED = 0.0833333333333333333333333333333*INV(dx); + CCTK_REAL /*const*/ p1o12dy CCTK_ATTRIBUTE_UNUSED = 0.0833333333333333333333333333333*INV(dy); + CCTK_REAL /*const*/ p1o12dz CCTK_ATTRIBUTE_UNUSED = 0.0833333333333333333333333333333*INV(dz); + CCTK_REAL /*const*/ p1o144dxdy CCTK_ATTRIBUTE_UNUSED = 0.00694444444444444444444444444444*INV(dx*dy); + CCTK_REAL /*const*/ p1o144dxdz CCTK_ATTRIBUTE_UNUSED = 0.00694444444444444444444444444444*INV(dx*dz); + CCTK_REAL /*const*/ p1o144dydz CCTK_ATTRIBUTE_UNUSED = 0.00694444444444444444444444444444*INV(dy*dz); + CCTK_REAL /*const*/ p1o2dx CCTK_ATTRIBUTE_UNUSED = 0.5*INV(dx); + CCTK_REAL /*const*/ p1o2dy CCTK_ATTRIBUTE_UNUSED = 0.5*INV(dy); + CCTK_REAL /*const*/ p1o2dz CCTK_ATTRIBUTE_UNUSED = 0.5*INV(dz); + CCTK_REAL /*const*/ p1o4dx2 CCTK_ATTRIBUTE_UNUSED = 0.25*INV(SQR(dx)); + CCTK_REAL /*const*/ p1o4dxdy CCTK_ATTRIBUTE_UNUSED = 0.25*INV(dx*dy); + CCTK_REAL /*const*/ p1o4dxdz CCTK_ATTRIBUTE_UNUSED = 0.25*INV(dx*dz); + CCTK_REAL /*const*/ p1o4dy2 CCTK_ATTRIBUTE_UNUSED = 0.25*INV(SQR(dy)); + CCTK_REAL /*const*/ p1o4dydz CCTK_ATTRIBUTE_UNUSED = 0.25*INV(dy*dz); + CCTK_REAL /*const*/ p1o4dz2 CCTK_ATTRIBUTE_UNUSED = 0.25*INV(SQR(dz)); + CCTK_REAL /*const*/ p1odx CCTK_ATTRIBUTE_UNUSED = INV(dx); + CCTK_REAL /*const*/ p1odx2 CCTK_ATTRIBUTE_UNUSED = INV(SQR(dx)); + CCTK_REAL /*const*/ p1odxdydz CCTK_ATTRIBUTE_UNUSED = INV(dx*dy*dz); + CCTK_REAL /*const*/ p1ody CCTK_ATTRIBUTE_UNUSED = INV(dy); + CCTK_REAL /*const*/ p1ody2 CCTK_ATTRIBUTE_UNUSED = INV(SQR(dy)); + CCTK_REAL /*const*/ p1odz CCTK_ATTRIBUTE_UNUSED = INV(dz); + CCTK_REAL /*const*/ p1odz2 CCTK_ATTRIBUTE_UNUSED = INV(SQR(dz)); + CCTK_REAL /*const*/ pm1o12dx2 CCTK_ATTRIBUTE_UNUSED = -0.0833333333333333333333333333333*INV(SQR(dx)); + CCTK_REAL /*const*/ pm1o12dy2 CCTK_ATTRIBUTE_UNUSED = -0.0833333333333333333333333333333*INV(SQR(dy)); + CCTK_REAL /*const*/ pm1o12dz2 CCTK_ATTRIBUTE_UNUSED = -0.0833333333333333333333333333333*INV(SQR(dz)); + CCTK_REAL /*const*/ pm1o2dx CCTK_ATTRIBUTE_UNUSED = -0.5*INV(dx); + CCTK_REAL /*const*/ pm1o2dy CCTK_ATTRIBUTE_UNUSED = -0.5*INV(dy); + CCTK_REAL /*const*/ pm1o2dz CCTK_ATTRIBUTE_UNUSED = -0.5*INV(dz); /* Assign local copies of arrays functions */ @@ -110,16 +110,16 @@ static void wave_boundary_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 phiL = phi[index]; - CCTK_REAL piL = pi[index]; - CCTK_REAL rL = r[index]; - CCTK_REAL xL = x[index]; - CCTK_REAL yL = y[index]; - CCTK_REAL zL = z[index]; + CCTK_REAL phiL CCTK_ATTRIBUTE_UNUSED = phi[index]; + CCTK_REAL piL CCTK_ATTRIBUTE_UNUSED = pi[index]; + CCTK_REAL rL CCTK_ATTRIBUTE_UNUSED = r[index]; + CCTK_REAL xL CCTK_ATTRIBUTE_UNUSED = x[index]; + CCTK_REAL yL CCTK_ATTRIBUTE_UNUSED = y[index]; + CCTK_REAL zL CCTK_ATTRIBUTE_UNUSED = z[index]; /* Include user supplied include files */ @@ -136,25 +136,26 @@ static void wave_boundary_Body(cGH const * restrict const cctkGH, int const dir, } /* Calculate temporaries and grid functions */ - CCTK_REAL norm1 = -(xL*INV(rL)); + CCTK_REAL CCTK_ATTRIBUTE_UNUSED norm1 = -1.*xL*INV(rL); - CCTK_REAL norm2 = -(yL*INV(rL)); + CCTK_REAL CCTK_ATTRIBUTE_UNUSED norm2 = -1.*yL*INV(rL); - CCTK_REAL norm3 = -(zL*INV(rL)); + CCTK_REAL CCTK_ATTRIBUTE_UNUSED norm3 = -1.*zL*INV(rL); - ptrdiff_t dir1 = isgn(norm1); + ptrdiff_t CCTK_ATTRIBUTE_UNUSED dir1 = isgn(norm1); - ptrdiff_t dir2 = isgn(norm2); + ptrdiff_t CCTK_ATTRIBUTE_UNUSED dir2 = isgn(norm2); - ptrdiff_t dir3 = isgn(norm3); + ptrdiff_t CCTK_ATTRIBUTE_UNUSED dir3 = isgn(norm3); - CCTK_REAL phirhsL = PDonesided2nd1(&phi[index])*norm1 + + CCTK_REAL CCTK_ATTRIBUTE_UNUSED phirhsL = + PDonesided2nd1(&phi[index])*norm1 + PDonesided2nd2(&phi[index])*norm2 + - PDonesided2nd3(&phi[index])*norm3 - phiL*INV(rL); + PDonesided2nd3(&phi[index])*norm3 - 1.*phiL*INV(rL); - CCTK_REAL pirhsL = PDonesided2nd1(&pi[index])*norm1 + - PDonesided2nd2(&pi[index])*norm2 + PDonesided2nd3(&pi[index])*norm3 - - piL*INV(rL); + CCTK_REAL CCTK_ATTRIBUTE_UNUSED pirhsL = + PDonesided2nd1(&pi[index])*norm1 + PDonesided2nd2(&pi[index])*norm2 + + PDonesided2nd3(&pi[index])*norm3 - 1.*piL*INV(rL); /* Copy local copies back to grid functions */ phirhs[index] = phirhsL; diff --git a/Examples/Wave/src/wave_calc_errors.cc b/Examples/Wave/src/wave_calc_errors.cc index 1547953..66eb7ab 100644 --- a/Examples/Wave/src/wave_calc_errors.cc +++ b/Examples/Wave/src/wave_calc_errors.cc @@ -31,58 +31,58 @@ static void wave_calc_errors_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 p1o12dx = 0.0833333333333333333333333333333*INV(dx); - CCTK_REAL const p1o12dy = 0.0833333333333333333333333333333*INV(dy); - CCTK_REAL const p1o12dz = 0.0833333333333333333333333333333*INV(dz); - CCTK_REAL const p1o144dxdy = 0.00694444444444444444444444444444*INV(dx*dy); - CCTK_REAL const p1o144dxdz = 0.00694444444444444444444444444444*INV(dx*dz); - CCTK_REAL const p1o144dydz = 0.00694444444444444444444444444444*INV(dy*dz); - CCTK_REAL const p1o2dx = 0.5*INV(dx); - CCTK_REAL const p1o2dy = 0.5*INV(dy); - CCTK_REAL const p1o2dz = 0.5*INV(dz); - CCTK_REAL const p1o4dx2 = 0.25*INV(SQR(dx)); - CCTK_REAL const p1o4dxdy = 0.25*INV(dx*dy); - CCTK_REAL const p1o4dxdz = 0.25*INV(dx*dz); - CCTK_REAL const p1o4dy2 = 0.25*INV(SQR(dy)); - CCTK_REAL const p1o4dydz = 0.25*INV(dy*dz); - CCTK_REAL const p1o4dz2 = 0.25*INV(SQR(dz)); - CCTK_REAL const p1odx = INV(dx); - CCTK_REAL const p1odx2 = INV(SQR(dx)); - CCTK_REAL const p1odxdydz = INV(dx*dy*dz); - CCTK_REAL const p1ody = INV(dy); - CCTK_REAL const p1ody2 = INV(SQR(dy)); - CCTK_REAL const p1odz = INV(dz); - CCTK_REAL const p1odz2 = INV(SQR(dz)); - CCTK_REAL const pm1o12dx2 = -0.0833333333333333333333333333333*INV(SQR(dx)); - CCTK_REAL const pm1o12dy2 = -0.0833333333333333333333333333333*INV(SQR(dy)); - CCTK_REAL const pm1o12dz2 = -0.0833333333333333333333333333333*INV(SQR(dz)); - CCTK_REAL const pm1o2dx = -0.5*INV(dx); - CCTK_REAL const pm1o2dy = -0.5*INV(dy); - CCTK_REAL const pm1o2dz = -0.5*INV(dz); + CCTK_REAL /*const*/ p1o12dx CCTK_ATTRIBUTE_UNUSED = 0.0833333333333333333333333333333*INV(dx); + CCTK_REAL /*const*/ p1o12dy CCTK_ATTRIBUTE_UNUSED = 0.0833333333333333333333333333333*INV(dy); + CCTK_REAL /*const*/ p1o12dz CCTK_ATTRIBUTE_UNUSED = 0.0833333333333333333333333333333*INV(dz); + CCTK_REAL /*const*/ p1o144dxdy CCTK_ATTRIBUTE_UNUSED = 0.00694444444444444444444444444444*INV(dx*dy); + CCTK_REAL /*const*/ p1o144dxdz CCTK_ATTRIBUTE_UNUSED = 0.00694444444444444444444444444444*INV(dx*dz); + CCTK_REAL /*const*/ p1o144dydz CCTK_ATTRIBUTE_UNUSED = 0.00694444444444444444444444444444*INV(dy*dz); + CCTK_REAL /*const*/ p1o2dx CCTK_ATTRIBUTE_UNUSED = 0.5*INV(dx); + CCTK_REAL /*const*/ p1o2dy CCTK_ATTRIBUTE_UNUSED = 0.5*INV(dy); + CCTK_REAL /*const*/ p1o2dz CCTK_ATTRIBUTE_UNUSED = 0.5*INV(dz); + CCTK_REAL /*const*/ p1o4dx2 CCTK_ATTRIBUTE_UNUSED = 0.25*INV(SQR(dx)); + CCTK_REAL /*const*/ p1o4dxdy CCTK_ATTRIBUTE_UNUSED = 0.25*INV(dx*dy); + CCTK_REAL /*const*/ p1o4dxdz CCTK_ATTRIBUTE_UNUSED = 0.25*INV(dx*dz); + CCTK_REAL /*const*/ p1o4dy2 CCTK_ATTRIBUTE_UNUSED = 0.25*INV(SQR(dy)); + CCTK_REAL /*const*/ p1o4dydz CCTK_ATTRIBUTE_UNUSED = 0.25*INV(dy*dz); + CCTK_REAL /*const*/ p1o4dz2 CCTK_ATTRIBUTE_UNUSED = 0.25*INV(SQR(dz)); + CCTK_REAL /*const*/ p1odx CCTK_ATTRIBUTE_UNUSED = INV(dx); + CCTK_REAL /*const*/ p1odx2 CCTK_ATTRIBUTE_UNUSED = INV(SQR(dx)); + CCTK_REAL /*const*/ p1odxdydz CCTK_ATTRIBUTE_UNUSED = INV(dx*dy*dz); + CCTK_REAL /*const*/ p1ody CCTK_ATTRIBUTE_UNUSED = INV(dy); + CCTK_REAL /*const*/ p1ody2 CCTK_ATTRIBUTE_UNUSED = INV(SQR(dy)); + CCTK_REAL /*const*/ p1odz CCTK_ATTRIBUTE_UNUSED = INV(dz); + CCTK_REAL /*const*/ p1odz2 CCTK_ATTRIBUTE_UNUSED = INV(SQR(dz)); + CCTK_REAL /*const*/ pm1o12dx2 CCTK_ATTRIBUTE_UNUSED = -0.0833333333333333333333333333333*INV(SQR(dx)); + CCTK_REAL /*const*/ pm1o12dy2 CCTK_ATTRIBUTE_UNUSED = -0.0833333333333333333333333333333*INV(SQR(dy)); + CCTK_REAL /*const*/ pm1o12dz2 CCTK_ATTRIBUTE_UNUSED = -0.0833333333333333333333333333333*INV(SQR(dz)); + CCTK_REAL /*const*/ pm1o2dx CCTK_ATTRIBUTE_UNUSED = -0.5*INV(dx); + CCTK_REAL /*const*/ pm1o2dy CCTK_ATTRIBUTE_UNUSED = -0.5*INV(dy); + CCTK_REAL /*const*/ pm1o2dz CCTK_ATTRIBUTE_UNUSED = -0.5*INV(dz); /* Assign local copies of arrays functions */ @@ -98,14 +98,14 @@ static void wave_calc_errors_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 */ - CCTK_REAL phiL = phi[index]; - CCTK_REAL phiExactL = phiExact[index]; - CCTK_REAL piL = pi[index]; - CCTK_REAL piExactL = piExact[index]; + CCTK_REAL phiL CCTK_ATTRIBUTE_UNUSED = phi[index]; + CCTK_REAL phiExactL CCTK_ATTRIBUTE_UNUSED = phiExact[index]; + CCTK_REAL piL CCTK_ATTRIBUTE_UNUSED = pi[index]; + CCTK_REAL piExactL CCTK_ATTRIBUTE_UNUSED = piExact[index]; /* Include user supplied include files */ @@ -122,9 +122,9 @@ static void wave_calc_errors_Body(cGH const * restrict const cctkGH, int const d } /* Calculate temporaries and grid functions */ - CCTK_REAL phiErrorL = -phiExactL + phiL; + CCTK_REAL CCTK_ATTRIBUTE_UNUSED phiErrorL = -1.*phiExactL + phiL; - CCTK_REAL piErrorL = -piExactL + piL; + CCTK_REAL CCTK_ATTRIBUTE_UNUSED piErrorL = -1.*piExactL + piL; /* Copy local copies back to grid functions */ phiError[index] = phiErrorL; diff --git a/Examples/Wave/src/wave_calc_norm.cc b/Examples/Wave/src/wave_calc_norm.cc index 58d8de3..d41ca04 100644 --- a/Examples/Wave/src/wave_calc_norm.cc +++ b/Examples/Wave/src/wave_calc_norm.cc @@ -27,7 +27,7 @@ extern "C" void wave_calc_norm_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 */, "Wave::norms","flat"); if (ierr < 0) CCTK_WARN(1, "Failed to register flat BC for Wave::norms."); @@ -43,58 +43,58 @@ static void wave_calc_norm_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 p1o12dx = 0.0833333333333333333333333333333*INV(dx); - CCTK_REAL const p1o12dy = 0.0833333333333333333333333333333*INV(dy); - CCTK_REAL const p1o12dz = 0.0833333333333333333333333333333*INV(dz); - CCTK_REAL const p1o144dxdy = 0.00694444444444444444444444444444*INV(dx*dy); - CCTK_REAL const p1o144dxdz = 0.00694444444444444444444444444444*INV(dx*dz); - CCTK_REAL const p1o144dydz = 0.00694444444444444444444444444444*INV(dy*dz); - CCTK_REAL const p1o2dx = 0.5*INV(dx); - CCTK_REAL const p1o2dy = 0.5*INV(dy); - CCTK_REAL const p1o2dz = 0.5*INV(dz); - CCTK_REAL const p1o4dx2 = 0.25*INV(SQR(dx)); - CCTK_REAL const p1o4dxdy = 0.25*INV(dx*dy); - CCTK_REAL const p1o4dxdz = 0.25*INV(dx*dz); - CCTK_REAL const p1o4dy2 = 0.25*INV(SQR(dy)); - CCTK_REAL const p1o4dydz = 0.25*INV(dy*dz); - CCTK_REAL const p1o4dz2 = 0.25*INV(SQR(dz)); - CCTK_REAL const p1odx = INV(dx); - CCTK_REAL const p1odx2 = INV(SQR(dx)); - CCTK_REAL const p1odxdydz = INV(dx*dy*dz); - CCTK_REAL const p1ody = INV(dy); - CCTK_REAL const p1ody2 = INV(SQR(dy)); - CCTK_REAL const p1odz = INV(dz); - CCTK_REAL const p1odz2 = INV(SQR(dz)); - CCTK_REAL const pm1o12dx2 = -0.0833333333333333333333333333333*INV(SQR(dx)); - CCTK_REAL const pm1o12dy2 = -0.0833333333333333333333333333333*INV(SQR(dy)); - CCTK_REAL const pm1o12dz2 = -0.0833333333333333333333333333333*INV(SQR(dz)); - CCTK_REAL const pm1o2dx = -0.5*INV(dx); - CCTK_REAL const pm1o2dy = -0.5*INV(dy); - CCTK_REAL const pm1o2dz = -0.5*INV(dz); + CCTK_REAL /*const*/ p1o12dx CCTK_ATTRIBUTE_UNUSED = 0.0833333333333333333333333333333*INV(dx); + CCTK_REAL /*const*/ p1o12dy CCTK_ATTRIBUTE_UNUSED = 0.0833333333333333333333333333333*INV(dy); + CCTK_REAL /*const*/ p1o12dz CCTK_ATTRIBUTE_UNUSED = 0.0833333333333333333333333333333*INV(dz); + CCTK_REAL /*const*/ p1o144dxdy CCTK_ATTRIBUTE_UNUSED = 0.00694444444444444444444444444444*INV(dx*dy); + CCTK_REAL /*const*/ p1o144dxdz CCTK_ATTRIBUTE_UNUSED = 0.00694444444444444444444444444444*INV(dx*dz); + CCTK_REAL /*const*/ p1o144dydz CCTK_ATTRIBUTE_UNUSED = 0.00694444444444444444444444444444*INV(dy*dz); + CCTK_REAL /*const*/ p1o2dx CCTK_ATTRIBUTE_UNUSED = 0.5*INV(dx); + CCTK_REAL /*const*/ p1o2dy CCTK_ATTRIBUTE_UNUSED = 0.5*INV(dy); + CCTK_REAL /*const*/ p1o2dz CCTK_ATTRIBUTE_UNUSED = 0.5*INV(dz); + CCTK_REAL /*const*/ p1o4dx2 CCTK_ATTRIBUTE_UNUSED = 0.25*INV(SQR(dx)); + CCTK_REAL /*const*/ p1o4dxdy CCTK_ATTRIBUTE_UNUSED = 0.25*INV(dx*dy); + CCTK_REAL /*const*/ p1o4dxdz CCTK_ATTRIBUTE_UNUSED = 0.25*INV(dx*dz); + CCTK_REAL /*const*/ p1o4dy2 CCTK_ATTRIBUTE_UNUSED = 0.25*INV(SQR(dy)); + CCTK_REAL /*const*/ p1o4dydz CCTK_ATTRIBUTE_UNUSED = 0.25*INV(dy*dz); + CCTK_REAL /*const*/ p1o4dz2 CCTK_ATTRIBUTE_UNUSED = 0.25*INV(SQR(dz)); + CCTK_REAL /*const*/ p1odx CCTK_ATTRIBUTE_UNUSED = INV(dx); + CCTK_REAL /*const*/ p1odx2 CCTK_ATTRIBUTE_UNUSED = INV(SQR(dx)); + CCTK_REAL /*const*/ p1odxdydz CCTK_ATTRIBUTE_UNUSED = INV(dx*dy*dz); + CCTK_REAL /*const*/ p1ody CCTK_ATTRIBUTE_UNUSED = INV(dy); + CCTK_REAL /*const*/ p1ody2 CCTK_ATTRIBUTE_UNUSED = INV(SQR(dy)); + CCTK_REAL /*const*/ p1odz CCTK_ATTRIBUTE_UNUSED = INV(dz); + CCTK_REAL /*const*/ p1odz2 CCTK_ATTRIBUTE_UNUSED = INV(SQR(dz)); + CCTK_REAL /*const*/ pm1o12dx2 CCTK_ATTRIBUTE_UNUSED = -0.0833333333333333333333333333333*INV(SQR(dx)); + CCTK_REAL /*const*/ pm1o12dy2 CCTK_ATTRIBUTE_UNUSED = -0.0833333333333333333333333333333*INV(SQR(dy)); + CCTK_REAL /*const*/ pm1o12dz2 CCTK_ATTRIBUTE_UNUSED = -0.0833333333333333333333333333333*INV(SQR(dz)); + CCTK_REAL /*const*/ pm1o2dx CCTK_ATTRIBUTE_UNUSED = -0.5*INV(dx); + CCTK_REAL /*const*/ pm1o2dy CCTK_ATTRIBUTE_UNUSED = -0.5*INV(dy); + CCTK_REAL /*const*/ pm1o2dz CCTK_ATTRIBUTE_UNUSED = -0.5*INV(dz); /* Assign local copies of arrays functions */ @@ -110,22 +110,22 @@ static void wave_calc_norm_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 phiL = phi[index]; - CCTK_REAL phiErrorL = phiError[index]; - CCTK_REAL piL = pi[index]; - CCTK_REAL piErrorL = piError[index]; + CCTK_REAL phiL CCTK_ATTRIBUTE_UNUSED = phi[index]; + CCTK_REAL phiErrorL CCTK_ATTRIBUTE_UNUSED = phiError[index]; + CCTK_REAL piL CCTK_ATTRIBUTE_UNUSED = pi[index]; + CCTK_REAL piErrorL CCTK_ATTRIBUTE_UNUSED = piError[index]; /* Include user supplied include files */ /* Precompute derivatives */ - CCTK_REAL PDplus1phi; - CCTK_REAL PDplus2phi; - CCTK_REAL PDplus3phi; + CCTK_REAL PDplus1phi CCTK_ATTRIBUTE_UNUSED ; + CCTK_REAL PDplus2phi CCTK_ATTRIBUTE_UNUSED ; + CCTK_REAL PDplus3phi CCTK_ATTRIBUTE_UNUSED ; switch(fdOrder) { @@ -143,18 +143,19 @@ static void wave_calc_norm_Body(cGH const * restrict const cctkGH, int const dir } /* Calculate temporaries and grid functions */ - CCTK_REAL VL2squared = SQR(phiL) + SQR(piL); + CCTK_REAL CCTK_ATTRIBUTE_UNUSED VL2squared = SQR(phiL) + SQR(piL); - CCTK_REAL VL2L = sqrt(VL2squared); + CCTK_REAL CCTK_ATTRIBUTE_UNUSED VL2L = sqrt(VL2squared); - CCTK_REAL VDPsquared = SQR(piL) + SQR(PDplus1phi) + SQR(PDplus2phi) - + SQR(PDplus3phi); + CCTK_REAL CCTK_ATTRIBUTE_UNUSED VDPsquared = SQR(piL) + + SQR(PDplus1phi) + SQR(PDplus2phi) + SQR(PDplus3phi); - CCTK_REAL VDPL = sqrt(VDPsquared); + CCTK_REAL CCTK_ATTRIBUTE_UNUSED VDPL = sqrt(VDPsquared); - CCTK_REAL EL2squared = SQR(phiErrorL) + SQR(piErrorL); + CCTK_REAL CCTK_ATTRIBUTE_UNUSED EL2squared = SQR(phiErrorL) + + SQR(piErrorL); - CCTK_REAL EL2L = sqrt(EL2squared); + CCTK_REAL CCTK_ATTRIBUTE_UNUSED EL2L = sqrt(EL2squared); /* Copy local copies back to grid functions */ EL2[index] = EL2L; diff --git a/Examples/Wave/src/wave_evolve.cc b/Examples/Wave/src/wave_evolve.cc index 5cee251..e0f4c15 100644 --- a/Examples/Wave/src/wave_evolve.cc +++ b/Examples/Wave/src/wave_evolve.cc @@ -27,7 +27,7 @@ extern "C" void wave_evolve_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 */, "Wave::evolvedrhs","flat"); if (ierr < 0) CCTK_WARN(1, "Failed to register flat BC for Wave::evolvedrhs."); @@ -43,58 +43,58 @@ static void wave_evolve_Body(cGH const * restrict const cctkGH, int const dir, i /* 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 p1o12dx = 0.0833333333333333333333333333333*INV(dx); - CCTK_REAL const p1o12dy = 0.0833333333333333333333333333333*INV(dy); - CCTK_REAL const p1o12dz = 0.0833333333333333333333333333333*INV(dz); - CCTK_REAL const p1o144dxdy = 0.00694444444444444444444444444444*INV(dx*dy); - CCTK_REAL const p1o144dxdz = 0.00694444444444444444444444444444*INV(dx*dz); - CCTK_REAL const p1o144dydz = 0.00694444444444444444444444444444*INV(dy*dz); - CCTK_REAL const p1o2dx = 0.5*INV(dx); - CCTK_REAL const p1o2dy = 0.5*INV(dy); - CCTK_REAL const p1o2dz = 0.5*INV(dz); - CCTK_REAL const p1o4dx2 = 0.25*INV(SQR(dx)); - CCTK_REAL const p1o4dxdy = 0.25*INV(dx*dy); - CCTK_REAL const p1o4dxdz = 0.25*INV(dx*dz); - CCTK_REAL const p1o4dy2 = 0.25*INV(SQR(dy)); - CCTK_REAL const p1o4dydz = 0.25*INV(dy*dz); - CCTK_REAL const p1o4dz2 = 0.25*INV(SQR(dz)); - CCTK_REAL const p1odx = INV(dx); - CCTK_REAL const p1odx2 = INV(SQR(dx)); - CCTK_REAL const p1odxdydz = INV(dx*dy*dz); - CCTK_REAL const p1ody = INV(dy); - CCTK_REAL const p1ody2 = INV(SQR(dy)); - CCTK_REAL const p1odz = INV(dz); - CCTK_REAL const p1odz2 = INV(SQR(dz)); - CCTK_REAL const pm1o12dx2 = -0.0833333333333333333333333333333*INV(SQR(dx)); - CCTK_REAL const pm1o12dy2 = -0.0833333333333333333333333333333*INV(SQR(dy)); - CCTK_REAL const pm1o12dz2 = -0.0833333333333333333333333333333*INV(SQR(dz)); - CCTK_REAL const pm1o2dx = -0.5*INV(dx); - CCTK_REAL const pm1o2dy = -0.5*INV(dy); - CCTK_REAL const pm1o2dz = -0.5*INV(dz); + CCTK_REAL /*const*/ p1o12dx CCTK_ATTRIBUTE_UNUSED = 0.0833333333333333333333333333333*INV(dx); + CCTK_REAL /*const*/ p1o12dy CCTK_ATTRIBUTE_UNUSED = 0.0833333333333333333333333333333*INV(dy); + CCTK_REAL /*const*/ p1o12dz CCTK_ATTRIBUTE_UNUSED = 0.0833333333333333333333333333333*INV(dz); + CCTK_REAL /*const*/ p1o144dxdy CCTK_ATTRIBUTE_UNUSED = 0.00694444444444444444444444444444*INV(dx*dy); + CCTK_REAL /*const*/ p1o144dxdz CCTK_ATTRIBUTE_UNUSED = 0.00694444444444444444444444444444*INV(dx*dz); + CCTK_REAL /*const*/ p1o144dydz CCTK_ATTRIBUTE_UNUSED = 0.00694444444444444444444444444444*INV(dy*dz); + CCTK_REAL /*const*/ p1o2dx CCTK_ATTRIBUTE_UNUSED = 0.5*INV(dx); + CCTK_REAL /*const*/ p1o2dy CCTK_ATTRIBUTE_UNUSED = 0.5*INV(dy); + CCTK_REAL /*const*/ p1o2dz CCTK_ATTRIBUTE_UNUSED = 0.5*INV(dz); + CCTK_REAL /*const*/ p1o4dx2 CCTK_ATTRIBUTE_UNUSED = 0.25*INV(SQR(dx)); + CCTK_REAL /*const*/ p1o4dxdy CCTK_ATTRIBUTE_UNUSED = 0.25*INV(dx*dy); + CCTK_REAL /*const*/ p1o4dxdz CCTK_ATTRIBUTE_UNUSED = 0.25*INV(dx*dz); + CCTK_REAL /*const*/ p1o4dy2 CCTK_ATTRIBUTE_UNUSED = 0.25*INV(SQR(dy)); + CCTK_REAL /*const*/ p1o4dydz CCTK_ATTRIBUTE_UNUSED = 0.25*INV(dy*dz); + CCTK_REAL /*const*/ p1o4dz2 CCTK_ATTRIBUTE_UNUSED = 0.25*INV(SQR(dz)); + CCTK_REAL /*const*/ p1odx CCTK_ATTRIBUTE_UNUSED = INV(dx); + CCTK_REAL /*const*/ p1odx2 CCTK_ATTRIBUTE_UNUSED = INV(SQR(dx)); + CCTK_REAL /*const*/ p1odxdydz CCTK_ATTRIBUTE_UNUSED = INV(dx*dy*dz); + CCTK_REAL /*const*/ p1ody CCTK_ATTRIBUTE_UNUSED = INV(dy); + CCTK_REAL /*const*/ p1ody2 CCTK_ATTRIBUTE_UNUSED = INV(SQR(dy)); + CCTK_REAL /*const*/ p1odz CCTK_ATTRIBUTE_UNUSED = INV(dz); + CCTK_REAL /*const*/ p1odz2 CCTK_ATTRIBUTE_UNUSED = INV(SQR(dz)); + CCTK_REAL /*const*/ pm1o12dx2 CCTK_ATTRIBUTE_UNUSED = -0.0833333333333333333333333333333*INV(SQR(dx)); + CCTK_REAL /*const*/ pm1o12dy2 CCTK_ATTRIBUTE_UNUSED = -0.0833333333333333333333333333333*INV(SQR(dy)); + CCTK_REAL /*const*/ pm1o12dz2 CCTK_ATTRIBUTE_UNUSED = -0.0833333333333333333333333333333*INV(SQR(dz)); + CCTK_REAL /*const*/ pm1o2dx CCTK_ATTRIBUTE_UNUSED = -0.5*INV(dx); + CCTK_REAL /*const*/ pm1o2dy CCTK_ATTRIBUTE_UNUSED = -0.5*INV(dy); + CCTK_REAL /*const*/ pm1o2dz CCTK_ATTRIBUTE_UNUSED = -0.5*INV(dz); /* Assign local copies of arrays functions */ @@ -110,20 +110,20 @@ static void wave_evolve_Body(cGH const * restrict const cctkGH, int const dir, i 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 phiL = phi[index]; - CCTK_REAL piL = pi[index]; + CCTK_REAL phiL CCTK_ATTRIBUTE_UNUSED = phi[index]; + CCTK_REAL piL CCTK_ATTRIBUTE_UNUSED = pi[index]; /* Include user supplied include files */ /* Precompute derivatives */ - CCTK_REAL PDstandard11phi; - CCTK_REAL PDstandard22phi; - CCTK_REAL PDstandard33phi; + CCTK_REAL PDstandard11phi CCTK_ATTRIBUTE_UNUSED ; + CCTK_REAL PDstandard22phi CCTK_ATTRIBUTE_UNUSED ; + CCTK_REAL PDstandard33phi CCTK_ATTRIBUTE_UNUSED ; switch(fdOrder) { @@ -141,10 +141,10 @@ static void wave_evolve_Body(cGH const * restrict const cctkGH, int const dir, i } /* Calculate temporaries and grid functions */ - CCTK_REAL phirhsL = piL; + CCTK_REAL CCTK_ATTRIBUTE_UNUSED phirhsL = piL; - CCTK_REAL pirhsL = PDstandard11phi + PDstandard22phi + - PDstandard33phi; + CCTK_REAL CCTK_ATTRIBUTE_UNUSED pirhsL = PDstandard11phi + + PDstandard22phi + PDstandard33phi; /* Copy local copies back to grid functions */ phirhs[index] = phirhsL; diff --git a/Examples/Wave/src/wave_exact_gaussian.cc b/Examples/Wave/src/wave_exact_gaussian.cc index 82c1bc7..9025cf6 100644 --- a/Examples/Wave/src/wave_exact_gaussian.cc +++ b/Examples/Wave/src/wave_exact_gaussian.cc @@ -31,58 +31,58 @@ static void wave_exact_gaussian_Body(cGH const * restrict const cctkGH, int cons /* 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 p1o12dx = 0.0833333333333333333333333333333*INV(dx); - CCTK_REAL const p1o12dy = 0.0833333333333333333333333333333*INV(dy); - CCTK_REAL const p1o12dz = 0.0833333333333333333333333333333*INV(dz); - CCTK_REAL const p1o144dxdy = 0.00694444444444444444444444444444*INV(dx*dy); - CCTK_REAL const p1o144dxdz = 0.00694444444444444444444444444444*INV(dx*dz); - CCTK_REAL const p1o144dydz = 0.00694444444444444444444444444444*INV(dy*dz); - CCTK_REAL const p1o2dx = 0.5*INV(dx); - CCTK_REAL const p1o2dy = 0.5*INV(dy); - CCTK_REAL const p1o2dz = 0.5*INV(dz); - CCTK_REAL const p1o4dx2 = 0.25*INV(SQR(dx)); - CCTK_REAL const p1o4dxdy = 0.25*INV(dx*dy); - CCTK_REAL const p1o4dxdz = 0.25*INV(dx*dz); - CCTK_REAL const p1o4dy2 = 0.25*INV(SQR(dy)); - CCTK_REAL const p1o4dydz = 0.25*INV(dy*dz); - CCTK_REAL const p1o4dz2 = 0.25*INV(SQR(dz)); - CCTK_REAL const p1odx = INV(dx); - CCTK_REAL const p1odx2 = INV(SQR(dx)); - CCTK_REAL const p1odxdydz = INV(dx*dy*dz); - CCTK_REAL const p1ody = INV(dy); - CCTK_REAL const p1ody2 = INV(SQR(dy)); - CCTK_REAL const p1odz = INV(dz); - CCTK_REAL const p1odz2 = INV(SQR(dz)); - CCTK_REAL const pm1o12dx2 = -0.0833333333333333333333333333333*INV(SQR(dx)); - CCTK_REAL const pm1o12dy2 = -0.0833333333333333333333333333333*INV(SQR(dy)); - CCTK_REAL const pm1o12dz2 = -0.0833333333333333333333333333333*INV(SQR(dz)); - CCTK_REAL const pm1o2dx = -0.5*INV(dx); - CCTK_REAL const pm1o2dy = -0.5*INV(dy); - CCTK_REAL const pm1o2dz = -0.5*INV(dz); + CCTK_REAL /*const*/ p1o12dx CCTK_ATTRIBUTE_UNUSED = 0.0833333333333333333333333333333*INV(dx); + CCTK_REAL /*const*/ p1o12dy CCTK_ATTRIBUTE_UNUSED = 0.0833333333333333333333333333333*INV(dy); + CCTK_REAL /*const*/ p1o12dz CCTK_ATTRIBUTE_UNUSED = 0.0833333333333333333333333333333*INV(dz); + CCTK_REAL /*const*/ p1o144dxdy CCTK_ATTRIBUTE_UNUSED = 0.00694444444444444444444444444444*INV(dx*dy); + CCTK_REAL /*const*/ p1o144dxdz CCTK_ATTRIBUTE_UNUSED = 0.00694444444444444444444444444444*INV(dx*dz); + CCTK_REAL /*const*/ p1o144dydz CCTK_ATTRIBUTE_UNUSED = 0.00694444444444444444444444444444*INV(dy*dz); + CCTK_REAL /*const*/ p1o2dx CCTK_ATTRIBUTE_UNUSED = 0.5*INV(dx); + CCTK_REAL /*const*/ p1o2dy CCTK_ATTRIBUTE_UNUSED = 0.5*INV(dy); + CCTK_REAL /*const*/ p1o2dz CCTK_ATTRIBUTE_UNUSED = 0.5*INV(dz); + CCTK_REAL /*const*/ p1o4dx2 CCTK_ATTRIBUTE_UNUSED = 0.25*INV(SQR(dx)); + CCTK_REAL /*const*/ p1o4dxdy CCTK_ATTRIBUTE_UNUSED = 0.25*INV(dx*dy); + CCTK_REAL /*const*/ p1o4dxdz CCTK_ATTRIBUTE_UNUSED = 0.25*INV(dx*dz); + CCTK_REAL /*const*/ p1o4dy2 CCTK_ATTRIBUTE_UNUSED = 0.25*INV(SQR(dy)); + CCTK_REAL /*const*/ p1o4dydz CCTK_ATTRIBUTE_UNUSED = 0.25*INV(dy*dz); + CCTK_REAL /*const*/ p1o4dz2 CCTK_ATTRIBUTE_UNUSED = 0.25*INV(SQR(dz)); + CCTK_REAL /*const*/ p1odx CCTK_ATTRIBUTE_UNUSED = INV(dx); + CCTK_REAL /*const*/ p1odx2 CCTK_ATTRIBUTE_UNUSED = INV(SQR(dx)); + CCTK_REAL /*const*/ p1odxdydz CCTK_ATTRIBUTE_UNUSED = INV(dx*dy*dz); + CCTK_REAL /*const*/ p1ody CCTK_ATTRIBUTE_UNUSED = INV(dy); + CCTK_REAL /*const*/ p1ody2 CCTK_ATTRIBUTE_UNUSED = INV(SQR(dy)); + CCTK_REAL /*const*/ p1odz CCTK_ATTRIBUTE_UNUSED = INV(dz); + CCTK_REAL /*const*/ p1odz2 CCTK_ATTRIBUTE_UNUSED = INV(SQR(dz)); + CCTK_REAL /*const*/ pm1o12dx2 CCTK_ATTRIBUTE_UNUSED = -0.0833333333333333333333333333333*INV(SQR(dx)); + CCTK_REAL /*const*/ pm1o12dy2 CCTK_ATTRIBUTE_UNUSED = -0.0833333333333333333333333333333*INV(SQR(dy)); + CCTK_REAL /*const*/ pm1o12dz2 CCTK_ATTRIBUTE_UNUSED = -0.0833333333333333333333333333333*INV(SQR(dz)); + CCTK_REAL /*const*/ pm1o2dx CCTK_ATTRIBUTE_UNUSED = -0.5*INV(dx); + CCTK_REAL /*const*/ pm1o2dy CCTK_ATTRIBUTE_UNUSED = -0.5*INV(dy); + CCTK_REAL /*const*/ pm1o2dz CCTK_ATTRIBUTE_UNUSED = -0.5*INV(dz); /* Assign local copies of arrays functions */ @@ -98,11 +98,11 @@ static void wave_exact_gaussian_Body(cGH const * restrict const cctkGH, int cons 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 rL = r[index]; + CCTK_REAL rL CCTK_ATTRIBUTE_UNUSED = r[index]; /* Include user supplied include files */ @@ -119,22 +119,23 @@ static void wave_exact_gaussian_Body(cGH const * restrict const cctkGH, int cons } /* Calculate temporaries and grid functions */ - CCTK_REAL rEps = pow(1.e-24 + QAD(rL),0.25); + CCTK_REAL CCTK_ATTRIBUTE_UNUSED rEps = pow(1.e-24 + QAD(rL),0.25); - CCTK_REAL phiExactL = (-(CUB(-rL + t + - ToReal(t0))*exp(-(INV(SQR(ToReal(nSigma)))*SQR(-rL + t + - ToReal(t0))))) + CUB(rL + t + - ToReal(t0))*exp(-(INV(SQR(ToReal(nSigma)))*SQR(rL + t + - ToReal(t0)))))*INV(rEps); + CCTK_REAL CCTK_ATTRIBUTE_UNUSED phiExactL = (-1.*CUB(-1.*rL + t + + ToReal(t0))*exp(-1.*INV(SQR(ToReal(nSigma)))*SQR(-1.*rL + t + + ToReal(t0))) + CUB(rL + t + + ToReal(t0))*exp(-1.*INV(SQR(ToReal(nSigma)))*SQR(rL + t + + ToReal(t0))))*INV(rEps); - CCTK_REAL piExactL = - INV(rEps)*(INV(SQR(ToReal(nSigma)))*(2*exp(-(INV(SQR(ToReal(nSigma)))*SQR(-rL - + t + ToReal(t0))))*QAD(-rL + t + ToReal(t0)) - - 2*exp(-(INV(SQR(ToReal(nSigma)))*SQR(rL + t + ToReal(t0))))*QAD(rL - + t + ToReal(t0))) - 3*exp(-(INV(SQR(ToReal(nSigma)))*SQR(-rL + t + - ToReal(t0))))*SQR(-rL + t + ToReal(t0)) + - 3*exp(-(INV(SQR(ToReal(nSigma)))*SQR(rL + t + ToReal(t0))))*SQR(rL - + t + ToReal(t0))); + CCTK_REAL CCTK_ATTRIBUTE_UNUSED piExactL = + INV(rEps)*(INV(SQR(ToReal(nSigma)))*(2.*exp(-1.*INV(SQR(ToReal(nSigma)))*SQR(-1.*rL + + t + ToReal(t0)))*QAD(-1.*rL + t + ToReal(t0)) - + 2.*exp(-1.*INV(SQR(ToReal(nSigma)))*SQR(rL + t + + ToReal(t0)))*QAD(rL + t + ToReal(t0))) - + 3.*exp(-1.*INV(SQR(ToReal(nSigma)))*SQR(-1.*rL + t + + ToReal(t0)))*SQR(-1.*rL + t + ToReal(t0)) + + 3.*exp(-1.*INV(SQR(ToReal(nSigma)))*SQR(rL + t + + ToReal(t0)))*SQR(rL + t + ToReal(t0))); /* Copy local copies back to grid functions */ phiExact[index] = phiExactL; diff --git a/Examples/Wave/src/wave_exact_sine.cc b/Examples/Wave/src/wave_exact_sine.cc index 3f6dc52..d04c310 100644 --- a/Examples/Wave/src/wave_exact_sine.cc +++ b/Examples/Wave/src/wave_exact_sine.cc @@ -31,58 +31,58 @@ static void wave_exact_sine_Body(cGH const * restrict const cctkGH, int const di /* 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 p1o12dx = 0.0833333333333333333333333333333*INV(dx); - CCTK_REAL const p1o12dy = 0.0833333333333333333333333333333*INV(dy); - CCTK_REAL const p1o12dz = 0.0833333333333333333333333333333*INV(dz); - CCTK_REAL const p1o144dxdy = 0.00694444444444444444444444444444*INV(dx*dy); - CCTK_REAL const p1o144dxdz = 0.00694444444444444444444444444444*INV(dx*dz); - CCTK_REAL const p1o144dydz = 0.00694444444444444444444444444444*INV(dy*dz); - CCTK_REAL const p1o2dx = 0.5*INV(dx); - CCTK_REAL const p1o2dy = 0.5*INV(dy); - CCTK_REAL const p1o2dz = 0.5*INV(dz); - CCTK_REAL const p1o4dx2 = 0.25*INV(SQR(dx)); - CCTK_REAL const p1o4dxdy = 0.25*INV(dx*dy); - CCTK_REAL const p1o4dxdz = 0.25*INV(dx*dz); - CCTK_REAL const p1o4dy2 = 0.25*INV(SQR(dy)); - CCTK_REAL const p1o4dydz = 0.25*INV(dy*dz); - CCTK_REAL const p1o4dz2 = 0.25*INV(SQR(dz)); - CCTK_REAL const p1odx = INV(dx); - CCTK_REAL const p1odx2 = INV(SQR(dx)); - CCTK_REAL const p1odxdydz = INV(dx*dy*dz); - CCTK_REAL const p1ody = INV(dy); - CCTK_REAL const p1ody2 = INV(SQR(dy)); - CCTK_REAL const p1odz = INV(dz); - CCTK_REAL const p1odz2 = INV(SQR(dz)); - CCTK_REAL const pm1o12dx2 = -0.0833333333333333333333333333333*INV(SQR(dx)); - CCTK_REAL const pm1o12dy2 = -0.0833333333333333333333333333333*INV(SQR(dy)); - CCTK_REAL const pm1o12dz2 = -0.0833333333333333333333333333333*INV(SQR(dz)); - CCTK_REAL const pm1o2dx = -0.5*INV(dx); - CCTK_REAL const pm1o2dy = -0.5*INV(dy); - CCTK_REAL const pm1o2dz = -0.5*INV(dz); + CCTK_REAL /*const*/ p1o12dx CCTK_ATTRIBUTE_UNUSED = 0.0833333333333333333333333333333*INV(dx); + CCTK_REAL /*const*/ p1o12dy CCTK_ATTRIBUTE_UNUSED = 0.0833333333333333333333333333333*INV(dy); + CCTK_REAL /*const*/ p1o12dz CCTK_ATTRIBUTE_UNUSED = 0.0833333333333333333333333333333*INV(dz); + CCTK_REAL /*const*/ p1o144dxdy CCTK_ATTRIBUTE_UNUSED = 0.00694444444444444444444444444444*INV(dx*dy); + CCTK_REAL /*const*/ p1o144dxdz CCTK_ATTRIBUTE_UNUSED = 0.00694444444444444444444444444444*INV(dx*dz); + CCTK_REAL /*const*/ p1o144dydz CCTK_ATTRIBUTE_UNUSED = 0.00694444444444444444444444444444*INV(dy*dz); + CCTK_REAL /*const*/ p1o2dx CCTK_ATTRIBUTE_UNUSED = 0.5*INV(dx); + CCTK_REAL /*const*/ p1o2dy CCTK_ATTRIBUTE_UNUSED = 0.5*INV(dy); + CCTK_REAL /*const*/ p1o2dz CCTK_ATTRIBUTE_UNUSED = 0.5*INV(dz); + CCTK_REAL /*const*/ p1o4dx2 CCTK_ATTRIBUTE_UNUSED = 0.25*INV(SQR(dx)); + CCTK_REAL /*const*/ p1o4dxdy CCTK_ATTRIBUTE_UNUSED = 0.25*INV(dx*dy); + CCTK_REAL /*const*/ p1o4dxdz CCTK_ATTRIBUTE_UNUSED = 0.25*INV(dx*dz); + CCTK_REAL /*const*/ p1o4dy2 CCTK_ATTRIBUTE_UNUSED = 0.25*INV(SQR(dy)); + CCTK_REAL /*const*/ p1o4dydz CCTK_ATTRIBUTE_UNUSED = 0.25*INV(dy*dz); + CCTK_REAL /*const*/ p1o4dz2 CCTK_ATTRIBUTE_UNUSED = 0.25*INV(SQR(dz)); + CCTK_REAL /*const*/ p1odx CCTK_ATTRIBUTE_UNUSED = INV(dx); + CCTK_REAL /*const*/ p1odx2 CCTK_ATTRIBUTE_UNUSED = INV(SQR(dx)); + CCTK_REAL /*const*/ p1odxdydz CCTK_ATTRIBUTE_UNUSED = INV(dx*dy*dz); + CCTK_REAL /*const*/ p1ody CCTK_ATTRIBUTE_UNUSED = INV(dy); + CCTK_REAL /*const*/ p1ody2 CCTK_ATTRIBUTE_UNUSED = INV(SQR(dy)); + CCTK_REAL /*const*/ p1odz CCTK_ATTRIBUTE_UNUSED = INV(dz); + CCTK_REAL /*const*/ p1odz2 CCTK_ATTRIBUTE_UNUSED = INV(SQR(dz)); + CCTK_REAL /*const*/ pm1o12dx2 CCTK_ATTRIBUTE_UNUSED = -0.0833333333333333333333333333333*INV(SQR(dx)); + CCTK_REAL /*const*/ pm1o12dy2 CCTK_ATTRIBUTE_UNUSED = -0.0833333333333333333333333333333*INV(SQR(dy)); + CCTK_REAL /*const*/ pm1o12dz2 CCTK_ATTRIBUTE_UNUSED = -0.0833333333333333333333333333333*INV(SQR(dz)); + CCTK_REAL /*const*/ pm1o2dx CCTK_ATTRIBUTE_UNUSED = -0.5*INV(dx); + CCTK_REAL /*const*/ pm1o2dy CCTK_ATTRIBUTE_UNUSED = -0.5*INV(dy); + CCTK_REAL /*const*/ pm1o2dz CCTK_ATTRIBUTE_UNUSED = -0.5*INV(dz); /* Assign local copies of arrays functions */ @@ -98,13 +98,13 @@ static void wave_exact_sine_Body(cGH const * restrict const cctkGH, int const di 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 yL = y[index]; - CCTK_REAL zL = z[index]; + CCTK_REAL xL CCTK_ATTRIBUTE_UNUSED = x[index]; + CCTK_REAL yL CCTK_ATTRIBUTE_UNUSED = y[index]; + CCTK_REAL zL CCTK_ATTRIBUTE_UNUSED = z[index]; /* Include user supplied include files */ @@ -121,16 +121,16 @@ static void wave_exact_sine_Body(cGH const * restrict const cctkGH, int const di } /* Calculate temporaries and grid functions */ - CCTK_REAL piconst = 3.1415926535897932385; + CCTK_REAL CCTK_ATTRIBUTE_UNUSED piconst = 3.1415926535897932385; - CCTK_REAL phiExactL = - sin(2*piconst*INV(ToReal(periodicity))*(-(t*sqrt(SQR(ToReal(n1)) + - SQR(ToReal(n2)) + SQR(ToReal(n3)))) + xL*ToReal(n1) + yL*ToReal(n2) + CCTK_REAL CCTK_ATTRIBUTE_UNUSED phiExactL = + sin(2.*piconst*INV(ToReal(periodicity))*(-1.*t*sqrt(SQR(ToReal(n1)) + + SQR(ToReal(n2)) + SQR(ToReal(n3))) + xL*ToReal(n1) + yL*ToReal(n2) + zL*ToReal(n3)))*ToReal(amplitude); - CCTK_REAL piExactL = - -2*piconst*cos(2*piconst*INV(ToReal(periodicity))*(-(t*sqrt(SQR(ToReal(n1)) - + SQR(ToReal(n2)) + SQR(ToReal(n3)))) + xL*ToReal(n1) + + CCTK_REAL CCTK_ATTRIBUTE_UNUSED piExactL = + -2.*piconst*cos(2.*piconst*INV(ToReal(periodicity))*(-1.*t*sqrt(SQR(ToReal(n1)) + + SQR(ToReal(n2)) + SQR(ToReal(n3))) + xL*ToReal(n1) + yL*ToReal(n2) + zL*ToReal(n3)))*INV(ToReal(periodicity))*sqrt(SQR(ToReal(n1)) + SQR(ToReal(n2)) + SQR(ToReal(n3)))*ToReal(amplitude); diff --git a/Examples/Wave/src/wave_import_exact.cc b/Examples/Wave/src/wave_import_exact.cc index 8f8b0f8..e9ece99 100644 --- a/Examples/Wave/src/wave_import_exact.cc +++ b/Examples/Wave/src/wave_import_exact.cc @@ -31,58 +31,58 @@ static void wave_import_exact_Body(cGH const * restrict const cctkGH, int const /* 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 p1o12dx = 0.0833333333333333333333333333333*INV(dx); - CCTK_REAL const p1o12dy = 0.0833333333333333333333333333333*INV(dy); - CCTK_REAL const p1o12dz = 0.0833333333333333333333333333333*INV(dz); - CCTK_REAL const p1o144dxdy = 0.00694444444444444444444444444444*INV(dx*dy); - CCTK_REAL const p1o144dxdz = 0.00694444444444444444444444444444*INV(dx*dz); - CCTK_REAL const p1o144dydz = 0.00694444444444444444444444444444*INV(dy*dz); - CCTK_REAL const p1o2dx = 0.5*INV(dx); - CCTK_REAL const p1o2dy = 0.5*INV(dy); - CCTK_REAL const p1o2dz = 0.5*INV(dz); - CCTK_REAL const p1o4dx2 = 0.25*INV(SQR(dx)); - CCTK_REAL const p1o4dxdy = 0.25*INV(dx*dy); - CCTK_REAL const p1o4dxdz = 0.25*INV(dx*dz); - CCTK_REAL const p1o4dy2 = 0.25*INV(SQR(dy)); - CCTK_REAL const p1o4dydz = 0.25*INV(dy*dz); - CCTK_REAL const p1o4dz2 = 0.25*INV(SQR(dz)); - CCTK_REAL const p1odx = INV(dx); - CCTK_REAL const p1odx2 = INV(SQR(dx)); - CCTK_REAL const p1odxdydz = INV(dx*dy*dz); - CCTK_REAL const p1ody = INV(dy); - CCTK_REAL const p1ody2 = INV(SQR(dy)); - CCTK_REAL const p1odz = INV(dz); - CCTK_REAL const p1odz2 = INV(SQR(dz)); - CCTK_REAL const pm1o12dx2 = -0.0833333333333333333333333333333*INV(SQR(dx)); - CCTK_REAL const pm1o12dy2 = -0.0833333333333333333333333333333*INV(SQR(dy)); - CCTK_REAL const pm1o12dz2 = -0.0833333333333333333333333333333*INV(SQR(dz)); - CCTK_REAL const pm1o2dx = -0.5*INV(dx); - CCTK_REAL const pm1o2dy = -0.5*INV(dy); - CCTK_REAL const pm1o2dz = -0.5*INV(dz); + CCTK_REAL /*const*/ p1o12dx CCTK_ATTRIBUTE_UNUSED = 0.0833333333333333333333333333333*INV(dx); + CCTK_REAL /*const*/ p1o12dy CCTK_ATTRIBUTE_UNUSED = 0.0833333333333333333333333333333*INV(dy); + CCTK_REAL /*const*/ p1o12dz CCTK_ATTRIBUTE_UNUSED = 0.0833333333333333333333333333333*INV(dz); + CCTK_REAL /*const*/ p1o144dxdy CCTK_ATTRIBUTE_UNUSED = 0.00694444444444444444444444444444*INV(dx*dy); + CCTK_REAL /*const*/ p1o144dxdz CCTK_ATTRIBUTE_UNUSED = 0.00694444444444444444444444444444*INV(dx*dz); + CCTK_REAL /*const*/ p1o144dydz CCTK_ATTRIBUTE_UNUSED = 0.00694444444444444444444444444444*INV(dy*dz); + CCTK_REAL /*const*/ p1o2dx CCTK_ATTRIBUTE_UNUSED = 0.5*INV(dx); + CCTK_REAL /*const*/ p1o2dy CCTK_ATTRIBUTE_UNUSED = 0.5*INV(dy); + CCTK_REAL /*const*/ p1o2dz CCTK_ATTRIBUTE_UNUSED = 0.5*INV(dz); + CCTK_REAL /*const*/ p1o4dx2 CCTK_ATTRIBUTE_UNUSED = 0.25*INV(SQR(dx)); + CCTK_REAL /*const*/ p1o4dxdy CCTK_ATTRIBUTE_UNUSED = 0.25*INV(dx*dy); + CCTK_REAL /*const*/ p1o4dxdz CCTK_ATTRIBUTE_UNUSED = 0.25*INV(dx*dz); + CCTK_REAL /*const*/ p1o4dy2 CCTK_ATTRIBUTE_UNUSED = 0.25*INV(SQR(dy)); + CCTK_REAL /*const*/ p1o4dydz CCTK_ATTRIBUTE_UNUSED = 0.25*INV(dy*dz); + CCTK_REAL /*const*/ p1o4dz2 CCTK_ATTRIBUTE_UNUSED = 0.25*INV(SQR(dz)); + CCTK_REAL /*const*/ p1odx CCTK_ATTRIBUTE_UNUSED = INV(dx); + CCTK_REAL /*const*/ p1odx2 CCTK_ATTRIBUTE_UNUSED = INV(SQR(dx)); + CCTK_REAL /*const*/ p1odxdydz CCTK_ATTRIBUTE_UNUSED = INV(dx*dy*dz); + CCTK_REAL /*const*/ p1ody CCTK_ATTRIBUTE_UNUSED = INV(dy); + CCTK_REAL /*const*/ p1ody2 CCTK_ATTRIBUTE_UNUSED = INV(SQR(dy)); + CCTK_REAL /*const*/ p1odz CCTK_ATTRIBUTE_UNUSED = INV(dz); + CCTK_REAL /*const*/ p1odz2 CCTK_ATTRIBUTE_UNUSED = INV(SQR(dz)); + CCTK_REAL /*const*/ pm1o12dx2 CCTK_ATTRIBUTE_UNUSED = -0.0833333333333333333333333333333*INV(SQR(dx)); + CCTK_REAL /*const*/ pm1o12dy2 CCTK_ATTRIBUTE_UNUSED = -0.0833333333333333333333333333333*INV(SQR(dy)); + CCTK_REAL /*const*/ pm1o12dz2 CCTK_ATTRIBUTE_UNUSED = -0.0833333333333333333333333333333*INV(SQR(dz)); + CCTK_REAL /*const*/ pm1o2dx CCTK_ATTRIBUTE_UNUSED = -0.5*INV(dx); + CCTK_REAL /*const*/ pm1o2dy CCTK_ATTRIBUTE_UNUSED = -0.5*INV(dy); + CCTK_REAL /*const*/ pm1o2dz CCTK_ATTRIBUTE_UNUSED = -0.5*INV(dz); /* Assign local copies of arrays functions */ @@ -98,12 +98,12 @@ static void wave_import_exact_Body(cGH const * restrict const cctkGH, int const 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 phiExactL = phiExact[index]; - CCTK_REAL piExactL = piExact[index]; + CCTK_REAL phiExactL CCTK_ATTRIBUTE_UNUSED = phiExact[index]; + CCTK_REAL piExactL CCTK_ATTRIBUTE_UNUSED = piExact[index]; /* Include user supplied include files */ @@ -120,9 +120,9 @@ static void wave_import_exact_Body(cGH const * restrict const cctkGH, int const } /* Calculate temporaries and grid functions */ - CCTK_REAL phiL = phiExactL; + CCTK_REAL CCTK_ATTRIBUTE_UNUSED phiL = phiExactL; - CCTK_REAL piL = piExactL; + CCTK_REAL CCTK_ATTRIBUTE_UNUSED piL = piExactL; /* Copy local copies back to grid functions */ phi[index] = phiL; |