#ifndef __MultiArrayRotation_h__
#define __MultiArrayRotation_h__

#include <iostream>
using namespace std;

#include "boost/multi_array/types.hpp"
#include "boost/array.hpp"
#include "boost/multi_array/algorithm.hpp"
#include <algorithm>
#include <cstddef>
#include <functional>
#include <numeric>
#include <vector>

namespace boost {

  class cyclic_storage_order
  {
    typedef detail::multi_array::size_type size_type;
  public:
	  cyclic_storage_order( const std::size_t shift )
		  :shift_(shift)
	  {}
    template <std::size_t NumDims>
    operator general_storage_order<NumDims>() const {
      boost::array<size_type,NumDims> ordering;
      boost::array<bool,NumDims> ascending;

      for (size_type i=0; i != NumDims; ++i) {
        ordering[i] = (i+shift_)%NumDims;
        ascending[i] = true;
      }
      return general_storage_order<NumDims>(ordering.begin(),
                                            ascending.begin());
    }
  private:
	  std::size_t shift_;
  };
}

#endif	//	__MultiArrayRotation_h__