seq/stencil.h

Go to the documentation of this file.

#ifndef GRPPI_SEQ_STENCIL_H
#define GRPPI_SEQ_STENCIL_H

#include "sequential_execution.h"

#include "../common/iterator.h"

namespace grppi {

template <typename InputIt, typename OutputIt, typename StencilTransformer, 
          typename Neighbourhood>
void stencil(sequential_execution & ex, 
             InputIt first, InputIt last, OutputIt out, 
             StencilTransformer transform_op, 
             Neighbourhood neighbour_op) 
{
  while (first!=last) {
    *out = transform_op(first, neighbour_op(first));
    first++;
    out++;
  }
}

template <typename InputIt, typename OutputIt, typename StencilTransformer, 
          typename Neighbourhood, typename ... OtherInputIts>
void stencil(sequential_execution & ex, 
             InputIt first, InputIt last, OutputIt out, 
             StencilTransformer transform_op, 
             Neighbourhood neighbour_op, OtherInputIts ... other_firsts) 
{
  while (first!=last) {
    *out = transform_op(first, neighbour_op(first, other_firsts ...));
    advance_iterators(other_firsts...);
    first++;
    out++;
  }
}

}

#endif