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

Merge pull request #1198 from dr-robertk/PR/fix-densewell-for-2p

Browse Source 
[bugfix] Fix invalid read in WellDenseStandard when used with 2P Model.
This commit is contained in:
Tor Harald Sandve 2017-06-07 13:57:44 +02:00 committed by GitHub
commit 2a66aa23a9
1 changed files with 12 additions and 5 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

17
opm/autodiff/StandardWellsDense_impl.hpp
View File

@ -700,11 +700,18 @@ namespace Opm {
setWellVariables(const WellState& xw)
{
const int nw = wells().number_of_wells;
for (int eqIdx = 0; eqIdx < numWellEq; ++eqIdx) {
// for two-phase numComp < numEq
const int numComp = numComponents();
for (int eqIdx = 0; eqIdx < numComp; ++eqIdx) {
for (int w = 0; w < nw; ++w) {
wellVariables_[w + nw*eqIdx] = 0.0;
wellVariables_[w + nw*eqIdx].setValue(xw.wellSolutions()[w + nw* eqIdx]);
wellVariables_[w + nw*eqIdx].setDerivative(numWellEq + eqIdx, 1.0);
const unsigned int idx = nw * eqIdx + w;
assert( idx < wellVariables_.size() );
assert( idx < xw.wellSolutions().size() );
EvalWell& eval = wellVariables_[ idx ];
eval = 0.0;
eval.setValue( xw.wellSolutions()[ idx ] );
eval.setDerivative(numWellEq + eqIdx, 1.0);
}
}
}
@ -1025,7 +1032,7 @@ namespace Opm {
{
bval/=global_nc_;
}
auto res = residual();
const int nw = res.size() / numComp;