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 issues in init().

Browse Source 
- Failed to set all rates to zero initially.
 - Could set out-of-bounds index if number of controls changed.
This commit is contained in:
Atgeirr Flø Rasmussen 2015-01-26 16:46:51 +01:00
parent 59f9b58d54
commit ce3519dadb
1 changed files with 8 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

10
opm/autodiff/WellStateFullyImplicitBlackoil.hpp
View File

@ -72,6 +72,8 @@ namespace Opm
// Initialize perfphaserates_, which must be done here. // Initialize perfphaserates_, which must be done here.
const int np = wells->number_of_phases; const int np = wells->number_of_phases;
const int nperf = wells->well_connpos[nw]; const int nperf = wells->well_connpos[nw];
// Ensure that we start out with zero rates by default.
perfphaserates_.clear();
perfphaserates_.resize(nperf * np, 0.0); perfphaserates_.resize(nperf * np, 0.0);
for (int w = 0; w < nw; ++w) { for (int w = 0; w < nw; ++w) {
assert((wells->type[w] == INJECTOR) || (wells->type[w] == PRODUCER)); assert((wells->type[w] == INJECTOR) || (wells->type[w] == PRODUCER));
@ -151,8 +153,12 @@ namespace Opm
} }
// currentControls // currentControls
// WARNING: This may be error prone if the number of controls change. const int old_control_index = prevState.currentControls()[ oldIndex ];
currentControls()[ newIndex ] = prevState.currentControls()[ oldIndex ]; if (old_control_index < well_controls_get_num(wells->ctrls[w])) {
// If the set of controls have changed, this may not be identical
// to the last control, but it must be a valid control.
currentControls()[ newIndex ] = old_control_index;
}
} }
} }