Function ca_statefill

Source/problems/Castro_bc_fill_nd.cpp:10–124

using namespace amrex;

void ca_statefill(Box const& bx, FArrayBox& data,
                  const int dcomp, const int numcomp,
                  Geometry const& geom, const Real time,
                  const Vector<BCRec>& bcr, const int bcomp,
                  const int scomp)
{
    // Here dcomp is the component in the destination array that we
    // are filling and bcr is a vector of length ncomp which are the
    // BC values corresponding to components dcomp to dcomp + ncomp -
    // 1

    // First, fill all the BC data using the default routines.
    // We replace inflow with outflow in the generic fill to ensure that
    // valid data is always present.

    Vector<BCRec> bcr_noinflow{bcr};
    for (auto & bc : bcr_noinflow) {
        for (int dir = 0; dir < AMREX_SPACEDIM; ++dir) {
            if (bc.lo(dir) == amrex::BCType::ext_dir) {
                bc.setLo(dir, amrex::BCType::foextrap);
            }
            if (bc.hi(dir) == amrex::BCType::ext_dir) {
                bc.setHi(dir, amrex::BCType::foextrap);
            }
        }
    }

    GpuBndryFuncFab<CastroGenericFill> gpu_bndry_func(CastroGenericFill{});
    gpu_bndry_func(bx, data, dcomp, numcomp, geom, time, bcr_noinflow, bcomp, scomp);

    // At this point, if we filling anything other than the full state (numcomp == NUM_STATE),
    // we assume that we are doing a derive or some other routine which is not actually setting
    // boundary conditions on the state data. So we immediately return. This means that the derive
    // data may not be exactly what the user wants at the physical boundary, but the impact of that
    // is usually negligible.

    if (numcomp != NUM_STATE) {
        return;
    }

    // Fill ambient BCs.

    ambient_fill(bx, data.array(dcomp), geom, bcr);

    // Now we consider external BCs (HSE).  Note, if we are at a
    // corner where two (or three) faces want to do HSE, we may run
    // into a situation that the data is not valid in the corner where
    // we start the integration.  We'll abort, for now, if we run into
    // this case.
    //
    // The future fix is to first call ext_fill on the ghost cells
    // that are not corners and then call it a second time on just the
    // corners.

#if AMREX_SPACEDIM == 2
    if ((bcr[URHO].lo(0) == amrex::BCType::ext_dir && bcr[URHO].lo(1) == amrex::BCType::ext_dir) ||
        (bcr[URHO].lo(0) == amrex::BCType::ext_dir && bcr[URHO].hi(1) == amrex::BCType::ext_dir) ||
        (bcr[URHO].hi(0) == amrex::BCType::ext_dir && bcr[URHO].lo(1) == amrex::BCType::ext_dir) ||