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

Do not require unnecessary high quadrature precision.

Browse Source
This commit is contained in:
Atgeirr Flø Rasmussen 2013-01-15 13:33:55 +01:00
parent 7f05b1bb0e
commit 81c8870ba8
1 changed files with 3 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

4
opm/core/transport/reorder/TransportModelTracerTofDiscGal.cpp
View File

@ -329,7 +329,9 @@ namespace Opm
// Compute cell jacobian contribution. We use Fortran ordering // Compute cell jacobian contribution. We use Fortran ordering
// for jac_, i.e. rows cycling fastest. // for jac_, i.e. rows cycling fastest.
{ {
const int deg_needed = use_cvi_ ? 2*degree_ : 2*degree_ - 1; // Even with ECVI velocity interpolation, degree of precision 1
// is sufficient for optimal convergence order for DG1.
const int deg_needed = 2*degree_ - 1;
CellQuadrature quad(grid_, cell, deg_needed); CellQuadrature quad(grid_, cell, deg_needed);
for (int quad_pt = 0; quad_pt < quad.numQuadPts(); ++quad_pt) { for (int quad_pt = 0; quad_pt < quad.numQuadPts(); ++quad_pt) {
// b_i (v \cdot \grad b_j) // b_i (v \cdot \grad b_j)