OilfieldToolsNet/opm-core
4
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

Prettified.

Browse Source
This commit is contained in:
Atgeirr Flø Rasmussen 2012-03-01 15:12:50 +01:00
parent c21299085e
commit 565430c2e9
1 changed files with 17 additions and 21 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

38
opm/core/ColumnExtract.hpp
View File

@ -2,28 +2,24 @@
#include <vector>
namespace Opm {
/**
* Extract each column of the grid.
* \note Assumes the pillars of the grid are all vertically aligned.
* \param grid The grid from which to extract the columns.
* \param columns will for each i + cartgrim[0]*j contain the k values contained in the column
* centered at (i, j).
*/
void extractColumn( const UnstructuredGrid& grid, std::vector<std::vector<int> >& columns ) {
columns.resize(grid.cartdims[0]*grid.cartdims[1]);
for(int i = 0; i < grid.cartdims[0]; i++) {
for(int j = 0; j < grid.cartdims[1]; j++) {
int plane_index = i + j*grid.cartdims[0];
for(int k = 0; k < grid.cartdims[2]; k++) {
columns[plane_index].push_back(plane_index + k*grid.cartdims[0]*grid.cartdims[1]);
}
}
/// Extract each column of the grid.
/// \note Assumes the pillars of the grid are all vertically aligned.
/// \param grid The grid from which to extract the columns.
/// \param columns will for each i + cartgrim[0]*j contain the k values contained in the column
/// centered at (i, j).
void extractColumn( const UnstructuredGrid& grid, std::vector<std::vector<int> >& columns )
{
const int* dims = grid.cartdims;
columns.resize(dims[0]*dims[1]);
for (int i = 0; i < dims[0]; ++i) {
for (int j = 0; j < dims[1]; ++j) {
int plane_index = i + j*dims[0];
for (int k = 0; k < dims[2]; ++k) {
columns[plane_index].push_back(plane_index + k*dims[0]*dims[1]);
}
}
}
}
}
} // namespace Opm