multiindex.hh

Go to the documentation of this file.

#ifndef DUNE_GRID_UTILITY_MULTIINDEX_HH
#define DUNE_GRID_UTILITY_MULTIINDEX_HH
 
#include<array>
 
namespace Dune
{
 namespace FactoryUtilities
 {
  template<std::size_t dim>
  class MultiIndex : public std::array<unsigned int,dim>
  {
    // The range of each component
    std::array<unsigned int,dim> limits_;
 
  public:
    MultiIndex(const std::array<unsigned int,dim>& limits) : limits_(limits)
    {
      std::fill(this->begin(), this->end(), 0);
    }
 
    MultiIndex<dim>& operator++()
    {
      for (std::size_t i=0; i<dim; i++)
      {
        // Augment digit
        (*this)[i]++;
 
        // If there is no carry-over we can stop here
        if ((*this)[i]<limits_[i])
          break;
 
        (*this)[i] = 0;
      }
      return *this;
    }
 
    size_t cycle() const
    {
      size_t result = 1;
      for (std::size_t i=0; i<dim; i++)
        result *= limits_[i];
      return result;
    }
  };
 }
}
 
#endif