#include <memory>
#include "ERF_Constants.H"
#include "AMReX_buildInfo.H"
#include "ERF.H"
#include "ERF_Utils.H"
#include "ERF_ProbCommon.H"
#include "ERF_SBMStateManager.H"
#include "ERF_SBMTransport.H"

Include dependency graph for ERF_MakeNewLevel.cpp:

Functions
static Vector< std::unique_ptr< MultiFab > >	save_qmoist_restart_vars (const Vector< MultiFab * > &qmoist_lev, const std::vector< int > &idx)

static void	fill_qmoist_restart_var (MultiFab &dst, const MultiFab old_mf, const MultiFab crse_mf, const Geometry &crse_geom, const Geometry &fine_geom, const IntVect &ratio)

Function Documentation

◆ fill_qmoist_restart_var()

static void fill_qmoist_restart_var	(	MultiFab &	dst,
		const MultiFab *	old_mf,
		const MultiFab *	crse_mf,
		const Geometry &	crse_geom,
		const Geometry &	fine_geom,
		const IntVect &	ratio
	)

static

 {
     const int ncomp = dst.nComp();
     if (crse_mf != nullptr) {
         MultiFab tmp(dst.boxArray(), dst.DistributionMap(), ncomp, 0);
         Vector<BCRec> bcs(ncomp); // unused by pc_interp, which needs no coarse ghost cells
         for (int n = 0; n < ncomp; ++n) {
             for (int dir = 0; dir < AMREX_SPACEDIM; ++dir) {
                 bcs[n].setLo(dir, BCType::foextrap);
                 bcs[n].setHi(dir, BCType::foextrap);
             }
         }
         InterpFromCoarseLevel(tmp, IntVect(0), IntVect(0), *crse_mf, 0, 0, ncomp,
                               crse_geom, fine_geom, ratio, &pc_interp, bcs, 0);
         const int klo = fine_geom.Domain().smallEnd(2);
         for (MFIter mfi(dst, TilingIfNotGPU()); mfi.isValid(); ++mfi) {
             Box bx = mfi.tilebox();
             if (bx.smallEnd(2) > klo) { continue; }
             bx.setBig(2, klo);
             const auto& dst_arr = dst.array(mfi);
             const auto& tmp_arr = tmp.const_array(mfi);
             ParallelFor(bx, ncomp, [=] AMREX_GPU_DEVICE (int i, int j, int k, int n) noexcept
             {
                 dst_arr(i,j,k,n) = tmp_arr(i,j,k,n);
             });
         }
     }
     if (old_mf != nullptr) { dst.ParallelCopy(*old_mf, 0, 0, ncomp); }
     dst.FillBoundary(fine_geom.periodicity());
 }

Referenced by ERF::MakeNewLevelFromCoarse(), and ERF::RemakeLevel().

Here is the call graph for this function:

Here is the caller graph for this function:

◆ save_qmoist_restart_vars()

static Vector<std::unique_ptr<MultiFab> > save_qmoist_restart_vars	(	const Vector< MultiFab * > &	qmoist_lev,
		const std::vector< int > &	idx
	)