diff options
Diffstat (limited to 'Source/FieldSolver/FiniteDifferenceSolver/FiniteDifferenceAlgorithms/YeeAlgorithm.H')
| -rw-r--r-- | Source/FieldSolver/FiniteDifferenceSolver/FiniteDifferenceAlgorithms/YeeAlgorithm.H | 24 |
1 files changed, 12 insertions, 12 deletions
diff --git a/Source/FieldSolver/FiniteDifferenceSolver/FiniteDifferenceAlgorithms/YeeAlgorithm.H b/Source/FieldSolver/FiniteDifferenceSolver/FiniteDifferenceAlgorithms/YeeAlgorithm.H index 0ad9c6345..75d191467 100644 --- a/Source/FieldSolver/FiniteDifferenceSolver/FiniteDifferenceAlgorithms/YeeAlgorithm.H +++ b/Source/FieldSolver/FiniteDifferenceSolver/FiniteDifferenceAlgorithms/YeeAlgorithm.H @@ -61,12 +61,12 @@ struct YeeAlgorithm { amrex::Real const* coefs_y, int const n_coefs_y, int const i, int const j, int const k ) { - #if defined WARPX_DIM_3D +#if defined WARPX_DIM_3D amrex::Real const inv_dy = coefs_y[0]; return inv_dy*( F(i,j+1,k) - F(i,j,k) ); - #elif (defined WARPX_DIM_XZ) +#elif (defined WARPX_DIM_XZ) return 0; // 2D Cartesian: derivative along y is 0 - #endif +#endif }; /** @@ -77,12 +77,12 @@ struct YeeAlgorithm { amrex::Real const* coefs_y, int const n_coefs_y, int const i, int const j, int const k ) { - #if defined WARPX_DIM_3D +#if defined WARPX_DIM_3D amrex::Real const inv_dy = coefs_y[0]; return inv_dy*( F(i,j,k) - F(i,j-1,k) ); - #elif (defined WARPX_DIM_XZ) +#elif (defined WARPX_DIM_XZ) return 0; // 2D Cartesian: derivative along y is 0 - #endif +#endif }; /** @@ -94,11 +94,11 @@ struct YeeAlgorithm { int const i, int const j, int const k ) { amrex::Real const inv_dz = coefs_z[0]; - #if defined WARPX_DIM_3D +#if defined WARPX_DIM_3D return inv_dz*( F(i,j,k+1) - F(i,j,k) ); - #elif (defined WARPX_DIM_XZ) +#elif (defined WARPX_DIM_XZ) return inv_dz*( F(i,j+1,k) - F(i,j,k) ); - #endif +#endif }; /** @@ -110,11 +110,11 @@ struct YeeAlgorithm { int const i, int const j, int const k ) { amrex::Real const inv_dz = coefs_z[0]; - #if defined WARPX_DIM_3D +#if defined WARPX_DIM_3D return inv_dz*( F(i,j,k) - F(i,j,k-1) ); - #elif (defined WARPX_DIM_XZ) +#elif (defined WARPX_DIM_XZ) return inv_dz*( F(i,j,k) - F(i,j-1,k) ); - #endif +#endif }; }; |