diff options
Diffstat (limited to 'Source/Evolve/WarpXEvolveEM.cpp')
| -rw-r--r-- | Source/Evolve/WarpXEvolveEM.cpp | 22 |
1 files changed, 1 insertions, 21 deletions
diff --git a/Source/Evolve/WarpXEvolveEM.cpp b/Source/Evolve/WarpXEvolveEM.cpp index 531cd6a56..49147338b 100644 --- a/Source/Evolve/WarpXEvolveEM.cpp +++ b/Source/Evolve/WarpXEvolveEM.cpp @@ -302,11 +302,7 @@ WarpX::OneStep_nosub (Real cur_time) FillBoundaryE(); FillBoundaryB(); #else - // amrex::Print()<< "WarpXEvolveEM.cpp : before CopyJinPMLs "<<std::endl; if (do_pml && pml_has_particles){ - // do current deposition in PMLs - // copy current computed on mother grid to PMLs - // amrex::Print()<< "WarpXEvolveEM.cpp : IN CopyJinPMLs "<<std::endl; for (int lev = 0; lev <= finest_level; ++lev) { if (pml[lev]->ok()){ @@ -321,14 +317,7 @@ WarpX::OneStep_nosub (Real cur_time) } if (do_pml && do_pml_j_damping){ - // amrex::Print()<< "WarpXEvolveEM.cpp : DampJ "<<std::endl; - // damp current in pmls - // amrex::Print() << "===== DAMPING IN PMLs =====" << std::endl; - - // amrex::Print()<< "===== DAMPING J ====="<< std::endl; - // amrex::Print()<< "[AV DAMP] max_Jx_pml = "<< pml[0]->Getj_fp()[0]->min(0) << std::endl; DampJPML(); - // amrex::Print()<< "[AP DAMP] max_Jx_pml = "<< pml[0]->Getj_fp()[0]->min(0) << std::endl; } EvolveF(0.5*dt[0], DtType::FirstHalf); @@ -340,32 +329,23 @@ WarpX::OneStep_nosub (Real cur_time) FillBoundaryE(); EvolveF(0.5*dt[0], DtType::SecondHalf); - // amrex::Print()<< "WarpXEvolveEM.cpp : before CopyJinReg "<<std::endl; + if (do_pml && pml_has_particles){ - // do current deposition in PMLs - // copy current computed on mother grid to PMLs - // amrex::Print()<< "WarpXEvolveEM.cpp : IN CopyJinReg "<<std::endl; for (int lev = 0; lev <= finest_level; ++lev) { if (pml[lev]->ok()){ - // amrex::Print()<< "[AV COPY] max_Jx = "<< current_fp[lev][0].get()->min(0) << std::endl; - // amrex::Print()<< "[AV COPY] max_Jx_pml = "<< pml[lev]->Getj_fp()[0]->min(0) << std::endl; - // amrex::Print()<< "===== Copy J from PML to Reg ====="<< std::endl; pml[lev]->CopyJinReg({ current_fp[lev][0].get(), current_fp[lev][1].get(), current_fp[lev][2].get() }, { current_cp[lev][0].get(), current_cp[lev][1].get(), current_cp[lev][2].get() }); - // amrex::Print()<< "[AP COPY] max_Jx = "<< current_fp[lev][0].get()->min(0) << std::endl; - } } } EvolveB(0.5*dt[0]); // We now have B^{n+1} if (do_pml) { - // amrex::Print()<< "WarpXEvolveEM.cpp : Damp "<<std::endl; DampPML(); FillBoundaryE(); } |