ERF
Energy Research and Forecasting: An Atmospheric Modeling Code
•All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Macros Pages
Public Member Functions | Private Attributes | Static Private Attributes | List of all members
WENO_Z7 Struct Reference

#include <ERF_Interpolation_WENO_Z.H>

Collaboration diagram for WENO_Z7:
Collaboration graph
[legend]

Public Member Functions

 WENO_Z7 (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 &sm3, const amrex::Real &sm2, const amrex::Real &sm1, const amrex::Real &s, const amrex::Real &sp1, const amrex::Real &sp2, const amrex::Real &sp3) const
 

Private Attributes

amrex::Array4< const amrex::Real > m_phi
 
const amrex::Real eps =1.0e-40
 
const amrex::Real tol =1.0e-12
 

Static Private Attributes

static constexpr amrex::Real g1 =( 1.0/64.0)
 
static constexpr amrex::Real g2 =(21.0/64.0)
 
static constexpr amrex::Real g3 =(35.0/64.0)
 
static constexpr amrex::Real g4 =( 7.0/64.0)
 

Detailed Description

Interpolation operators used for WENO_Z-7 scheme

Constructor & Destructor Documentation

◆ WENO_Z7()

WENO_Z7::WENO_Z7 ( const amrex::Array4< const amrex::Real > &  phi,
const amrex::Real   
)
inline
391  : m_phi(phi) {}
WENO_Z7::m_phi
amrex::Array4< const amrex::Real > m_phi
Definition: ERF_Interpolation_WENO_Z.H:554

Member Function Documentation

◆ Evaluate()

AMREX_GPU_DEVICE AMREX_FORCE_INLINE amrex::Real WENO_Z7::Evaluate ( const amrex::Real &  sm3,
const amrex::Real &  sm2,
const amrex::Real &  sm1,
const amrex::Real &  s,
const amrex::Real &  sp1,
const amrex::Real &  sp2,
const amrex::Real &  sp3 
) const
inline
490  {
491  // Smoothing factors
492  amrex::Real b1 = ( sm3 * sm3 * 6649./2880.0
493  - sm3 * sm2 * 2623./160.0
494  + sm3 * sm1 * 9449./480.0
495  - sm3 * s * 11389./1440.0
496  + sm2 * sm2 * 28547./960.0
497  - sm2 * sm1 * 35047./480.0
498  + sm2 * s * 14369./480.0
499  + sm1 * sm1 * 44747./960.0
500  - sm1 * s * 6383./160.0
501  + s * s * 25729./2880.0 );
502  amrex::Real b2 = ( sm2 * sm2 * 3169/2880.0
503  - sm2 * sm1 * 3229/480.0
504  + sm2 * s * 3169/480.0
505  - sm2 * sp1 * 2989/1440.0
506  + sm1 * sm1 * 11147/960.0
507  - sm1 * s * 11767/480.0
508  + sm1 * sp1 * 1283/160.0
509  + s * s * 13667/960.0
510  - s * sp1 * 5069/480.0
511  + sp1 * sp1 * 6649/2880.0 );
512  amrex::Real b3 = ( sm1 * sm1 * 6649./2880.0
513  - sm1 * s * 5069./480.0
514  + sm1 * sp1 * 1283./160.0
515  - sm1 * sp2 * 2989./1440.0
516  + s * s * 13667./960.0
517  - s * sp1 * 11767./480.0
518  + s * sp2 * 3169./480.0
519  + sp1 * sp1 * 11147./960.0
520  - sp1 * sp2 * 3229./480.0
521  + sp2 * sp2 * 3169./2880.0 );
522  amrex::Real b4 = ( s * s * 25729./2880.
523  - s * sp1 * 6383./160.
524  + s * sp2 * 14369./480.
525  - s * sp3 * 11389./1440.
526  + sp1 * sp1 * 44747./960.
527  - sp1 * sp2 * 35047./480.
528  + sp1 * sp3 * 9449./480.
529  + sp2 * sp2 * 28547./960.
530  - sp2 * sp3 * 2623./160.
531  + sp3 * sp3 * 6649./2880. );
532 
533  // Weight factors
534  amrex::Real t5 = std::abs(b1 - b2 - b3 + b4);
535  amrex::Real w1 = g1 * ( 1.0 + (t5*t5) / ((eps + b1) * (eps + b1)) );
536  amrex::Real w2 = g2 * ( 1.0 + (t5*t5) / ((eps + b2) * (eps + b2)) );
537  amrex::Real w3 = g3 * ( 1.0 + (t5*t5) / ((eps + b3) * (eps + b3)) );
538  amrex::Real w4 = g4 * ( 1.0 + (t5*t5) / ((eps + b3) * (eps + b3)) );
539 
540  // Weight factor norm
541  amrex::Real wsum = w1 + w2 + w3 + w4;
542 
543  // Taylor expansions
544  amrex::Real v1 = (-0.3125)*sm3 + ( 1.3125)*sm2 + (-2.1875)*sm1 + ( 2.1875)*s;
545  amrex::Real v2 = ( 0.0625)*sm2 + (-0.3125)*sm1 + ( 0.9375)*s + ( 0.3125)*sp1;
546  amrex::Real v3 = (-0.0625)*sm1 + ( 0.5625)*s + ( 0.5625)*sp1 + (-0.0625)*sp2;
547  amrex::Real v4 = ( 0.3125)*s + ( 0.9375)*sp1 + (-0.3125)*sp2 + ( 0.0625)*sp3;
548 
549  // Interpolated value
550  return ( (w1 * v1 + w2 * v2 + w3 * v3 + w4 * v4) / (wsum) );
551  }
WENO_Z7::g1
static constexpr amrex::Real g1
Definition: ERF_Interpolation_WENO_Z.H:557
WENO_Z7::eps
const amrex::Real eps
Definition: ERF_Interpolation_WENO_Z.H:555
WENO_Z7::g4
static constexpr amrex::Real g4
Definition: ERF_Interpolation_WENO_Z.H:560
WENO_Z7::g2
static constexpr amrex::Real g2
Definition: ERF_Interpolation_WENO_Z.H:558
WENO_Z7::g3
static constexpr amrex::Real g3
Definition: ERF_Interpolation_WENO_Z.H:559

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

Here is the caller graph for this function:

◆ InterpolateInX()

AMREX_GPU_DEVICE AMREX_FORCE_INLINE void WENO_Z7::InterpolateInX ( const int &  i,
const int &  j,
const int &  k,
const int &  qty_index,
amrex::Real &  val_lo,
amrex::Real  upw_lo 
) const
inline
402  {
403  // Data to interpolate on
404  amrex::Real sp3 = m_phi(i+3, j , k , qty_index);
405  amrex::Real sp2 = m_phi(i+2, j , k , qty_index);
406  amrex::Real sp1 = m_phi(i+1, j , k , qty_index);
407  amrex::Real s = m_phi(i , j , k , qty_index);
408  amrex::Real sm1 = m_phi(i-1, j , k , qty_index);
409  amrex::Real sm2 = m_phi(i-2, j , k , qty_index);
410  amrex::Real sm3 = m_phi(i-3, j , k , qty_index);
411  amrex::Real sm4 = m_phi(i-4, j , k , qty_index);
412 
413  if (upw_lo > tol) {
414  val_lo = Evaluate(sm4,sm3,sm2,sm1,s ,sp1,sp2);
415  } else if (upw_lo < -tol) {
416  val_lo = Evaluate(sp3,sp2,sp1,s ,sm1,sm2,sm3);
417  } else {
418  val_lo = 0.5 * (s + sm1);
419  }
420  }
WENO_Z7::Evaluate
AMREX_GPU_DEVICE AMREX_FORCE_INLINE amrex::Real Evaluate(const amrex::Real &sm3, const amrex::Real &sm2, const amrex::Real &sm1, const amrex::Real &s, const amrex::Real &sp1, const amrex::Real &sp2, const amrex::Real &sp3) const
Definition: ERF_Interpolation_WENO_Z.H:483
WENO_Z7::tol
const amrex::Real tol
Definition: ERF_Interpolation_WENO_Z.H:556
Here is the call graph for this function:

◆ InterpolateInY()

AMREX_GPU_DEVICE AMREX_FORCE_INLINE void WENO_Z7::InterpolateInY ( const int &  i,
const int &  j,
const int &  k,
const int &  qty_index,
amrex::Real &  val_lo,
amrex::Real  upw_lo 
) const
inline
431  {
432  // Data to interpolate on
433  amrex::Real sp3 = m_phi(i , j+3, k , qty_index);
434  amrex::Real sp2 = m_phi(i , j+2, k , qty_index);
435  amrex::Real sp1 = m_phi(i , j+1, k , qty_index);
436  amrex::Real s = m_phi(i , j , k , qty_index);
437  amrex::Real sm1 = m_phi(i , j-1, k , qty_index);
438  amrex::Real sm2 = m_phi(i , j-2, k , qty_index);
439  amrex::Real sm3 = m_phi(i , j-3, k , qty_index);
440  amrex::Real sm4 = m_phi(i , j-4, k , qty_index);
441 
442  if (upw_lo > tol) {
443  val_lo = Evaluate(sm4,sm3,sm2,sm1,s ,sp1,sp2);
444  } else if (upw_lo < -tol) {
445  val_lo = Evaluate(sp3,sp2,sp1,s ,sm1,sm2,sm3);
446  } else {
447  val_lo = 0.5 * (s + sm1);
448  }
449  }
Here is the call graph for this function:

◆ InterpolateInZ()

AMREX_GPU_DEVICE AMREX_FORCE_INLINE void WENO_Z7::InterpolateInZ ( const int &  i,
const int &  j,
const int &  k,
const int &  qty_index,
amrex::Real &  val_lo,
amrex::Real  upw_lo 
) const
inline
460  {
461  // Data to interpolate on
462  amrex::Real sp3 = m_phi(i , j , k+2, qty_index);
463  amrex::Real sp2 = m_phi(i , j , k+2, qty_index);
464  amrex::Real sp1 = m_phi(i , j , k+1, qty_index);
465  amrex::Real s = m_phi(i , j , k , qty_index);
466  amrex::Real sm1 = m_phi(i , j , k-1, qty_index);
467  amrex::Real sm2 = m_phi(i , j , k-2, qty_index);
468  amrex::Real sm3 = m_phi(i , j , k-3, qty_index);
469  amrex::Real sm4 = m_phi(i , j , k-3, qty_index);
470 
471  if (upw_lo > tol) {
472  val_lo = Evaluate(sm4,sm3,sm2,sm1,s ,sp1,sp2);
473  } else if (upw_lo < -tol) {
474  val_lo = Evaluate(sp3,sp2,sp1,s ,sm1,sm2,sm3);
475  } else {
476  val_lo = 0.5 * (s + sm1);
477  }
478  }
Here is the call graph for this function:

Member Data Documentation

◆ eps

const amrex::Real WENO_Z7::eps =1.0e-40
private

Referenced by Evaluate().

◆ g1

constexpr amrex::Real WENO_Z7::g1 =( 1.0/64.0)
staticconstexprprivate

Referenced by Evaluate().

◆ g2

constexpr amrex::Real WENO_Z7::g2 =(21.0/64.0)
staticconstexprprivate

Referenced by Evaluate().

◆ g3

constexpr amrex::Real WENO_Z7::g3 =(35.0/64.0)
staticconstexprprivate

Referenced by Evaluate().

◆ g4

constexpr amrex::Real WENO_Z7::g4 =( 7.0/64.0)
staticconstexprprivate

Referenced by Evaluate().

◆ m_phi

amrex::Array4<const amrex::Real> WENO_Z7::m_phi
private

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

◆ tol

const amrex::Real WENO_Z7::tol =1.0e-12
private

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

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