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

Fixed namespacing issue

Browse Source
This commit is contained in:
Kjetil Olsen Lye 2012-04-10 14:47:29 +02:00
parent 88d6a626e3
commit 7db1027935
1 changed files with 17 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

19
examples/wells_example.cpp
View File

@ -3,8 +3,11 @@
#include <opm/core/utility/parameters/ParameterGroup.hpp> #include <opm/core/utility/parameters/ParameterGroup.hpp>
#include <opm/core/WellsManager.hpp> #include <opm/core/WellsManager.hpp>
#include <opm/core/GridManager.hpp> #include <opm/core/GridManager.hpp>
#include <opm/core/pressure/IncompTpfa.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"
int main(int argc, char** argv) { int main(int argc, char** argv) {
@ -12,7 +15,7 @@ int main(int argc, char** argv) {
using namespace Opm; using namespace Opm;
ParameterGroup parameters( argc, argv, false ); ParameterGroup parameters( argc, argv, false );
std::string file_name = parameters.getDefault<std::string>("inputdeck", "data.data"); std::string file_name = parameters.getDefault<std::string>("inputdeck", "data.data");
// Read input file // Read input file
EclipseGridParser parser(file_name); EclipseGridParser parser(file_name);
std::cout << "Done!" << std::endl; std::cout << "Done!" << std::endl;
@ -22,6 +25,18 @@ int main(int argc, char** argv) {
// Finally handle the wells // Finally handle the wells
WellsManager wells(parser, *grid.c_grid(), NULL); WellsManager wells(parser, *grid.c_grid(), NULL);
std::vector<int> global_cells(grid.c_grid()->global_cell, grid.c_grid()->global_cell + grid.c_grid()->number_of_cells);
std::cout << "ahoi" << std::endl;
double gravity[3] = {0.0, 0.0, parameters.getDefault<double>("gravity", 0.0)};
IncompPropertiesFromDeck incomp_properties(parser, global_cells);
std::cout << "there" << std::endl;
LinearSolverUmfpack umfpack_solver;
std::cout << "here" << std::endl;
IncompTpfa pressure_solver(*grid.c_grid(), incomp_properties.permeability(),
gravity, umfpack_solver, wells.c_wells());
return 0; return 0;
} }