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

Follow change to IncompTpfa (adding BCs).

Browse Source
This commit is contained in:
Atgeirr Flø Rasmussen 2012-03-07 09:38:09 +01:00
parent d12d444411
commit 93a747403e
1 changed files with 5 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

6
examples/polymer_reorder.cpp
View File

@ -21,6 +21,7 @@
#include "config.h"
#include <opm/core/pressure/IncompTpfa.hpp>
#include <opm/core/pressure/FlowBCManager.hpp>
#include <opm/core/grid.h>
#include <opm/core/GridManager.hpp>
@ -394,6 +395,9 @@ main(int argc, char** argv)
Opm::TransportModelPolymer tmodel(*grid->c_grid(), props->porosity(), &porevol[0], *props, polydata,
method, nltol, maxit);
// Boundary conditions.
Opm::FlowBCManager bcs;
// State-related and source-related variables init.
int num_cells = grid->c_grid()->number_of_cells;
std::vector<double> totmob;
@ -466,7 +470,7 @@ main(int argc, char** argv)
totmob);
}
pressure_timer.start();
psolver.solve(totmob, omega, src, state.pressure(), state.faceflux());
psolver.solve(totmob, omega, src, bcs.c_bcs(), state.pressure(), state.faceflux());
pressure_timer.stop();
double pt = pressure_timer.secsSinceStart();
std::cout << "Pressure solver took: " << pt << " seconds." << std::endl;