intersectioniterator.hh

// -*- tab-width: 4; indent-tabs-mode: nil; c-basic-offset: 2 -*-
// vi: set et ts=4 sw=2 sts=2:
#ifndef DUNE_GEOGRID_INTERSECTIONITERATOR_HH
#define DUNE_GEOGRID_INTERSECTIONITERATOR_HH
 
#include <dune/grid/geometrygrid/entitypointer.hh>
#include <dune/grid/geometrygrid/intersection.hh>
 
namespace Dune
{
 
  namespace GeoGrid
  {
 
    // IntersectionIterator
    // --------------------
 
    template< class Grid, class HostIntersectionIterator >
    class IntersectionIterator
    {
      typedef typename remove_const< Grid >::type::Traits Traits;
 
      typedef GeoGrid::Intersection< Grid, typename HostIntersectionIterator::Intersection > IntersectionImpl;
 
      typedef typename Traits::template Codim< 0 >::EntityPointerImpl EntityPointerImpl;
      typedef typename Traits::template Codim< 0 >::Geometry ElementGeometry;
 
    public:
      typedef Dune::Intersection< Grid, IntersectionImpl > Intersection;
 
      typedef typename Traits::template Codim< 0 >::EntityPointer EntityPointer;
 
      template< class Entity >
      IntersectionIterator ( const Entity &inside,
                             const HostIntersectionIterator &hostIterator )
        : hostIterator_( hostIterator ),
          intersection_( IntersectionImpl( inside.geometry() ) )
      {}
 
      IntersectionIterator ( const IntersectionIterator &other )
        : hostIterator_( other.hostIterator_ ),
          intersection_( IntersectionImpl( Grid::getRealImplementation( other.intersection_ ) ) )
      {}
 
      IntersectionIterator &operator= ( const IntersectionIterator &other )
      {
        hostIterator_ = other.hostIterator_;
        Grid::getRealImplementation( intersection_ ) = Grid::getRealImplementation( other.intersection_ );
        return *this;
      }
 
      bool equals ( const IntersectionIterator &other ) const
      {
        return (hostIterator_ == other.hostIterator_);
      }
 
      void increment ()
      {
        ++hostIterator_;
        intersectionImpl().invalidate();
      }
 
      const Intersection &dereference () const
      {
        if( !intersectionImpl() )
          intersectionImpl().initialize( *hostIterator_ );
        return intersection_;
      }
 
    private:
      IntersectionImpl &intersectionImpl () const
      {
        return Grid::getRealImplementation( intersection_ );
      }
 
      HostIntersectionIterator hostIterator_;
      mutable Intersection intersection_;
    };
 
  } // namespace GeoGrid
 
} // namespace Dune
 
#endif // #ifndef DUNE_GEOGRID_INTERSECTIONITERATOR_HH