23#ifndef OPM_REGIONATTRIBUTEHELPERS_HPP_HEADER_INCLUDED
24#define OPM_REGIONATTRIBUTEHELPERS_HPP_HEADER_INCLUDED
26#include <opm/core/props/BlackoilPhases.hpp>
27#include <opm/grid/utility/RegionMapping.hpp>
29#include <dune/grid/common/gridenums.hh>
35#include <unordered_map>
40 namespace RegionAttributeHelpers {
46 template <
class RegionID,
bool>
50 typename std::remove_reference<RegionID>::type &;
53 template <
class RegionID>
56 using type = RegionID;
66 template<
bool is_parallel>
80 std::tuple<double, double, double, double, int>
82 const std::vector<double>& temperature,
83 const std::vector<double>& rs,
84 const std::vector<double>& rv,
85 const std::vector<double>& ownership,
87 if ( ownership[cell] )
89 return std::make_tuple(pressure[cell],
97 return std::make_tuple(0, 0, 0, 0, 0);
104 std::tuple<double, double, double, double, int>
105 operator()(
const std::vector<double>& pressure,
106 const std::vector<double>& temperature,
107 const std::vector<double>& rs,
108 const std::vector<double>& rv,
109 const std::vector<double>&,
111 return std::make_tuple(pressure[cell],
130 template <
typename RegionId,
class Attributes>
140 <RegionId, std::is_integral<RegionId>::value>::type;
143 typename std::remove_reference<RegionId>::type;
150 Value(
const Attributes& attr)
160 std::unordered_map<ID, std::unique_ptr<Value>>;
176 template <
class RMap>
178 const Attributes& attr)
180 using VT =
typename AttributeMap::value_type;
182 for (
const auto& r : rmap.activeRegions()) {
183 auto v = std::make_unique<Value>(attr);
185 const auto stat = attr_.insert(VT(r, std::move(v)));
189 const auto& cells = rmap.cells(r);
191 assert (! cells.empty());
194 stat.first->second->cell_ = cells[0];
208 return this->find(reg).cell_;
213 return this->attr_.find(reg) != this->attr_.end();
216 void insert(
const RegionID r,
const Attributes& attr)
218 auto [pos, inserted] = this->attr_.try_emplace(r, std::make_unique<Value>(attr));
220 pos->second->cell_ = -1;
246 return this->find(reg).attr_;
259 return this->find(reg).attr_;
269 const Value& find(
const RegionID reg)
const
271 const auto& i = attr_.find(reg);
273 if (i == attr_.end()) {
274 throw std::invalid_argument(
"Unknown region ID");
285 const auto& i = attr_.find(reg);
287 if (i == attr_.end()) {
288 throw std::invalid_argument(
"Unknown region ID");
299 namespace PhaseUsed {
310 return pu.phase_used[ BlackoilPhases::Aqua ] != 0;
323 return pu.phase_used[ BlackoilPhases::Liquid ] != 0;
336 return pu.phase_used[ BlackoilPhases::Vapour ] != 0;
359 p = pu.phase_pos[ BlackoilPhases::Aqua ];
379 p = pu.phase_pos[ BlackoilPhases::Liquid ];
399 p = pu.phase_pos[ BlackoilPhases::Vapour ];
Provide mapping from Region IDs to user-specified collection of per-region attributes.
Definition: RegionAttributeHelpers.hpp:132
int cell(const RegionID reg) const
Retrieve representative cell in region.
Definition: RegionAttributeHelpers.hpp:206
const AttributeMap & attributes() const
Request read-only access to region's attributes.
Definition: RegionAttributeHelpers.hpp:230
typename Select::RegionIDParameter< RegionId, std::is_integral< RegionId >::value >::type RegionID
Expose RegionId as a vocabulary type for use in query methods.
Definition: RegionAttributeHelpers.hpp:140
const Attributes & attributes(const RegionID reg) const
Request read-only access to region's attributes.
Definition: RegionAttributeHelpers.hpp:244
RegionAttributes(const RMap &rmap, const Attributes &attr)
Constructor.
Definition: RegionAttributeHelpers.hpp:177
Attributes & attributes(const RegionID reg)
Request modifiable access to region's attributes.
Definition: RegionAttributeHelpers.hpp:257
int gas(const PhaseUsage &pu)
Numerical ID of active gas phase.
Definition: RegionAttributeHelpers.hpp:394
int oil(const PhaseUsage &pu)
Numerical ID of active oil phase.
Definition: RegionAttributeHelpers.hpp:374
int water(const PhaseUsage &pu)
Numerical ID of active water phase.
Definition: RegionAttributeHelpers.hpp:354
bool water(const PhaseUsage &pu)
Active water predicate.
Definition: RegionAttributeHelpers.hpp:308
bool oil(const PhaseUsage &pu)
Active oil predicate.
Definition: RegionAttributeHelpers.hpp:321
bool gas(const PhaseUsage &pu)
Active gas predicate.
Definition: RegionAttributeHelpers.hpp:334
This file contains a set of helper functions used by VFPProd / VFPInj.
Definition: BlackoilPhases.hpp:27
Definition: BlackoilPhases.hpp:46
Computes the temperature, pressure, and counter increment.
Definition: RegionAttributeHelpers.hpp:68
std::tuple< double, double, double, double, int > operator()(const std::vector< double > &pressure, const std::vector< double > &temperature, const std::vector< double > &rs, const std::vector< double > &rv, const std::vector< double > &ownership, std::size_t cell)
Computes the temperature, pressure, and counter increment.
Definition: RegionAttributeHelpers.hpp:81
Aggregate per-region attributes along with region's representative cell.
Definition: RegionAttributeHelpers.hpp:149
Definition: RegionAttributeHelpers.hpp:48