UPWIND3 Struct Reference

#include <ERF_Interpolation_UPW.H>

Collaboration diagram for UPWIND3:

Public Member Functions
	UPWIND3 (const amrex::Array4< const amrex::Real > &phi, const amrex::Real upw_frac)
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 &sp1, const amrex::Real &s, const amrex::Real &sm1, const amrex::Real &sm2, const amrex::Real &upw) const
int	GetUpwindCellNumber () const
void	SetUpwinding (amrex::Real upw_frac)

Private Attributes
amrex::Array4< const amrex::Real >	m_phi
amrex::Real	m_upw_frac

Static Private Attributes
static constexpr amrex::Real	g1 =(7.0/12.0)
static constexpr amrex::Real	g2 =(1.0/12.0)

Detailed Description

Interpolation operators used for 3rd order upwind scheme

Constructor & Destructor Documentation

UPWIND3()

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

inline

        : m_phi(phi)
        , m_upw_frac(upw_frac)
    {}

Member Function Documentation

Evaluate()

AMREX_GPU_DEVICE AMREX_FORCE_INLINE amrex::Real UPWIND3::Evaluate ( const amrex::Real &  sp1, const amrex::Real &  s, const amrex::Real &  sm1, const amrex::Real &  sm2, const amrex::Real &  upw  ) const

inline

    {
        // Averages and diffs
        amrex::Real a1 = (s   + sm1);
        amrex::Real d1 = (s   - sm1);
        amrex::Real a2 = (sp1 + sm2);
        amrex::Real d2 = (sp1 - sm2);
 
        // Interpolated value
        return ( g1*a1 - g2*a2 + upw * g2 * (d2 - 3.0*d1) );
    }

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

Here is the caller graph for this function:

GetUpwindCellNumber()

int UPWIND3::GetUpwindCellNumber ( ) const

inline

{return 2;}

InterpolateInX()

AMREX_GPU_DEVICE AMREX_FORCE_INLINE void UPWIND3::InterpolateInX ( const int &  i, const int &  j, const int &  k, const int &  qty_index, amrex::Real &  val_lo, amrex::Real  upw_lo  ) 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);
 
        // Upwinding flags
        if (upw_lo != 0.) upw_lo = (upw_lo > 0) ? 1. : -1.;
 
        // Add blending
        upw_lo *= m_upw_frac;
 
        // Interpolate lo
        val_lo = Evaluate(sp1,s,sm1,sm2,upw_lo);
    }

Referenced by EBAdvectionSrcForScalarsWrapper().

Here is the call graph for this function:

Here is the caller graph for this function:

InterpolateInY()

AMREX_GPU_DEVICE AMREX_FORCE_INLINE void UPWIND3::InterpolateInY ( const int &  i, const int &  j, const int &  k, const int &  qty_index, amrex::Real &  val_lo, amrex::Real  upw_lo  ) 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);
 
        // Upwinding flags
        if (upw_lo != 0.) upw_lo = (upw_lo > 0) ? 1. : -1.;
 
        // Add blending
        upw_lo *= m_upw_frac;
 
        // Interpolate lo
        val_lo = Evaluate(sp1,s,sm1,sm2,upw_lo);
    }

Referenced by EBAdvectionSrcForScalarsWrapper().

Here is the call graph for this function:

Here is the caller graph for this function:

InterpolateInZ()

AMREX_GPU_DEVICE AMREX_FORCE_INLINE void UPWIND3::InterpolateInZ ( const int &  i, const int &  j, const int &  k, const int &  qty_index, amrex::Real &  val_lo, amrex::Real  upw_lo  ) 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);
 
        // Upwinding flags
        if (upw_lo != 0.) upw_lo = (upw_lo > 0) ? 1. : -1.;
 
        // Add blending
        upw_lo *= m_upw_frac;
 
        // Interpolate lo
        val_lo = Evaluate(sp1,s,sm1,sm2,upw_lo);
    }

Referenced by EBAdvectionSrcForScalarsWrapper().

Here is the call graph for this function:

Here is the caller graph for this function:

SetUpwinding()

void UPWIND3::SetUpwinding ( amrex::Real  upw_frac )

inline

{m_upw_frac = upw_frac;}

Referenced by EBAdvectionSrcForScalarsWrapper().

Here is the caller graph for this function:

Member Data Documentation

g1

constexpr amrex::Real UPWIND3::g1 =(7.0/12.0)

staticconstexprprivate

Referenced by Evaluate().

g2

constexpr amrex::Real UPWIND3::g2 =(1.0/12.0)

staticconstexprprivate

Referenced by Evaluate().

m_phi

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

private

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

m_upw_frac

amrex::Real UPWIND3::m_upw_frac

private

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

---

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

• Source/Utils/ERF_Interpolation_UPW.H