stencils.H

Go to the documentation of this file.

#ifndef STENCILS_H
#define STENCILS_H

 
#include "AMReX_REAL.H"
#include "AMReX_Box.H"
#include "AMReX_Orientation.H"
#include "AMReX_Geometry.H"
 
using namespace amrex::literals;
 
namespace kynema_sgf::fvm::stencil {
 
namespace impl {

AMREX_INLINE amrex::Box
box_lo(const amrex::Box& bx, const amrex::Geometry& geom, const int idir)
{
    if (geom.isPeriodic(idir)) {
        return {};
    }
    if (bx.smallEnd(idir) != geom.Domain().smallEnd(idir)) {
        return {};
    }
 
    amrex::IntVect low(bx.smallEnd());
    amrex::IntVect hi(bx.bigEnd());
    int sm = low[idir];
    low.setVal(idir, sm);
    hi.setVal(idir, sm);
    return {low, hi};
}

AMREX_INLINE amrex::Box
box_hi(const amrex::Box& bx, const amrex::Geometry& geom, const int idir)
{
    if (geom.isPeriodic(idir)) {
        return {};
    }
    if (bx.bigEnd(idir) != geom.Domain().bigEnd(idir)) {
        return {};
    }
 
    amrex::IntVect low(bx.smallEnd());
    amrex::IntVect hi(bx.bigEnd());
    int sm = hi[idir];
    low.setVal(idir, sm);
    hi.setVal(idir, sm);
    return {low, hi};
}
} // namespace impl

struct StencilInterior
{
    // First derivatives
    static constexpr amrex::Real c00 = 0.5_rt;
    static constexpr amrex::Real c01 = 0.0_rt;
    static constexpr amrex::Real c02 = -0.5_rt;
    static constexpr amrex::Real c10 = c00;
    static constexpr amrex::Real c11 = c01;
    static constexpr amrex::Real c12 = c02;
    static constexpr amrex::Real c20 = c00;
    static constexpr amrex::Real c21 = c01;
    static constexpr amrex::Real c22 = c02;
 
    // Second derivatives
    static constexpr amrex::Real s00 = 1.0_rt;
    static constexpr amrex::Real s01 = -2.0_rt;
    static constexpr amrex::Real s02 = 1.0_rt;
    static constexpr amrex::Real s10 = s00;
    static constexpr amrex::Real s11 = s01;
    static constexpr amrex::Real s12 = s02;
    static constexpr amrex::Real s20 = s00;
    static constexpr amrex::Real s21 = s01;
    static constexpr amrex::Real s22 = s02;
 
    // Filter
    static constexpr amrex::Real f00 = 0.25_rt;
    static constexpr amrex::Real f01 = 0.5_rt;
    static constexpr amrex::Real f02 = 0.25_rt;
    static constexpr amrex::Real f10 = f00;
    static constexpr amrex::Real f11 = f01;
    static constexpr amrex::Real f12 = f02;
    static constexpr amrex::Real f20 = f00;
    static constexpr amrex::Real f21 = f01;
    static constexpr amrex::Real f22 = f02;
 
    AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE static bool applies_to(
        int /*i*/,
        int /*j*/,
        int /*k*/,
        amrex::IntVect const& /*dlo*/,
        amrex::IntVect const& /*dhi*/) noexcept
    {
        return true;
    }
 
    static amrex::Box
    box(const amrex::Box& bx, const amrex::Geometry& /*unused*/)
    {
        return bx;
    }
};
 
struct StencilILO
{
    static constexpr int direction = 0;
    static constexpr int side = amrex::Orientation::Side::low;
 
    static constexpr amrex::Real c00 = 1.0_rt / 3.0_rt;
    static constexpr amrex::Real c01 = 1.0_rt;
    static constexpr amrex::Real c02 = -4.0_rt / 3.0_rt;
    static constexpr amrex::Real c10 = 0.5_rt;
    static constexpr amrex::Real c11 = 0.0_rt;
    static constexpr amrex::Real c12 = -0.5_rt;
    static constexpr amrex::Real c20 = c10;
    static constexpr amrex::Real c21 = c11;
    static constexpr amrex::Real c22 = c12;
 
    static constexpr amrex::Real s00 = 4.0_rt / 3.0_rt; // axx
    static constexpr amrex::Real s01 = -4.0_rt;         // bxx
    static constexpr amrex::Real s02 = 8.0_rt / 3.0_rt; // cxx
    static constexpr amrex::Real s10 = 1.0_rt;          // ayy
    static constexpr amrex::Real s11 = -2.0_rt;         // byy
    static constexpr amrex::Real s12 = 1.0_rt;          // cyy
    static constexpr amrex::Real s20 = s10;             // azz
    static constexpr amrex::Real s21 = s11;             // bzz
    static constexpr amrex::Real s22 = s12;             // czz
 
    static constexpr amrex::Real f00 = 0.125_rt;
    static constexpr amrex::Real f01 = 0.625_rt;
    static constexpr amrex::Real f02 = 0.25_rt;
    static constexpr amrex::Real f10 = 0.25_rt;
    static constexpr amrex::Real f11 = 0.5_rt;
    static constexpr amrex::Real f12 = 0.25_rt;
    static constexpr amrex::Real f20 = f10;
    static constexpr amrex::Real f21 = f11;
    static constexpr amrex::Real f22 = f12;
 
    AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE static bool applies_to(
        int i,
        int /*j*/,
        int /*k*/,
        amrex::IntVect const& dlo,
        amrex::IntVect const& /*dhi*/) noexcept
    {
        return i == dlo[0];
    }
 
    static amrex::Box box(const amrex::Box& bx, const amrex::Geometry& geom)
    {
        return impl::box_lo(bx, geom, StencilILO::direction);
    }
};
 
struct StencilJLO
{
    static constexpr int direction = 1;
    static constexpr int side = amrex::Orientation::Side::low;
 
    static constexpr amrex::Real c00 = 0.5_rt;
    static constexpr amrex::Real c01 = 0.0_rt;
    static constexpr amrex::Real c02 = -0.5_rt;
    static constexpr amrex::Real c10 = 1.0_rt / 3.0_rt;
    static constexpr amrex::Real c11 = 1.0_rt;
    static constexpr amrex::Real c12 = -4.0_rt / 3.0_rt;
    static constexpr amrex::Real c20 = c00;
    static constexpr amrex::Real c21 = c01;
    static constexpr amrex::Real c22 = c02;
 
    static constexpr amrex::Real s00 = 1.0_rt;          // axx
    static constexpr amrex::Real s01 = -2.0_rt;         // bxx
    static constexpr amrex::Real s02 = 1.0_rt;          // cxx
    static constexpr amrex::Real s10 = 4.0_rt / 3.0_rt; // ayy
    static constexpr amrex::Real s11 = -4.0_rt;         // byy
    static constexpr amrex::Real s12 = 8.0_rt / 3.0_rt; // cyy
    static constexpr amrex::Real s20 = s00;             // azz
    static constexpr amrex::Real s21 = s01;             // bzz
    static constexpr amrex::Real s22 = s02;             // czz
 
    static constexpr amrex::Real f00 = 0.25_rt;
    static constexpr amrex::Real f01 = 0.5_rt;
    static constexpr amrex::Real f02 = 0.25_rt;
    static constexpr amrex::Real f10 = 0.125_rt;
    static constexpr amrex::Real f11 = 0.625_rt;
    static constexpr amrex::Real f12 = 0.25_rt;
    static constexpr amrex::Real f20 = f00;
    static constexpr amrex::Real f21 = f01;
    static constexpr amrex::Real f22 = f02;
 
    AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE static bool applies_to(
        int /*i*/,
        int j,
        int /*k*/,
        amrex::IntVect const& dlo,
        amrex::IntVect const& /*dhi*/) noexcept
    {
        return j == dlo[1];
    }
 
    static amrex::Box box(const amrex::Box& bx, const amrex::Geometry& geom)
    {
        return impl::box_lo(bx, geom, StencilJLO::direction);
    }
};
 
struct StencilKLO
{
    static constexpr int direction = 2;
    static constexpr int side = amrex::Orientation::Side::low;
 
    static constexpr amrex::Real c00 = 0.5_rt;
    static constexpr amrex::Real c01 = 0.0_rt;
    static constexpr amrex::Real c02 = -0.5_rt;
    static constexpr amrex::Real c10 = c00;
    static constexpr amrex::Real c11 = c01;
    static constexpr amrex::Real c12 = c02;
    static constexpr amrex::Real c20 = 1.0_rt / 3.0_rt;
    static constexpr amrex::Real c21 = 1.0_rt;
    static constexpr amrex::Real c22 = -4.0_rt / 3.0_rt;
 
    static constexpr amrex::Real s00 = 1.0_rt;          // axx
    static constexpr amrex::Real s01 = -2.0_rt;         // bxx
    static constexpr amrex::Real s02 = 1.0_rt;          // cxx
    static constexpr amrex::Real s10 = s00;             // ayy
    static constexpr amrex::Real s11 = s01;             // byy
    static constexpr amrex::Real s12 = s02;             // cyy
    static constexpr amrex::Real s20 = 4.0_rt / 3.0_rt; // azz
    static constexpr amrex::Real s21 = -4.0_rt;         // bzz
    static constexpr amrex::Real s22 = 8.0_rt / 3.0_rt; // czz
 
    static constexpr amrex::Real f00 = 0.25_rt;
    static constexpr amrex::Real f01 = 0.5_rt;
    static constexpr amrex::Real f02 = 0.25_rt;
    static constexpr amrex::Real f10 = f00;
    static constexpr amrex::Real f11 = f01;
    static constexpr amrex::Real f12 = f02;
    static constexpr amrex::Real f20 = 0.125_rt;
    static constexpr amrex::Real f21 = 0.625_rt;
    static constexpr amrex::Real f22 = 0.25_rt;
 
    AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE static bool applies_to(
        int /*i*/,
        int /*j*/,
        int k,
        amrex::IntVect const& dlo,
        amrex::IntVect const& /*dhi*/) noexcept
    {
        return k == dlo[2];
    }
 
    static amrex::Box box(const amrex::Box& bx, const amrex::Geometry& geom)
    {
        return impl::box_lo(bx, geom, StencilKLO::direction);
    }
};
 
struct StencilIHI
{
    static constexpr int direction = 0;
    static constexpr int side = amrex::Orientation::Side::high;
 
    static constexpr amrex::Real c00 = 4.0_rt / 3.0_rt;
    static constexpr amrex::Real c01 = -1.0_rt;
    static constexpr amrex::Real c02 = -1.0_rt / 3.0_rt;
    static constexpr amrex::Real c10 = 0.5_rt;
    static constexpr amrex::Real c11 = 0.0_rt;
    static constexpr amrex::Real c12 = -0.5_rt;
    static constexpr amrex::Real c20 = c10;
    static constexpr amrex::Real c21 = c11;
    static constexpr amrex::Real c22 = c12;
 
    static constexpr amrex::Real s00 = 8.0_rt / 3.0_rt; // axx
    static constexpr amrex::Real s01 = -4.0_rt;         // bxx
    static constexpr amrex::Real s02 = 4.0_rt / 3.0_rt; // cxx
    static constexpr amrex::Real s10 = 1.0_rt;          // ayy
    static constexpr amrex::Real s11 = -2.0_rt;         // byy
    static constexpr amrex::Real s12 = 1.0_rt;          // cyy
    static constexpr amrex::Real s20 = s10;             // azz
    static constexpr amrex::Real s21 = s11;             // bzz
    static constexpr amrex::Real s22 = s12;             // czz
 
    static constexpr amrex::Real f00 = 0.25_rt;
    static constexpr amrex::Real f01 = 0.625_rt;
    static constexpr amrex::Real f02 = 0.125_rt;
    static constexpr amrex::Real f10 = 0.25_rt;
    static constexpr amrex::Real f11 = 0.5_rt;
    static constexpr amrex::Real f12 = 0.25_rt;
    static constexpr amrex::Real f20 = f10;
    static constexpr amrex::Real f21 = f11;
    static constexpr amrex::Real f22 = f12;
 
    AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE static bool applies_to(
        int i,
        int /*j*/,
        int /*k*/,
        amrex::IntVect const& /*dlo*/,
        amrex::IntVect const& dhi) noexcept
    {
        return i == dhi[0];
    }
 
    static amrex::Box box(const amrex::Box& bx, const amrex::Geometry& geom)
    {
        return impl::box_hi(bx, geom, StencilIHI::direction);
    }
};
 
struct StencilJHI
{
    static constexpr int direction = 1;
    static constexpr int side = amrex::Orientation::Side::high;
 
    static constexpr amrex::Real c00 = 0.5_rt;
    static constexpr amrex::Real c01 = 0.0_rt;
    static constexpr amrex::Real c02 = -0.5_rt;
    static constexpr amrex::Real c10 = 4.0_rt / 3.0_rt;
    static constexpr amrex::Real c11 = -1.0_rt;
    static constexpr amrex::Real c12 = -1.0_rt / 3.0_rt;
    static constexpr amrex::Real c20 = c00;
    static constexpr amrex::Real c21 = c01;
    static constexpr amrex::Real c22 = c02;
 
    static constexpr amrex::Real s00 = 1.0_rt;          // axx
    static constexpr amrex::Real s01 = -2.0_rt;         // bxx
    static constexpr amrex::Real s02 = 1.0_rt;          // cxx
    static constexpr amrex::Real s10 = 8.0_rt / 3.0_rt; // ayy
    static constexpr amrex::Real s11 = -4.0_rt;         // byy
    static constexpr amrex::Real s12 = 4.0_rt / 3.0_rt; // cyy
    static constexpr amrex::Real s20 = s00;             // azz
    static constexpr amrex::Real s21 = s01;             // bzz
    static constexpr amrex::Real s22 = s02;             // czz
 
    static constexpr amrex::Real f00 = 0.25_rt;
    static constexpr amrex::Real f01 = 0.5_rt;
    static constexpr amrex::Real f02 = 0.25_rt;
    static constexpr amrex::Real f10 = 0.25_rt;
    static constexpr amrex::Real f11 = 0.625_rt;
    static constexpr amrex::Real f12 = 0.125_rt;
    static constexpr amrex::Real f20 = f00;
    static constexpr amrex::Real f21 = f01;
    static constexpr amrex::Real f22 = f02;
 
    AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE static bool applies_to(
        int /*i*/,
        int j,
        int /*k*/,
        amrex::IntVect const& /*dlo*/,
        amrex::IntVect const& dhi) noexcept
    {
        return j == dhi[1];
    }
 
    static amrex::Box box(const amrex::Box& bx, const amrex::Geometry& geom)
    {
        return impl::box_hi(bx, geom, StencilJHI::direction);
    }
};
 
struct StencilKHI
{
    static constexpr int direction = 2;
    static constexpr int side = amrex::Orientation::Side::high;
 
    static constexpr amrex::Real c00 = 0.5_rt;
    static constexpr amrex::Real c01 = 0.0_rt;
    static constexpr amrex::Real c02 = -0.5_rt;
    static constexpr amrex::Real c10 = c00;
    static constexpr amrex::Real c11 = c01;
    static constexpr amrex::Real c12 = c02;
    static constexpr amrex::Real c20 = 4.0_rt / 3.0_rt;
    static constexpr amrex::Real c21 = -1.0_rt;
    static constexpr amrex::Real c22 = -1.0_rt / 3.0_rt;
 
    static constexpr amrex::Real s00 = 1.0_rt;          // axx
    static constexpr amrex::Real s01 = -2.0_rt;         // bxx
    static constexpr amrex::Real s02 = 1.0_rt;          // cxx
    static constexpr amrex::Real s10 = s00;             // ayy
    static constexpr amrex::Real s11 = s01;             // byy
    static constexpr amrex::Real s12 = s02;             // cyy
    static constexpr amrex::Real s20 = 8.0_rt / 3.0_rt; // azz
    static constexpr amrex::Real s21 = -4.0_rt;         // bzz
    static constexpr amrex::Real s22 = 4.0_rt / 3.0_rt; // czz
 
    static constexpr amrex::Real f00 = 0.25_rt;
    static constexpr amrex::Real f01 = 0.5_rt;
    static constexpr amrex::Real f02 = 0.25_rt;
    static constexpr amrex::Real f10 = f00;
    static constexpr amrex::Real f11 = f01;
    static constexpr amrex::Real f12 = f02;
    static constexpr amrex::Real f20 = 0.25_rt;
    static constexpr amrex::Real f21 = 0.625_rt;
    static constexpr amrex::Real f22 = 0.125_rt;
 
    AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE static bool applies_to(
        int /*i*/,
        int /*j*/,
        int k,
        amrex::IntVect const& /*dlo*/,
        amrex::IntVect const& dhi) noexcept
    {
        return k == dhi[2];
    }
 
    static amrex::Box box(const amrex::Box& bx, const amrex::Geometry& geom)
    {
        return impl::box_hi(bx, geom, StencilKHI::direction);
    }
};
 
struct StencilIloJlo
{
    static constexpr int direction = 10;
 
    static constexpr amrex::Real c00 = 1.0_rt / 3.0_rt;
    static constexpr amrex::Real c01 = 1.0_rt;
    static constexpr amrex::Real c02 = -4.0_rt / 3.0_rt;
    static constexpr amrex::Real c10 = c00;
    static constexpr amrex::Real c11 = c01;
    static constexpr amrex::Real c12 = c02;
    static constexpr amrex::Real c20 = 0.5_rt;
    static constexpr amrex::Real c21 = 0.0_rt;
    static constexpr amrex::Real c22 = -0.5_rt;
 
    static constexpr amrex::Real s00 = 4.0_rt / 3.0_rt; // axx
    static constexpr amrex::Real s01 = -4.0_rt;         // bxx
    static constexpr amrex::Real s02 = 8.0_rt / 3.0_rt; // cxx
    static constexpr amrex::Real s10 = s00;             // ayy
    static constexpr amrex::Real s11 = s01;             // byy
    static constexpr amrex::Real s12 = s02;             // cyy
    static constexpr amrex::Real s20 = 1.0_rt;          // azz
    static constexpr amrex::Real s21 = -2.0_rt;         // bzz
    static constexpr amrex::Real s22 = 1.0_rt;          // czz
 
    static constexpr amrex::Real f00 = 0.125_rt;
    static constexpr amrex::Real f01 = 0.625_rt;
    static constexpr amrex::Real f02 = 0.25_rt;
    static constexpr amrex::Real f10 = f00;
    static constexpr amrex::Real f11 = f01;
    static constexpr amrex::Real f12 = f02;
    static constexpr amrex::Real f20 = 0.25_rt;
    static constexpr amrex::Real f21 = 0.5_rt;
    static constexpr amrex::Real f22 = 0.25_rt;
 
    AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE static bool applies_to(
        int i,
        int j,
        int /*k*/,
        amrex::IntVect const& dlo,
        amrex::IntVect const& /*dhi*/) noexcept
    {
        return i == dlo[0] && j == dlo[1];
    }
 
    static amrex::Box box(const amrex::Box& bx, const amrex::Geometry& geom)
    {
        return StencilILO::box(bx, geom) & StencilJLO::box(bx, geom);
    }
};
 
struct StencilJloKlo
{
    static constexpr int direction = 11;
 
    static constexpr amrex::Real c00 = 0.5_rt;
    static constexpr amrex::Real c01 = 0.0_rt;
    static constexpr amrex::Real c02 = -0.5_rt;
    static constexpr amrex::Real c10 = 1.0_rt / 3.0_rt;
    static constexpr amrex::Real c11 = 1.0_rt;
    static constexpr amrex::Real c12 = -4.0_rt / 3.0_rt;
    static constexpr amrex::Real c20 = c10;
    static constexpr amrex::Real c21 = c11;
    static constexpr amrex::Real c22 = c12;
 
    static constexpr amrex::Real s00 = 1.0_rt;          // axx
    static constexpr amrex::Real s01 = -2.0_rt;         // bxx
    static constexpr amrex::Real s02 = 1.0_rt;          // cxx
    static constexpr amrex::Real s10 = 4.0_rt / 3.0_rt; // ayy
    static constexpr amrex::Real s11 = -4.0_rt;         // byy
    static constexpr amrex::Real s12 = 8.0_rt / 3.0_rt; // cyy
    static constexpr amrex::Real s20 = s10;             // azz
    static constexpr amrex::Real s21 = s11;             // bzz
    static constexpr amrex::Real s22 = s12;             // czz
 
    static constexpr amrex::Real f00 = 0.25_rt;
    static constexpr amrex::Real f01 = 0.5_rt;
    static constexpr amrex::Real f02 = 0.25_rt;
    static constexpr amrex::Real f10 = 0.125_rt;
    static constexpr amrex::Real f11 = 0.625_rt;
    static constexpr amrex::Real f12 = 0.25_rt;
    static constexpr amrex::Real f20 = f10;
    static constexpr amrex::Real f21 = f11;
    static constexpr amrex::Real f22 = f12;
 
    AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE static bool applies_to(
        int /*i*/,
        int j,
        int k,
        amrex::IntVect const& dlo,
        amrex::IntVect const& /*dhi*/) noexcept
    {
        return j == dlo[1] && k == dlo[2];
    }
 
    static amrex::Box box(const amrex::Box& bx, const amrex::Geometry& geom)
    {
        return StencilJLO::box(bx, geom) & StencilKLO::box(bx, geom);
    }
};
 
struct StencilIloKlo
{
    static constexpr int direction = 12;
 
    static constexpr amrex::Real c00 = 1.0_rt / 3.0_rt;
    static constexpr amrex::Real c01 = 1.0_rt;
    static constexpr amrex::Real c02 = -4.0_rt / 3.0_rt;
    static constexpr amrex::Real c10 = 0.5_rt;
    static constexpr amrex::Real c11 = 0.0_rt;
    static constexpr amrex::Real c12 = -0.5_rt;
    static constexpr amrex::Real c20 = c00;
    static constexpr amrex::Real c21 = c01;
    static constexpr amrex::Real c22 = c02;
 
    static constexpr amrex::Real s00 = 4.0_rt / 3.0_rt; // axx
    static constexpr amrex::Real s01 = -4.0_rt;         // bxx
    static constexpr amrex::Real s02 = 8.0_rt / 3.0_rt; // cxx
    static constexpr amrex::Real s10 = 1.0_rt;          // ayy
    static constexpr amrex::Real s11 = -2.0_rt;         // byy
    static constexpr amrex::Real s12 = 1.0_rt;          // cyy
    static constexpr amrex::Real s20 = s00;             // azz
    static constexpr amrex::Real s21 = s01;             // bzz
    static constexpr amrex::Real s22 = s02;             // czz
 
    static constexpr amrex::Real f00 = 0.125_rt;
    static constexpr amrex::Real f01 = 0.625_rt;
    static constexpr amrex::Real f02 = 0.25_rt;
    static constexpr amrex::Real f10 = 0.25_rt;
    static constexpr amrex::Real f11 = 0.5_rt;
    static constexpr amrex::Real f12 = 0.25_rt;
    static constexpr amrex::Real f20 = f00;
    static constexpr amrex::Real f21 = f01;
    static constexpr amrex::Real f22 = f02;
 
    AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE static bool applies_to(
        int i,
        int /*j*/,
        int k,
        amrex::IntVect const& dlo,
        amrex::IntVect const& /*dhi*/) noexcept
    {
        return i == dlo[0] && k == dlo[2];
    }
 
    static amrex::Box box(const amrex::Box& bx, const amrex::Geometry& geom)
    {
        return StencilILO::box(bx, geom) & StencilKLO::box(bx, geom);
    }
};
 
struct StencilIhiJhi
{
    static constexpr int direction = 13;
 
    static constexpr amrex::Real c00 = 4.0_rt / 3.0_rt;
    static constexpr amrex::Real c01 = -1.0_rt;
    static constexpr amrex::Real c02 = -1.0_rt / 3.0_rt;
    static constexpr amrex::Real c10 = c00;
    static constexpr amrex::Real c11 = c01;
    static constexpr amrex::Real c12 = c02;
    static constexpr amrex::Real c20 = 0.5_rt;
    static constexpr amrex::Real c21 = 0.0_rt;
    static constexpr amrex::Real c22 = -0.5_rt;
 
    static constexpr amrex::Real s00 = 8.0_rt / 3.0_rt; // axx
    static constexpr amrex::Real s01 = -4.0_rt;         // bxx
    static constexpr amrex::Real s02 = 4.0_rt / 3.0_rt; // cxx
    static constexpr amrex::Real s10 = s00;             // ayy
    static constexpr amrex::Real s11 = s01;             // byy
    static constexpr amrex::Real s12 = s02;             // cyy
    static constexpr amrex::Real s20 = 1.0_rt;          // azz
    static constexpr amrex::Real s21 = -2.0_rt;         // bzz
    static constexpr amrex::Real s22 = 1.0_rt;          // czz
 
    static constexpr amrex::Real f00 = 0.25_rt;
    static constexpr amrex::Real f01 = 0.625_rt;
    static constexpr amrex::Real f02 = 0.125_rt;
    static constexpr amrex::Real f10 = f00;
    static constexpr amrex::Real f11 = f01;
    static constexpr amrex::Real f12 = f02;
    static constexpr amrex::Real f20 = 0.25_rt;
    static constexpr amrex::Real f21 = 0.5_rt;
    static constexpr amrex::Real f22 = 0.25_rt;
 
    AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE static bool applies_to(
        int i,
        int j,
        int /*k*/,
        amrex::IntVect const& /*dlo*/,
        amrex::IntVect const& dhi) noexcept
    {
        return i == dhi[0] && j == dhi[1];
    }
 
    static amrex::Box box(const amrex::Box& bx, const amrex::Geometry& geom)
    {
        return StencilIHI::box(bx, geom) & StencilJHI::box(bx, geom);
    }
};
 
struct StencilJhiKhi
{
    static constexpr int direction = 14;
 
    static constexpr amrex::Real c00 = 0.5_rt;
    static constexpr amrex::Real c01 = 0.0_rt;
    static constexpr amrex::Real c02 = -0.5_rt;
    static constexpr amrex::Real c10 = 4.0_rt / 3.0_rt;
    static constexpr amrex::Real c11 = -1.0_rt;
    static constexpr amrex::Real c12 = -1.0_rt / 3.0_rt;
    static constexpr amrex::Real c20 = c10;
    static constexpr amrex::Real c21 = c11;
    static constexpr amrex::Real c22 = c12;
 
    static constexpr amrex::Real s00 = 1.0_rt;          // axx
    static constexpr amrex::Real s01 = -2.0_rt;         // bxx
    static constexpr amrex::Real s02 = 1.0_rt;          // cxx
    static constexpr amrex::Real s10 = 8.0_rt / 3.0_rt; // ayy
    static constexpr amrex::Real s11 = -4.0_rt;         // byy
    static constexpr amrex::Real s12 = 4.0_rt / 3.0_rt; // cyy
    static constexpr amrex::Real s20 = s10;             // azz
    static constexpr amrex::Real s21 = s11;             // bzz
    static constexpr amrex::Real s22 = s12;             // czz
 
    static constexpr amrex::Real f00 = 0.25_rt;
    static constexpr amrex::Real f01 = 0.5_rt;
    static constexpr amrex::Real f02 = 0.25_rt;
    static constexpr amrex::Real f10 = 0.25_rt;
    static constexpr amrex::Real f11 = 0.625_rt;
    static constexpr amrex::Real f12 = 0.125_rt;
    static constexpr amrex::Real f20 = f10;
    static constexpr amrex::Real f21 = f11;
    static constexpr amrex::Real f22 = f12;
 
    AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE static bool applies_to(
        int /*i*/,
        int j,
        int k,
        amrex::IntVect const& /*dlo*/,
        amrex::IntVect const& dhi) noexcept
    {
        return j == dhi[1] && k == dhi[2];
    }
 
    static amrex::Box box(const amrex::Box& bx, const amrex::Geometry& geom)
    {
        return StencilJHI::box(bx, geom) & StencilKHI::box(bx, geom);
    }
};
 
struct StencilIhiKhi
{
    static constexpr int direction = 15;
 
    static constexpr amrex::Real c00 = 4.0_rt / 3.0_rt;
    static constexpr amrex::Real c01 = -1.0_rt;
    static constexpr amrex::Real c02 = -1.0_rt / 3.0_rt;
    static constexpr amrex::Real c10 = 0.5_rt;
    static constexpr amrex::Real c11 = 0.0_rt;
    static constexpr amrex::Real c12 = -0.5_rt;
    static constexpr amrex::Real c20 = c00;
    static constexpr amrex::Real c21 = c01;
    static constexpr amrex::Real c22 = c02;
 
    static constexpr amrex::Real s00 = 8.0_rt / 3.0_rt; // axx
    static constexpr amrex::Real s01 = -4.0_rt;         // bxx
    static constexpr amrex::Real s02 = 4.0_rt / 3.0_rt; // cxx
    static constexpr amrex::Real s10 = 1.0_rt;          // ayy
    static constexpr amrex::Real s11 = -2.0_rt;         // byy
    static constexpr amrex::Real s12 = 1.0_rt;          // cyy
    static constexpr amrex::Real s20 = s00;             // azz
    static constexpr amrex::Real s21 = s01;             // bzz
    static constexpr amrex::Real s22 = s02;             // czz
 
    static constexpr amrex::Real f00 = 0.25_rt;
    static constexpr amrex::Real f01 = 0.625_rt;
    static constexpr amrex::Real f02 = 0.125_rt;
    static constexpr amrex::Real f10 = 0.25_rt;
    static constexpr amrex::Real f11 = 0.5_rt;
    static constexpr amrex::Real f12 = 0.25_rt;
    static constexpr amrex::Real f20 = f00;
    static constexpr amrex::Real f21 = f01;
    static constexpr amrex::Real f22 = f02;
 
    AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE static bool applies_to(
        int i,
        int /*j*/,
        int k,
        amrex::IntVect const& /*dlo*/,
        amrex::IntVect const& dhi) noexcept
    {
        return i == dhi[0] && k == dhi[2];
    }
 
    static amrex::Box box(const amrex::Box& bx, const amrex::Geometry& geom)
    {
        return StencilIHI::box(bx, geom) & StencilKHI::box(bx, geom);
    }
};
 
struct StencilIloJhi
{
    static constexpr int direction = 16;
 
    static constexpr amrex::Real c00 = 1.0_rt / 3.0_rt;
    static constexpr amrex::Real c01 = 1.0_rt;
    static constexpr amrex::Real c02 = -4.0_rt / 3.0_rt;
    static constexpr amrex::Real c10 = 4.0_rt / 3.0_rt;
    static constexpr amrex::Real c11 = -1.0_rt;
    static constexpr amrex::Real c12 = -1.0_rt / 3.0_rt;
    static constexpr amrex::Real c20 = 0.5_rt;
    static constexpr amrex::Real c21 = 0.0_rt;
    static constexpr amrex::Real c22 = -0.5_rt;
 
    static constexpr amrex::Real s00 = 4.0_rt / 3.0_rt; // axx
    static constexpr amrex::Real s01 = -4.0_rt;         // bxx
    static constexpr amrex::Real s02 = 8.0_rt / 3.0_rt; // cxx
    static constexpr amrex::Real s10 = 8.0_rt / 3.0_rt; // ayy
    static constexpr amrex::Real s11 = -4.0_rt;         // byy
    static constexpr amrex::Real s12 = 4.0_rt / 3.0_rt; // cyy
    static constexpr amrex::Real s20 = 1.0_rt;          // azz
    static constexpr amrex::Real s21 = -2.0_rt;         // bzz
    static constexpr amrex::Real s22 = 1.0_rt;          // czz
 
    static constexpr amrex::Real f00 = 0.125_rt;
    static constexpr amrex::Real f01 = 0.625_rt;
    static constexpr amrex::Real f02 = 0.25_rt;
    static constexpr amrex::Real f10 = 0.25_rt;
    static constexpr amrex::Real f11 = 0.625_rt;
    static constexpr amrex::Real f12 = 0.125_rt;
    static constexpr amrex::Real f20 = 0.25_rt;
    static constexpr amrex::Real f21 = 0.5_rt;
    static constexpr amrex::Real f22 = 0.25_rt;
 
    AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE static bool applies_to(
        int i,
        int j,
        int /*k*/,
        amrex::IntVect const& dlo,
        amrex::IntVect const& dhi) noexcept
    {
        return i == dlo[0] && j == dhi[1];
    }
 
    static amrex::Box box(const amrex::Box& bx, const amrex::Geometry& geom)
    {
        return StencilILO::box(bx, geom) & StencilJHI::box(bx, geom);
    }
};
 
struct StencilJloKhi
{
    static constexpr int direction = 17;
 
    static constexpr amrex::Real c00 = 0.5_rt;
    static constexpr amrex::Real c01 = 0.0_rt;
    static constexpr amrex::Real c02 = -0.5_rt;
    static constexpr amrex::Real c10 = 1.0_rt / 3.0_rt;
    static constexpr amrex::Real c11 = 1.0_rt;
    static constexpr amrex::Real c12 = -4.0_rt / 3.0_rt;
    static constexpr amrex::Real c20 = 4.0_rt / 3.0_rt;
    static constexpr amrex::Real c21 = -1.0_rt;
    static constexpr amrex::Real c22 = -1.0_rt / 3.0_rt;
 
    static constexpr amrex::Real s00 = 1.0_rt;          // axx
    static constexpr amrex::Real s01 = -2.0_rt;         // bxx
    static constexpr amrex::Real s02 = 1.0_rt;          // cxx
    static constexpr amrex::Real s10 = 4.0_rt / 3.0_rt; // ayy
    static constexpr amrex::Real s11 = -4.0_rt;         // byy
    static constexpr amrex::Real s12 = 8.0_rt / 3.0_rt; // cyy
    static constexpr amrex::Real s20 = 8.0_rt / 3.0_rt; // azz
    static constexpr amrex::Real s21 = -4.0_rt;         // bzz
    static constexpr amrex::Real s22 = 4.0_rt / 3.0_rt; // czz
 
    static constexpr amrex::Real f00 = 0.25_rt;
    static constexpr amrex::Real f01 = 0.5_rt;
    static constexpr amrex::Real f02 = 0.25_rt;
    static constexpr amrex::Real f10 = 0.125_rt;
    static constexpr amrex::Real f11 = 0.625_rt;
    static constexpr amrex::Real f12 = 0.25_rt;
    static constexpr amrex::Real f20 = 0.25_rt;
    static constexpr amrex::Real f21 = 0.625_rt;
    static constexpr amrex::Real f22 = 0.125_rt;
 
    AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE static bool applies_to(
        int /*i*/,
        int j,
        int k,
        amrex::IntVect const& dlo,
        amrex::IntVect const& dhi) noexcept
    {
        return j == dlo[1] && k == dhi[2];
    }
 
    static amrex::Box box(const amrex::Box& bx, const amrex::Geometry& geom)
    {
        return StencilJLO::box(bx, geom) & StencilKHI::box(bx, geom);
    }
};
 
struct StencilIloKhi
{
    static constexpr int direction = 18;
 
    static constexpr amrex::Real c00 = 1.0_rt / 3.0_rt;
    static constexpr amrex::Real c01 = 1.0_rt;
    static constexpr amrex::Real c02 = -4.0_rt / 3.0_rt;
    static constexpr amrex::Real c10 = 0.5_rt;
    static constexpr amrex::Real c11 = 0.0_rt;
    static constexpr amrex::Real c12 = -0.5_rt;
    static constexpr amrex::Real c20 = 4.0_rt / 3.0_rt;
    static constexpr amrex::Real c21 = -1.0_rt;
    static constexpr amrex::Real c22 = -1.0_rt / 3.0_rt;
 
    static constexpr amrex::Real s00 = 4.0_rt / 3.0_rt; // axx
    static constexpr amrex::Real s01 = -4.0_rt;         // bxx
    static constexpr amrex::Real s02 = 8.0_rt / 3.0_rt; // cxx
    static constexpr amrex::Real s10 = 1;               // ayy
    static constexpr amrex::Real s11 = -2;              // byy
    static constexpr amrex::Real s12 = 1;               // cyy
    static constexpr amrex::Real s20 = 8.0_rt / 3.0_rt; // azz
    static constexpr amrex::Real s21 = -4.0_rt;         // bzz
    static constexpr amrex::Real s22 = 4.0_rt / 3.0_rt; // czz
 
    static constexpr amrex::Real f00 = 0.125_rt;
    static constexpr amrex::Real f01 = 0.625_rt;
    static constexpr amrex::Real f02 = 0.25_rt;
    static constexpr amrex::Real f10 = 0.25_rt;
    static constexpr amrex::Real f11 = 0.5_rt;
    static constexpr amrex::Real f12 = 0.25_rt;
    static constexpr amrex::Real f20 = 0.25_rt;
    static constexpr amrex::Real f21 = 0.625_rt;
    static constexpr amrex::Real f22 = 0.125_rt;
 
    AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE static bool applies_to(
        int i,
        int /*j*/,
        int k,
        amrex::IntVect const& dlo,
        amrex::IntVect const& dhi) noexcept
    {
        return i == dlo[0] && k == dhi[2];
    }
 
    static amrex::Box box(const amrex::Box& bx, const amrex::Geometry& geom)
    {
        return StencilILO::box(bx, geom) & StencilKHI::box(bx, geom);
    }
};
 
struct StencilIhiJlo
{
    static constexpr int direction = 19;
 
    static constexpr amrex::Real c00 = 4.0_rt / 3.0_rt;
    static constexpr amrex::Real c01 = -1.0_rt;
    static constexpr amrex::Real c02 = -1.0_rt / 3.0_rt;
    static constexpr amrex::Real c10 = 1.0_rt / 3.0_rt;
    static constexpr amrex::Real c11 = 1;
    static constexpr amrex::Real c12 = -4.0_rt / 3.0_rt;
    static constexpr amrex::Real c20 = 0.5_rt;
    static constexpr amrex::Real c21 = 0.0_rt;
    static constexpr amrex::Real c22 = -0.5_rt;
 
    static constexpr amrex::Real s00 = 8.0_rt / 3.0_rt; // axx
    static constexpr amrex::Real s01 = -4.0_rt;         // bxx
    static constexpr amrex::Real s02 = 4.0_rt / 3.0_rt; // cxx
    static constexpr amrex::Real s10 = 4.0_rt / 3.0_rt; // ayy
    static constexpr amrex::Real s11 = -4.0_rt;         // byy
    static constexpr amrex::Real s12 = 8.0_rt / 3.0_rt; // cyy
    static constexpr amrex::Real s20 = 1.0_rt;          // azz
    static constexpr amrex::Real s21 = -2.0_rt;         // bzz
    static constexpr amrex::Real s22 = 1.0_rt;          // czz
 
    static constexpr amrex::Real f00 = 0.25_rt;
    static constexpr amrex::Real f01 = 0.625_rt;
    static constexpr amrex::Real f02 = 0.125_rt;
    static constexpr amrex::Real f10 = 0.125_rt;
    static constexpr amrex::Real f11 = 0.625_rt;
    static constexpr amrex::Real f12 = 0.25_rt;
    static constexpr amrex::Real f20 = 0.25_rt;
    static constexpr amrex::Real f21 = 0.5_rt;
    static constexpr amrex::Real f22 = 0.25_rt;
 
    AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE static bool applies_to(
        int i,
        int j,
        int /*k*/,
        amrex::IntVect const& dlo,
        amrex::IntVect const& dhi) noexcept
    {
        return i == dhi[0] && j == dlo[1];
    }
 
    static amrex::Box box(const amrex::Box& bx, const amrex::Geometry& geom)
    {
        return StencilIHI::box(bx, geom) & StencilJLO::box(bx, geom);
    }
};
 
struct StencilJhiKlo
{
    static constexpr int direction = 20;
 
    static constexpr amrex::Real c00 = 0.5_rt;
    static constexpr amrex::Real c01 = 0.0_rt;
    static constexpr amrex::Real c02 = -0.5_rt;
    static constexpr amrex::Real c10 = 4.0_rt / 3.0_rt;
    static constexpr amrex::Real c11 = -1.0_rt;
    static constexpr amrex::Real c12 = -1.0_rt / 3.0_rt;
    static constexpr amrex::Real c20 = 1.0_rt / 3.0_rt;
    static constexpr amrex::Real c21 = 1.0_rt;
    static constexpr amrex::Real c22 = -4.0_rt / 3.0_rt;
 
    static constexpr amrex::Real s00 = 1.0_rt;          // axx
    static constexpr amrex::Real s01 = -2.0_rt;         // bxx
    static constexpr amrex::Real s02 = 1.0_rt;          // cxx
    static constexpr amrex::Real s10 = 8.0_rt / 3.0_rt; // ayy
    static constexpr amrex::Real s11 = -4.0_rt;         // byy
    static constexpr amrex::Real s12 = 4.0_rt / 3.0_rt; // cyy
    static constexpr amrex::Real s20 = 4.0_rt / 3.0_rt; // azz
    static constexpr amrex::Real s21 = -4.0_rt;         // bzz
    static constexpr amrex::Real s22 = 8.0_rt / 3.0_rt; // czz
 
    static constexpr amrex::Real f00 = 0.25_rt;
    static constexpr amrex::Real f01 = 0.5_rt;
    static constexpr amrex::Real f02 = 0.25_rt;
    static constexpr amrex::Real f10 = 0.25_rt;
    static constexpr amrex::Real f11 = 0.625_rt;
    static constexpr amrex::Real f12 = 0.125_rt;
    static constexpr amrex::Real f20 = 0.125_rt;
    static constexpr amrex::Real f21 = 0.625_rt;
    static constexpr amrex::Real f22 = 0.25_rt;
 
    AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE static bool applies_to(
        int /*i*/,
        int j,
        int k,
        amrex::IntVect const& dlo,
        amrex::IntVect const& dhi) noexcept
    {
        return j == dhi[1] && k == dlo[2];
    }
 
    static amrex::Box box(const amrex::Box& bx, const amrex::Geometry& geom)
    {
        return StencilJHI::box(bx, geom) & StencilKLO::box(bx, geom);
    }
};
 
struct StencilIhiKlo
{
    static constexpr int direction = 21;
 
    static constexpr amrex::Real c00 = 4.0_rt / 3.0_rt;
    static constexpr amrex::Real c01 = -1.0_rt;
    static constexpr amrex::Real c02 = -1.0_rt / 3.0_rt;
    static constexpr amrex::Real c10 = 0.5_rt;
    static constexpr amrex::Real c11 = 0.0_rt;
    static constexpr amrex::Real c12 = -0.5_rt;
    static constexpr amrex::Real c20 = 1.0_rt / 3.0_rt;
    static constexpr amrex::Real c21 = 1.0_rt;
    static constexpr amrex::Real c22 = -4.0_rt / 3.0_rt;
 
    static constexpr amrex::Real s00 = 8.0_rt / 3.0_rt; // axx
    static constexpr amrex::Real s01 = -4.0_rt;         // bxx
    static constexpr amrex::Real s02 = 4.0_rt / 3.0_rt; // cxx
    static constexpr amrex::Real s10 = 1.0_rt;          // ayy
    static constexpr amrex::Real s11 = -2.0_rt;         // byy
    static constexpr amrex::Real s12 = 1.0_rt;          // cyy
    static constexpr amrex::Real s20 = 4.0_rt / 3.0_rt; // azz
    static constexpr amrex::Real s21 = -4.0_rt;         // bzz
    static constexpr amrex::Real s22 = 8.0_rt / 3.0_rt; // czz
 
    static constexpr amrex::Real f00 = 0.25_rt;
    static constexpr amrex::Real f01 = 0.625_rt;
    static constexpr amrex::Real f02 = 0.125_rt;
    static constexpr amrex::Real f10 = 0.25_rt;
    static constexpr amrex::Real f11 = 0.5_rt;
    static constexpr amrex::Real f12 = 0.25_rt;
    static constexpr amrex::Real f20 = 0.125_rt;
    static constexpr amrex::Real f21 = 0.625_rt;
    static constexpr amrex::Real f22 = 0.25_rt;
 
    AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE static bool applies_to(
        int i,
        int /*j*/,
        int k,
        amrex::IntVect const& dlo,
        amrex::IntVect const& dhi) noexcept
    {
        return i == dhi[0] && k == dlo[2];
    }
 
    static amrex::Box box(const amrex::Box& bx, const amrex::Geometry& geom)
    {
        return StencilIHI::box(bx, geom) & StencilKLO::box(bx, geom);
    }
};
 
struct StencilIloJloKlo
{
    static constexpr int direction = 22;
 
    static constexpr amrex::Real c00 = 1.0_rt / 3.0_rt;
    static constexpr amrex::Real c01 = 1.0_rt;
    static constexpr amrex::Real c02 = -4.0_rt / 3.0_rt;
    static constexpr amrex::Real c10 = c00;
    static constexpr amrex::Real c11 = c01;
    static constexpr amrex::Real c12 = c02;
    static constexpr amrex::Real c20 = c00;
    static constexpr amrex::Real c21 = c01;
    static constexpr amrex::Real c22 = c02;
 
    static constexpr amrex::Real s00 = 4.0_rt / 3.0_rt; // axx
    static constexpr amrex::Real s01 = -4.0_rt;         // bxx
    static constexpr amrex::Real s02 = 8.0_rt / 3.0_rt; // cxx
    static constexpr amrex::Real s10 = s00;             // ayy
    static constexpr amrex::Real s11 = s01;             // byy
    static constexpr amrex::Real s12 = s02;             // cyy
    static constexpr amrex::Real s20 = s00;             // azz
    static constexpr amrex::Real s21 = s01;             // bzz
    static constexpr amrex::Real s22 = s02;             // czz
 
    static constexpr amrex::Real f00 = 0.125_rt;
    static constexpr amrex::Real f01 = 0.625_rt;
    static constexpr amrex::Real f02 = 0.25_rt;
    static constexpr amrex::Real f10 = f00;
    static constexpr amrex::Real f11 = f01;
    static constexpr amrex::Real f12 = f02;
    static constexpr amrex::Real f20 = f00;
    static constexpr amrex::Real f21 = f01;
    static constexpr amrex::Real f22 = f02;
 
    AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE static bool applies_to(
        int i,
        int j,
        int k,
        amrex::IntVect const& dlo,
        amrex::IntVect const& /*dhi*/) noexcept
    {
        return i == dlo[0] && j == dlo[1] && k == dlo[2];
    }
 
    static amrex::Box box(const amrex::Box& bx, const amrex::Geometry& geom)
    {
        return StencilILO::box(bx, geom) & StencilJLO::box(bx, geom) &
               StencilKLO::box(bx, geom);
    }
};
 
struct StencilIloJhiKlo
{
    static constexpr int direction = 23;
 
    static constexpr amrex::Real c00 = 1.0_rt / 3.0_rt;
    static constexpr amrex::Real c01 = 1.0_rt;
    static constexpr amrex::Real c02 = -4.0_rt / 3.0_rt;
    static constexpr amrex::Real c10 = 4.0_rt / 3.0_rt;
    static constexpr amrex::Real c11 = -1.0_rt;
    static constexpr amrex::Real c12 = -1.0_rt / 3.0_rt;
    static constexpr amrex::Real c20 = c00;
    static constexpr amrex::Real c21 = c01;
    static constexpr amrex::Real c22 = c02;
 
    static constexpr amrex::Real s00 = 4.0_rt / 3.0_rt; // axx
    static constexpr amrex::Real s01 = -4.0_rt;         // bxx
    static constexpr amrex::Real s02 = 8.0_rt / 3.0_rt; // cxx
    static constexpr amrex::Real s10 = 8.0_rt / 3.0_rt; // ayy
    static constexpr amrex::Real s11 = -4.0_rt;         // byy
    static constexpr amrex::Real s12 = 4.0_rt / 3.0_rt; // cyy
    static constexpr amrex::Real s20 = s00;             // azz
    static constexpr amrex::Real s21 = s01;             // bzz
    static constexpr amrex::Real s22 = s02;             // czz
 
    static constexpr amrex::Real f00 = 0.125_rt;
    static constexpr amrex::Real f01 = 0.625_rt;
    static constexpr amrex::Real f02 = 0.25_rt;
    static constexpr amrex::Real f10 = 0.25_rt;
    static constexpr amrex::Real f11 = 0.625_rt;
    static constexpr amrex::Real f12 = 0.125_rt;
    static constexpr amrex::Real f20 = f00;
    static constexpr amrex::Real f21 = f01;
    static constexpr amrex::Real f22 = f02;
 
    AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE static bool applies_to(
        int i,
        int j,
        int k,
        amrex::IntVect const& dlo,
        amrex::IntVect const& dhi) noexcept
    {
        return i == dlo[0] && j == dhi[1] && k == dlo[2];
    }
 
    static amrex::Box box(const amrex::Box& bx, const amrex::Geometry& geom)
    {
        return StencilILO::box(bx, geom) & StencilJHI::box(bx, geom) &
               StencilKLO::box(bx, geom);
    }
};
 
struct StencilIhiJloKlo
{
    static constexpr int direction = 24;
 
    static constexpr amrex::Real c00 = 4.0_rt / 3.0_rt;
    static constexpr amrex::Real c01 = -1.0_rt;
    static constexpr amrex::Real c02 = -1.0_rt / 3.0_rt;
    static constexpr amrex::Real c10 = 1.0_rt / 3.0_rt;
    static constexpr amrex::Real c11 = 1.0_rt;
    static constexpr amrex::Real c12 = -4.0_rt / 3.0_rt;
    static constexpr amrex::Real c20 = c10;
    static constexpr amrex::Real c21 = c11;
    static constexpr amrex::Real c22 = c12;
 
    static constexpr amrex::Real s00 = 8.0_rt / 3.0_rt; // axx
    static constexpr amrex::Real s01 = -4.0_rt;         // bxx
    static constexpr amrex::Real s02 = 4.0_rt / 3.0_rt; // cxx
    static constexpr amrex::Real s10 = 4.0_rt / 3.0_rt; // ayy
    static constexpr amrex::Real s11 = -4.0_rt;         // byy
    static constexpr amrex::Real s12 = 8.0_rt / 3.0_rt; // cyy
    static constexpr amrex::Real s20 = s10;             // azz
    static constexpr amrex::Real s21 = s11;             // bzz
    static constexpr amrex::Real s22 = s12;             // czz
 
    static constexpr amrex::Real f00 = 0.25_rt;
    static constexpr amrex::Real f01 = 0.625_rt;
    static constexpr amrex::Real f02 = 0.125_rt;
    static constexpr amrex::Real f10 = 0.125_rt;
    static constexpr amrex::Real f11 = 0.625_rt;
    static constexpr amrex::Real f12 = 0.25_rt;
    static constexpr amrex::Real f20 = f10;
    static constexpr amrex::Real f21 = f11;
    static constexpr amrex::Real f22 = f12;
 
    AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE static bool applies_to(
        int i,
        int j,
        int k,
        amrex::IntVect const& dlo,
        amrex::IntVect const& dhi) noexcept
    {
        return i == dhi[0] && j == dlo[1] && k == dlo[2];
    }
 
    static amrex::Box box(const amrex::Box& bx, const amrex::Geometry& geom)
    {
        return StencilIHI::box(bx, geom) & StencilJLO::box(bx, geom) &
               StencilKLO::box(bx, geom);
    }
};
 
struct StencilIhiJhiKlo
{
    static constexpr int direction = 25;
 
    static constexpr amrex::Real c00 = 4.0_rt / 3.0_rt;
    static constexpr amrex::Real c01 = -1.0_rt;
    static constexpr amrex::Real c02 = -1.0_rt / 3.0_rt;
    static constexpr amrex::Real c10 = 4.0_rt / 3.0_rt;
    static constexpr amrex::Real c11 = -1.0_rt;
    static constexpr amrex::Real c12 = -1.0_rt / 3.0_rt;
    static constexpr amrex::Real c20 = 1.0_rt / 3.0_rt;
    static constexpr amrex::Real c21 = 1.0_rt;
    static constexpr amrex::Real c22 = -4.0_rt / 3.0_rt;
 
    static constexpr amrex::Real s00 = 8.0_rt / 3.0_rt; // axx
    static constexpr amrex::Real s01 = -4.0_rt;         // bxx
    static constexpr amrex::Real s02 = 4.0_rt / 3.0_rt; // cxx
    static constexpr amrex::Real s10 = 8.0_rt / 3.0_rt; // ayy
    static constexpr amrex::Real s11 = -4.0_rt;         // byy
    static constexpr amrex::Real s12 = 4.0_rt / 3.0_rt; // cyy
    static constexpr amrex::Real s20 = 4.0_rt / 3.0_rt; // azz
    static constexpr amrex::Real s21 = -4.0_rt;         // bzz
    static constexpr amrex::Real s22 = 8.0_rt / 3.0_rt; // czz
 
    static constexpr amrex::Real f00 = 0.25_rt;
    static constexpr amrex::Real f01 = 0.625_rt;
    static constexpr amrex::Real f02 = 0.125_rt;
    static constexpr amrex::Real f10 = 0.25_rt;
    static constexpr amrex::Real f11 = 0.625_rt;
    static constexpr amrex::Real f12 = 0.125_rt;
    static constexpr amrex::Real f20 = 0.125_rt;
    static constexpr amrex::Real f21 = 0.625_rt;
    static constexpr amrex::Real f22 = 0.25_rt;
 
    AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE static bool applies_to(
        int i,
        int j,
        int k,
        amrex::IntVect const& dlo,
        amrex::IntVect const& dhi) noexcept
    {
        return i == dhi[0] && j == dhi[1] && k == dlo[2];
    }
 
    static amrex::Box box(const amrex::Box& bx, const amrex::Geometry& geom)
    {
        return StencilIHI::box(bx, geom) & StencilJHI::box(bx, geom) &
               StencilKLO::box(bx, geom);
    }
};
 
struct StencilIloJloKhi
{
    static constexpr int direction = 26;
 
    static constexpr amrex::Real c00 = 1.0_rt / 3.0_rt;
    static constexpr amrex::Real c01 = 1.0_rt;
    static constexpr amrex::Real c02 = -4.0_rt / 3.0_rt;
    static constexpr amrex::Real c10 = c00;
    static constexpr amrex::Real c11 = c01;
    static constexpr amrex::Real c12 = c02;
    static constexpr amrex::Real c20 = 4.0_rt / 3.0_rt;
    static constexpr amrex::Real c21 = -1.0_rt;
    static constexpr amrex::Real c22 = -1.0_rt / 3.0_rt;
 
    static constexpr amrex::Real s00 = 4.0_rt / 3.0_rt; // axx
    static constexpr amrex::Real s01 = -4.0_rt;         // bxx
    static constexpr amrex::Real s02 = 8.0_rt / 3.0_rt; // cxx
    static constexpr amrex::Real s10 = s00;             // ayy
    static constexpr amrex::Real s11 = s01;             // byy
    static constexpr amrex::Real s12 = s02;             // cyy
    static constexpr amrex::Real s20 = 8.0_rt / 3.0_rt; // azz
    static constexpr amrex::Real s21 = -4.0_rt;         // bzz
    static constexpr amrex::Real s22 = 4.0_rt / 3.0_rt; // czz
 
    static constexpr amrex::Real f00 = 0.125_rt;
    static constexpr amrex::Real f01 = 0.625_rt;
    static constexpr amrex::Real f02 = 0.25_rt;
    static constexpr amrex::Real f10 = f00;
    static constexpr amrex::Real f11 = f01;
    static constexpr amrex::Real f12 = f02;
    static constexpr amrex::Real f20 = 0.25_rt;
    static constexpr amrex::Real f21 = 0.625_rt;
    static constexpr amrex::Real f22 = 0.125_rt;
 
    AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE static bool applies_to(
        int i,
        int j,
        int k,
        amrex::IntVect const& dlo,
        amrex::IntVect const& dhi) noexcept
    {
        return i == dlo[0] && j == dlo[1] && k == dhi[2];
    }
 
    static amrex::Box box(const amrex::Box& bx, const amrex::Geometry& geom)
    {
        return StencilILO::box(bx, geom) & StencilJLO::box(bx, geom) &
               StencilKHI::box(bx, geom);
    }
};
 
struct StencilIloJhiKhi
{
    static constexpr int direction = 27;
 
    static constexpr amrex::Real c00 = 1.0_rt / 3.0_rt;
    static constexpr amrex::Real c01 = 1.0_rt;
    static constexpr amrex::Real c02 = -4.0_rt / 3.0_rt;
    static constexpr amrex::Real c10 = 4.0_rt / 3.0_rt;
    static constexpr amrex::Real c11 = -1.0_rt;
    static constexpr amrex::Real c12 = -1.0_rt / 3.0_rt;
    static constexpr amrex::Real c20 = 4.0_rt / 3.0_rt;
    static constexpr amrex::Real c21 = -1.0_rt;
    static constexpr amrex::Real c22 = -1.0_rt / 3.0_rt;
 
    static constexpr amrex::Real s00 = 4.0_rt / 3.0_rt; // axx
    static constexpr amrex::Real s01 = -4.0_rt;         // bxx
    static constexpr amrex::Real s02 = 8.0_rt / 3.0_rt; // cxx
    static constexpr amrex::Real s10 = 8.0_rt / 3.0_rt; // ayy
    static constexpr amrex::Real s11 = -4.0_rt;         // byy
    static constexpr amrex::Real s12 = 4.0_rt / 3.0_rt; // cyy
    static constexpr amrex::Real s20 = 8.0_rt / 3.0_rt; // azz
    static constexpr amrex::Real s21 = -4.0_rt;         // bzz
    static constexpr amrex::Real s22 = 4.0_rt / 3.0_rt; // czz
 
    static constexpr amrex::Real f00 = 0.125_rt;
    static constexpr amrex::Real f01 = 0.625_rt;
    static constexpr amrex::Real f02 = 0.25_rt;
    static constexpr amrex::Real f10 = 0.25_rt;
    static constexpr amrex::Real f11 = 0.625_rt;
    static constexpr amrex::Real f12 = 0.125_rt;
    static constexpr amrex::Real f20 = 0.25_rt;
    static constexpr amrex::Real f21 = 0.625_rt;
    static constexpr amrex::Real f22 = 0.125_rt;
 
    AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE static bool applies_to(
        int i,
        int j,
        int k,
        amrex::IntVect const& dlo,
        amrex::IntVect const& dhi) noexcept
    {
        return i == dlo[0] && j == dhi[1] && k == dhi[2];
    }
 
    static amrex::Box box(const amrex::Box& bx, const amrex::Geometry& geom)
    {
        return StencilILO::box(bx, geom) & StencilJHI::box(bx, geom) &
               StencilKHI::box(bx, geom);
    }
};
 
struct StencilIhiJloKhi
{
    static constexpr int direction = 28;
 
    static constexpr amrex::Real c00 = 4.0_rt / 3.0_rt;
    static constexpr amrex::Real c01 = -1.0_rt;
    static constexpr amrex::Real c02 = -1.0_rt / 3.0_rt;
    static constexpr amrex::Real c10 = 1.0_rt / 3.0_rt;
    static constexpr amrex::Real c11 = 1.0_rt;
    static constexpr amrex::Real c12 = -4.0_rt / 3.0_rt;
    static constexpr amrex::Real c20 = 4.0_rt / 3.0_rt;
    static constexpr amrex::Real c21 = -1.0_rt;
    static constexpr amrex::Real c22 = -1.0_rt / 3.0_rt;
 
    static constexpr amrex::Real s00 = 8.0_rt / 3.0_rt; // axx
    static constexpr amrex::Real s01 = -4.0_rt;         // bxx
    static constexpr amrex::Real s02 = 4.0_rt / 3.0_rt; // cxx
    static constexpr amrex::Real s10 = 4.0_rt / 3.0_rt; // ayy
    static constexpr amrex::Real s11 = -4.0_rt;         // byy
    static constexpr amrex::Real s12 = 8.0_rt / 3.0_rt; // cyy
    static constexpr amrex::Real s20 = 8.0_rt / 3.0_rt; // azz
    static constexpr amrex::Real s21 = -4.0_rt;         // bzz
    static constexpr amrex::Real s22 = 4.0_rt / 3.0_rt; // czz
 
    static constexpr amrex::Real f00 = 0.25_rt;
    static constexpr amrex::Real f01 = 0.625_rt;
    static constexpr amrex::Real f02 = 0.125_rt;
    static constexpr amrex::Real f10 = 0.125_rt;
    static constexpr amrex::Real f11 = 0.625_rt;
    static constexpr amrex::Real f12 = 0.25_rt;
    static constexpr amrex::Real f20 = 0.25_rt;
    static constexpr amrex::Real f21 = 0.625_rt;
    static constexpr amrex::Real f22 = 0.125_rt;
 
    AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE static bool applies_to(
        int i,
        int j,
        int k,
        amrex::IntVect const& dlo,
        amrex::IntVect const& dhi) noexcept
    {
        return i == dhi[0] && j == dlo[1] && k == dhi[2];
    }
 
    static amrex::Box box(const amrex::Box& bx, const amrex::Geometry& geom)
    {
        return StencilIHI::box(bx, geom) & StencilJLO::box(bx, geom) &
               StencilKHI::box(bx, geom);
    }
};
 
struct StencilIhiJhiKhi
{
    static constexpr int direction = 29;
 
    static constexpr amrex::Real c00 = 4.0_rt / 3.0_rt;
    static constexpr amrex::Real c01 = -1.0_rt;
    static constexpr amrex::Real c02 = -1.0_rt / 3.0_rt;
    static constexpr amrex::Real c10 = 4.0_rt / 3.0_rt;
    static constexpr amrex::Real c11 = -1.0_rt;
    static constexpr amrex::Real c12 = -1.0_rt / 3.0_rt;
    static constexpr amrex::Real c20 = 4.0_rt / 3.0_rt;
    static constexpr amrex::Real c21 = -1.0_rt;
    static constexpr amrex::Real c22 = -1.0_rt / 3.0_rt;
 
    static constexpr amrex::Real s00 = 8.0_rt / 3.0_rt; // axx
    static constexpr amrex::Real s01 = -4.0_rt;         // bxx
    static constexpr amrex::Real s02 = 4.0_rt / 3.0_rt; // cxx
    static constexpr amrex::Real s10 = 8.0_rt / 3.0_rt; // ayy
    static constexpr amrex::Real s11 = -4.0_rt;         // byy
    static constexpr amrex::Real s12 = 4.0_rt / 3.0_rt; // cyy
    static constexpr amrex::Real s20 = 8.0_rt / 3.0_rt; // azz
    static constexpr amrex::Real s21 = -4.0_rt;         // bzz
    static constexpr amrex::Real s22 = 4.0_rt / 3.0_rt; // czz
 
    static constexpr amrex::Real f00 = 0.25_rt;
    static constexpr amrex::Real f01 = 0.625_rt;
    static constexpr amrex::Real f02 = 0.125_rt;
    static constexpr amrex::Real f10 = 0.25_rt;
    static constexpr amrex::Real f11 = 0.625_rt;
    static constexpr amrex::Real f12 = 0.125_rt;
    static constexpr amrex::Real f20 = 0.25_rt;
    static constexpr amrex::Real f21 = 0.625_rt;
    static constexpr amrex::Real f22 = 0.125_rt;
 
    AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE static bool applies_to(
        int i,
        int j,
        int k,
        amrex::IntVect const& /*dlo*/,
        amrex::IntVect const& dhi) noexcept
    {
        return i == dhi[0] && j == dhi[1] && k == dhi[2];
    }
 
    static amrex::Box box(const amrex::Box& bx, const amrex::Geometry& geom)
    {
        return StencilIHI::box(bx, geom) & StencilJHI::box(bx, geom) &
               StencilKHI::box(bx, geom);
    }
};
 
} // namespace kynema_sgf::fvm::stencil
 
#endif /* STENCILS_H */