CD_MFInterfaceFABImplem.H

Go to the documentation of this file.

/*
 * SPDX-FileCopyrightText: 2021-2026 SINTEF Energy Research
 *
 * SPDX-License-Identifier: GPL-3.0-or-later
 */
 
#ifndef CD_MFINTERFACEFABIMPLEM_H
#define CD_MFINTERFACEFABIMPLEM_H
 
// Our includes
#include <CD_MFInterfaceFAB.H>
#include <CD_NamespaceHeader.H>
 
template <class T>
MFInterfaceFAB<T>::MFInterfaceFAB()
{}
 
template <class T>
MFInterfaceFAB<T>::MFInterfaceFAB(const MFLevelGrid& a_mflg, const DataIndex& a_dit)
{
  this->define(a_mflg, a_dit);
}
 
template <class T>
MFInterfaceFAB<T>::~MFInterfaceFAB()
{
  for (int i = 0; i < m_phase.size(); i++) {
    delete m_phase[i];
  }
}
 
template <class T>
void
MFInterfaceFAB<T>::define(const MFLevelGrid& a_mflg, const DataIndex& a_dit)
{
 
  const DisjointBoxLayout&                   dbl                  = a_mflg.getGrids();
  const ProblemDomain&                       domain               = a_mflg.getDomain();
  const RefCountedPtr<MultiFluidIndexSpace>& MultiFluidIndexSpace = a_mflg.getMfIndexSpace();
  const int                                  numPhases            = MultiFluidIndexSpace->numPhases();
 
  m_ivs = a_mflg.interfaceRegion(dbl.get(a_dit), a_dit);
  m_phase.resize(numPhases, NULL);
 
  for (int i = 0; i < numPhases; i++) {
    const EBLevelGrid& eblg  = a_mflg.getEBLevelGrid(i);
    const EBISLayout&  ebisl = eblg.getEBISL();
    m_phase[i]               = new BaseIVFAB<T>(m_ivs, ebisl[a_dit].getEBGraph(), 1);
  }
}
 
template <class T>
const IntVectSet&
MFInterfaceFAB<T>::getIVS() const
{
  return m_ivs;
}
 
template <class T>
BaseIVFAB<T>&
MFInterfaceFAB<T>::getIVFAB(const int a_phase)
{
  CH_assert(a_phase >= 0 && a_phase < m_phase.size());
  CH_assert(m_phase[a_phase] != nullptr);
  return *m_phase[a_phase];
}
 
template <class T>
const BaseIVFAB<T>&
MFInterfaceFAB<T>::getIVFAB(const int a_phase) const
{
  CH_assert(a_phase >= 0 && a_phase < m_phase.size());
  CH_assert(m_phase[a_phase] != nullptr);
  return *m_phase[a_phase];
}
 
template <class T>
BaseIVFAB<T>*
MFInterfaceFAB<T>::getPhasePtr(int a_phase)
{
  CH_assert(a_phase >= 0 && a_phase < m_phase.size());
  return m_phase[a_phase];
}
 
template <class T>
int
MFInterfaceFAB<T>::numPhases() const
{
  return m_phase.size();
}
 
template <class T>
void
MFInterfaceFAB<T>::setVal(const T& a_value)
{
  for (int i = 0; i < m_phase.size(); i++) {
    m_phase[i]->setVal(a_value);
  }
}
 
#include <CD_NamespaceFooter.H>
 
#endif