/usr/include/dune/grid/onedgrid/onedgridleafiterator.hh is in libdune-grid-dev 2.5.1-1.
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 | // -*- tab-width: 4; indent-tabs-mode: nil; c-basic-offset: 2 -*-
// vi: set et ts=4 sw=2 sts=2:
#ifndef DUNE_ONE_D_GRID_LEAFITERATOR_HH
#define DUNE_ONE_D_GRID_LEAFITERATOR_HH
/** \file
* \brief The OneDGridLeafIterator class
*/
namespace Dune {
/** \brief Iterator over all entities of a given codimension and level of a grid.
* \ingroup OneDGrid
*/
template<int codim, PartitionIteratorType pitype, class GridImp>
class OneDGridLeafIterator
{
enum {dim = GridImp::dimension};
friend class OneDGridEntity<codim,dim,GridImp>;
public:
typedef typename GridImp::template Codim<codim>::Entity Entity;
enum {codimension = codim};
OneDGridLeafIterator(const GridImp& grid) : grid_(&grid) {
/** \todo Can a make the fullRefineLevel work somehow? */
const int fullRefineLevel = 0;
GridImp::getRealImplementation(this->virtualEntity_).setToTarget((OneDEntityImp<1-codim>*) std::get<1-codim>(grid_->entityImps_[fullRefineLevel]).begin());
if (!GridImp::getRealImplementation(this->virtualEntity_).target_->isLeaf())
increment();
}
//! Constructor
OneDGridLeafIterator()
{
GridImp::getRealImplementation(this->virtualEntity_).setToTarget(OneDGridNullIteratorFactory<1-codim>::null());
}
//! prefix increment
void increment() {
// Increment until you find a leaf entity
do {
globalIncrement();
} while (GridImp::getRealImplementation(this->virtualEntity_).target_
&& !GridImp::getRealImplementation(this->virtualEntity_).target_->isLeaf());
}
//! dereferencing
const Entity& dereference() const {return virtualEntity_;}
//! equality
bool equals(const OneDGridLeafIterator<codim,pitype,GridImp>& other) const {
return virtualEntity_ == other.virtualEntity_;
}
private:
/** \brief This increment makes the iterator wander over all entities on all levels */
void globalIncrement() {
// Backup current level because it may not be accessible anymore after
// setting the pointer to the next entity.
const int oldLevel = this->virtualEntity_.level();
// Increment on this level
GridImp::getRealImplementation(this->virtualEntity_).setToTarget(GridImp::getRealImplementation(this->virtualEntity_).target_->succ_);
// If beyond the end of this level set to first of next level
if (!GridImp::getRealImplementation(this->virtualEntity_).target_ && oldLevel < grid_->maxLevel()) {
GridImp::getRealImplementation(this->virtualEntity_).setToTarget(const_cast<OneDEntityImp<dim-codim>*>(std::get<1-codim>(grid_->entityImps_[oldLevel+1]).begin()));
}
}
// /////////////////////////////////////
// Data members
// /////////////////////////////////////
const GridImp* grid_;
//! The entity that the iterator is pointing to
Entity virtualEntity_;
};
} // namespace Dune
#endif
|