diff options
| author |  Remi Lehe <remi.lehe@normalesup.org>
| 2022-02-16 10:58:26 -0800 |
|---|---|---|
| committer |  GitHub <noreply@github.com>
| 2022-02-16 18:58:26 +0000 |
| commit | 64e8c51a26202def1d9aec73ec43ed07cb82a9e2 (patch) | |
| tree | cdba16427b780824939c879cb3f349e670535cc0 /Source/FieldSolver/ElectrostaticSolver.cpp | |
| parent | 455c8cdfda44962e18a0476afa52e9337cfbbc99 (diff) | |
| download | WarpX-64e8c51a26202def1d9aec73ec43ed07cb82a9e2.tar.gz WarpX-64e8c51a26202def1d9aec73ec43ed07cb82a9e2.tar.zst WarpX-64e8c51a26202def1d9aec73ec43ed07cb82a9e2.zip | |
Correct typo in the relativistic Poisson solver (#2853)
* Correct typo in the relativistic Poisson solver
* Fix unused variable
* Update benchmark
Diffstat (limited to 'Source/FieldSolver/ElectrostaticSolver.cpp')
| -rw-r--r-- | Source/FieldSolver/ElectrostaticSolver.cpp | 9 |
1 files changed, 5 insertions, 4 deletions
diff --git a/Source/FieldSolver/ElectrostaticSolver.cpp b/Source/FieldSolver/ElectrostaticSolver.cpp index e1d70be16..d0922b029 100644 --- a/Source/FieldSolver/ElectrostaticSolver.cpp +++ b/Source/FieldSolver/ElectrostaticSolver.cpp @@ -511,7 +511,7 @@ WarpX::computeE (amrex::Vector<std::array<std::unique_ptr<amrex::MultiFab>, 3> > + phi_arr(i+1,j,k+1)-phi_arr(i-1,j,k+1)) +beta_z*beta_y*0.25_rt*inv_dy*(phi_arr(i,j+1,k )-phi_arr(i,j-1,k ) + phi_arr(i,j+1,k+1)-phi_arr(i,j-1,k+1)) - +(beta_y*beta_z-1)*inv_dz*( phi_arr(i,j,k+1)-phi_arr(i,j,k) ); + +(beta_z*beta_z-1)*inv_dz*( phi_arr(i,j,k+1)-phi_arr(i,j,k) ); } ); #elif defined(WARPX_DIM_XZ) || defined(WARPX_DIM_RZ) @@ -526,17 +526,18 @@ WarpX::computeE (amrex::Vector<std::array<std::unique_ptr<amrex::MultiFab>, 3> > Ez_arr(i,j,k) += +beta_z*beta_x*0.25_rt*inv_dx*(phi_arr(i+1,j ,k)-phi_arr(i-1,j ,k) + phi_arr(i+1,j+1,k)-phi_arr(i-1,j+1,k)) - +(beta_y*beta_z-1)*inv_dz*( phi_arr(i,j+1,k)-phi_arr(i,j,k) ); + +(beta_z*beta_z-1)*inv_dz*( phi_arr(i,j+1,k)-phi_arr(i,j,k) ); } ); + ignore_unused(beta_y); #else amrex::ParallelFor( tbz, [=] AMREX_GPU_DEVICE (int i, int j, int k) { Ez_arr(i,j,k) += - +(beta_y*beta_z-1)*inv_dz*( phi_arr(i+1,j,k)-phi_arr(i,j,k) ); + +(beta_z*beta_z-1)*inv_dz*( phi_arr(i+1,j,k)-phi_arr(i,j,k) ); } ); - ignore_unused(beta_x); + ignore_unused(beta_x,beta_y); #endif } } |