CD_ScanShop.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_SCANSHOP_H
#define CD_SCANSHOP_H
 
// Chombo includes
#include <GeometryShop.H>
#include <GeometryService.H>
 
// Our includes
#include <CD_Timer.H>
#include <CD_BoxSorting.H>
#include <CD_NamespaceHeader.H>
 
class ScanShop : public GeometryShop
{
public:
  ScanShop() = delete;
 
  ScanShop(const BaseIF&        a_localGeom,
           int                  a_verbosity,
           Real                 a_dx,
           const RealVect&      a_probLo,
           const ProblemDomain& a_finestDomain,
           const ProblemDomain& a_scanLevel,
           int                  a_ebGhost   = 4,
           Real                 a_thrshdVoF = 1.0e-16);
 
  ~ScanShop() override;
 
  ScanShop(const ScanShop&) = delete;
 
  ScanShop&
  operator=(const ScanShop&) = delete;
 
  ScanShop(ScanShop&&) = delete;
 
  ScanShop&
  operator=(ScanShop&&) = delete;
 
  using GeometryShop::InsideOutside;
  using GeometryShop::isCovered;
  using GeometryShop::isRegular;
 
  void
  setProfileFileName(std::string a_fileName);
 
  void
  makeGrids(const ProblemDomain& a_domain,
            DisjointBoxLayout&   a_grids,
            const int&           a_maxGridSize,
            const int&           a_maxIrregGridSize) override;
 
  GeometryService::InOut
  InsideOutside(const Box&           a_region,
                const ProblemDomain& a_domain,
                const RealVect&      a_probLo,
                const Real&          a_dx,
                const DataIndex&     a_dit) const override;
 
  void
  fillGraph(BaseFab<int>&        a_regIrregCovered,
            Vector<IrregNode>&   a_nodes,
            const Box&           a_validRegion,
            const Box&           a_ghostRegion,
            const ProblemDomain& a_domain,
            const RealVect&      a_probLo,
            const Real&          a_dx,
            const DataIndex&     a_di) const override;
 
protected:
  std::string m_fileName;
 
  BoxSorting m_boxSorting;
 
  mutable Timer m_timer;
 
  bool m_profile;
 
  int m_scanLevel;
 
  int m_ebGhost;
 
  RealVect m_probLo;
 
  const BaseIF* m_baseIF;
 
  bool m_hasScanLevel;
 
  Vector<Real> m_dx;
 
  Vector<ProblemDomain> m_domains;
 
  Vector<DisjointBoxLayout> m_grids;
 
  Vector<RefCountedPtr<LayoutData<GeometryService::InOut>>> m_boxMap;
 
  std::vector<bool> m_hasThisLevel;
 
  void
  makeDomains(Real                 a_dx,
              const RealVect&      a_probLo,
              const ProblemDomain& a_finestDomain,
              const ProblemDomain& a_scanLevel);
 
  void
  buildFinerLevels(int a_coarserLevel, int a_maxGridSize);
 
  void
  buildCoarseLevel(int a_level, int a_maxGridSize);
 
  inline bool
  isRegular(Box a_box, const RealVect& a_probLo, Real a_dx) const;
 
  inline bool
  isCovered(Box a_box, const RealVect& a_probLo, Real a_dx) const;
 
  inline std::vector<std::pair<Box, int>>
  getSortedBoxesAndTypes(const Vector<Box>& a_boxes, const Vector<int>& a_types) const;
 
  void
  defineLevel(Vector<Box>& coveredBoxes, Vector<Box>& regularBoxes, Vector<Box>& cutCellBoxes, int a_level);
};
 
#include <CD_NamespaceFooter.H>
 
#include <CD_ScanShopImplem.H>
 
#endif