CD_EBHelmholtzRobinDomainBCFactory.H

/*
 * SPDX-FileCopyrightText: 2021-2026 SINTEF Energy Research
 *
 * SPDX-License-Identifier: GPL-3.0-or-later
 */
 
/*
  @file   CD_EBHelmholtzRobinDomainBCFactory.H
  @brief  Declaration of a factory class for making Robin domain boundary conditions with EBHelmholtzOp
  @author Robert Marskar
*/
 
#ifndef CD_EBHELMHOLTZROBINDOMAINBCFACTORY_H
#define CD_EBHELMHOLTZROBINDOMAINBCFACTORY_H
 
// Std includes
#include <functional>
 
// Our includes
#include <CD_EBHelmholtzDomainBCFactory.H>
#include <CD_NamespaceHeader.H>
 
class EBHelmholtzRobinDomainBCFactory : public EBHelmholtzDomainBCFactory
{
public:
  EBHelmholtzRobinDomainBCFactory();
 
  EBHelmholtzRobinDomainBCFactory(const Real a_A, const Real a_B, const Real a_C);
 
  EBHelmholtzRobinDomainBCFactory(const std::function<Real(const RealVect& a_pos)>& a_A,
                                  const std::function<Real(const RealVect& a_pos)>& a_B,
                                  const std::function<Real(const RealVect& a_pos)>& a_C);
 
  virtual ~EBHelmholtzRobinDomainBCFactory();
 
  void
  setCoefficients(const Real a_A, const Real a_B, const Real a_C);
 
  void
  setCoefficients(const std::function<Real(const RealVect& a_pos)>& a_A,
                  const std::function<Real(const RealVect& a_pos)>& a_B,
                  const std::function<Real(const RealVect& a_pos)>& a_C);
 
  virtual RefCountedPtr<EBHelmholtzDomainBC>
  create() const override;
 
protected:
  /*
    @brief Use constant for BC
  */
  bool m_useConstant;
 
  bool m_useFunction;
 
  Real m_constantA;
 
  Real m_constantB;
 
  Real m_constantC;
 
  std::function<Real(const RealVect& a_pos)> m_functionA;
 
  std::function<Real(const RealVect& a_pos)> m_functionB;
 
  std::function<Real(const RealVect& a_pos)> m_functionC;
};
 
#include <CD_NamespaceFooter.H>
 
#endif