IpolDelimiter.hpp

//
// $Id: IpolDelimiter.hpp,v 1.3 2008/01/14 19:05:50 werner Exp $
//
#ifndef __IPOL_INTERPOLATION_TRAIT_HPP
#define __IPOL_INTERPOLATION_TRAIT_HPP "Created 11.06.2004 22:28:21 by bzfbenge"

/*
 Do NOT include math.h, because this is incompatible
 with the Intel's math header file that is used by VecAl.hpp!
#include <math.h>
*/

#include "Index.hpp"

namespace Fiber
{

template <class Type>
        class   NoDelimiter
{
public:

static  const Type&limit(const Type&T)
        {
                return T;
        }
};


#if 0

class   ConstantBoundary
{
public:

static  index_t left(index_t i, index_t size)
        {
                if (i<0 || i>=size) return 0;
                return i;
        }

static  index_t right(index_t i, index_t size)
        {
                if (i<0 || i>=size) return size;
                return i;
        }
};

class   PeriodicBoundary
{
public:
static  index_t left(index_t i, index_t size)
        {
                if (i>=size) return i % size;
                if (i<0)     return i % size + size;
                return i;
        }

static  index_t right(index_t i, index_t size)
        {
                return left(i,size);
        }
};

class   ReflectedBoundary
{
public:
static  index_t left(index_t i, index_t size)
        {
                if (i>=2*size) i %= 2*size;
                if (i>=size) return size-i;
                if (i<size ) i %= 2*size;
                if (i<0)     return size+i;
                return i;
        }

static  index_t right(index_t i, index_t size)
        {
                return left(i,size);
        }
};
#endif


} /* namespace Fiber */ 

#endif /* __IPOL_INTERPOLATION_TRAIT_HPP */