diff options
Diffstat (limited to 'Source/FieldSolver')
| -rw-r--r-- | Source/FieldSolver/FiniteDifferenceSolver/FiniteDifferenceAlgorithms/CKCAlgorithm.H | 8 | ||||
| -rw-r--r-- | Source/FieldSolver/FiniteDifferenceSolver/FiniteDifferenceSolver.H | 7 |
2 files changed, 8 insertions, 7 deletions
diff --git a/Source/FieldSolver/FiniteDifferenceSolver/FiniteDifferenceAlgorithms/CKCAlgorithm.H b/Source/FieldSolver/FiniteDifferenceSolver/FiniteDifferenceAlgorithms/CKCAlgorithm.H index d31693a71..5b0d5e718 100644 --- a/Source/FieldSolver/FiniteDifferenceSolver/FiniteDifferenceAlgorithms/CKCAlgorithm.H +++ b/Source/FieldSolver/FiniteDifferenceSolver/FiniteDifferenceAlgorithms/CKCAlgorithm.H @@ -176,17 +176,17 @@ struct CKCAlgorithm { amrex::Real gammaz = coefs_z[4]; #if defined WARPX_DIM_3D return alphaz * (F(i ,j ,k+1) - F(i ,j ,k )) - - betazx * (F(i+1,j ,k+1) - F(i+1,j ,k ) + + betazx * (F(i+1,j ,k+1) - F(i+1,j ,k ) + F(i-1,j ,k+1) - F(i-1,j ,k )) - - betazy * (F(i ,j+1,k+1) - F(i ,j+1,k ) + + betazy * (F(i ,j+1,k+1) - F(i ,j+1,k ) + F(i ,j-1,k+1) - F(i ,j-1,k )) - - gammaz * (F(i+1,j+1,k+1) - F(i+1,j+1,k ) + + gammaz * (F(i+1,j+1,k+1) - F(i+1,j+1,k ) + F(i-1,j+1,k+1) - F(i-1,j+1,k ) + F(i+1,j-1,k+1) - F(i+1,j-1,k ) + F(i-1,j-1,k+1) - F(i-1,j-1,k )); #elif (defined WARPX_DIM_XZ) return alphaz * (F(i ,j+1,k ) - F(i ,j ,k )) - - betazx * (F(i+1,j+1,k ) - F(i+1,j ,k ) + + betazx * (F(i+1,j+1,k ) - F(i+1,j ,k ) + F(i-1,j+1,k ) - F(i-1,j ,k )); #endif }; diff --git a/Source/FieldSolver/FiniteDifferenceSolver/FiniteDifferenceSolver.H b/Source/FieldSolver/FiniteDifferenceSolver/FiniteDifferenceSolver.H index 8653ce62b..4c0e1955d 100644 --- a/Source/FieldSolver/FiniteDifferenceSolver/FiniteDifferenceSolver.H +++ b/Source/FieldSolver/FiniteDifferenceSolver/FiniteDifferenceSolver.H @@ -3,6 +3,7 @@ #include "WarpXAlgorithmSelection.H" #include "FiniteDifferenceAlgorithms/YeeAlgorithm.H" +#include "FiniteDifferenceAlgorithms/CKCAlgorithm.H" #include <AMReX_MultiFab.H> /** @@ -27,9 +28,9 @@ class FiniteDifferenceSolver if (fdtd_algo == MaxwellSolverAlgo::Yee){ YeeAlgorithm::InitializeStencilCoefficients( cell_size, stencil_coefs_x, stencil_coefs_y, stencil_coefs_z ); -// } else if (fdtd_algo == MaxwellSolverAlgo::CKC) { -// CKCAlgorithm::InitializeStencilCoefficients( cell_size, -// stencil_coefs_x, stencil_coefs_y, stencil_coefs_z ); + } else if (fdtd_algo == MaxwellSolverAlgo::CKC) { + CKCAlgorithm::InitializeStencilCoefficients( cell_size, + stencil_coefs_x, stencil_coefs_y, stencil_coefs_z ); } else { amrex::Abort("Unknown algorithm"); } |