WENO_MZQ3 Struct Reference

#include <ERF_Interpolation_WENO_Z.H>

Collaboration diagram for WENO_MZQ3:

Public Member Functions
	WENO_MZQ3 (const amrex::Array4< const amrex::Real > &phi)
AMREX_GPU_DEVICE AMREX_FORCE_INLINE void	InterpolateInX (const int &i, const int &j, const int &k, const int &qty_index, amrex::Real &val_lo, amrex::Real upw_lo, const amrex::Real) const
AMREX_GPU_DEVICE AMREX_FORCE_INLINE void	InterpolateInY (const int &i, const int &j, const int &k, const int &qty_index, amrex::Real &val_lo, amrex::Real upw_lo, const amrex::Real) const
AMREX_GPU_DEVICE AMREX_FORCE_INLINE void	InterpolateInZ (const int &i, const int &j, const int &k, const int &qty_index, amrex::Real &val_lo, amrex::Real upw_lo, const amrex::Real) const
AMREX_GPU_DEVICE AMREX_FORCE_INLINE amrex::Real	Evaluate (const amrex::Real &sm1, const amrex::Real &s, const amrex::Real &sp1) const

Private Attributes
amrex::Array4< const amrex::Real >	m_phi
const amrex::Real	eps =1.0e-6
const amrex::Real	tol =1.0e-12

Static Private Attributes
static constexpr amrex::Real	g1 =(1.0/3.0)
static constexpr amrex::Real	g2 =(1.0/3.0)
static constexpr amrex::Real	g3 =(1.0/3.0)

Detailed Description

Interpolation operators used for WENO_MZQ3 scheme

Constructor & Destructor Documentation

WENO_MZQ3()

WENO_MZQ3::WENO_MZQ3 ( const amrex::Array4< const amrex::Real > &  phi )

inline

        : m_phi(phi) {}

Member Function Documentation

Evaluate()

AMREX_GPU_DEVICE AMREX_FORCE_INLINE amrex::Real WENO_MZQ3::Evaluate ( const amrex::Real &  sm1, const amrex::Real &  s, const amrex::Real &  sp1  ) const

inline

    {
        // Smoothing factors
        amrex::Real b1 = (s - sm1) * (s - sm1);
        amrex::Real b2 = (sp1 - s) * (sp1 - s);
        amrex::Real b3 = ( (13.0 / 12.0) * ((sm1 - 2.0*s + sp1)*(sm1 - 2.0*s + sp1)) ) + ( ((sp1 - sm1)*(sp1 - sm1)) / 4.0 );
 
        // Weight factors
        amrex::Real t5 = ( std::abs(b3 - b1) + std::abs(b3 - b2) ) / 32.0;
        amrex::Real a1 = g1 * ( 1.0 + (t5*t5) / ((eps + b1) * (eps + b1)) );
        amrex::Real a2 = g2 * ( 1.0 + (t5*t5) / ((eps + b2) * (eps + b2)) );
        amrex::Real a3 = g3 * ( 1.0 + (t5*t5) / ((eps + b3) * (eps + b3)) );
        amrex::Real asum = a1 + a2 + a3;
        amrex::Real w1 = a1 / asum;
        amrex::Real w2 = a2 / asum;
        amrex::Real w3 = a3 / asum;
 
        // Taylor expansions
        amrex::Real v1 = (-sm1 + 3.0 * s) / 2.0;
        amrex::Real v2 = (  s  +  sp1) / 2.0;
        amrex::Real v3 = (-sm1 + 5.0 * s + 2.0 * sp1) / 6.0;
 
        // Interpolated value
        return ( (w3 / g3) * (v3 - g1 * v1 - g2 * v2) + w1 * v1 + w2 * v2 );
    }

Referenced by InterpolateInX(), InterpolateInY(), and InterpolateInZ().

Here is the caller graph for this function:

InterpolateInX()

AMREX_GPU_DEVICE AMREX_FORCE_INLINE void WENO_MZQ3::InterpolateInX ( const int &  i, const int &  j, const int &  k, const int &  qty_index, amrex::Real &  val_lo, amrex::Real  upw_lo, const amrex::Real    ) const

inline

    {
        // Data to interpolate on
        amrex::Real sp1 = m_phi(i+1, j  , k  , qty_index);
        amrex::Real s   = m_phi(i  , j  , k  , qty_index);
        amrex::Real sm1 = m_phi(i-1, j  , k  , qty_index);
        amrex::Real sm2 = m_phi(i-2, j  , k  , qty_index);
 
        if (upw_lo > tol) {
            val_lo = Evaluate(sm2,sm1,s  );
        } else if (upw_lo < -tol) {
            val_lo = Evaluate(sp1,s  ,sm1);
        } else {
            val_lo = 0.5 * (s + sm1);
        }
    }

Here is the call graph for this function:

InterpolateInY()

AMREX_GPU_DEVICE AMREX_FORCE_INLINE void WENO_MZQ3::InterpolateInY ( const int &  i, const int &  j, const int &  k, const int &  qty_index, amrex::Real &  val_lo, amrex::Real  upw_lo, const amrex::Real    ) const

inline

    {
        // Data to interpolate on
        amrex::Real sp1 = m_phi(i  , j+1, k  , qty_index);
        amrex::Real s   = m_phi(i  , j  , k  , qty_index);
        amrex::Real sm1 = m_phi(i  , j-1, k  , qty_index);
        amrex::Real sm2 = m_phi(i  , j-2, k  , qty_index);
 
        if (upw_lo > tol) {
            val_lo = Evaluate(sm2,sm1,s  );
        } else if (upw_lo < -tol) {
            val_lo = Evaluate(sp1,s  ,sm1);
        } else {
            val_lo = 0.5 * (s + sm1);
        }
    }

Here is the call graph for this function:

InterpolateInZ()

AMREX_GPU_DEVICE AMREX_FORCE_INLINE void WENO_MZQ3::InterpolateInZ ( const int &  i, const int &  j, const int &  k, const int &  qty_index, amrex::Real &  val_lo, amrex::Real  upw_lo, const amrex::Real    ) const

inline

    {
        // Data to interpolate on
        amrex::Real sp1 = m_phi(i  , j  , k+1, qty_index);
        amrex::Real s   = m_phi(i  , j  , k  , qty_index);
        amrex::Real sm1 = m_phi(i  , j  , k-1, qty_index);
        amrex::Real sm2 = m_phi(i  , j  , k-2, qty_index);
 
        if (upw_lo > tol) {
            val_lo = Evaluate(sm2,sm1,s  );
        } else if (upw_lo < -tol) {
            val_lo = Evaluate(sp1,s  ,sm1);
        } else {
            val_lo = 0.5 * (s + sm1);
        }
    }

Here is the call graph for this function:

Member Data Documentation

eps

const amrex::Real WENO_MZQ3::eps =1.0e-6

private

Referenced by Evaluate().

g1

constexpr amrex::Real WENO_MZQ3::g1 =(1.0/3.0)

staticconstexprprivate

Referenced by Evaluate().

g2

constexpr amrex::Real WENO_MZQ3::g2 =(1.0/3.0)

staticconstexprprivate

Referenced by Evaluate().

g3

constexpr amrex::Real WENO_MZQ3::g3 =(1.0/3.0)

staticconstexprprivate

Referenced by Evaluate().

m_phi

amrex::Array4<const amrex::Real> WENO_MZQ3::m_phi

private

Referenced by InterpolateInX(), InterpolateInY(), and InterpolateInZ().

tol

const amrex::Real WENO_MZQ3::tol =1.0e-12

private

Referenced by InterpolateInX(), InterpolateInY(), and InterpolateInZ().

---

The documentation for this struct was generated from the following file:

• Source/Utils/ERF_Interpolation_WENO_Z.H