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

Print face fluxes after solve.

Browse Source 
Also define the HACK_INCOMPRESSIBLE_GRAVITY symbol, since test case
is meaningless without it.
This commit is contained in:
Atgeirr Flø Rasmussen 2013-05-16 14:49:12 +02:00
parent f6c6ee0fa0
commit 37677fe032
1 changed files with 7 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

7
examples/test_impestpfa_ad.cpp
View File

@ -20,6 +20,8 @@
#include <config.h> #include <config.h>
#define HACK_INCOMPRESSIBLE_GRAVITY 1
#include <opm/autodiff/ImpesTPFAAD.hpp> #include <opm/autodiff/ImpesTPFAAD.hpp>
#include <opm/core/grid.h> #include <opm/core/grid.h>
@ -140,8 +142,13 @@ main(int argc, char* argv[])
ps.solve(1.0, state, well_state); ps.solve(1.0, state, well_state);
std::cout << "Cell pressure:" << std::endl;
std::copy(state.pressure().begin(), state.pressure().end(), std::ostream_iterator<double>(std::cout, " ")); std::copy(state.pressure().begin(), state.pressure().end(), std::ostream_iterator<double>(std::cout, " "));
std::cout << std::endl; std::cout << std::endl;
std::cout << "Face flux:" << std::endl;
std::copy(state.faceflux().begin(), state.faceflux().end(), std::ostream_iterator<double>(std::cout, " "));
std::cout << std::endl;
std::cout << "Well bhp pressure:" << std::endl;
std::copy(well_state.bhp().begin(), well_state.bhp().end(), std::ostream_iterator<double>(std::cout, " ")); std::copy(well_state.bhp().begin(), well_state.bhp().end(), std::ostream_iterator<double>(std::cout, " "));
std::cout << std::endl; std::cout << std::endl;