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

Merge pull request #1339 from blattms/write-ranks-without-nonactive

Browse Source 
Disregard non-active cells when writing MPI_RANK in writeInitial.
This commit is contained in:
Atgeirr Flø Rasmussen
2017-11-27 07:22:34 +01:00
committed by GitHub
parent 401aa28c08 f32743cf31
commit f62275cb06
2 changed files with 5 additions and 9 deletions
Download Patch File Download Diff File Expand all files Collapse all files

6
opm/autodiff/FlowMainEbos.hpp
View File

@@ -813,11 +813,9 @@ namespace Opm
using ElementMapper = Dune::MultipleCodimMultipleGeomTypeMapper<GridView, Dune::MCMGElementLayout>;
using Handle = CellOwnerDataHandle<ElementMapper>;
ElementMapper globalMapper(this->globalGrid().leafGridView());
const auto* dims = UgGridHelpers::cartDims(grid());
const auto globalSize = dims[0]*dims[1]*dims[2];
const auto globalSize = this->globalGrid().size(0);
std::vector<int> ranks(globalSize, -1);
Handle handle(globalMapper, ranks,
this->globalGrid().globalCell());
Handle handle(globalMapper, ranks);
this->grid().gatherData(handle);
integerVectors.emplace("MPI_RANK", ranks);
}

8
opm/autodiff/RedistributeDataHandles.hpp
View File

@@ -126,9 +126,8 @@ class CellOwnerDataHandle
public:
using DataType = int;
CellOwnerDataHandle(const Mapper& globalMapper, std::vector<int>& globalData,
const std::vector<int>& globalCell)
: globalMapper_(globalMapper), globalData_(globalData), globalCell_(globalCell)
CellOwnerDataHandle(const Mapper& globalMapper, std::vector<int>& globalData)
: globalMapper_(globalMapper), globalData_(globalData)
{
int argc = 0;
char** argv = nullptr;
@@ -158,7 +157,7 @@ public:
template<class B, class T>
void scatter(B& buffer, const T& e, std::size_t /* size */)
{
const auto& index = globalCell_[globalMapper_.index(e)];
const auto& index = globalMapper_.index(e);
buffer.read(globalData_[index]);
}
bool contains(int dim, int codim)
@@ -170,7 +169,6 @@ private:
int my_rank_;
const Mapper& globalMapper_;
std::vector<int>& globalData_;
const std::vector<int>& globalCell_;
};
#if HAVE_OPM_GRID && HAVE_MPI