CD_CoarseInterpQuadCF.H

Go to the documentation of this file.

/* chombo-discharge
 * Copyright © 2023 SINTEF Energy Research.
 * Please refer to Copyright.txt and LICENSE in the chombo-discharge root directory.
 */
 
#ifndef CD_CoarseInterpQuadCF_H
#define CD_CoarseInterpQuadCF_H
 
// Chombo includes
#include <IntVect.H>
#include <FArrayBox.H>
#include <LoHiSide.H>
#include <DerivStencil.H>
#include <DisjointBoxLayout.H>
 
// Our includes
#include <CD_NamespaceHeader.H>
 
class CoarseInterpQuadCF
{
public:
  CoarseInterpQuadCF() noexcept;
 
  virtual ~CoarseInterpQuadCF() noexcept;
 
  virtual void
  define(const DisjointBoxLayout& a_dblFine,
         const ProblemDomain&     a_domainCoar,
         const DataIndex&         a_dit,
         const Box&               a_fineGhostCells,
         const int                a_refRat,
         const int                a_ignoreDir) noexcept;
 
  virtual Real
  computeFirstDeriv(const FArrayBox& a_coarPhi,
                    const IntVect&   a_ivCoar,
                    const int        a_dir,
                    const int        a_coarVar) const noexcept;
 
  virtual Real
  computeSecondDeriv(const FArrayBox& a_coarPhi,
                     const IntVect&   a_ivCoar,
                     const int        a_dir,
                     const int        a_coarVar) const noexcept;
 
  virtual Real
  computeMixedDeriv(const FArrayBox& a_coarPhi, const IntVect& a_ivCoar, const int a_coarVar) const noexcept;
 
  virtual DerivStencil
  getFirstDerivStencil(const IntVect& a_ivCoar, const int a_dir) const noexcept;
 
  virtual DerivStencil
  getSecondDerivStencil(const IntVect& a_ivCoar, const int a_dir) const noexcept;
 
  virtual DerivStencil
  getMixedDerivStencil(const IntVect& a_ivCoar) const noexcept;
 
protected:
  enum class FirstDerivStencil : unsigned short
  {
    Centered2,
    Backward2,
    Forward2,
    Backward1,
    Forward1,
  };
 
  enum class SecondDerivStencil : unsigned short
  {
    Centered2,
    Backward1,
    Forward1
  };
 
  enum class MixedDerivStencil : unsigned short
  {
    Standard,
    Explicit
  };
 
  bool m_isDefined;
 
  DisjointBoxLayout m_dblFine;
 
  ProblemDomain m_domainCoar;
 
  DataIndex m_dit;
 
  int m_ignoreDir;
 
  int m_tanDir1;
 
  int m_tanDir2;
 
  int m_refRat;
 
  Box m_stencilBox;
 
  BaseFab<FirstDerivStencil> m_firstDerivStencils[SpaceDim];
 
  BaseFab<SecondDerivStencil> m_secondDerivStencils[SpaceDim];
 
  BaseFab<MixedDerivStencil> m_mixedDerivStencils;
 
  BaseFab<DerivStencil> m_explicitMixedDerivStencils;
 
  virtual void
  defineStencils() noexcept;
 
  virtual void
  defineMixedDerivStencils() noexcept;
};
 
#include <CD_NamespaceFooter.H>
 
#endif