compute sigma and sigma_star

1 files changed, 26 insertions, 3 deletions

diff --git a/Source/WarpXEvolve.cpp b/Source/WarpXEvolve.cpp
index f73d1f32b..0530c159d 100644
--- a/Source/WarpXEvolve.cpp
+++ b/Source/WarpXEvolve.cpp
@@ -158,9 +158,9 @@ WarpX::EvolveB (int lev, Real dt)
 #endif
     for ( MFIter mfi(*Bfield[lev][0],true); mfi.isValid(); ++mfi )
     {
-      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 Box& tbx  = mfi.tilebox(Bx_nodal_flag);
+        const Box& tby  = mfi.tilebox(By_nodal_flag);
+        const Box& tbz  = mfi.tilebox(Bz_nodal_flag);
 
 	// Call picsar routine for each tile
 	WRPX_PXR_PUSH_BVEC(
@@ -176,6 +176,29 @@ WarpX::EvolveB (int lev, Real dt)
 	    &dtsdx[0], &dtsdx[1], &dtsdx[2],
 	    &norder);
     }
+
+    if (do_pml && lev == 0)
+    {
+#if 0
+        const Geometry& gm = Geom(0);
+        const Box& domain = gm.Domain();
+        const int* dlo = domain.loVect();
+        const int* dhi = domain.hiVect();
+
+#ifdef _OPENMP
+#pragma omp parallel
+#endif
+        for ( MFIter mfi(*pmlB[0]); mfi.isValid(); ++mfi )
+        {
+            const Box& ccbx = amrex::enclosedCells(mfi.validbox());
+            const int* cclo = ccbx.loVect();
+            const int* cchi = ccbx.hiVect();
+            Box Bx_bx = amrex::convert(ccbx, Bx_nodal_flag); 
+            Box By_bx = amrex::convert(ccbx, By_nodal_flag); 
+            Box Bz_bx = amrex::convert(ccbx, Bz_nodal_flag);
+        }
+#endif
+    }
 }
 
 void