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

Use UMFPACK solver if available in sim_simple.cpp.

Browse Source 
Fall back to Eigen's BiCGSTAB solver if UMFPACK is not available.  The
BiCGSTAB is built into Eigen and consequently always available when
Eigen is found.
This commit is contained in:
Bård Skaflestad 2013-07-09 14:39:44 +02:00
parent e9ad8f0b5e
commit a9b90b34f8
2 changed files with 12 additions and 7 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

6
CMakeLists.txt
View File

@ -132,12 +132,6 @@ opm_out_dirs ()
# tests files in tests/, examples in examples/
opm_sources (${project})
if (NOT SuiteSparse_FOUND)
list(REMOVE_ITEM examples_SOURCES
${PROJECT_SOURCE_DIR}/examples/sim_simple.cpp
)
endif (NOT SuiteSparse_FOUND)
# Create configuration header which describes available features
# necessary to compile this library. Singular version is the names that
# is required by this project alone, plural version transitively

13
examples/sim_simple.cpp
View File

@ -27,7 +27,12 @@
#include <opm/core/utility/Units.hpp>
#include <opm/core/utility/StopWatch.hpp>
#include <opm/core/pressure/tpfa/trans_tpfa.h>
#if HAVE_SUITESPARSE_UMFPACK_H
#include <Eigen/UmfPackSupport>
#else
#include <Eigen/IterativeLinearSolvers>
#endif
#include <iostream>
#include <cstdlib>
@ -202,7 +207,13 @@ int main()
// Where R(p0) and J(p0) are contained in residual.value() and
// residual.derived()[0].
Eigen::BiCGSTAB<M> solver;
#if HAVE_SUITESPARSE_UMFPACK_H
typedef Eigen::UmfPackLU<M> LinSolver;
#else
typedef Eigen::BiCGSTAB<M> LinSolver;
#endif // HAVE_SUITESPARSE_UMFPACK_H
LinSolver solver;
M pmatr = residual.derivative()[0];
pmatr.coeffRef(0,0) *= 2.0;
pmatr.makeCompressed();