OilfieldToolsNet/opm-simulators
4
0
Fork 0
mirror of https://github.com/OPM/opm-simulators.git synced 2025-01-18 17:52:58 -06:00
Code Issues Packages Projects Releases Wiki Activity

Support variable number of physical grid dimensions (typically 2 or 3).

Browse Source
This commit is contained in:
Bård Skaflestad 2012-03-09 08:56:28 +01:00
parent 29fa476715
commit 4fc45dea82
1 changed files with 18 additions and 5 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

23
opm/core/pressure/FlowBCManager.cpp
View File

@ -169,15 +169,21 @@ namespace Opm
void cartCoord(const int log_cart_coord, void cartCoord(const int ndims,
const int log_cart_coord,
const int* dims, const int* dims,
int* ijk) int* ijk)
{ {
int ix = log_cart_coord; int ix = log_cart_coord;
ijk[0] = ix % dims[0]; ix /= dims[0]; for (int dim = 0; dim < ndims; ++dim) {
ijk[1] = ix % dims[1]; ix /= dims[1]; ijk[dim] = ix % dims[dim];
ijk[2] = ix % dims[2]; ix /= dims[dim];
}
ASSERT2 (ix == 0,
"Lexicographic index is not consistent "
"with grid dimensions.");
} }
@ -195,13 +201,20 @@ namespace Opm
THROW("Faces not tagged - cannot extract " << sideString(side) << " faces."); THROW("Faces not tagged - cannot extract " << sideString(side) << " faces.");
} }
ASSERT2 (grid.dimensions <= 3,
"Grid must have three dimensions or less.");
ASSERT2 (side < 2 * grid.dimensions,
"Boundary condition side not consistent with "
"number of physical grid dimensions.");
// Get all boundary faces with the correct tag and with // Get all boundary faces with the correct tag and with
// min/max i/j/k (depending on side). // min/max i/j/k (depending on side).
const int correct_ijk = (side % 2) ? grid.cartdims[side/2] : 0; const int correct_ijk = (side % 2) ? grid.cartdims[side/2] : 0;
for (int c = 0; c < grid.number_of_cells; ++c) { for (int c = 0; c < grid.number_of_cells; ++c) {
int ijk[3] = { -1, -1, -1 }; int ijk[3] = { -1, -1, -1 };
int gc = (grid.global_cell != 0) ? grid.global_cell[c] : c; int gc = (grid.global_cell != 0) ? grid.global_cell[c] : c;
cartCoord(gc, grid.cartdims, ijk); cartCoord(grid.dimensions, gc, grid.cartdims, ijk);
if (ijk[side/2] != correct_ijk) { if (ijk[side/2] != correct_ijk) {
continue; continue;
} }