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

Re-introduced tempoary fix for singular matrix.

Browse Source
This commit is contained in:
osae
2014-07-03 12:16:00 +02:00
parent 3512a7b2f0
commit bb12bdd1fd
1 changed files with 12 additions and 0 deletions
Download Patch File Download Diff File Expand all files Collapse all files

12
opm/autodiff/FullyImplicitBlackoilSolver_impl.hpp
View File

@@ -233,6 +233,18 @@ namespace {
std::vector<std::vector<double>> residual_history; std::vector<std::vector<double>> residual_history;
// do a zero update once to switch to the right primary
// variables. this is a bit hacky, but it is by far the
// easiest thing to do...
int numCells = Opm::AutoDiffGrid::numCells(grid_);
int numPhases = x.numPhases();
int numWells = xw.basicWellState().bhp().size();
int numDofs = numPhases * (numCells + numWells) + numWells;
V zero = V::Zero(numDofs);
updateState(zero, x, xw);
assemble(pvdt, x, xw); assemble(pvdt, x, xw);