Use compressedIndexForInterior() for mapping cartesian indices.

This commit is contained in:
Atgeirr Flø Rasmussen 2022-09-05 12:42:54 +02:00
parent 69b32f7004
commit e9cba9e0c1
3 changed files with 1 additions and 7 deletions

View File

@ -316,9 +316,6 @@ namespace Opm {
// a vector of all the wells. // a vector of all the wells.
std::vector<WellInterfacePtr > well_container_{}; std::vector<WellInterfacePtr > well_container_{};
// map from logically cartesian cell indices to compressed ones
std::vector<int> cartesian_to_compressed_;
std::vector<bool> is_cell_perforated_{}; std::vector<bool> is_cell_perforated_{};

View File

@ -439,8 +439,6 @@ protected:
std::vector<std::reference_wrapper<ParallelWellInfo>> local_parallel_well_info_; std::vector<std::reference_wrapper<ParallelWellInfo>> local_parallel_well_info_;
std::vector<WellProdIndexCalculator> prod_index_calc_; std::vector<WellProdIndexCalculator> prod_index_calc_;
std::vector<int> cartesian_to_compressed_;
std::vector<int> pvt_region_idx_; std::vector<int> pvt_region_idx_;

View File

@ -128,8 +128,7 @@ namespace Opm {
for ( size_t c=0; c < connectionSet.size(); c++ ) for ( size_t c=0; c < connectionSet.size(); c++ )
{ {
const auto& connection = connectionSet.get(c); const auto& connection = connectionSet.get(c);
int compressed_idx = cartesian_to_compressed_.at(connection.global_index()); int compressed_idx = compressedIndexForInterior(connection.global_index());
if ( compressed_idx >= 0 ) { // Ignore connections in inactive/remote cells. if ( compressed_idx >= 0 ) { // Ignore connections in inactive/remote cells.
wellCells.push_back(compressed_idx); wellCells.push_back(compressed_idx);
} }