ERF_ShocDriver.cpp File Reference

#include "ERF_ShocDriver.H"
#include "ERF_ShocImplicit.H"
#include "ERF_Constants.H"
#include "ERF_IndexDefines.H"
#include <AMReX_BLProfiler.H>
#include <AMReX_Gpu.H>
#include <algorithm>
#include <cstring>
#include <iomanip>
#include <limits>
#include <sstream>
#include <vector>

Include dependency graph for ERF_ShocDriver.cpp:

Functions
bool	shoc_boxarray_spans_full_height (const BoxArray &ba, const Box &domain)
void	shoc_fill_physical_boundary_ghosts (MultiFab &mf, const Geometry &geom, int comp, int ncomp)

Function Documentation

shoc_boxarray_spans_full_height()

bool shoc_boxarray_spans_full_height	(	const BoxArray &	ba,
		const Box &	domain
	)

{
    const int dom_klo = domain.smallEnd(2);
    const int dom_khi = domain.bigEnd(2);
    for (int ibox = 0; ibox < ba.size(); ++ibox) {
        const Box& bx = ba[ibox];
        if (bx.smallEnd(2) != dom_klo || bx.bigEnd(2) != dom_khi) {
            return false;
        }
    }
    return true;
}

shoc_fill_physical_boundary_ghosts()

void shoc_fill_physical_boundary_ghosts	(	MultiFab &	mf,
		const Geometry &	geom,
		int	comp,
		int	ncomp
	)

{
    mf.FillBoundary(comp, ncomp, geom.periodicity());
 
    const Box domain = geom.Domain();
    const bool nonperiodic_x = !geom.isPeriodic(0);
    const bool nonperiodic_y = !geom.isPeriodic(1);
    const bool nonperiodic_z = !geom.isPeriodic(2);
 
    for (MFIter mfi(mf, false); mfi.isValid(); ++mfi) {
        const Box gbx = mfi.growntilebox();
        auto data = mf.array(mfi);
        ParallelFor(gbx, ncomp,
                    [=] AMREX_GPU_DEVICE (int i, int j, int k, int n) noexcept {
            const int ii = nonperiodic_x
                ? ((i < domain.smallEnd(0)) ? domain.smallEnd(0)
                   : ((i > domain.bigEnd(0)) ? domain.bigEnd(0) : i))
                : i;
            const int jj = nonperiodic_y
                ? ((j < domain.smallEnd(1)) ? domain.smallEnd(1)
                   : ((j > domain.bigEnd(1)) ? domain.bigEnd(1) : j))
                : j;
            const int kk = nonperiodic_z
                ? ((k < domain.smallEnd(2)) ? domain.smallEnd(2)
                   : ((k > domain.bigEnd(2)) ? domain.bigEnd(2) : k))
                : k;
 
            if (ii != i || jj != j || kk != k) {
                data(i,j,k,comp+n) = data(ii,jj,kk,comp+n);
            }
        });
    }
}

Referenced by ShocDriver::set_eddy_diffs().

Here is the call graph for this function:

Here is the caller graph for this function: