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

drop using statement for ipr_b_

Browse Source 
rather qualify member variable with this->
This commit is contained in:
Arne Morten Kvarving 2021-09-06 12:58:16 +02:00
parent 870cb04e68
commit ae212ff883
2 changed files with 5 additions and 9 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

4
opm/simulators/wells/StandardWell.hpp
View File

@ -261,10 +261,6 @@ namespace Opm
protected: protected:
// protected member variables from the Base class
using Base::ipr_b_;
Eval getPerfCellPressure(const FluidState& fs) const; Eval getPerfCellPressure(const FluidState& fs) const;
// xw = inv(D)*(rw - C*x) // xw = inv(D)*(rw - C*x)

10
opm/simulators/wells/StandardWell_impl.hpp
View File

@ -935,7 +935,7 @@ namespace Opm
// initialize all the values to be zero to begin with // initialize all the values to be zero to begin with
std::fill(this->ipr_a_.begin(), this->ipr_a_.end(), 0.); std::fill(this->ipr_a_.begin(), this->ipr_a_.end(), 0.);
std::fill(ipr_b_.begin(), ipr_b_.end(), 0.); std::fill(this->ipr_b_.begin(), this->ipr_b_.end(), 0.);
for (int perf = 0; perf < this->number_of_perforations_; ++perf) { for (int perf = 0; perf < this->number_of_perforations_; ++perf) {
std::vector<EvalWell> mob(this->num_components_, {this->numWellEq_ + Indices::numEq, 0.0}); std::vector<EvalWell> mob(this->num_components_, {this->numWellEq_ + Indices::numEq, 0.0});
@ -978,7 +978,7 @@ namespace Opm
// phases vs components // phases vs components
// ipr values for the perforation // ipr values for the perforation
std::vector<double> ipr_a_perf(this->ipr_a_.size()); std::vector<double> ipr_a_perf(this->ipr_a_.size());
std::vector<double> ipr_b_perf(ipr_b_.size()); std::vector<double> ipr_b_perf(this->ipr_b_.size());
for (int p = 0; p < this->number_of_phases_; ++p) { for (int p = 0; p < this->number_of_phases_; ++p) {
const double tw_mob = tw_perf * mob[p].value() * b_perf[p]; const double tw_mob = tw_perf * mob[p].value() * b_perf[p];
ipr_a_perf[p] += tw_mob * pressure_diff; ipr_a_perf[p] += tw_mob * pressure_diff;
@ -1008,11 +1008,11 @@ namespace Opm
for (int p = 0; p < this->number_of_phases_; ++p) { for (int p = 0; p < this->number_of_phases_; ++p) {
// TODO: double check the indices here // TODO: double check the indices here
this->ipr_a_[this->ebosCompIdxToFlowCompIdx(p)] += ipr_a_perf[p]; this->ipr_a_[this->ebosCompIdxToFlowCompIdx(p)] += ipr_a_perf[p];
ipr_b_[this->ebosCompIdxToFlowCompIdx(p)] += ipr_b_perf[p]; this->ipr_b_[this->ebosCompIdxToFlowCompIdx(p)] += ipr_b_perf[p];
} }
} }
this->parallel_well_info_.communication().sum(this->ipr_a_.data(), this->ipr_a_.size()); this->parallel_well_info_.communication().sum(this->ipr_a_.data(), this->ipr_a_.size());
this->parallel_well_info_.communication().sum(ipr_b_.data(), ipr_b_.size()); this->parallel_well_info_.communication().sum(this->ipr_b_.data(), this->ipr_b_.size());
} }
@ -1031,7 +1031,7 @@ namespace Opm
// we need to check the BHP limit // we need to check the BHP limit
for (int p = 0; p < this->number_of_phases_; ++p) { for (int p = 0; p < this->number_of_phases_; ++p) {
const double temp = this->ipr_a_[p] - ipr_b_[p] * bhp_limit; const double temp = this->ipr_a_[p] - this->ipr_b_[p] * bhp_limit;
if (temp < 0.) { if (temp < 0.) {
this->operability_status_.operable_under_only_bhp_limit = false; this->operability_status_.operable_under_only_bhp_limit = false;
break; break;