permute.h

#ifndef _MOBY_PERMUTE_
#define _MOBY_PERMUTE_

#include <iterator>

namespace Moby {

template <class ForwardIterator, class RandomAccessIterator1, class RandomAccessIterator2>
void permute(ForwardIterator mapping_begin, ForwardIterator mapping_end, RandomAccessIterator1 input_begin, RandomAccessIterator2 output_begin)
{
  // setup the target index
  unsigned target_index = 0;

  for (ForwardIterator i = mapping_begin; i != mapping_end; i++)
    output_begin[target_index++] = input_begin[*i]; 
}

} // end namespace

#endif