Loading...
Searching...
No Matches
elementborderlistfromgrid.hh
Go to the documentation of this file.
// -*- mode: C++; tab-width: 4; indent-tabs-mode: nil; c-basic-offset: 4 -*-
// vi: set et ts=4 sw=4 sts=4:
/*
This file is part of the Open Porous Media project (OPM).
OPM is free software: you can redistribute it and/or modify
it under the terms of the GNU General Public License as published by
the Free Software Foundation, either version 2 of the License, or
(at your option) any later version.
OPM is distributed in the hope that it will be useful,
but WITHOUT ANY WARRANTY; without even the implied warranty of
MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
GNU General Public License for more details.
You should have received a copy of the GNU General Public License
along with OPM. If not, see <http://www.gnu.org/licenses/>.
Consult the COPYING file in the top-level source directory of this
module for the precise wording of the license and the list of
copyright holders.
*/
#ifndef EWOMS_ELEMENT_BORDER_LIST_FROM_GRID_HH
#define EWOMS_ELEMENT_BORDER_LIST_FROM_GRID_HH
#include "overlaptypes.hh"
#include "blacklist.hh"
#include <dune/grid/common/datahandleif.hh>
#include <dune/grid/common/gridenums.hh>
#include <dune/grid/common/partitionset.hh>
#include <dune/common/version.hh>
#include <set>
namespace Opm {
namespace Linear {
template <class GridView, class ElementMapper>
class ElementBorderListFromGrid
{
using PeerBlackListedEntry = BlackList::PeerBlackListedEntry;
using PeerBlackList = BlackList::PeerBlackList;
using PeerBlackLists = BlackList::PeerBlackLists;
using Element = typename GridView::template Codim<0>::Entity;
class BorderListHandle_
: public Dune::CommDataHandleIF<BorderListHandle_, ProcessRank>
{
public:
BorderListHandle_(const GridView& gridView,
const ElementMapper& map,
BlackList& blackList,
BorderList& borderList)
: gridView_(gridView)
, map_(map)
, blackList_(blackList)
, borderList_(borderList)
{
for (const auto& elem : elements(gridView_)) {
if (elem.partitionType() != Dune::InteriorEntity) {
Index elemIdx = static_cast<Index>(map_.index(elem));
blackList_.addIndex(elemIdx);
}
}
}
// data handle methods
bool contains(int, int codim) const
{ return codim == 0; }
bool fixedSize(int, int) const
{ return true; }
template <class EntityType>
size_t size(const EntityType&) const
{ return 2; }
template <class MessageBufferImp, class EntityType>
void gather(MessageBufferImp& buff, const EntityType& e) const
{
unsigned myIdx = static_cast<unsigned>(map_.index(e));
buff.write(static_cast<unsigned>(gridView_.comm().rank()));
buff.write(myIdx);
}
template <class MessageBufferImp>
void scatter(MessageBufferImp& buff,
const Element& e,
size_t)
{
// discard the index if it is not on the process boundary
bool isInteriorNeighbor = false;
auto isIt = gridView_.ibegin(e);
const auto& isEndIt = gridView_.iend(e);
for (; isIt != isEndIt; ++isIt) {
if (!isIt->neighbor())
continue;
else if (isIt->outside().partitionType() == Dune::InteriorEntity) {
isInteriorNeighbor = true;
break;
}
}
if (!isInteriorNeighbor)
return;
BorderIndex bIdx;
bIdx.localIdx = static_cast<Index>(map_.index(e));
{
ProcessRank tmp;
buff.read(tmp);
bIdx.peerRank = tmp;
peerSet_.insert(tmp);
}
{
unsigned tmp;
buff.read(tmp);
bIdx.peerIdx = static_cast<Index>(tmp);
}
bIdx.borderDistance = 1;
borderList_.push_back(bIdx);
}
// this template method is needed because the above one only works for codim-0
// entities (i.e., elements) but the dune grid uses some code which causes the
// compiler to invoke the scatter method for every codim...
template <class MessageBufferImp, class EntityType>
void scatter(MessageBufferImp&,
const EntityType&,
size_t)
{ }
const std::set<ProcessRank>& peerSet() const
{ return peerSet_; }
private:
GridView gridView_;
const ElementMapper& map_;
std::set<ProcessRank> peerSet_;
BlackList& blackList_;
BorderList& borderList_;
};
class PeerBlackListHandle_
: public Dune::CommDataHandleIF<PeerBlackListHandle_, int>
{
public:
PeerBlackListHandle_(const GridView& gridView,
const ElementMapper& map,
PeerBlackLists& peerBlackLists)
: gridView_(gridView)
, map_(map)
, peerBlackLists_(peerBlackLists)
{}
// data handle methods
bool contains(int, int codim) const
{ return codim == 0; }
bool fixedSize(int, int) const
{ return true; }
template <class EntityType>
size_t size(const EntityType&) const
{ return 2; }
template <class MessageBufferImp, class EntityType>
void gather(MessageBufferImp& buff, const EntityType& e) const
{
buff.write(static_cast<int>(gridView_.comm().rank()));
buff.write(static_cast<int>(map_.index(e)));
}
template <class MessageBufferImp, class EntityType>
void scatter(MessageBufferImp& buff, const EntityType& e, size_t)
{
int peerRank;
int peerIdx;
Index localIdx;
buff.read(peerRank);
buff.read(peerIdx);
localIdx = static_cast<Index>(map_.index(e));
PeerBlackListedEntry pIdx;
pIdx.nativeIndexOfPeer = static_cast<Index>(peerIdx);
pIdx.myOwnNativeIndex = static_cast<Index>(localIdx);
peerBlackLists_[static_cast<ProcessRank>(peerRank)].push_back(pIdx);
}
const PeerBlackList& peerBlackList(ProcessRank peerRank) const
{ return peerBlackLists_.at(peerRank); }
private:
GridView gridView_;
const ElementMapper& map_;
PeerBlackLists peerBlackLists_;
};
public:
ElementBorderListFromGrid(const GridView& gridView, const ElementMapper& map)
: gridView_(gridView)
, map_(map)
{
BorderListHandle_ blh(gridView, map, blackList_, borderList_);
gridView.communicate(blh,
Dune::InteriorBorder_All_Interface,
Dune::BackwardCommunication);
PeerBlackListHandle_ pblh(gridView, map, peerBlackLists_);
gridView.communicate(pblh,
Dune::InteriorBorder_All_Interface,
Dune::BackwardCommunication);
auto peerIt = blh.peerSet().begin();
const auto& peerEndIt = blh.peerSet().end();
for (; peerIt != peerEndIt; ++peerIt) {
blackList_.setPeerList(*peerIt, pblh.peerBlackList(*peerIt));
}
}
// Access to the border list.
const BorderList& borderList() const
{ return borderList_; }
// Access to the black-list indices.
const BlackList& blackList() const
{ return blackList_; }
private:
const GridView gridView_;
const ElementMapper& map_;
BorderList borderList_;
BlackList blackList_;
PeerBlackLists peerBlackLists_;
};
} // namespace Linear
} // namespace Opm
#endif
blacklist.hh
Expresses which degrees of freedom are blacklisted for the parallel linear solvers and which domestic...
Opm::Linear::BlackList
Expresses which degrees of freedom are blacklisted for the parallel linear solvers and which domestic...
Definition blacklist.hh:49
Opm
This file contains a set of helper functions used by VFPProd / VFPInj.
Definition blackoilbioeffectsmodules.hh:43
overlaptypes.hh
This files provides several data structures for storing tuples of indices of remote and/or local proc...
Opm::Linear::ProcessRank
unsigned ProcessRank
The type of the rank of a process.
Definition overlaptypes.hh:49
Opm::Linear::Index
int Index
The type of an index of a degree of freedom.
Definition overlaptypes.hh:44
Opm::Linear::BorderList
std::list< BorderIndex > BorderList
This class managages a list of indices which are on the border of a process' partition of the grid.
Definition overlaptypes.hh:120
Opm::Linear::BlackList::PeerBlackListedEntry
Definition blacklist.hh:51
Opm::Linear::BorderIndex
A single index intersecting with the process boundary.
Definition overlaptypes.hh:102
Opm::Linear::BorderIndex::localIdx
Index localIdx
Index of the entity for the local process.
Definition overlaptypes.hh:104
Opm::Linear::BorderIndex::borderDistance
BorderDistance borderDistance
Distance to the process border for the peer (in hops).
Definition overlaptypes.hh:113
Opm::Linear::BorderIndex::peerRank
ProcessRank peerRank
Rank of the peer process.
Definition overlaptypes.hh:110
Opm::Linear::BorderIndex::peerIdx
Index peerIdx
Index of the entity for the peer process.
Definition overlaptypes.hh:107