stencil.h

//
// Created by milinda on 3/31/17.
//

#ifndef SFCSORTBENCH_STENCIL_H
#define SFCSORTBENCH_STENCIL_H

#include <iostream>
#include <array>
#include "assert.h"

enum StencilDirection{STENCIL_DIR_X,STENCIL_DIR_Y,STENCIL_DIR_Z,STENCIL_DIR_XY,STENCIL_DIR_YZ,STENCIL_DIR_XZ,STENCIL_DIR_XYZ};

template  <typename T,unsigned int length, unsigned int offset>
struct Stencil
{

private:
    std::array<T,length> m_uiS;
    StencilDirection m_uiDirection;
    unsigned int m_uiOffset=offset;


public:
    Stencil(){};

    Stencil(T* coeff, int n, StencilDirection pDir)
    {
        assert(n==length);
        for(unsigned int k=0;k<length;k++)
            m_uiS[k]=coeff[k];

        m_uiDirection=pDir;
    }

    const T& operator[](unsigned int i) const
    {
        assert(i<length);
        return m_uiS[i];
    }

    inline StencilDirection getStencilDirection() const {return m_uiDirection;}
    inline unsigned int getOffset() const {return m_uiOffset;}
    inline unsigned int getStencilLength() const {return m_uiS.size();}


};


#endif //SFCSORTBENCH_STENCIL_H