GridPathSegmenter.h

#ifndef IPPL_GRIDPATH_SEGMENTER_H
#define IPPL_GRIDPATH_SEGMENTER_H
#include <array>
 
namespace ippl {
 
template<unsigned Dim, typename T>
struct Segment {
  Vector<T,Dim> p0, p1;
};
 
struct DefaultCellCrossingRule {};
 
template<unsigned Dim, typename T, typename Rule = DefaultCellCrossingRule>
struct GridPathSegmenter {
 
  static KOKKOS_INLINE_FUNCTION
  Kokkos::Array<Segment<Dim,T>, Dim+1>
  split(const Vector<T, Dim>& A,
        const Vector<T, Dim>& B,
        const Vector<T, Dim>& origin,
        const Vector<T, Dim>& h);
};
 
} // namespace ippl 
 
#include "GridPathSegmenter.hpp"
 
#endif