CD_MFHelmholtzRobinDomainBCFactory.H

/*
 * SPDX-FileCopyrightText: 2021-2026 SINTEF Energy Research
 *
 * SPDX-License-Identifier: GPL-3.0-or-later
 */
 
/*
  @file   CD_MFHelmholtzRobinDomainBCFactory.H
  @brief  Declaration of a factory class for making Robin domain boundary conditions with MFHelmholtzOp
  @author Robert Marskar
*/
 
#ifndef CD_MFHELMHOLTZROBINDOMAINBCFACTORY_H
#define CD_MFHELMHOLTZROBINDOMAINBCFACTORY_H
 
// Std includes
#include <functional>
 
// Our includes
#include <CD_MFHelmholtzDomainBCFactory.H>
#include <CD_NamespaceHeader.H>
 
class MFHelmholtzRobinDomainBCFactory : public MFHelmholtzDomainBCFactory
{
public:
  MFHelmholtzRobinDomainBCFactory();
 
  MFHelmholtzRobinDomainBCFactory(const Real a_A, const Real a_B, const Real a_C);
 
  MFHelmholtzRobinDomainBCFactory(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 ~MFHelmholtzRobinDomainBCFactory();
 
  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);
 
  RefCountedPtr<EBHelmholtzDomainBC>
  create(const int a_iphase) 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