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

Added limiter for DG1, parameter 'use_limiter'.

Browse Source 
The limiter is experimental and unfinished, untested work in progress.
Limiter is therefore inactive by default.

Also fixed a minor bug: use_cvi_ was not initialized.
This commit is contained in:
Atgeirr Flø Rasmussen
2012-12-18 14:15:31 +01:00
parent 94b128d6b6
commit 3e723bc965
3 changed files with 108 additions and 7 deletions
Download Patch File Download Diff File Expand all files Collapse all files

4
examples/compute_tof_from_files.cpp
View File

@@ -124,10 +124,12 @@ main(int argc, char** argv)
bool use_dg = param.getDefault("use_dg", false);
int dg_degree = -1;
bool use_cvi = false;
bool use_limiter = false;
bool use_multidim_upwind = false;
if (use_dg) {
dg_degree = param.getDefault("dg_degree", 0);
use_cvi = param.getDefault("use_cvi", false);
use_limiter = param.getDefault("use_limiter", false);
} else {
use_multidim_upwind = param.getDefault("use_multidim_upwind", false);
}
@@ -157,7 +159,7 @@ main(int argc, char** argv)
transport_timer.start();
std::vector<double> tof;
if (use_dg) {
Opm::TransportModelTracerTofDiscGal tofsolver(grid, use_cvi);
Opm::TransportModelTracerTofDiscGal tofsolver(grid, use_cvi, use_limiter);
tofsolver.solveTof(&flux[0], &porevol[0], &src[0], dg_degree, tof);
} else {
Opm::TransportModelTracerTof tofsolver(grid, use_multidim_upwind);