CD_EBHelmholtzDomainBC.H

/* chombo-discharge
 * Copyright © 2021 SINTEF Energy Research.
 * Please refer to Copyright.txt and LICENSE in the chombo-discharge root directory.
 */
 
/*
  @file   CD_EBHelmholtzDomainBC.H
  @brief  Declaration of a domain boundary condition class for EBHelmholtzOp.
  @author Robert Marskar
*/
 
#ifndef CD_EBHelmholtzDomainBC_H
#define CD_EBHelmholtzDomainBC_H
 
// Chombo includes
#include <EBLevelGrid.H>
#include <RefCountedPtr.H>
#include <LevelData.H>
#include <EBFluxFAB.H>
 
// Our includes
#include <CD_Location.H>
#include <CD_NamespaceHeader.H>
 
class EBHelmholtzDomainBC
{
public:
  EBHelmholtzDomainBC();
 
  virtual ~EBHelmholtzDomainBC();
 
  EBHelmholtzDomainBC(const EBHelmholtzDomainBC& a_other) = delete;
 
  EBHelmholtzDomainBC(const EBHelmholtzDomainBC&& a_other) = delete;
 
  EBHelmholtzDomainBC&
  operator=(const EBHelmholtzDomainBC& a_other) = delete;
 
  EBHelmholtzDomainBC&
  operator=(const EBHelmholtzDomainBC&& a_other) = delete;
 
  virtual void
  define(const Location::Cell a_dataLocation, const EBLevelGrid& a_eblg, const RealVect& a_probLo, const Real a_dx);
 
  virtual void
  getFaceFlux(BaseFab<Real>&        a_faceFlux,
              const BaseFab<Real>&  a_phi,
              const BaseFab<Real>&  a_Bcoef,
              const int&            a_dir,
              const Side::LoHiSide& a_side,
              const DataIndex&      a_dit,
              const bool            a_useHomogeneous) const = 0;
 
  virtual Real
  getFaceFlux(const VolIndex&       a_vof,
              const EBCellFAB&      a_phi,
              const EBFaceFAB&      a_Bcoef,
              const int&            a_dir,
              const Side::LoHiSide& a_side,
              const DataIndex&      a_dit,
              const bool            a_useHomogeneous) const = 0;
 
protected:
  constexpr static int m_comp = 0;
 
  constexpr static int m_nComp = 1;
 
  Location::Cell m_dataLocation;
 
  Real m_dx;
 
  RealVect m_probLo;
 
  EBLevelGrid m_eblg;
 
  inline RealVect
  getBoundaryPosition(const IntVect& a_iv, const int& a_dir, const Side::LoHiSide& a_side) const;
 
  void
  multiplyByBcoef(BaseFab<Real>&       a_flux,
                  const BaseFab<Real>& a_bco,
                  const int            a_dir,
                  const Side::LoHiSide a_side) const;
};
 
#include <CD_NamespaceFooter.H>
 
#include <CD_EBHelmholtzDomainBCImplem.H>
 
#endif