OilfieldToolsNet/opm-simulators
4
0
Fork 0
mirror of https://github.com/OPM/opm-simulators.git synced 2025-02-10 02:55:33 -06:00
Code Issues Packages Projects Releases Wiki Activity

Fix bugs relating to misused face_cells variable.

Browse Source
This commit is contained in:
Atgeirr Flø Rasmussen 2014-04-15 20:49:47 +02:00
parent e4a6ab1c86
commit e9d6a0e8e0
1 changed files with 2 additions and 4 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

6
opm/autodiff/ImpesTPFAAD.cpp
View File

@ -61,7 +61,7 @@ namespace {
std::vector<int> f2hf(2 * numFaces(grid), -1); std::vector<int> f2hf(2 * numFaces(grid), -1);
Eigen::Array<int, Eigen::Dynamic, 2, Eigen::RowMajor> Eigen::Array<int, Eigen::Dynamic, 2, Eigen::RowMajor>
face_cells; face_cells = faceCells(grid);
typedef typename Opm::UgGridHelpers::Cell2FacesTraits<UnstructuredGrid>::Type typedef typename Opm::UgGridHelpers::Cell2FacesTraits<UnstructuredGrid>::Type
Cell2Faces; Cell2Faces;
@ -71,7 +71,7 @@ namespace {
cell_faces = c2f[c]; cell_faces = c2f[c];
typedef typename Cell2Faces::row_type::iterator Iter; typedef typename Cell2Faces::row_type::iterator Iter;
for (Iter f=cell_faces.begin(), end=cell_faces.end(); for (Iter f=cell_faces.begin(), end=cell_faces.end();
f!=end; ++end) { f!=end; ++f) {
const int p = 0 + (face_cells(*f,0) != c); const int p = 0 + (face_cells(*f,0) != c);
f2hf[2*(*f) + p] = f-c2f[0].begin(); f2hf[2*(*f) + p] = f-c2f[0].begin();
} }
@ -89,8 +89,6 @@ namespace {
std::vector<Tri> grav; grav.reserve(2 * ni); std::vector<Tri> grav; grav.reserve(2 * ni);
for (HelperOps::IFaces::Index i = 0; i < ni; ++i) { for (HelperOps::IFaces::Index i = 0; i < ni; ++i) {
const int f = ops.internal_faces[ i ]; const int f = ops.internal_faces[ i ];
Eigen::Array<int, Eigen::Dynamic, 2, Eigen::RowMajor>
face_cells=faceCells(grid);
const int c1 = face_cells(f,0); const int c1 = face_cells(f,0);
const int c2 = face_cells(f,1); const int c2 = face_cells(f,1);