OilfieldToolsNet/opm-core
4
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

Made use of the new LinearSolverFactory in wells_example

Browse Source
This commit is contained in:
Kjetil Olsen Lye 2012-04-16 12:18:37 +02:00
parent 6877cc8cd1
commit 11f0a14e61
1 changed files with 4 additions and 10 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

14
examples/wells_example.cpp
View File

@ -8,16 +8,12 @@
#include <opm/core/GridManager.hpp> #include <opm/core/GridManager.hpp>
#include <opm/core/pressure/IncompTpfa.hpp> #include <opm/core/pressure/IncompTpfa.hpp>
#include <opm/core/fluid/IncompPropertiesFromDeck.hpp> #include <opm/core/fluid/IncompPropertiesFromDeck.hpp>
#include <opm/core/linalg/LinearSolverUmfpack.hpp>
#include <opm/core/newwells.h> #include <opm/core/newwells.h>
#include <opm/core/grid.h> #include <opm/core/grid.h>
#include <opm/core/utility/miscUtilities.hpp> #include <opm/core/utility/miscUtilities.hpp>
#include <opm/core/TwophaseState.hpp> #include <opm/core/TwophaseState.hpp>
#include <opm/core/pressure/FlowBCManager.hpp> #include <opm/core/pressure/FlowBCManager.hpp>
#include <opm/core/linalg/LinearSolverFactory.hpp>
#ifdef EXPERIMENT_ISTL
#include <opm/core/linalg/LinearSolverIstl.hpp>
#endif
int main(int argc, char** argv) { int main(int argc, char** argv) {
using namespace Opm::parameter; using namespace Opm::parameter;
@ -39,11 +35,9 @@ int main(int argc, char** argv) {
double gravity[3] = {0.0, 0.0, parameters.getDefault<double>("gravity", 0.0)}; double gravity[3] = {0.0, 0.0, parameters.getDefault<double>("gravity", 0.0)};
IncompPropertiesFromDeck incomp_properties(parser, global_cells); IncompPropertiesFromDeck incomp_properties(parser, global_cells);
#ifdef EXPERIMENT_ISTL Opm::LinearSolverFactory linsolver(parameters);
Opm::LinearSolverIstl linsolver(parameters);
#else // EXPERIMENT_ISTL
Opm::LinearSolverUmfpack linsolver;
#endif // EXPERIMENT_ISTL
IncompTpfa pressure_solver(*grid.c_grid(), incomp_properties.permeability(), IncompTpfa pressure_solver(*grid.c_grid(), incomp_properties.permeability(),
gravity, linsolver, wells.c_wells()); gravity, linsolver, wells.c_wells());