OilfieldToolsNet/opm-simulators
4
0
Fork 0
mirror of https://github.com/OPM/opm-simulators.git synced 2025-02-25 18:55:30 -06:00
Code Issues Packages Projects Releases Wiki Activity

Merge pull request #5715 from aritorto/missedCommentPR5712

Browse Source 
LevelCartesianIndexMapper constructor for AluGrid taking only a grid as an argument
This commit is contained in:
Lisa Julia Nebel
2024-12-20 07:32:02 +01:00
committed by GitHub
parent 4f5c86d72c 423b6596fb
commit 78708e1213
3 changed files with 16 additions and 8 deletions
Download Patch File Download Diff File Expand all files Collapse all files

18
opm/simulators/flow/AluGridLevelCartesianIndexMapper.hpp
View File

@@ -43,6 +43,15 @@
namespace Opm {
// Interface class to access the local Cartesian grid of each level grid (when refinement).
// Further documentation in opm/grid/common/LevelCartesianIndexMapper.hpp
//
// Adapter Design Pattern: In this case, LevelCartesianIndexMapper uses the Object Adapter variant, where it holds an instance
// (here, a std::unique_ptr) of CartesianIndexMapper, the wrapped type. The goal is to provide a standardized interface, allowing
// incompatible functionality (such as Cartesian indexing in the context of refinement that may not be supported - yet -for all
// grid types, like CpGrid) to integrate smoothly within the existing conventions.
//
// Specialization for AluGrid
template<>
class LevelCartesianIndexMapper<Dune::ALUGrid<3, 3, Dune::cube, Dune::nonconforming>>
{
@@ -56,12 +65,10 @@ class LevelCartesianIndexMapper<Dune::ALUGrid<3, 3, Dune::cube, Dune::nonconform
public:
static constexpr int dimension = 3 ;
LevelCartesianIndexMapper(const Grid& grid,
const Dune::CartesianIndexMapper<Grid>& cartesianIndexMapper)
: grid_{&grid}
, cartesianIndexMapper_{std::make_unique<Dune::CartesianIndexMapper<Grid>>(cartesianIndexMapper)}
explicit LevelCartesianIndexMapper(const Dune::CartesianIndexMapper<Grid>& cartesianIndexMapper)
: cartesianIndexMapper_{std::make_unique<Dune::CartesianIndexMapper<Grid>>(cartesianIndexMapper)}
{}
const std::array<int,3>& cartesianDimensions(int level) const
{
throwIfLevelPositive(level);
@@ -93,7 +100,6 @@ class LevelCartesianIndexMapper<Dune::ALUGrid<3, 3, Dune::cube, Dune::nonconform
}
private:
[[maybe_unused]] const Grid* grid_;
std::unique_ptr<Dune::CartesianIndexMapper<Grid>> cartesianIndexMapper_;
void throwIfLevelPositive(int level) const

3
opm/simulators/flow/AluGridVanguard.hpp
View File

@@ -238,9 +238,10 @@ public:
/*!
* \brief Returns the object which maps a global element index of the simulation grid
* to the corresponding element index of the level logically Cartesian index.
* No refinement is supported for AluGrid so it coincides with CartesianIndexMapper.
*/
const LevelCartesianIndexMapper levelCartesianIndexMapper() const
{ return LevelCartesianIndexMapper(*grid_, *cartesianIndexMapper_); }
{ return LevelCartesianIndexMapper(*cartesianIndexMapper_); }
/*!
* \brief Returns mapper from compressed to cartesian indices for the EQUIL grid

3
opm/simulators/flow/PolyhedralGridVanguard.hpp
View File

@@ -174,9 +174,10 @@ public:
/*!
* \brief Returns the object which maps a global element index of the simulation grid
* to the corresponding element index of the level logically Cartesian index.
* No refinement is supported for AluGrid so it coincides with CartesianIndexMapper.
*/
const LevelCartesianIndexMapper levelCartesianIndexMapper() const
{ return LevelCartesianIndexMapper(*grid_); }
{ return LevelCartesianIndexMapper(*cartesianIndexMapper_); }
/*!
* \brief Returns mapper from compressed to cartesian indices for the EQUIL grid