diff options
Diffstat (limited to 'Source/FieldSolver/WarpXPushFieldsEM.cpp')
| -rw-r--r-- | Source/FieldSolver/WarpXPushFieldsEM.cpp | 6 |
1 files changed, 5 insertions, 1 deletions
diff --git a/Source/FieldSolver/WarpXPushFieldsEM.cpp b/Source/FieldSolver/WarpXPushFieldsEM.cpp index 6d588ae86..97204cc17 100644 --- a/Source/FieldSolver/WarpXPushFieldsEM.cpp +++ b/Source/FieldSolver/WarpXPushFieldsEM.cpp @@ -76,6 +76,7 @@ WarpX::EvolveB (int lev, PatchType patch_type, amrex::Real dt) const Box& tbx = mfi.tilebox(Bx_nodal_flag); const Box& tby = mfi.tilebox(By_nodal_flag); const Box& tbz = mfi.tilebox(Bz_nodal_flag); + const Real xmin = mfi.tilebox().smallEnd(0)*dx[0]; if (do_nodal) { auto const& Bxfab = Bx->array(mfi); @@ -112,6 +113,7 @@ WarpX::EvolveB (int lev, PatchType patch_type, amrex::Real dt) BL_TO_FORTRAN_3D((*By)[mfi]), BL_TO_FORTRAN_3D((*Bz)[mfi]), &dtsdx, &dtsdy, &dtsdz, + &xmin, &dx[0], &WarpX::maxwell_fdtd_solver_id); } @@ -230,6 +232,7 @@ WarpX::EvolveE (int lev, PatchType patch_type, amrex::Real dt) const Box& tex = mfi.tilebox(Ex_nodal_flag); const Box& tey = mfi.tilebox(Ey_nodal_flag); const Box& tez = mfi.tilebox(Ez_nodal_flag); + const Real xmin = mfi.tilebox().smallEnd(0)*dx[0]; if (do_nodal) { auto const& Exfab = Ex->array(mfi); @@ -272,7 +275,8 @@ WarpX::EvolveE (int lev, PatchType patch_type, amrex::Real dt) BL_TO_FORTRAN_3D((*jy)[mfi]), BL_TO_FORTRAN_3D((*jz)[mfi]), &mu_c2_dt, - &dtsdx_c2, &dtsdy_c2, &dtsdz_c2); + &dtsdx_c2, &dtsdy_c2, &dtsdz_c2, + &xmin, &dx[0]); } if (F) |