CD_EBHelmholtzDomainBC.H

/*
 * SPDX-FileCopyrightText: 2021-2026 SINTEF Energy Research
 *
 * SPDX-License-Identifier: GPL-3.0-or-later
 */
 
/*
  @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;
 
  virtual Real
  getDiagWeight(const int a_dir, const Side::LoHiSide a_side) 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;
 
  static void
  multiplyByBcoef(BaseFab<Real>& a_flux, const BaseFab<Real>& a_bco, const int a_dir, const Side::LoHiSide a_side);
};
 
#include <CD_NamespaceFooter.H>
 
#include <CD_EBHelmholtzDomainBCImplem.H>
 
#endif