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

do use eigen's umfpack support anymore

Browse Source 
use BiCGStab instead. it might be slower, but at least it is always
available...
This commit is contained in:
Andreas Lauser
2013-07-09 13:33:57 +02:00
parent df2964729a
commit 7139f4addf
1 changed files with 2 additions and 2 deletions
Download Patch File Download Diff File Expand all files Collapse all files

4
examples/sim_simple.cpp
View File

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