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

Merge pull request #1317 from GitPaean/throw_when_no_umfpack

Browse Source 
giving more information when umfpack is not available.
This commit is contained in:
Atgeirr Flø Rasmussen
2017-11-03 14:43:18 +01:00
committed by GitHub
parent 63415ebd96 d9ed5a36b5
commit 1d0282c6fe
3 changed files with 19 additions and 19 deletions
Download Patch File Download Diff File Expand all files Collapse all files

7
opm/autodiff/MSWellHelpers.hpp
View File

@@ -33,12 +33,12 @@ namespace Opm {
namespace mswellhelpers
{
#if HAVE_UMFPACK
// obtain y = D^-1 * x with a direct solver
template <typename MatrixType, typename VectorType>
VectorType
invDXDirect(const MatrixType& D, VectorType x)
{
#if HAVE_UMFPACK
VectorType y(x.size());
y = 0.;
@@ -61,8 +61,11 @@ namespace mswellhelpers
}
return y;
}
#else
OPM_THROW(std::runtime_error, "Cannot use invDXDirect() without UMFPACK. "
"Reconfigure opm-simulator with SuiteSparse/UMFPACK support and recompile.");
#endif // HAVE_UMFPACK
}

8
opm/autodiff/MultisegmentWell_impl.hpp
View File

@@ -494,11 +494,7 @@ namespace Opm
duneB_.mv(x, Bx);
// invDBx = duneD^-1 * Bx_
#if HAVE_UMFPACK
const BVectorWell invDBx = mswellhelpers::invDXDirect(duneD_, Bx);
#else
const BVectorWell invDBx = mswellhelpers::invDX(duneD_, Bx);
#endif
// Ax = Ax - duneC_^T * invDBx
duneC_.mmtv(invDBx,Ax);
@@ -514,11 +510,7 @@ namespace Opm
apply(BVector& r) const
{
// invDrw_ = duneD^-1 * resWell_
#if HAVE_UMFPACK
const BVectorWell invDrw = mswellhelpers::invDXDirect(duneD_, resWell_);
#else
const BVectorWell invDrw = mswellhelpers::invDX(duneD_, resWell_);
#endif // HAVE_UMFPACK
// r = r - duneC_^T * invDrw
duneC_.mmtv(invDrw, r);
}