aboutsummaryrefslogtreecommitdiff
path: root/Source/Diagnostics/WarpXIO.cpp
diff options
context:
space:
mode:
Diffstat (limited to 'Source/Diagnostics/WarpXIO.cpp')
-rw-r--r--Source/Diagnostics/WarpXIO.cpp40
1 files changed, 0 insertions, 40 deletions
diff --git a/Source/Diagnostics/WarpXIO.cpp b/Source/Diagnostics/WarpXIO.cpp
index 548c95cd8..1785a04c1 100644
--- a/Source/Diagnostics/WarpXIO.cpp
+++ b/Source/Diagnostics/WarpXIO.cpp
@@ -352,43 +352,3 @@ WarpX::InitFromCheckpoint ()
mypc->Restart(restart_chkfile);
}
-
-
-std::unique_ptr<MultiFab>
-WarpX::GetCellCenteredData() {
-
- WARPX_PROFILE("WarpX::GetCellCenteredData()");
-
- const amrex::IntVect ng(1);
- const int nc = 10;
-
- Vector<std::unique_ptr<MultiFab> > cc(finest_level+1);
-
- for (int lev = 0; lev <= finest_level; ++lev)
- {
- cc[lev] = std::make_unique<MultiFab>(grids[lev], dmap[lev], nc, ng );
-
- int dcomp = 0;
- // first the electric field
- AverageAndPackVectorField( *cc[lev], Efield_aux[lev], dmap[lev], dcomp, ng );
- dcomp += 3;
- // then the magnetic field
- AverageAndPackVectorField( *cc[lev], Bfield_aux[lev], dmap[lev], dcomp, ng );
- dcomp += 3;
- // then the current density
- AverageAndPackVectorField( *cc[lev], current_fp[lev], dmap[lev], dcomp, ng );
- dcomp += 3;
- // then the charge density
- const std::unique_ptr<MultiFab>& charge_density = mypc->GetChargeDensity(lev);
- AverageAndPackScalarField( *cc[lev], *charge_density, dmap[lev], dcomp, ng );
-
- ablastr::utils::communication::FillBoundary(*cc[lev], WarpX::do_single_precision_comms, geom[lev].periodicity());
- }
-
- for (int lev = finest_level; lev > 0; --lev)
- {
- CoarsenIO::Coarsen( *cc[lev-1], *cc[lev], 0, 0, nc, 0, refRatio(lev-1) );
- }
-
- return std::move(cc[0]);
-}