From f7131d6d8c6b576531de0796214ef715a7716b16 Mon Sep 17 00:00:00 2001 From: Prabhat Kumar Date: Fri, 26 Apr 2024 12:41:40 -0700 Subject: [PATCH] fillBoundary J --- Source/Evolve/WarpXEvolve.cpp | 5 ++--- Source/FieldSolver/FerroE/FerroE.cpp | 3 --- 2 files changed, 2 insertions(+), 6 deletions(-) diff --git a/Source/Evolve/WarpXEvolve.cpp b/Source/Evolve/WarpXEvolve.cpp index 153be1c03..31241551c 100644 --- a/Source/Evolve/WarpXEvolve.cpp +++ b/Source/Evolve/WarpXEvolve.cpp @@ -156,7 +156,7 @@ WarpX::Evolve (int numsteps) m_ferroe->EvolveP(-0.5_rt*dt[0]); // P^(n) to P^(n-1/2) using E^(n) FillBoundaryP(guard_cells.ng_alloc_EB); m_ferroe->EvolveFerroEJ(-0.5_rt*dt[0]); // J^(n) to J^(n-1/2) using E^(n) - //FillBoundaryJ(guard_cells.ng_alloc_EB); + FillBoundaryJ(guard_cells.ng_alloc_EB); } #endif is_synchronized = false; @@ -459,8 +459,7 @@ WarpX::OneStep_nosub (Real cur_time) m_ferroe->EvolveP(dt[0]); // P^(n-1/2) to P^(n+1/2) using E^(n) FillBoundaryP(guard_cells.ng_alloc_EB); m_ferroe->EvolveFerroEJ(dt[0]); // J^(n-1/2) to J^(n+1/2) using E^(n) - //EvolveB(0.5_rt * dt[0], DtType::FirstHalf); // We now have B^{n+1/2} - //FillBoundaryJ(guard_cells.ng_alloc_EB); + FillBoundaryJ(guard_cells.ng_alloc_EB); // fill boundary here } #endif diff --git a/Source/FieldSolver/FerroE/FerroE.cpp b/Source/FieldSolver/FerroE/FerroE.cpp index 2efcb85cf..69911bf54 100644 --- a/Source/FieldSolver/FerroE/FerroE.cpp +++ b/Source/FieldSolver/FerroE/FerroE.cpp @@ -191,8 +191,6 @@ FerroE::EvolveP (amrex::Real dt) { amrex::Print() << " evolve P \n"; auto & warpx = WarpX::GetInstance(); -// int include_Landau = ferroe.include_Landau; -// int include_grad = ferroe.include_grad; const int lev = 0; const amrex::GpuArray dx = warpx.Geom(lev).CellSizeArray(); @@ -218,7 +216,6 @@ FerroE::EvolveP (amrex::Real dt) amrex::Box const& tpy = mfi.tilebox(Py->ixType().toIntVect()); amrex::Box const& tpz = mfi.tilebox(Pz->ixType().toIntVect()); - amrex::Print() << "include_Landau = " << include_Landau << "\n"; amrex::ParallelFor(tpx, tpy, tpz, [=] AMREX_GPU_DEVICE (int i, int j, int k) { if (fe_arr(i,j,k)==1 and fe_arr(i+1,j,k)==1) {