CD_IrregStencil.H

Go to the documentation of this file.

/* chombo-discharge
 * Copyright © 2021 SINTEF Energy Research.
 * Please refer to Copyright.txt and LICENSE in the chombo-discharge root directory.
 */
 
#ifndef CD_IrregStencil_H
#define CD_IrregStencil_H
 
// Chombo includes
#include <DisjointBoxLayout.H>
#include <EBISLayout.H>
#include <BaseIVFAB.H>
#include <EBCellFAB.H>
#include <RefCountedPtr.H>
#include <Stencils.H>
#include <LayoutData.H>
#include <VoFIterator.H>
 
// Our includes
#include <CD_NamespaceHeader.H>
 
class IrregStencil
{
public:
  enum class StencilType
  {
    NoType,
    Linear,
    TaylorExtrapolation,
    LeastSquares,
    PiecewiseLinear
  };
 
  IrregStencil();
 
  IrregStencil(const DisjointBoxLayout&        a_dbl,
               const EBISLayout&               a_ebisl,
               const ProblemDomain&            a_domain,
               const Real&                     a_dx,
               const int                       a_order,
               const int                       a_radius,
               const IrregStencil::StencilType a_type);
 
  virtual ~IrregStencil();
 
  const BaseIVFAB<VoFStencil>&
  operator[](const DataIndex& a_dit) const;
 
  BaseIVFAB<VoFStencil>&
  operator[](const DataIndex& a_dit);
 
  virtual void
  apply(EBCellFAB& a_dst, const EBCellFAB& a_src, const DataIndex& a_dit) const;
 
  virtual void
  apply(BaseIVFAB<Real>& a_dst, const EBCellFAB& a_src, const DataIndex& a_dit) const;
 
protected:
  static constexpr int m_defaultStenComp = 0;
 
  static constexpr int m_defaultNumSten = 1;
 
  LayoutData<RefCountedPtr<BaseIVFAB<VoFStencil>>> m_stencils;
 
  mutable LayoutData<VoFIterator> m_vofIter;
 
  DisjointBoxLayout m_dbl;
 
  EBISLayout m_ebisl;
 
  Real m_dx;
 
  int m_radius;
 
  int m_order;
 
  IrregStencil::StencilType m_stencilType;
 
  ProblemDomain m_domain;
 
  virtual void
  define(const DisjointBoxLayout&        a_dbl,
         const EBISLayout&               a_ebisl,
         const ProblemDomain&            a_domain,
         const Real&                     a_dx,
         const int                       a_order,
         const int                       a_radius,
         const IrregStencil::StencilType a_type);
 
  virtual void
  buildStencil(VoFStencil&              a_sten,
               const VolIndex&          a_vof,
               const DisjointBoxLayout& a_dbl,
               const ProblemDomain&     a_domain,
               const EBISBox&           a_ebisbox,
               const Box&               a_box,
               const Real&              a_dx,
               const IntVectSet&        a_cfivs) = 0;
};
 
#include <CD_NamespaceFooter.H>
 
#endif