CD_MFHelmholtzJumpBCImplem.H

/*
 * SPDX-FileCopyrightText: 2021-2026 SINTEF Energy Research
 *
 * SPDX-License-Identifier: GPL-3.0-or-later
 */
 
/*
  @file   CD_MFHelmholtzJumpBCImplem.H
  @brief  Implementation of CD_MFHelmholtzJumpBC.H
  @author Robert Marskar
*/
 
#ifndef CD_MFHELMHOLTZJUMPBCIMPLEM_H
#define CD_MFHELMHOLTZJUMPBCIMPLEM_H
 
// Chombo includes
#include <NeighborIterator.H>
 
// Our includes
#include <CD_MFHelmholtzJumpBC.H>
#include <CD_NamespaceHeader.H>
 
inline bool
MFHelmholtzJumpBC::isStencilValidCF(const VoFStencil& a_stencil, const DataIndex& a_dit) const
{
  const DisjointBoxLayout& dbl    = m_mflg.getGrids();
  const ProblemDomain&     domain = m_mflg.getDomain();
 
  // Construct boxes that contain all the valid cells for this stencil.
  std::vector<Box> validBoxes;
 
  Box curBox = dbl[a_dit];
  curBox.grow(m_ghostCF);
  curBox &= domain;
  validBoxes.emplace_back(curBox);
 
  NeighborIterator nit(dbl);
  for (nit.begin(a_dit); nit.ok(); ++nit) {
    Box neighBox = dbl[nit()];
    neighBox.grow(m_ghostCF);
    neighBox &= domain;
 
    validBoxes.emplace_back(neighBox);
  }
 
  // Now check that the stencil if the stencil. We set valid = false
  // if any of the stencil points reaches out of the ghosted boxes.
  bool valid = true;
 
  for (int i = 0; i < a_stencil.size(); i++) {
    const VolIndex& vof = a_stencil.vof(i);
    const IntVect&  iv  = vof.gridIndex();
 
    bool insideOneBox = false;
    for (const auto& b : validBoxes) {
      if (b.contains(iv))
        insideOneBox = true;
    }
 
    if (!insideOneBox) {
      valid = false;
      break;
    }
  }
 
  return valid;
}
 
inline Real
MFHelmholtzJumpBC::applyStencil(const VoFStencil& a_stencil, const EBCellFAB& a_phi) const
{
  Real contrib = 0.0;
  for (int i = 0; i < a_stencil.size(); i++) {
    contrib += a_stencil.weight(i) * a_phi(a_stencil.vof(i), m_comp);
  }
 
  return contrib;
}
 
#include <CD_NamespaceFooter.H>
 
#endif