OilfieldToolsNet/opm-simulators
4
0
Fork 0
mirror of https://github.com/OPM/opm-simulators.git synced 2024-11-29 04:23:48 -06:00
Code Issues Packages Projects Releases Wiki Activity

Make Darcy flux actually constant.

Browse Source 
The definition from commit 6f22d10 included a non-zero derivative
which interfered with the automatic differentiation solver.  Using a
fully constant Darcy flux enables using dt=150 while still computing
the correct solution.

While here, output the solution (s1) using more precision/digits for
comparison with 'sim_simple.m'.
This commit is contained in:
Bård Skaflestad 2013-05-05 22:02:40 +02:00
parent 6f22d10b06
commit 056b7ae292
1 changed files with 6 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

7
sim_simple.cpp
View File

@ -409,7 +409,12 @@ int main()
const std::vector<ADB> pmobc0 = phaseMobility<ADB>(props, allcells, s00.value());
const std::vector<ADB> pmobf0 = upws.select(p1, pmobc0);
const std::vector<ADB::M> null = { ADB::M(transi.size(), nc) };
const ADB dflux = (pmobf0[0] + pmobf0[1]) * ADB::function(nkdp.value(), null);
const ADB dflux = (ADB::function((pmobf0[0] + pmobf0[1]).value(), null) *
ADB::function(nkdp.value() , null));
std::cout.setf(std::ios::scientific);
std::cout.precision(16);
do {
const std::vector<int>& bp = block_pattern;
ADB s = ADB::variable(0, s1, bp);