aboutsummaryrefslogtreecommitdiff
path: root/Source/FieldSolver/FiniteDifferenceSolver/FiniteDifferenceAlgorithms
diff options
context:
space:
mode:
authorGravatar Remi Lehe <remi.lehe@normalesup.org> 2020-01-24 22:00:56 -0800
committerGravatar Remi Lehe <remi.lehe@normalesup.org> 2020-01-24 22:00:56 -0800
commitbfd3d3bf383bae468237f556a910afc83bacd6a9 (patch)
treef6b2064b40846b18f1c15a9b1943714e6a988379 /Source/FieldSolver/FiniteDifferenceSolver/FiniteDifferenceAlgorithms
parent2fa99b8e54678226fb19c40343ef00f4a68a851c (diff)
downloadWarpX-bfd3d3bf383bae468237f556a910afc83bacd6a9.tar.gz
WarpX-bfd3d3bf383bae468237f556a910afc83bacd6a9.tar.zst
WarpX-bfd3d3bf383bae468237f556a910afc83bacd6a9.zip
Bug fixes for CKC
Diffstat (limited to 'Source/FieldSolver/FiniteDifferenceSolver/FiniteDifferenceAlgorithms')
-rw-r--r--Source/FieldSolver/FiniteDifferenceSolver/FiniteDifferenceAlgorithms/CKCAlgorithm.H8
1 files changed, 4 insertions, 4 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
};