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

Bugfix: use sufficient quadrature order.

Browse Source
This commit is contained in:
Atgeirr Flø Rasmussen 2012-12-18 14:12:56 +01:00
parent 5e04ec1448
commit 94b128d6b6
1 changed files with 1 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

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

@ -242,7 +242,7 @@ namespace Opm
// \int_{\partial K} u_h^{ext} (v(x) \cdot n) b_j ds // \int_{\partial K} u_h^{ext} (v(x) \cdot n) b_j ds
// (where u_h^{ext} is the upstream unknown (tof)). // (where u_h^{ext} is the upstream unknown (tof)).
const double normal_velocity = flux / grid_.face_areas[face]; const double normal_velocity = flux / grid_.face_areas[face];
FaceQuadrature quad(grid_, face, degree_); FaceQuadrature quad(grid_, face, 2*degree_);
for (int quad_pt = 0; quad_pt < quad.numQuadPts(); ++quad_pt) { for (int quad_pt = 0; quad_pt < quad.numQuadPts(); ++quad_pt) {
quad.quadPtCoord(quad_pt, &coord_[0]); quad.quadPtCoord(quad_pt, &coord_[0]);
DGBasis::eval(grid_, cell, degree_, &coord_[0], &basis_[0]); DGBasis::eval(grid_, cell, degree_, &coord_[0], &basis_[0]);