20 #ifndef OPM_GROUPSTATE_HEADER_INCLUDED
21 #define OPM_GROUPSTATE_HEADER_INCLUDED
26 #include <opm/core/props/BlackoilPhases.hpp>
27 #include <opm/parser/eclipse/EclipseState/Schedule/Group/Group.hpp>
28 #include <opm/parser/eclipse/EclipseState/Runspec.hpp>
29 #include <opm/parser/eclipse/EclipseState/Schedule/Group/GPMaint.hpp>
30 #include <opm/simulators/wells/WellContainer.hpp>
37 bool operator==(
const GroupState& other)
const;
39 bool has_production_rates(
const std::string& gname)
const;
40 void update_production_rates(
const std::string& gname,
const std::vector<double>& rates);
41 const std::vector<double>& production_rates(
const std::string& gname)
const;
43 bool has_production_reduction_rates(
const std::string& gname)
const;
44 void update_production_reduction_rates(
const std::string& gname,
const std::vector<double>& rates);
45 const std::vector<double>& production_reduction_rates(
const std::string& gname)
const;
47 bool has_injection_reduction_rates(
const std::string& gname)
const;
48 void update_injection_reduction_rates(
const std::string& gname,
const std::vector<double>& rates);
49 const std::vector<double>& injection_reduction_rates(
const std::string& gname)
const;
51 bool has_injection_reservoir_rates(
const std::string& gname)
const;
52 void update_injection_reservoir_rates(
const std::string& gname,
const std::vector<double>& rates);
53 const std::vector<double>& injection_reservoir_rates(
const std::string& gname)
const;
55 bool has_injection_surface_rates(
const std::string& gname)
const;
56 void update_injection_surface_rates(
const std::string& gname,
const std::vector<double>& rates);
57 const std::vector<double>& injection_surface_rates(
const std::string& gname)
const;
60 void update_injection_rein_rates(
const std::string& gname,
const std::vector<double>& rates);
61 const std::vector<double>& injection_rein_rates(
const std::string& gname)
const;
63 void update_injection_vrep_rate(
const std::string& gname,
double rate);
64 double injection_vrep_rate(
const std::string& gname)
const;
66 void update_grat_sales_target(
const std::string& gname,
double target);
67 double grat_sales_target(
const std::string& gname)
const;
68 bool has_grat_sales_target(
const std::string& gname)
const;
70 void update_gpmaint_target(
const std::string& gname,
double target);
71 double gpmaint_target(
const std::string& gname)
const;
72 bool has_gpmaint_target(
const std::string& gname)
const;
74 bool has_production_control(
const std::string& gname)
const;
75 void production_control(
const std::string& gname, Group::ProductionCMode cmode);
76 Group::ProductionCMode production_control(
const std::string& gname)
const;
78 bool has_injection_control(
const std::string& gname, Phase phase)
const;
79 void injection_control(
const std::string& gname, Phase phase, Group::InjectionCMode cmode);
80 Group::InjectionCMode injection_control(
const std::string& gname, Phase phase)
const;
82 std::size_t data_size()
const;
83 std::size_t collect(
double * data)
const;
84 std::size_t distribute(
const double * data);
86 GPMaint::State& gpmaint(
const std::string& gname);
90 void communicate_rates(
const Comm& comm)
99 auto iterateContainer = [](
auto& container,
auto& func) {
100 for (
auto& x : container) {
106 auto forAllGroupData = [&](
auto& func) {
107 iterateContainer(m_production_rates, func);
108 iterateContainer(prod_red_rates, func);
109 iterateContainer(inj_red_rates, func);
110 iterateContainer(inj_resv_rates, func);
111 iterateContainer(inj_rein_rates, func);
112 iterateContainer(inj_surface_rates, func);
117 auto computeSize = [&sz](
const auto& v) {
120 forAllGroupData(computeSize);
121 sz += this->inj_vrep_rate.size();
124 std::vector<double> data(sz);
131 auto collect = [&data, &pos](
auto& v) {
137 forAllGroupData(collect);
138 for (
const auto& x : this->inj_vrep_rate) {
139 data[pos++] = x.second;
142 throw std::logic_error(
"Internal size mismatch when collecting groupData");
145 comm.sum(data.data(), data.size());
149 auto distribute = [&data, &pos](
auto& v) {
154 forAllGroupData(distribute);
155 for (
auto& x : this->inj_vrep_rate) {
156 x.second = data[pos++];
159 throw std::logic_error(
"Internal size mismatch when distributing groupData");
162 std::string dump()
const;
166 std::size_t num_phases;
167 std::map<std::string, std::vector<double>> m_production_rates;
168 std::map<std::string, Group::ProductionCMode> production_controls;
169 std::map<std::string, std::vector<double>> prod_red_rates;
170 std::map<std::string, std::vector<double>> inj_red_rates;
171 std::map<std::string, std::vector<double>> inj_surface_rates;
172 std::map<std::string, std::vector<double>> inj_resv_rates;
173 std::map<std::string, std::vector<double>> inj_rein_rates;
174 std::map<std::string, double> inj_vrep_rate;
175 std::map<std::string, double> m_grat_sales_target;
176 std::map<std::string, double> m_gpmaint_target;
179 std::map<std::pair<Phase, std::string>, Group::InjectionCMode> injection_controls;
Definition: GroupState.hpp:34
This file contains a set of helper functions used by VFPProd / VFPInj.
Definition: BlackoilPhases.hpp:26