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"
 
namespace amr_wind::fvm::stencil {
 
namespace impl {
 
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};
}
 
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;
    static constexpr amrex::Real c01 = 0.0;
    static constexpr amrex::Real c02 = -0.5;
    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;
    static constexpr amrex::Real s01 = -2.0;
    static constexpr amrex::Real s02 = 1.0;
    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;
    static constexpr amrex::Real f01 = 0.5;
    static constexpr amrex::Real f02 = 0.25;
    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;
 
    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 / 3.0;
    static constexpr amrex::Real c01 = 1.0;
    static constexpr amrex::Real c02 = -4.0 / 3.0;
    static constexpr amrex::Real c10 = 0.5;
    static constexpr amrex::Real c11 = 0.0;
    static constexpr amrex::Real c12 = -0.5;
    static constexpr amrex::Real c20 = c10;
    static constexpr amrex::Real c21 = c11;
    static constexpr amrex::Real c22 = c12;
 
    static constexpr amrex::Real s00 = 4. / 3.; // axx
    static constexpr amrex::Real s01 = -4;      // bxx
    static constexpr amrex::Real s02 = 8. / 3.; // 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 = s10;     // azz
    static constexpr amrex::Real s21 = s11;     // bzz
    static constexpr amrex::Real s22 = s12;     // czz
 
    static constexpr amrex::Real f00 = 0.125;
    static constexpr amrex::Real f01 = 0.625;
    static constexpr amrex::Real f02 = 0.25;
    static constexpr amrex::Real f10 = 0.25;
    static constexpr amrex::Real f11 = 0.5;
    static constexpr amrex::Real f12 = 0.25;
    static constexpr amrex::Real f20 = f10;
    static constexpr amrex::Real f21 = f11;
    static constexpr amrex::Real f22 = f12;
 
    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;
    static constexpr amrex::Real c01 = 0.0;
    static constexpr amrex::Real c02 = -0.5;
    static constexpr amrex::Real c10 = 1.0 / 3.0;
    static constexpr amrex::Real c11 = 1.0;
    static constexpr amrex::Real c12 = -4.0 / 3.0;
    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;     // axx
    static constexpr amrex::Real s01 = -2.0;    // bxx
    static constexpr amrex::Real s02 = 1.0;     // cxx
    static constexpr amrex::Real s10 = 4. / 3.; // ayy
    static constexpr amrex::Real s11 = -4.;     // byy
    static constexpr amrex::Real s12 = 8. / 3.; // 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;
    static constexpr amrex::Real f01 = 0.5;
    static constexpr amrex::Real f02 = 0.25;
    static constexpr amrex::Real f10 = 0.125;
    static constexpr amrex::Real f11 = 0.625;
    static constexpr amrex::Real f12 = 0.25;
    static constexpr amrex::Real f20 = f00;
    static constexpr amrex::Real f21 = f01;
    static constexpr amrex::Real f22 = f02;
 
    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;
    static constexpr amrex::Real c01 = 0.0;
    static constexpr amrex::Real c02 = -0.5;
    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 / 3.0;
    static constexpr amrex::Real c21 = 1.0;
    static constexpr amrex::Real c22 = -4.0 / 3.0;
 
    static constexpr amrex::Real s00 = 1.0;     // axx
    static constexpr amrex::Real s01 = -2.0;    // bxx
    static constexpr amrex::Real s02 = 1.0;     // 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. / 3.; // azz
    static constexpr amrex::Real s21 = -4.;     // bzz
    static constexpr amrex::Real s22 = 8. / 3.; // czz
 
    static constexpr amrex::Real f00 = 0.25;
    static constexpr amrex::Real f01 = 0.5;
    static constexpr amrex::Real f02 = 0.25;
    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;
    static constexpr amrex::Real f21 = 0.625;
    static constexpr amrex::Real f22 = 0.25;
 
    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 / 3.0;
    static constexpr amrex::Real c01 = -1.0;
    static constexpr amrex::Real c02 = -1.0 / 3.0;
    static constexpr amrex::Real c10 = 0.5;
    static constexpr amrex::Real c11 = 0.0;
    static constexpr amrex::Real c12 = -0.5;
    static constexpr amrex::Real c20 = c10;
    static constexpr amrex::Real c21 = c11;
    static constexpr amrex::Real c22 = c12;
 
    static constexpr amrex::Real s00 = 8. / 3.; // axx
    static constexpr amrex::Real s01 = -4.;     // bxx
    static constexpr amrex::Real s02 = 4. / 3.; // 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 = s10;     // azz
    static constexpr amrex::Real s21 = s11;     // bzz
    static constexpr amrex::Real s22 = s12;     // czz
 
    static constexpr amrex::Real f00 = 0.25;
    static constexpr amrex::Real f01 = 0.625;
    static constexpr amrex::Real f02 = 0.125;
    static constexpr amrex::Real f10 = 0.25;
    static constexpr amrex::Real f11 = 0.5;
    static constexpr amrex::Real f12 = 0.25;
    static constexpr amrex::Real f20 = f10;
    static constexpr amrex::Real f21 = f11;
    static constexpr amrex::Real f22 = f12;
 
    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;
    static constexpr amrex::Real c01 = 0.0;
    static constexpr amrex::Real c02 = -0.5;
    static constexpr amrex::Real c10 = 4.0 / 3.0;
    static constexpr amrex::Real c11 = -1.0;
    static constexpr amrex::Real c12 = -1.0 / 3.0;
    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;     // axx
    static constexpr amrex::Real s01 = -2.0;    // bxx
    static constexpr amrex::Real s02 = 1.0;     // cxx
    static constexpr amrex::Real s10 = 8. / 3.; // ayy
    static constexpr amrex::Real s11 = -4.;     // byy
    static constexpr amrex::Real s12 = 4. / 3.; // 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;
    static constexpr amrex::Real f01 = 0.5;
    static constexpr amrex::Real f02 = 0.25;
    static constexpr amrex::Real f10 = 0.25;
    static constexpr amrex::Real f11 = 0.625;
    static constexpr amrex::Real f12 = 0.125;
    static constexpr amrex::Real f20 = f00;
    static constexpr amrex::Real f21 = f01;
    static constexpr amrex::Real f22 = f02;
 
    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;
    static constexpr amrex::Real c01 = 0.0;
    static constexpr amrex::Real c02 = -0.5;
    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 / 3.0;
    static constexpr amrex::Real c21 = -1.0;
    static constexpr amrex::Real c22 = -1.0 / 3.0;
 
    static constexpr amrex::Real s00 = 1.0;     // axx
    static constexpr amrex::Real s01 = -2.0;    // bxx
    static constexpr amrex::Real s02 = 1.0;     // 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. / 3.; // azz
    static constexpr amrex::Real s21 = -4.;     // bzz
    static constexpr amrex::Real s22 = 4. / 3.; // czz
 
    static constexpr amrex::Real f00 = 0.25;
    static constexpr amrex::Real f01 = 0.5;
    static constexpr amrex::Real f02 = 0.25;
    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;
    static constexpr amrex::Real f21 = 0.625;
    static constexpr amrex::Real f22 = 0.125;
 
    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 / 3;
    static constexpr amrex::Real c01 = 1.0;
    static constexpr amrex::Real c02 = -4.0 / 3;
    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;
    static constexpr amrex::Real c21 = 0;
    static constexpr amrex::Real c22 = -0.5;
 
    static constexpr amrex::Real s00 = 4. / 3.; // axx
    static constexpr amrex::Real s01 = -4;      // bxx
    static constexpr amrex::Real s02 = 8. / 3.; // 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;       // azz
    static constexpr amrex::Real s21 = -2;      // bzz
    static constexpr amrex::Real s22 = 1;       // czz
 
    static constexpr amrex::Real f00 = 0.125;
    static constexpr amrex::Real f01 = 0.625;
    static constexpr amrex::Real f02 = 0.25;
    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;
    static constexpr amrex::Real f21 = 0.5;
    static constexpr amrex::Real f22 = 0.25;
 
    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;
    static constexpr amrex::Real c01 = 0;
    static constexpr amrex::Real c02 = -0.5;
    static constexpr amrex::Real c10 = 1.0 / 3;
    static constexpr amrex::Real c11 = 1.0;
    static constexpr amrex::Real c12 = -4.0 / 3;
    static constexpr amrex::Real c20 = c10;
    static constexpr amrex::Real c21 = c11;
    static constexpr amrex::Real c22 = c12;
 
    static constexpr amrex::Real s00 = 1;       // axx
    static constexpr amrex::Real s01 = -2;      // bxx
    static constexpr amrex::Real s02 = 1;       // cxx
    static constexpr amrex::Real s10 = 4. / 3.; // ayy
    static constexpr amrex::Real s11 = -4;      // byy
    static constexpr amrex::Real s12 = 8. / 3.; // 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;
    static constexpr amrex::Real f01 = 0.5;
    static constexpr amrex::Real f02 = 0.25;
    static constexpr amrex::Real f10 = 0.125;
    static constexpr amrex::Real f11 = 0.625;
    static constexpr amrex::Real f12 = 0.25;
    static constexpr amrex::Real f20 = f10;
    static constexpr amrex::Real f21 = f11;
    static constexpr amrex::Real f22 = f12;
 
    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 / 3;
    static constexpr amrex::Real c01 = 1.0;
    static constexpr amrex::Real c02 = -4.0 / 3;
    static constexpr amrex::Real c10 = 0.5;
    static constexpr amrex::Real c11 = 0;
    static constexpr amrex::Real c12 = -0.5;
    static constexpr amrex::Real c20 = c00;
    static constexpr amrex::Real c21 = c01;
    static constexpr amrex::Real c22 = c02;
 
    static constexpr amrex::Real s00 = 4. / 3.; // axx
    static constexpr amrex::Real s01 = -4;      // bxx
    static constexpr amrex::Real s02 = 8. / 3.; // 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 = s00;     // azz
    static constexpr amrex::Real s21 = s01;     // bzz
    static constexpr amrex::Real s22 = s02;     // czz
 
    static constexpr amrex::Real f00 = 0.125;
    static constexpr amrex::Real f01 = 0.625;
    static constexpr amrex::Real f02 = 0.25;
    static constexpr amrex::Real f10 = 0.25;
    static constexpr amrex::Real f11 = 0.5;
    static constexpr amrex::Real f12 = 0.25;
    static constexpr amrex::Real f20 = f00;
    static constexpr amrex::Real f21 = f01;
    static constexpr amrex::Real f22 = f02;
 
    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 / 3;
    static constexpr amrex::Real c01 = -1;
    static constexpr amrex::Real c02 = -1 / 3.;
    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;
    static constexpr amrex::Real c21 = 0.0;
    static constexpr amrex::Real c22 = -0.5;
 
    static constexpr amrex::Real s00 = 8. / 3.; // axx
    static constexpr amrex::Real s01 = -4.;     // bxx
    static constexpr amrex::Real s02 = 4. / 3.; // 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.;      // azz
    static constexpr amrex::Real s21 = -2.;     // bzz
    static constexpr amrex::Real s22 = 1.;      // czz
 
    static constexpr amrex::Real f00 = 0.25;
    static constexpr amrex::Real f01 = 0.625;
    static constexpr amrex::Real f02 = 0.125;
    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;
    static constexpr amrex::Real f21 = 0.5;
    static constexpr amrex::Real f22 = 0.25;
 
    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;
    static constexpr amrex::Real c01 = 0.0;
    static constexpr amrex::Real c02 = -0.5;
    static constexpr amrex::Real c10 = 4. / 3.;
    static constexpr amrex::Real c11 = -1.;
    static constexpr amrex::Real c12 = -1. / 3.;
    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;     // axx
    static constexpr amrex::Real s01 = -2.0;    // bxx
    static constexpr amrex::Real s02 = 1.0;     // cxx
    static constexpr amrex::Real s10 = 8. / 3.; // ayy
    static constexpr amrex::Real s11 = -4.;     // byy
    static constexpr amrex::Real s12 = 4. / 3.; // 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;
    static constexpr amrex::Real f01 = 0.5;
    static constexpr amrex::Real f02 = 0.25;
    static constexpr amrex::Real f10 = 0.25;
    static constexpr amrex::Real f11 = 0.625;
    static constexpr amrex::Real f12 = 0.125;
    static constexpr amrex::Real f20 = f10;
    static constexpr amrex::Real f21 = f11;
    static constexpr amrex::Real f22 = f12;
 
    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 / 3;
    static constexpr amrex::Real c01 = -1;
    static constexpr amrex::Real c02 = -1 / 3.;
    static constexpr amrex::Real c10 = 0.5;
    static constexpr amrex::Real c11 = 0.0;
    static constexpr amrex::Real c12 = -0.5;
    static constexpr amrex::Real c20 = c00;
    static constexpr amrex::Real c21 = c01;
    static constexpr amrex::Real c22 = c02;
 
    static constexpr amrex::Real s00 = 8. / 3.; // axx
    static constexpr amrex::Real s01 = -4.;     // bxx
    static constexpr amrex::Real s02 = 4. / 3.; // 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 = s00;     // azz
    static constexpr amrex::Real s21 = s01;     // bzz
    static constexpr amrex::Real s22 = s02;     // czz
 
    static constexpr amrex::Real f00 = 0.25;
    static constexpr amrex::Real f01 = 0.625;
    static constexpr amrex::Real f02 = 0.125;
    static constexpr amrex::Real f10 = 0.25;
    static constexpr amrex::Real f11 = 0.5;
    static constexpr amrex::Real f12 = 0.25;
    static constexpr amrex::Real f20 = f00;
    static constexpr amrex::Real f21 = f01;
    static constexpr amrex::Real f22 = f02;
 
    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 / 3;
    static constexpr amrex::Real c01 = 1.0;
    static constexpr amrex::Real c02 = -4.0 / 3;
    static constexpr amrex::Real c10 = 4. / 3;
    static constexpr amrex::Real c11 = -1;
    static constexpr amrex::Real c12 = -1. / 3;
    static constexpr amrex::Real c20 = 0.5;
    static constexpr amrex::Real c21 = 0;
    static constexpr amrex::Real c22 = -0.5;
 
    static constexpr amrex::Real s00 = 4. / 3.; // axx
    static constexpr amrex::Real s01 = -4;      // bxx
    static constexpr amrex::Real s02 = 8. / 3.; // cxx
    static constexpr amrex::Real s10 = 8. / 3.; // ayy
    static constexpr amrex::Real s11 = -4.;     // byy
    static constexpr amrex::Real s12 = 4. / 3.; // cyy
    static constexpr amrex::Real s20 = 1;       // azz
    static constexpr amrex::Real s21 = -2;      // bzz
    static constexpr amrex::Real s22 = 1;       // czz
 
    static constexpr amrex::Real f00 = 0.125;
    static constexpr amrex::Real f01 = 0.625;
    static constexpr amrex::Real f02 = 0.25;
    static constexpr amrex::Real f10 = 0.25;
    static constexpr amrex::Real f11 = 0.625;
    static constexpr amrex::Real f12 = 0.125;
    static constexpr amrex::Real f20 = 0.25;
    static constexpr amrex::Real f21 = 0.5;
    static constexpr amrex::Real f22 = 0.25;
 
    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;
    static constexpr amrex::Real c01 = 0;
    static constexpr amrex::Real c02 = -0.5;
    static constexpr amrex::Real c10 = 1.0 / 3;
    static constexpr amrex::Real c11 = 1.0;
    static constexpr amrex::Real c12 = -4.0 / 3;
    static constexpr amrex::Real c20 = 4. / 3.;
    static constexpr amrex::Real c21 = -1.;
    static constexpr amrex::Real c22 = -1. / 3.;
 
    static constexpr amrex::Real s00 = 1;       // axx
    static constexpr amrex::Real s01 = -2;      // bxx
    static constexpr amrex::Real s02 = 1;       // cxx
    static constexpr amrex::Real s10 = 4. / 3.; // ayy
    static constexpr amrex::Real s11 = -4;      // byy
    static constexpr amrex::Real s12 = 8. / 3.; // cyy
    static constexpr amrex::Real s20 = 8. / 3.; // azz
    static constexpr amrex::Real s21 = -4.;     // bzz
    static constexpr amrex::Real s22 = 4. / 3.; // czz
 
    static constexpr amrex::Real f00 = 0.25;
    static constexpr amrex::Real f01 = 0.5;
    static constexpr amrex::Real f02 = 0.25;
    static constexpr amrex::Real f10 = 0.125;
    static constexpr amrex::Real f11 = 0.625;
    static constexpr amrex::Real f12 = 0.25;
    static constexpr amrex::Real f20 = 0.25;
    static constexpr amrex::Real f21 = 0.625;
    static constexpr amrex::Real f22 = 0.125;
 
    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 / 3;
    static constexpr amrex::Real c01 = 1.0;
    static constexpr amrex::Real c02 = -4.0 / 3;
    static constexpr amrex::Real c10 = 0.5;
    static constexpr amrex::Real c11 = 0;
    static constexpr amrex::Real c12 = -0.5;
    static constexpr amrex::Real c20 = 4. / 3.;
    static constexpr amrex::Real c21 = -1.;
    static constexpr amrex::Real c22 = -1. / 3.;
 
    static constexpr amrex::Real s00 = 4. / 3.; // axx
    static constexpr amrex::Real s01 = -4;      // bxx
    static constexpr amrex::Real s02 = 8. / 3.; // 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. / 3.; // azz
    static constexpr amrex::Real s21 = -4.;     // bzz
    static constexpr amrex::Real s22 = 4. / 3.; // czz
 
    static constexpr amrex::Real f00 = 0.125;
    static constexpr amrex::Real f01 = 0.625;
    static constexpr amrex::Real f02 = 0.25;
    static constexpr amrex::Real f10 = 0.25;
    static constexpr amrex::Real f11 = 0.5;
    static constexpr amrex::Real f12 = 0.25;
    static constexpr amrex::Real f20 = 0.25;
    static constexpr amrex::Real f21 = 0.625;
    static constexpr amrex::Real f22 = 0.125;
 
    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 / 3;
    static constexpr amrex::Real c01 = -1;
    static constexpr amrex::Real c02 = -1 / 3.;
    static constexpr amrex::Real c10 = 1. / 3;
    static constexpr amrex::Real c11 = 1;
    static constexpr amrex::Real c12 = -4. / 3.;
    static constexpr amrex::Real c20 = 0.5;
    static constexpr amrex::Real c21 = 0.0;
    static constexpr amrex::Real c22 = -0.5;
 
    static constexpr amrex::Real s00 = 8. / 3.; // axx
    static constexpr amrex::Real s01 = -4.;     // bxx
    static constexpr amrex::Real s02 = 4. / 3.; // cxx
    static constexpr amrex::Real s10 = 4. / 3.; // ayy
    static constexpr amrex::Real s11 = -4.;     // byy
    static constexpr amrex::Real s12 = 8. / 3.; // cyy
    static constexpr amrex::Real s20 = 1.;      // azz
    static constexpr amrex::Real s21 = -2.;     // bzz
    static constexpr amrex::Real s22 = 1.;      // czz
 
    static constexpr amrex::Real f00 = 0.25;
    static constexpr amrex::Real f01 = 0.625;
    static constexpr amrex::Real f02 = 0.125;
    static constexpr amrex::Real f10 = 0.125;
    static constexpr amrex::Real f11 = 0.625;
    static constexpr amrex::Real f12 = 0.25;
    static constexpr amrex::Real f20 = 0.25;
    static constexpr amrex::Real f21 = 0.5;
    static constexpr amrex::Real f22 = 0.25;
 
    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;
    static constexpr amrex::Real c01 = 0.0;
    static constexpr amrex::Real c02 = -0.5;
    static constexpr amrex::Real c10 = 4. / 3.;
    static constexpr amrex::Real c11 = -1.;
    static constexpr amrex::Real c12 = -1. / 3.;
    static constexpr amrex::Real c20 = 1. / 3.;
    static constexpr amrex::Real c21 = 1.;
    static constexpr amrex::Real c22 = -4. / 3.;
 
    static constexpr amrex::Real s00 = 1.0;     // axx
    static constexpr amrex::Real s01 = -2.0;    // bxx
    static constexpr amrex::Real s02 = 1.0;     // cxx
    static constexpr amrex::Real s10 = 8. / 3.; // ayy
    static constexpr amrex::Real s11 = -4.;     // byy
    static constexpr amrex::Real s12 = 4. / 3.; // cyy
    static constexpr amrex::Real s20 = 4. / 3.; // azz
    static constexpr amrex::Real s21 = -4.;     // bzz
    static constexpr amrex::Real s22 = 8. / 3.; // czz
 
    static constexpr amrex::Real f00 = 0.25;
    static constexpr amrex::Real f01 = 0.5;
    static constexpr amrex::Real f02 = 0.25;
    static constexpr amrex::Real f10 = 0.25;
    static constexpr amrex::Real f11 = 0.625;
    static constexpr amrex::Real f12 = 0.125;
    static constexpr amrex::Real f20 = 0.125;
    static constexpr amrex::Real f21 = 0.625;
    static constexpr amrex::Real f22 = 0.25;
 
    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 / 3.0;
    static constexpr amrex::Real c01 = -1.0;
    static constexpr amrex::Real c02 = -1.0 / 3.0;
    static constexpr amrex::Real c10 = 0.5;
    static constexpr amrex::Real c11 = 0.0;
    static constexpr amrex::Real c12 = -0.5;
    static constexpr amrex::Real c20 = 1. / 3.;
    static constexpr amrex::Real c21 = 1.;
    static constexpr amrex::Real c22 = -4. / 3.;
 
    static constexpr amrex::Real s00 = 8. / 3.;   // axx
    static constexpr amrex::Real s01 = -4.;       // bxx
    static constexpr amrex::Real s02 = 4. / 3.;   // 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 = 4.0 / 3.0; // azz
    static constexpr amrex::Real s21 = -4.0;      // bzz
    static constexpr amrex::Real s22 = 8.0 / 3.0; // czz
 
    static constexpr amrex::Real f00 = 0.25;
    static constexpr amrex::Real f01 = 0.625;
    static constexpr amrex::Real f02 = 0.125;
    static constexpr amrex::Real f10 = 0.25;
    static constexpr amrex::Real f11 = 0.5;
    static constexpr amrex::Real f12 = 0.25;
    static constexpr amrex::Real f20 = 0.125;
    static constexpr amrex::Real f21 = 0.625;
    static constexpr amrex::Real f22 = 0.25;
 
    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 / 3;
    static constexpr amrex::Real c01 = 1.0;
    static constexpr amrex::Real c02 = -4.0 / 3;
    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. / 3.; // axx
    static constexpr amrex::Real s01 = -4;      // bxx
    static constexpr amrex::Real s02 = 8. / 3.; // 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;
    static constexpr amrex::Real f01 = 0.625;
    static constexpr amrex::Real f02 = 0.25;
    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;
 
    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 / 3;
    static constexpr amrex::Real c01 = 1.0;
    static constexpr amrex::Real c02 = -4.0 / 3;
    static constexpr amrex::Real c10 = 4. / 3.;
    static constexpr amrex::Real c11 = -1.;
    static constexpr amrex::Real c12 = -1. / 3.;
    static constexpr amrex::Real c20 = c00;
    static constexpr amrex::Real c21 = c01;
    static constexpr amrex::Real c22 = c02;
 
    static constexpr amrex::Real s00 = 4. / 3.; // axx
    static constexpr amrex::Real s01 = -4;      // bxx
    static constexpr amrex::Real s02 = 8. / 3.; // cxx
    static constexpr amrex::Real s10 = 8. / 3.; // ayy
    static constexpr amrex::Real s11 = -4.;     // byy
    static constexpr amrex::Real s12 = 4. / 3.; // 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;
    static constexpr amrex::Real f01 = 0.625;
    static constexpr amrex::Real f02 = 0.25;
    static constexpr amrex::Real f10 = 0.25;
    static constexpr amrex::Real f11 = 0.625;
    static constexpr amrex::Real f12 = 0.125;
    static constexpr amrex::Real f20 = f00;
    static constexpr amrex::Real f21 = f01;
    static constexpr amrex::Real f22 = f02;
 
    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. / 3.;
    static constexpr amrex::Real c01 = -1;
    static constexpr amrex::Real c02 = -1. / 3.;
    static constexpr amrex::Real c10 = 1. / 3;
    static constexpr amrex::Real c11 = 1;
    static constexpr amrex::Real c12 = -4. / 3.;
    static constexpr amrex::Real c20 = c10;
    static constexpr amrex::Real c21 = c11;
    static constexpr amrex::Real c22 = c12;
 
    static constexpr amrex::Real s00 = 8. / 3.; // axx
    static constexpr amrex::Real s01 = -4.;     // bxx
    static constexpr amrex::Real s02 = 4. / 3.; // cxx
    static constexpr amrex::Real s10 = 4. / 3.; // ayy
    static constexpr amrex::Real s11 = -4.;     // byy
    static constexpr amrex::Real s12 = 8. / 3.; // 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;
    static constexpr amrex::Real f01 = 0.625;
    static constexpr amrex::Real f02 = 0.125;
    static constexpr amrex::Real f10 = 0.125;
    static constexpr amrex::Real f11 = 0.625;
    static constexpr amrex::Real f12 = 0.25;
    static constexpr amrex::Real f20 = f10;
    static constexpr amrex::Real f21 = f11;
    static constexpr amrex::Real f22 = f12;
 
    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. / 3.;
    static constexpr amrex::Real c01 = -1;
    static constexpr amrex::Real c02 = -1. / 3.;
    static constexpr amrex::Real c10 = 4. / 3.;
    static constexpr amrex::Real c11 = -1.;
    static constexpr amrex::Real c12 = -1. / 3.;
    static constexpr amrex::Real c20 = 1.0 / 3.0;
    static constexpr amrex::Real c21 = 1.0;
    static constexpr amrex::Real c22 = -4.0 / 3.;
 
    static constexpr amrex::Real s00 = 8. / 3.;   // axx
    static constexpr amrex::Real s01 = -4.;       // bxx
    static constexpr amrex::Real s02 = 4. / 3.;   // cxx
    static constexpr amrex::Real s10 = 8. / 3.;   // ayy
    static constexpr amrex::Real s11 = -4.;       // byy
    static constexpr amrex::Real s12 = 4. / 3.;   // cyy
    static constexpr amrex::Real s20 = 4.0 / 3.0; // azz
    static constexpr amrex::Real s21 = -4.0;      // bzz
    static constexpr amrex::Real s22 = 8.0 / 3.0; // czz
 
    static constexpr amrex::Real f00 = 0.25;
    static constexpr amrex::Real f01 = 0.625;
    static constexpr amrex::Real f02 = 0.125;
    static constexpr amrex::Real f10 = 0.25;
    static constexpr amrex::Real f11 = 0.625;
    static constexpr amrex::Real f12 = 0.125;
    static constexpr amrex::Real f20 = 0.125;
    static constexpr amrex::Real f21 = 0.625;
    static constexpr amrex::Real f22 = 0.25;
 
    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 / 3;
    static constexpr amrex::Real c01 = 1.0;
    static constexpr amrex::Real c02 = -4.0 / 3;
    static constexpr amrex::Real c10 = c00;
    static constexpr amrex::Real c11 = c01;
    static constexpr amrex::Real c12 = c02;
    static constexpr amrex::Real c20 = 4. / 3.;
    static constexpr amrex::Real c21 = -1.;
    static constexpr amrex::Real c22 = -1. / 3.;
 
    static constexpr amrex::Real s00 = 4. / 3.; // axx
    static constexpr amrex::Real s01 = -4;      // bxx
    static constexpr amrex::Real s02 = 8. / 3.; // 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. / 3.; // azz
    static constexpr amrex::Real s21 = -4.;     // bzz
    static constexpr amrex::Real s22 = 4. / 3.; // czz
 
    static constexpr amrex::Real f00 = 0.125;
    static constexpr amrex::Real f01 = 0.625;
    static constexpr amrex::Real f02 = 0.25;
    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;
    static constexpr amrex::Real f21 = 0.625;
    static constexpr amrex::Real f22 = 0.125;
 
    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 / 3;
    static constexpr amrex::Real c01 = 1.0;
    static constexpr amrex::Real c02 = -4.0 / 3;
    static constexpr amrex::Real c10 = 4. / 3.;
    static constexpr amrex::Real c11 = -1.;
    static constexpr amrex::Real c12 = -1. / 3.;
    static constexpr amrex::Real c20 = 4. / 3.;
    static constexpr amrex::Real c21 = -1.;
    static constexpr amrex::Real c22 = -1. / 3.;
 
    static constexpr amrex::Real s00 = 4. / 3.; // axx
    static constexpr amrex::Real s01 = -4;      // bxx
    static constexpr amrex::Real s02 = 8. / 3.; // cxx
    static constexpr amrex::Real s10 = 8. / 3.; // ayy
    static constexpr amrex::Real s11 = -4.;     // byy
    static constexpr amrex::Real s12 = 4. / 3.; // cyy
    static constexpr amrex::Real s20 = 8. / 3.; // azz
    static constexpr amrex::Real s21 = -4.;     // bzz
    static constexpr amrex::Real s22 = 4. / 3.; // czz
 
    static constexpr amrex::Real f00 = 0.125;
    static constexpr amrex::Real f01 = 0.625;
    static constexpr amrex::Real f02 = 0.25;
    static constexpr amrex::Real f10 = 0.25;
    static constexpr amrex::Real f11 = 0.625;
    static constexpr amrex::Real f12 = 0.125;
    static constexpr amrex::Real f20 = 0.25;
    static constexpr amrex::Real f21 = 0.625;
    static constexpr amrex::Real f22 = 0.125;
 
    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. / 3.;
    static constexpr amrex::Real c01 = -1;
    static constexpr amrex::Real c02 = -1. / 3.;
    static constexpr amrex::Real c10 = 1. / 3;
    static constexpr amrex::Real c11 = 1;
    static constexpr amrex::Real c12 = -4. / 3.;
    static constexpr amrex::Real c20 = 4. / 3.;
    static constexpr amrex::Real c21 = -1.;
    static constexpr amrex::Real c22 = -1. / 3.;
 
    static constexpr amrex::Real s00 = 8. / 3.; // axx
    static constexpr amrex::Real s01 = -4.;     // bxx
    static constexpr amrex::Real s02 = 4. / 3.; // cxx
    static constexpr amrex::Real s10 = 4. / 3.; // ayy
    static constexpr amrex::Real s11 = -4.;     // byy
    static constexpr amrex::Real s12 = 8. / 3.; // cyy
    static constexpr amrex::Real s20 = 8. / 3.; // azz
    static constexpr amrex::Real s21 = -4.;     // bzz
    static constexpr amrex::Real s22 = 4. / 3.; // czz
 
    static constexpr amrex::Real f00 = 0.25;
    static constexpr amrex::Real f01 = 0.625;
    static constexpr amrex::Real f02 = 0.125;
    static constexpr amrex::Real f10 = 0.125;
    static constexpr amrex::Real f11 = 0.625;
    static constexpr amrex::Real f12 = 0.25;
    static constexpr amrex::Real f20 = 0.25;
    static constexpr amrex::Real f21 = 0.625;
    static constexpr amrex::Real f22 = 0.125;
 
    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. / 3.;
    static constexpr amrex::Real c01 = -1;
    static constexpr amrex::Real c02 = -1. / 3.;
    static constexpr amrex::Real c10 = 4. / 3.;
    static constexpr amrex::Real c11 = -1.;
    static constexpr amrex::Real c12 = -1. / 3.;
    static constexpr amrex::Real c20 = 4. / 3.;
    static constexpr amrex::Real c21 = -1.;
    static constexpr amrex::Real c22 = -1. / 3.;
 
    static constexpr amrex::Real s00 = 8. / 3.; // axx
    static constexpr amrex::Real s01 = -4.;     // bxx
    static constexpr amrex::Real s02 = 4. / 3.; // cxx
    static constexpr amrex::Real s10 = 8. / 3.; // ayy
    static constexpr amrex::Real s11 = -4.;     // byy
    static constexpr amrex::Real s12 = 4. / 3.; // cyy
    static constexpr amrex::Real s20 = 8. / 3.; // azz
    static constexpr amrex::Real s21 = -4.;     // bzz
    static constexpr amrex::Real s22 = 4. / 3.; // czz
 
    static constexpr amrex::Real f00 = 0.25;
    static constexpr amrex::Real f01 = 0.625;
    static constexpr amrex::Real f02 = 0.125;
    static constexpr amrex::Real f10 = 0.25;
    static constexpr amrex::Real f11 = 0.625;
    static constexpr amrex::Real f12 = 0.125;
    static constexpr amrex::Real f20 = 0.25;
    static constexpr amrex::Real f21 = 0.625;
    static constexpr amrex::Real f22 = 0.125;
 
    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 amr_wind::fvm::stencil
 
#endif /* STENCILS_H */