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

[bugfix] use SparseLU when UMFPack was not found.

Browse Source
This commit is contained in:
Robert K 2014-11-19 11:22:57 +01:00
parent 636cb3c52c
commit 9cffb51543
1 changed files with 12 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

12
opm/autodiff/NewtonIterationBlackoilCPR.cpp
View File

@ -44,7 +44,11 @@
#include <opm/core/utility/platform_dependent/reenable_warnings.h>
#if HAVE_UMFPACK
#include <Eigen/UmfPackSupport>
#else
#include <Eigen/SparseLU>
#endif
namespace Opm
{
@ -249,7 +253,11 @@ namespace Opm
const std::vector<M>& Jn = eqs[n].derivative();
// Use sparse LU to solve the block submatrices i.e compute inv(D)
#if HAVE_UMFPACK
const Eigen::UmfPackLU< M > solver(Jn[n]);
#else
const Eigen::SparseLU< M > solver(Jn[n]);
#endif
M id(Jn[n].rows(), Jn[n].cols());
id.setIdentity();
const M Di = solver.solve(id);
@ -323,7 +331,11 @@ namespace Opm
const M& C = eq_coll.derivative()[0];
// Use sparse LU to solve the block submatrices
#if HAVE_UMFPACK
const Eigen::UmfPackLU< M > solver(D);
#else
const Eigen::SparseLU< M > solver(D);
#endif
// Compute value of eliminated variable.
const Eigen::VectorXd b = (equation.value().matrix() - C * partial_solution.matrix());