SubsetIterator.h

// Copyright (C) 2010 Anders Logg
//
// This file is part of DOLFIN.
//
// DOLFIN is free software: you can redistribute it and/or modify
// it under the terms of the GNU Lesser General Public License as published by
// the Free Software Foundation, either version 3 of the License, or
// (at your option) any later version.
//
// DOLFIN 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 Lesser General Public License for more details.
//
// You should have received a copy of the GNU Lesser General Public License
// along with DOLFIN. If not, see <http://www.gnu.org/licenses/>.
//
// Modified by Andre Massing 2011.
//
// First added:  2010-11-17
// Last changed: 2011-11-10

#ifndef __SUBSET_ITERATOR_H
#define __SUBSET_ITERATOR_H

#include <memory>
#include <vector>

#include "Mesh.h"
#include "MeshEntity.h"
#include "MeshEntityIterator.h"
#include "MeshFunction.h"

namespace dolfin
{


  class SubsetIterator
  {
  public:

    SubsetIterator(const MeshFunction<std::size_t>& labels, std::size_t label)
      : _entity(*labels.mesh(), labels.dim(), 0),
        _subset(new std::vector<std::size_t>()), subset(*_subset)
    {
      // Extract subset
      subset.clear();
      for (MeshEntityIterator entity(*labels.mesh(), labels.dim());
           !entity.end(); ++entity)
      {
        if (labels[*entity] == label)
          subset.push_back(entity->index());
      }
      info("Iterating over subset, found %d entities out of %d.",
           subset.size(), labels.size());

      // Set iterator
      it = subset.begin();
    }

  SubsetIterator(const SubsetIterator& subset_iter)
    : _entity(subset_iter._entity), _subset(subset_iter._subset),
      subset(*_subset), it(subset_iter.it) {}

    virtual ~SubsetIterator() {}

    SubsetIterator& operator++()
    {
      ++it;
      return *this;
    }

    SubsetIterator& operator--()
    {
      --it;
      return *this;
    }

    bool operator==(const SubsetIterator& sub_iter) const
    {
      return ((const_cast<SubsetIterator *>(this))->operator*()
              == (const_cast<SubsetIterator *>(&sub_iter))->operator*()
              && it == sub_iter.it && &subset == &sub_iter.subset);
    }

    bool operator!=(const SubsetIterator & sub_iter) const
    { return !operator==(sub_iter); }

    MeshEntity& operator*()
    { return *operator->(); }

    MeshEntity* operator->()
    { _entity._local_index = *it; return &_entity; }

    bool end() const
    { return it == subset.end(); }

    SubsetIterator end_iterator()
    {
      SubsetIterator sg(*this);
      sg.set_end();
      return sg;
    }

  private:

    // Set pos to end position. To create a kind of mesh.end()
    // iterator.
    void set_end()
    { it = subset.end(); }

    // Mesh entity
    MeshEntity _entity;

    // Subset in shared data form
    std::shared_ptr<std::vector<std::size_t>> _subset;

    //Subset reference for convenience / speed
    std::vector<std::size_t>& subset;

    // Iterator
    std::vector<std::size_t>::iterator it;

  };

}

#endif