OilfieldToolsNet/opm-simulators
4
0
Fork 0
mirror of https://github.com/OPM/opm-simulators.git synced 2024-12-26 01:01:00 -06:00
Code Issues Packages Projects Releases Wiki Activity

lowering the singular limit for the invertMatrix function.

Browse Source 
to avoid failure due to singularity detection.
This commit is contained in:
Kai Bao 2018-05-28 15:54:16 +02:00
parent 97c51a7d02
commit 53ab8bc132
3 changed files with 1 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

1
opm/autodiff/ISTLSolver.hpp
View File

@ -214,6 +214,7 @@ static inline void invertMatrix (FieldMatrix<K,4,4> &matrix)
template <typename K, int n>
static inline void invertMatrix (FieldMatrix<K,n,n> &matrix)
{
Dune::FMatrixPrecision<K>::set_singular_limit(1.e-20);
matrix.invert();
}

1
opm/autodiff/StandardWell.hpp
View File

@ -27,7 +27,6 @@
#include <opm/autodiff/WellInterface.hpp>
#include <opm/autodiff/ISTLSolver.hpp>
#include <opm/autodiff/RateConverter.hpp>
#include <opm/autodiff/ISTLSolver.hpp>
namespace Opm
{

1
opm/autodiff/StandardWell_impl.hpp
View File

@ -1073,7 +1073,6 @@ namespace Opm
// If it is under THP control, it will be set to be the target value.
// TODO: a better standard is probably whether we have the table to calculate the THP value
// TODO: it is something we need to check the output to decide.
const PhaseUsage& pu = phaseUsage();
const WellControls* wc = well_controls_;
// TODO: we should only maintain one current control either from the well_state or from well_controls struct.
// Either one can be more favored depending on the final strategy for the initilzation of the well control