UPWIND3SL Struct Reference

#include <ERF_Interpolation_UPW.H>

Collaboration diagram for UPWIND3SL:

Public Member Functions
	UPWIND3SL (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 &s, const amrex::Real &sm1, const amrex::Real &sm2) const
int	GetUpwindCellNumber () const

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

Static Private Attributes
static constexpr amrex::Real	eps =(amrex::Real(1.0e-16))
static constexpr amrex::Real	l1 =(amrex::Real(1)/amrex::Real(6.0))
static constexpr amrex::Real	l2 =(amrex::Real(1)/three)

Detailed Description

Interpolation operators used for 3rd order upwind scheme

Constructor & Destructor Documentation

UPWIND3SL()

UPWIND3SL::UPWIND3SL ( 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 UPWIND3SL::Evaluate ( const amrex::Real &  s, const amrex::Real &  sm1, const amrex::Real &  sm2  ) const

inline

    {
        amrex::Real rmh = (s - sm1) / ((sm1 - sm2) + eps);
        amrex::Real K   = l2 + amrex::Real(2)*l2*rmh;
        amrex::Real phi = std::max(amrex::Real(0),std::min(amrex::Real(2)*rmh,std::min(amrex::Real(2),K)));
        return (sm1 + myhalf * phi * (sm1 - sm2));
    }

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

Here is the caller graph for this function:

GetUpwindCellNumber()

int UPWIND3SL::GetUpwindCellNumber ( ) const

inline

{return 2;}

InterpolateInX()

AMREX_GPU_DEVICE AMREX_FORCE_INLINE void UPWIND3SL::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 != amrex::Real(0)) upw_lo = (upw_lo > 0) ? amrex::Real(1) : -amrex::Real(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(s  ,sm1,sm2);
        amrex::Real Flo = Evaluate(sm1,s  ,sp1);
 
        // Numerical flux
        val_lo = (amrex::Real(1) + upw_lo)/amrex::Real(2) * Fhi + (amrex::Real(1) - upw_lo)/amrex::Real(2) * Flo;
    }

Here is the call graph for this function:

InterpolateInY()

AMREX_GPU_DEVICE AMREX_FORCE_INLINE void UPWIND3SL::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 != amrex::Real(0)) upw_lo = (upw_lo > 0) ? amrex::Real(1) : -amrex::Real(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(s  ,sm1,sm2);
        amrex::Real Flo = Evaluate(sm1,s  ,sp1);
 
        // Numerical flux
        val_lo = (amrex::Real(1) + upw_lo)/amrex::Real(2) * Fhi + (amrex::Real(1) - upw_lo)/amrex::Real(2) * Flo;
    }

Here is the call graph for this function:

InterpolateInZ()

AMREX_GPU_DEVICE AMREX_FORCE_INLINE void UPWIND3SL::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 != amrex::Real(0)) upw_lo = (upw_lo > 0) ? amrex::Real(1) : -amrex::Real(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(s  ,sm1,sm2);
        amrex::Real Flo = Evaluate(sm1,s  ,sp1);
 
        // Numerical flux
        val_lo = (amrex::Real(1) + upw_lo)/amrex::Real(2) * Fhi + (amrex::Real(1) - upw_lo)/amrex::Real(2) * Flo;
    }

Here is the call graph for this function:

Member Data Documentation

eps

constexpr amrex::Real UPWIND3SL::eps =(amrex::Real(1.0e-16))

staticconstexprprivate

Referenced by Evaluate().

l1

constexpr amrex::Real UPWIND3SL::l1 =(amrex::Real(1)/amrex::Real(6.0))

staticconstexprprivate

l2

constexpr amrex::Real UPWIND3SL::l2 =(amrex::Real(1)/three)

staticconstexprprivate

Referenced by Evaluate().

m_phi

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

private

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

---

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

• Source/Utils/ERF_Interpolation_UPW.H