diff options
Diffstat (limited to 'Source/FieldSolver')
| -rw-r--r-- | Source/FieldSolver/WarpXFFT.cpp | 14 |
1 files changed, 14 insertions, 0 deletions
diff --git a/Source/FieldSolver/WarpXFFT.cpp b/Source/FieldSolver/WarpXFFT.cpp index e8b1442bb..ffeb237c0 100644 --- a/Source/FieldSolver/WarpXFFT.cpp +++ b/Source/FieldSolver/WarpXFFT.cpp @@ -405,6 +405,20 @@ WarpX::PushPSATD (int lev, amrex::Real /* dt */) rho_fp_fft[lev]->ParallelCopy(*rho_fp[lev], 0, 0, 2, 0, 0, period_fp); BL_PROFILE_VAR_STOP(blp_copy); + + for (MFIter mfi(*rho_fp_fft[lev]); mfi.isValid(); ++mfi) + { + MultiFab &mf = *rho_fp_fft[lev]; + Box realspace_bx = mf[mfi].box(); // Copy the box + Array4<const Real> mf_arr = mf[mfi].array(); + amrex::Print() << " at begin Push PSATD rho " << mf_arr(0,0,0,0) ; + amrex::Print() << " new rho " << mf_arr(0,0,0,1) << "\n"; + MultiFab &mf_orig = *rho_fp[lev]; + Array4<const Real> mforig_arr = mf_orig[mfi].array(); + amrex::Print() << " at begin Push PSATD rho " << mforig_arr(0,0,0,0) ; + amrex::Print() << " new rho " << mforig_arr(0,0,0,1) << "\n"; + } + BL_PROFILE_VAR_START(blp_push_eb); if (fft_hybrid_mpi_decomposition){ if (Efield_fp_fft[lev][0]->local_size() == 1) |