CD_CellCentroidInterpolation.H

Go to the documentation of this file.

/* chombo-discharge
 * Copyright © 2025 SINTEF Energy Research.
 * Please refer to Copyright.txt and LICENSE in the chombo-discharge root directory.
 */
 
#ifndef CD_CellCentroidInterpolation_H
#define CD_CellCentroidInterpolation_H
 
// Chombo includes
#include <EBLevelGrid.H>
#include <LevelData.H>
#include <LayoutData.H>
#include <VoFIterator.H>
#include <Stencils.H>
#include <BaseIVFAB.H>
#include <EBCellFAB.H>
 
// Our includes
#include <CD_NamespaceHeader.H>
 
class CellCentroidInterpolation
{
public:
  enum class Type
  {
    Constant,           // Uses cell-centered value
    Linear,             // Uses bi/tri-linear interpolation
    Taylor,             // Uses Taylor-extrapolation
    LeastSquares,       // Uses least-squares reconstruction
    PiecewiseLinear,    // Uses piecewise-linear interpolation
    MinMod,             // Uses slope limiter with minmod slope
    MonotonizedCentral, // Uses slope limiter with van Leer limiter
    Superbee            // Uses slope limiter with superbee limiter
  };
 
  CellCentroidInterpolation() noexcept;
 
  CellCentroidInterpolation(const CellCentroidInterpolation& a_other) = delete;
 
  CellCentroidInterpolation(const EBLevelGrid& a_eblg, const Real& a_dx, const Type& a_interpolationType) noexcept;
 
  virtual ~CellCentroidInterpolation() noexcept;
 
  virtual void
  define(const EBLevelGrid& a_eblg, const Real& a_dx, const Type& a_interpolationType) noexcept;
 
  virtual void
  interpolate(LevelData<BaseIVFAB<Real>>& a_centroidData, const LevelData<EBCellFAB>& a_cellData) const noexcept;
 
  virtual void
  interpolate(LevelData<EBCellFAB>& a_centroidData, const LevelData<EBCellFAB>& a_cellData) const noexcept;
 
  virtual void
  interpolate(LevelData<EBCellFAB>& a_data) const noexcept;
 
  template <typename T>
  void
  interpolate(T& a_centroidData, const EBCellFAB& a_cellData, const DataIndex& a_din) const noexcept;
 
protected:
  bool m_isDefined;
 
  EBLevelGrid m_eblg;
 
  Type m_interpolationType;
 
  Real m_dx;
 
  mutable LayoutData<VoFIterator> m_vofIterator;
 
  LayoutData<BaseIVFAB<VoFStencil>> m_interpStencils;
 
  virtual bool
  getLinearStencil(VoFStencil&          a_stencil,
                   const VolIndex&      a_vof,
                   const EBISBox&       a_ebisBox,
                   const ProblemDomain& a_domain) const noexcept;
 
  virtual bool
  getTaylorExtrapolationStencil(VoFStencil&          a_stencil,
                                const VolIndex&      a_vof,
                                const EBISBox&       a_ebisBox,
                                const ProblemDomain& a_domain) const noexcept;
 
  virtual bool
  getLeastSquaresStencil(VoFStencil&          a_stencil,
                         const VolIndex&      a_vof,
                         const EBISBox&       a_ebisBox,
                         const ProblemDomain& a_domain) const noexcept;
 
  virtual bool
  getPiecewiseLinearStencil(VoFStencil&          a_stencil,
                            const VolIndex&      a_vof,
                            const EBISBox&       a_ebisBox,
                            const ProblemDomain& a_domain) const noexcept;
 
  inline Real
  MinMod(const Real& a_dwl, const Real& a_dwr) const noexcept;
 
  inline Real
  MonotonizedCentral(const Real& a_dwl, const Real& a_dwr) const noexcept;
 
  inline Real
  Superbee(const Real& a_dwl, const Real& a_dwr) const noexcept;
};
 
#include <CD_NamespaceFooter.H>
 
#include <CD_CellCentroidInterpolationImplem.H>
 
#endif