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 */