WENO_Z3 Struct Reference

#include <ERF_Interpolation_WENO_Z.H>

Collaboration diagram for WENO_Z3:

Public Member Functions
	WENO_Z3 (const amrex::Array4< const amrex::Real > &phi, const amrex::Real)
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_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_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_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-40

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

Detailed Description

Interpolation operators used for WENO_Z-3 scheme

Constructor & Destructor Documentation

WENO_Z3()

WENO_Z3::WENO_Z3 ( const amrex::Array4< const amrex::Real > &  phi, const amrex::Real    )

inline

        : m_phi(phi) {}

Member Function Documentation

Evaluate()

AMREX_GPU_DEVICE AMREX_FORCE_INLINE amrex::Real WENO_Z3::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);
 
        // Weight factors
        amrex::Real t5 = std::abs(b2 - b1);
        amrex::Real w1 = g1 * ( 1.0 + (t5*t5) / ((eps + b1) * (eps + b1)) );
        amrex::Real w2 = g2 * ( 1.0 + (t5*t5) / ((eps + b2) * (eps + b2)) );
 
        // Weight factor norm
        amrex::Real wsum = w1 + w2;
 
        // Taylor expansions
        amrex::Real v1 = -sm1 + 3.0 * s;
        amrex::Real v2 =  s   + sp1;
 
        // Interpolated value
        return ( (w1 * v1 + w2 * v2) / (2.0 * wsum) );
    }

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

Here is the caller graph for this function:

InterpolateInX()

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

inline

    {
        // Upwinding flags
        if (upw_lo != 0.) upw_lo = (upw_lo > 0) ? 1. : -1.;
 
        // 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);
 
        // Left and right fluxes
        amrex::Real Fhi = Evaluate(sm2,sm1,s  );
        amrex::Real Flo = Evaluate(sp1,s  ,sm1);
 
        // Numerical flux
        val_lo = (1.0 + upw_lo)/2.0 * Fhi + (1.0 - upw_lo)/2.0 * Flo;
    }

Here is the call graph for this function:

InterpolateInY()

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

inline

    {
        // Upwinding flags
        if (upw_lo != 0.) upw_lo = (upw_lo > 0) ? 1. : -1.;
 
        // 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);
 
        // Left and right fluxes
        amrex::Real Fhi = Evaluate(sm2,sm1,s  );
        amrex::Real Flo = Evaluate(sp1,s  ,sm1);
 
        // Numerical flux
        val_lo = (1.0 + upw_lo)/2.0 * Fhi + (1.0 - upw_lo)/2.0 * Flo;
    }

Here is the call graph for this function:

InterpolateInZ()

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

inline

    {
        // Upwinding flags
        if (upw_lo != 0.) upw_lo = (upw_lo > 0) ? 1. : -1.;
 
        // 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);
 
        // Left and right fluxes
        amrex::Real Fhi = Evaluate(sm2,sm1,s  );
        amrex::Real Flo = Evaluate(sp1,s  ,sm1);
 
        // Numerical flux
        val_lo = (1.0 + upw_lo)/2.0 * Fhi + (1.0 - upw_lo)/2.0 * Flo;
    }

Here is the call graph for this function:

Member Data Documentation

eps

const amrex::Real WENO_Z3::eps =1.0e-40

private

Referenced by Evaluate().

g1

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

staticconstexprprivate

Referenced by Evaluate().

g2

constexpr amrex::Real WENO_Z3::g2 =(2.0/3.0)

staticconstexprprivate

Referenced by Evaluate().

m_phi

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

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