diff options
Diffstat (limited to 'Source/Diagnostics/WarpXIO.cpp')
-rw-r--r-- | Source/Diagnostics/WarpXIO.cpp | 40 |
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]); -} |