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

Fix bug in access order for state.qs.

Browse Source 
State.qs is called well_phase_flow_rate in this function.

It is stored with wells running fastest, not phases.
This commit is contained in:
Atgeirr Flø Rasmussen 2014-03-25 17:42:14 +01:00
parent cb50728ac2
commit dcb59a2fab
1 changed files with 4 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

5
opm/autodiff/FullyImplicitBlackoilSolver.cpp
View File

@ -945,9 +945,12 @@ namespace {
const int num_phases, const int num_phases,
const double* distr) const double* distr)
{ {
const int num_wells = well_phase_flow_rate.size() / num_phases;
double rate = 0.0; double rate = 0.0;
for (int phase = 0; phase < num_phases; ++phase) { for (int phase = 0; phase < num_phases; ++phase) {
rate += well_phase_flow_rate.value()[well*num_phases + phase] * distr[phase]; // Important: well_phase_flow_rate is ordered with all rates for first
// phase coming first, then all for second phase etc.
rate += well_phase_flow_rate.value()[well + phase*num_wells] * distr[phase];
} }
return rate; return rate;
} }