fillBoundary J

prkkumar · Apr 26, 2024 · f7131d6 · f7131d6
1 parent 48d96b5
commit f7131d6
Show file tree

Hide file tree

Showing 2 changed files with 2 additions and 6 deletions.
diff --git 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
@@ -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<amrex::Real, AMREX_SPACEDIM> 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) {