/usr/include/dune/grid/uggrid/uggridleveliterator.hh is in libdune-grid-dev 2.2.1-2.
This file is owned by root:root, with mode 0o644.
The actual contents of the file can be viewed below.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 | #ifndef DUNE_UGGRIDLEVELITERATOR_HH
#define DUNE_UGGRIDLEVELITERATOR_HH
#include <dune/grid/uggrid/uggridentitypointer.hh>
/** \file
* \brief The UGGridLevelIterator class
*/
namespace Dune {
//**********************************************************************
//
// --UGGridLevelIterator
// --LevelIterator
/** \brief Iterator over all entities of a given codimension and level of a grid.
* \ingroup UGGrid
*/
template<int codim, PartitionIteratorType pitype, class GridImp>
class UGGridLevelIterator :
public Dune::UGGridEntityPointer <codim,GridImp>
{
enum {dim = GridImp::dimension};
friend class UGGridEntity<codim,GridImp::dimension,GridImp>;
friend class UGGridEntity<0, GridImp::dimension,GridImp>;
// The type of the UG entity we're pointing to
typedef typename UG_NS<dim>::template Entity<codim>::T UGEntity;
public:
typedef typename GridImp::template Codim<codim>::Entity Entity;
//! Constructor
explicit UGGridLevelIterator()
{
this->virtualEntity_.setToTarget(nullptr,nullptr);
}
//! Constructor
explicit UGGridLevelIterator(const GridImp& gridImp, int level) : gridImp_(&gridImp)
{
typename UG_NS<dim>::Grid *theGrid = const_cast<typename UG_NS<dim>::Grid* >(gridImp_->multigrid_->grids[level]);
assert(theGrid);
if (codim==dim) {
if (pitype==All_Partition || pitype==Ghost_Partition)
this->virtualEntity_.setToTarget((UGEntity*)UG_NS<dim>::PFirstNode(theGrid),gridImp_);
else if (pitype == Dune::Interior_Partition || pitype == Dune::InteriorBorder_Partition)
this->virtualEntity_.setToTarget((UGEntity*)UG_NS<dim>::FirstNode(theGrid),gridImp_);
else // overlap and overlap-front
this->virtualEntity_.setToTarget(0,nullptr);
}
else if (codim==0) {
if (pitype==All_Partition || pitype==Ghost_Partition)
this->virtualEntity_.setToTarget((UGEntity*)UG_NS<dim>::PFirstElement(theGrid),gridImp_);
else if (pitype == Dune::Interior_Partition || pitype == Dune::InteriorBorder_Partition)
this->virtualEntity_.setToTarget((UGEntity*)UG_NS<dim>::FirstElement(theGrid),gridImp_);
else // overlap and overlap-front
this->virtualEntity_.setToTarget(0,nullptr);
}
else
DUNE_THROW(NotImplemented, "UGGrid level iterators for codimension " << codim);
if (this->virtualEntity_.getTarget() && !entityOK_())
increment();
}
//! prefix increment
void increment()
{
assert(this->level() == UG_NS<dim>::myLevel(this->virtualEntity_.getTarget()));
// Increment
do {
this->virtualEntity_.setToTarget(UG_NS<dim>::succ(this->virtualEntity_.getTarget()),gridImp_);
}
while (this->virtualEntity_.getTarget() && !entityOK_());
}
private:
/**
* \brief Return true iff the current entity is within the right
* partition.
*/
bool entityOK_()
{
if (pitype == All_Partition)
return true;
Dune::PartitionType entityPIType = this->virtualEntity_.partitionType();
if (pitype == Ghost_Partition && entityPIType == GhostEntity)
return true;
else if (pitype == Interior_Partition && entityPIType == InteriorEntity)
return true;
else if (pitype == InteriorBorder_Partition &&
(entityPIType == BorderEntity ||
entityPIType == InteriorEntity))
return true;
return false;
}
// /////////////////////////////////////
// Data members
// /////////////////////////////////////
const GridImp* gridImp_;
};
} // namespace Dune
#endif
|