RefinementRequest.hpp

Go to the documentation of this file.

/*
  Copyright 2026 SINTEF Digital.
 
  This file is part of the Open Porous Media project (OPM).
 
  OPM is free software: you can redistribute it and/or modify
  it under the terms of the GNU General Public License as published by
  the Free Software Foundation, either version 3 of the License, or
  (at your option) any later version.
 
  OPM is distributed in the hope that it will be useful,
  but WITHOUT ANY WARRANTY; without even the implied warranty of
  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
  GNU General Public License for more details.
 
  You should have received a copy of the GNU General Public License
  along with OPM.  If not, see <http://www.gnu.org/licenses/>.
*/
#ifndef OPM_GRID_REFINEMENT_REQUEST_HEADER_INCLUDED
#define OPM_GRID_REFINEMENT_REQUEST_HEADER_INCLUDED
 
#include <array>
#include <cstddef>
#include <string>
#include <vector>
 
namespace Opm
{
namespace Refinement
{
 
struct AxisSubdivision
{
    std::vector<int> parentOffset{};
    std::vector<double> fracLo{};
    std::vector<double> fracHi{};
 
    bool empty() const { return parentOffset.empty(); }
    std::size_t size() const { return parentOffset.size(); }
};
 
struct BlockRefinement
{
    std::string name;
    std::string parentGridName{"GLOBAL"};
    std::array<int,3> cellsPerDim{};
    std::array<int,3> startIJK{};
    std::array<int,3> endIJK{};
    std::array<AxisSubdivision,3> subdivision{};
    std::vector<int> minpvRemoved{};
    bool pillarsFromBoxLayer{false};
};
 
AxisSubdivision axisSubdivision(const BlockRefinement& request, int dim);
 
std::array<int,3> refinedDims(const BlockRefinement& request);
 
struct AxisPositions
{
    std::vector<int> subIndex{};
    std::vector<int> parentCount{};
    std::vector<int> firstColumn{};
};
 
AxisPositions axisPositions(const AxisSubdivision& sub);
 
std::vector<BlockRefinement>
blockRefinements(const std::vector<std::array<int,3>>& cellsPerDim,
                 const std::vector<std::array<int,3>>& startIJK,
                 const std::vector<std::array<int,3>>& endIJK,
                 const std::vector<std::string>& names,
                 const std::vector<std::string>& parentGridNames);
 
void validateBlockRefinements(const std::vector<BlockRefinement>& requests);
 
} // namespace Refinement
} // namespace Opm
 
#endif // OPM_GRID_REFINEMENT_REQUEST_HEADER_INCLUDED