ERF_PoissonWallDist.cpp File Reference

#include "ERF.H"
#include "ERF_Constants.H"
#include "ERF_Utils.H"
#include "ERF_TerrainPoisson_3D_K.H"
#include "ERF_TerrainMetrics.H"
#include "AMReX_MLMG.H"
#include "AMReX_MLABecLaplacian.H"

Include dependency graph for ERF_PoissonWallDist.cpp:

Functions
static MultiFab	surface_slab (const MultiFab &z_nd, int klo)
static void	assert_surface_gathered (const MultiFab &znd_surf)

Function Documentation

assert_surface_gathered()

static void assert_surface_gathered ( const MultiFab &  znd_surf )

static

Stop unless every surface node, valid and x/y ghost, was gathered. The slab has no z ghosts, so the reduction runs over its own ghost vector.

Parameters

znd_surf

Surface slab filled by fill_surface_nodes

{
    Real znd_max = ReduceMax(znd_surf, znd_surf.nGrowVect(),
        [=] AMREX_GPU_HOST_DEVICE (Box const& bx, Array4<Real const> const& a) -> Real
        {
            Real m = std::numeric_limits<Real>::lowest();
            amrex::Loop(bx, [&] (int i, int j, int k) { m = amrex::max(m, a(i,j,k)); });
            return m;
        });
    ParallelAllReduce::Max(znd_max, ParallelContext::CommunicatorSub());
    // Same threshold as fill_surface_nodes: a node interpolated from a stencil that still
    // held bogus values is a fraction of bogus_large_value, not equal to it
    AMREX_ALWAYS_ASSERT_WITH_MESSAGE(znd_max < Real(0.5) * bogus_large_value,
        "poisson_wall_dist: the surface nodes were not gathered onto every box");
}

Referenced by ERF::poisson_wall_dist().

Here is the call graph for this function:

Here is the caller graph for this function:

surface_slab()

static MultiFab surface_slab ( const MultiFab &  z_nd, int  klo  )

static

Nodal slab at k = klo on the 2D footprint of every box of z_nd, so that boxes that do not touch the surface (a BoxArray split in z, a refined level that stops above the ground) still hold the surface nodes once fill_surface_nodes has run.

Parameters

z_nd	Nodal physical heights of the level
klo	Lowest k index of the level's domain

{
    BoxList bl_surf = z_nd.boxArray().boxList();
    for (auto& b : bl_surf) { b.setRange(2,klo); }
    BoxArray ba_surf(std::move(bl_surf));
    IntVect ng_surf = z_nd.nGrowVect(); ng_surf[2] = 0;
    return MultiFab(ba_surf, z_nd.DistributionMap(), 1, ng_surf);
}

Referenced by ERF::poisson_wall_dist().

Here is the caller graph for this function: