CD_ScanShop.H

Go to the documentation of this file.

/* chombo-discharge
 * Copyright © 2021 SINTEF Energy Research.
 * Please refer to Copyright.txt and LICENSE in the chombo-discharge root directory.
 */
 
#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_geometry,
           const int           a_verbosity,
           const Real          a_finestDx,
           const RealVect      a_probLo,
           const ProblemDomain a_finestDomain,
           const ProblemDomain a_scanLevel,
           const int           a_ebGhost   = 4,
           const Real          a_thrshdVoF = 1.0e-16);
 
  ~ScanShop();
 
  void
  setProfileFileName(const 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;
 
  virtual 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(const Real          a_dx,
              const RealVect      a_probLo,
              const ProblemDomain a_finestDomain,
              const ProblemDomain a_scanLevel);
 
  void
  buildFinerLevels(const int a_coarserLevel, const int a_maxGridSize);
 
  void
  buildCoarseLevel(const int a_finerLevel, const int a_maxGridSize);
 
  inline bool
  isRegular(const Box a_box, const RealVect a_probLo, const Real a_dx) const;
 
  inline bool
  isCovered(const Box a_box, const RealVect a_probLo, const 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, const int a_level);
};
 
#include <CD_NamespaceFooter.H>
 
#include <CD_ScanShopImplem.H>
 
#endif