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

NewtonIterationBlackoilInterleaved: single precision can be used whether UMFpack is

Browse Source 
there or not.
This commit is contained in:
Robert Kloefkorn 2016-02-09 16:07:56 +01:00
parent 4111797822
commit e3b578d685
1 changed files with 1 additions and 7 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

8
opm/autodiff/NewtonIterationBlackoilInterleaved.cpp
View File

@ -541,15 +541,9 @@ namespace Opm
{ {
// get np and call appropriate template method // get np and call appropriate template method
const int np = residual.material_balance_eq.size(); const int np = residual.material_balance_eq.size();
#if ! HAVE_UMFPACK const NewtonIterationBlackoilInterface& newtonIncrement = residual.singlePrecision ?
const bool singlePrecision = residual.singlePrecision ;
const NewtonIterationBlackoilInterface& newtonIncrement = singlePrecision ?
detail::NewtonIncrement< maxNumberEquations_, float > :: get( newtonIncrementSinglePrecision_, parameters_, parallelInformation_, np ) : detail::NewtonIncrement< maxNumberEquations_, float > :: get( newtonIncrementSinglePrecision_, parameters_, parallelInformation_, np ) :
detail::NewtonIncrement< maxNumberEquations_, double > :: get( newtonIncrementDoublePrecision_, parameters_, parallelInformation_, np ); detail::NewtonIncrement< maxNumberEquations_, double > :: get( newtonIncrementDoublePrecision_, parameters_, parallelInformation_, np );
#else
const NewtonIterationBlackoilInterface& newtonIncrement =
detail::NewtonIncrement< maxNumberEquations_, double > :: get( newtonIncrementDoublePrecision_, parameters_, parallelInformation_, np );
#endif
// compute newton increment // compute newton increment
SolutionVector dx = newtonIncrement.computeNewtonIncrement( residual ); SolutionVector dx = newtonIncrement.computeNewtonIncrement( residual );