OilfieldToolsNet/opm-core
4
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

Fill semantic hole in initIteration(). Store the current sat-change.

Browse Source 
This is needed to implement the accumulation term of the residual.
This commit is contained in:
Bård Skaflestad
2011-10-06 14:57:46 +02:00
parent 255a0cf6ee
commit b867283f05
1 changed files with 4 additions and 1 deletions
Download Patch File Download Diff File Expand all files Collapse all files

5
src/SinglePointUpwindTwoPhase.hpp
View File

@@ -270,10 +270,13 @@ namespace Opm {
std::array<double, 2 > mob ;
std::array<double, 2*2> dmob;
const ::std::vector<double>& x = sys.vector().solution();
const ::std::vector<double>& sat = state.saturation();
for (int c = 0; c < g.number_of_cells; ++c) {
s[0] = sat[c*2 + 0] + sys.vector().solution()[c];
store_.ds(c) = x[c]; // Store sat-change for accumulation().
s[0] = sat[c*2 + 0] + x[c];
s[1] = 1 - s[0];
fluid_.mobility(c, s, mob, dmob);