CD_TiledMeshRefine.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_TILEDMESHREFINE_H
#define CD_TILEDMESHREFINE_H
 
// Std includes
#include <cstdint>
#include <unordered_map>
#include <unordered_set>
#include <vector>
 
// Chombo includes
#include <IntVectSet.H>
#include <ProblemDomain.H>
#include <Box.H>
#include <BoxIterator.H>
 
// Our includes
#include <CD_NamespaceHeader.H>
 
class TiledMeshRefine
{
public:
  TiledMeshRefine() = delete;
 
  TiledMeshRefine(const ProblemDomain& a_coarsestDomain,
                  const Vector<int>&   a_refRatios,
                  const IntVect&       a_tileSize,
                  const IntVect&       a_maxBlockSize) noexcept;
 
  virtual ~TiledMeshRefine() noexcept;
 
  virtual int
  regrid(Vector<Vector<Box>>& a_newGrids, const Vector<IntVectSet>& a_tagsLevel) const noexcept;
 
protected:
  struct SuperTiles
  {
    std::unordered_set<std::uint64_t> m_full;
 
    std::unordered_map<std::uint64_t, std::vector<std::uint64_t>> m_partial;
  };
 
  Vector<ProblemDomain> m_amrDomains;
 
  Vector<int> m_refRatios;
 
  IntVect m_tileSize;
 
  IntVect m_maxBlockSize;
 
  IntVect m_superFactor;
 
  long long m_superVol;
 
  int m_bitmaskWords;
 
  std::uint64_t
  encodeSuper(const IntVect& a_super) const noexcept;
 
  IntVect
  decodeSuper(const std::uint64_t a_key) const noexcept;
 
  int
  subIndex(const IntVect& a_sub) const noexcept;
 
  void
  addFineTile(SuperTiles& a_tiles, const IntVect& a_fineTile) const noexcept;
 
  void
  classify(SuperTiles& a_tiles) const noexcept;
 
  void
  gatherSuperTiles(SuperTiles& a_tiles) const noexcept;
 
  void
  nestFrom(SuperTiles&       a_tiles,
           const SuperTiles& a_finer,
           const int         a_refToFine,
           const Box&        a_thisTileBox) const noexcept;
 
  virtual void
  makeLevelTiles(SuperTiles&          a_tiles,
                 const SuperTiles&    a_fineTiles,
                 const IntVectSet&    a_coarTags,
                 const ProblemDomain& a_domain,
                 const int            a_refToFine,
                 const int            a_refToCoar) const noexcept;
 
  virtual void
  makeBoxesFromTiles(Vector<Box>& a_boxes, const SuperTiles& a_tiles, const ProblemDomain& a_domain) const noexcept;
 
  void
  packTiles(std::vector<IntVect>& a_tiles,
            const IntVect&        a_maxTile,
            const IntVect&        a_probLo,
            Vector<Box>&          a_boxes) const noexcept;
};
 
#include <CD_NamespaceFooter.H>
 
#endif