My Project
|
Struct that hods all the data needed to represent a Cpgrid. More...
#include <CpGridData.hpp>
Public Types | |
enum | { MAX_DATA_PER_CELL = 2 } |
using | MPICommunicator = CpGridDataTraits::MPICommunicator |
The type of the mpi communicator. | |
using | Communication = CpGridDataTraits::Communication |
The type of the collective communication. | |
using | CollectiveCommunication = CpGridDataTraits::CollectiveCommunication |
Public Member Functions | |
CpGridData (MPIHelper::MPICommunicator comm, std::vector< std::shared_ptr< CpGridData > > &data) | |
Constructor for parallel grid data. | |
CpGridData (std::vector< std::shared_ptr< CpGridData > > &data) | |
Constructor. | |
~CpGridData () | |
Destructor. | |
int | size (int codim) const |
number of leaf entities per codim in this process | |
int | size (GeometryType type) const |
number of leaf entities per geometry type in this process | |
void | readSintefLegacyFormat (const std::string &grid_prefix) |
Read the Sintef legacy grid format ('topogeom'). | |
void | writeSintefLegacyFormat (const std::string &grid_prefix) const |
Write the Sintef legacy grid format ('topogeom'). | |
void | readEclipseFormat (const std::string &filename, bool periodic_extension, bool turn_normals=false) |
Read the Eclipse grid format ('grdecl'). | |
void | processEclipseFormat (const grdecl &input_data, std::array< std::set< std::pair< int, int > >, 2 > &nnc, bool remove_ij_boundary, bool turn_normals, bool pinchActive, double tolerance_unique_points) |
Read the Eclipse grid format ('grdecl'). | |
void | getIJK (int c, std::array< int, 3 > &ijk) const |
Extract Cartesian index triplet (i,j,k) of an active cell. | |
std::tuple< const std::shared_ptr< CpGridData >, const std::vector< std::array< int, 2 > >, const std::vector< std::tuple< int, std::vector< int > > >, const std::tuple< int, std::vector< int > >, const std::vector< std::array< int, 2 > >, const std::vector< std::array< int, 2 > > > | refineSingleCell (const std::array< int, 3 > &cells_per_dim, const int &parent_idx) const |
Refine a single cell and return a shared pointer of CpGridData type. | |
std::tuple< std::shared_ptr< CpGridData >, const std::vector< std::array< int, 2 > >, const std::vector< std::tuple< int, std::vector< int > > >, const std::vector< std::tuple< int, std::vector< int > > >, const std::vector< std::tuple< int, std::vector< int > > >, const std::vector< std::array< int, 2 > >, const std::vector< std::array< int, 2 > > > | refinePatch (const std::array< int, 3 > &cells_per_dim, const std::array< int, 3 > &startIJK, const std::array< int, 3 > &endIJK) const |
Refine a (connected block-shaped) patch of cells. | |
std::array< double, 3 > | computeEclCentroid (const int idx) const |
std::array< double, 3 > | computeEclCentroid (const Entity< 0 > &elem) const |
void | computeUniqueBoundaryIds () |
bool | uniqueBoundaryIds () const |
Is the grid currently using unique boundary ids? | |
void | setUniqueBoundaryIds (bool uids) |
Set whether we want to have unique boundary ids. | |
const std::vector< double > & | zcornData () const |
Return the internalized zcorn copy from the grid processing, if no cells were adjusted during the minpvprocessing this can be and empty vector. | |
const IndexSet & | indexSet () const |
Get the index set. | |
const std::array< int, 3 > & | logicalCartesianSize () const |
The logical cartesian size of the grid. | |
void | distributeGlobalGrid (CpGrid &grid, const CpGridData &view_data, const std::vector< int > &cell_part) |
Redistribute a global grid. | |
template<class DataHandle > | |
void | communicate (DataHandle &data, InterfaceType iftype, CommunicationDirection dir) |
communicate objects for all codims on a given level | |
const std::vector< int > & | sortedNumAquiferCells () const |
Get sorted active cell indices of numerical aquifer. | |
Struct that hods all the data needed to represent a Cpgrid.
anonymous enum |
|
explicit |
Constructor for parallel grid data.
comm | The MPI communicator Default constructor. |
void Dune::cpgrid::CpGridData::communicate | ( | DataHandle & | data, |
InterfaceType | iftype, | ||
CommunicationDirection | dir | ||
) |
communicate objects for all codims on a given level
data | The data handle describing the data. Has to adhere to the Dune::DataHandleIF interface. |
iftype | The interface to use for the communication. |
dir | The direction of the communication along the interface (forward or backward). |
void Dune::cpgrid::CpGridData::distributeGlobalGrid | ( | CpGrid & | grid, |
const CpGridData & | view_data, | ||
const std::vector< int > & | cell_part | ||
) |
Redistribute a global grid.
The whole grid must be available on all processors.
|
inline |
Extract Cartesian index triplet (i,j,k) of an active cell.
[in] | c | Active cell index. |
[out] | ijk | Cartesian index triplet |
|
inline |
Get the index set.
This is the lead as well as th level index set.
|
inline |
The logical cartesian size of the grid.
This function is not part of the Dune grid interface, and should be used with caution.
void Dune::cpgrid::CpGridData::processEclipseFormat | ( | const grdecl & | input_data, |
std::array< std::set< std::pair< int, int > >, 2 > & | nnc, | ||
bool | remove_ij_boundary, | ||
bool | turn_normals, | ||
bool | pinchActive, | ||
double | tolerance_unique_points | ||
) |
Read the Eclipse grid format ('grdecl').
Read the Eclipse grid format ('.grdecl').
input_data | the data in grdecl format, declared in preprocess.h. |
ecl_state | the object from opm-parser provide information regarding to pore volume, NNC, aquifer information when ecl_state is available. NNC and aquifer connection information will also be updated during the function call when available and necessary. |
nnc | is the non-neighboring connections |
remove_ij_boundary | if true, will remove (i, j) boundaries. Used internally. |
pinchActive | If true, we will add faces between vertical cells that have only inactive cells or cells with zero volume between them. If false these cells will not be connected. |
tolerance_unique_points | Tolerance used to identify points based on their cooridinate |
void Dune::cpgrid::CpGridData::readEclipseFormat | ( | const std::string & | filename, |
bool | periodic_extension, | ||
bool | turn_normals = false |
||
) |
Read the Eclipse grid format ('grdecl').
filename | the name of the file to read. |
periodic_extension | if true, the grid will be (possibly) refined, so that intersections/faces along i and j boundaries will match those on the other side. That is, i- faces will match i+ faces etc. |
void Dune::cpgrid::CpGridData::readSintefLegacyFormat | ( | const std::string & | grid_prefix | ) |
Read the Sintef legacy grid format ('topogeom').
grid_prefix | the grid name, such that topology is found in <grid_prefix>-topo.dat etc. |
std::tuple< std::shared_ptr< CpGridData >, const std::vector< std::array< int, 2 > >, const std::vector< std::tuple< int, std::vector< int > > >, const std::vector< std::tuple< int, std::vector< int > > >, const std::vector< std::tuple< int, std::vector< int > > >, const std::vector< std::array< int, 2 > >, const std::vector< std::array< int, 2 > > > Dune::cpgrid::CpGridData::refinePatch | ( | const std::array< int, 3 > & | cells_per_dim, |
const std::array< int, 3 > & | startIJK, | ||
const std::array< int, 3 > & | endIJK | ||
) | const |
Refine a (connected block-shaped) patch of cells.
Based on the patch, a Geometry<3,3> object is created and refined.
[in] | cells_per_dim | Number of (refined) cells in each direction that each parent cell should be refined to. |
[in] | startIJK | Cartesian triplet index where the patch starts. |
[in] | endIJK | Cartesian triplet index where the patch ends. Last cell part of the lgr will be {endijk[0]-1, ... endIJK[2]-1}. |
std::tuple< const std::shared_ptr< CpGridData >, const std::vector< std::array< int, 2 > >, const std::vector< std::tuple< int, std::vector< int > > >, const std::tuple< int, std::vector< int > >, const std::vector< std::array< int, 2 > >, const std::vector< std::array< int, 2 > > > Dune::cpgrid::CpGridData::refineSingleCell | ( | const std::array< int, 3 > & | cells_per_dim, |
const int & | parent_idx | ||
) | const |
Refine a single cell and return a shared pointer of CpGridData type.
refineSingleCell() takes a cell and refines it in a chosen amount of cells (per direction); creating the geometries, topological relations, etc. Stored in a CpGridData object. Additionally, containers for parent-to-new-born entities are buil, as well as, new-born-to-parent. Maps(<int,bool>) to detect parent faces or cells are also provided. (Cell with 6 faces required).
[in] | cells_per_dim | Number of (refined) cells in each direction that each parent cell should be refined to. |
[in] | parent_idx | Parent cell index, cell to be refined. |
|
inline |
Set whether we want to have unique boundary ids.
uids | if true, each boundary intersection will have a unique boundary id. |
|
inline |
Is the grid currently using unique boundary ids?
void Dune::cpgrid::CpGridData::writeSintefLegacyFormat | ( | const std::string & | grid_prefix | ) | const |
Write the Sintef legacy grid format ('topogeom').
Read the Sintef legacy grid format ('topogeom').
grid_prefix | the grid name, such that topology will be found in <grid_prefix>-topo.dat etc. |