PolynomReg.h

Go to the documentation of this file.

// PolynomReg.h:
// -------------------
//

#ifndef __IPSDKUTIL_INSTRUCTIONSET_ARITHMETIC_DETAIL_COMMON_POLYNOMREG_H__
#define __IPSDKUTIL_INSTRUCTIONSET_ARITHMETIC_DETAIL_COMMON_POLYNOMREG_H__

#include <IPSDKUtil/IPSDKUtilExports.h>
#include <IPSDKUtil/InstructionSet/Arithmetic/detail/PolynomReg.h>
#include <IPSDKUtil/InstructionSet/Arithmetic/detail/FmaddReg.h>
#include <IPSDKUtil/InstructionSet/Arithmetic/detail/MulReg.h>
#include <IPSDKUtil/InstructionSet/detail/AssignReg.h>
#include <IPSDKUtil/InstructionSet/RegType.h>
#include <IPSDKUtil/Tools/ForceInline.h>

#include <boost/preprocessor/repetition/enum_params.hpp>

namespace ipsdk {
namespace simd {
    namespace detail {


#define DEFINE_POLYNOMREG(N)    \
template <eInstructionSet::domain IS, eInstructionSet::domain ISFma, typename T>    \
struct PolynomReg <IS, ISFma, N, T> \
{   \
    static IPSDK_FORCEINLINE    \
    typename RegType<IS, T>::Type   \
    act(const typename RegType<IS, T>::Type& in,    \
        BOOST_PP_ENUM_PARAMS(BOOST_PP_INC(N), const T& c))  \
    {   \
        typename RegType<IS, T>::Type out;  \
        act(in, BOOST_PP_ENUM_PARAMS(BOOST_PP_INC(N), c), out); \
        return out; \
    }   \
    \
    static IPSDK_FORCEINLINE    \
    void    \
    act(const typename RegType<IS, T>::Type& x, \
        BOOST_PP_ENUM_PARAMS(BOOST_PP_INC(N), const T& c),  \
        typename RegType<IS, T>::Type& out) \
    {   \
        typedef typename RegType<IS, T>::Type Reg;  \
        Reg tmp;    \
        PolynomReg<IS, ISFma, BOOST_PP_DEC(N), T>::act(x, BOOST_PP_ENUM_SHIFTED_PARAMS(BOOST_PP_INC(N), c), tmp);  \
        Reg r0; \
        AssignReg<IS, T>::act(r0, c0); \
        FmaddReg<ISFma, T>::act(x, tmp, r0, out);   \
    }   \
};

DEFINE_POLYNOMREG(20)
DEFINE_POLYNOMREG(19)
DEFINE_POLYNOMREG(18)
DEFINE_POLYNOMREG(17)
DEFINE_POLYNOMREG(16)
DEFINE_POLYNOMREG(15)
DEFINE_POLYNOMREG(14)
DEFINE_POLYNOMREG(13)
DEFINE_POLYNOMREG(12)
DEFINE_POLYNOMREG(11)
DEFINE_POLYNOMREG(10)
DEFINE_POLYNOMREG(9)
DEFINE_POLYNOMREG(8)
DEFINE_POLYNOMREG(7)
DEFINE_POLYNOMREG(6)
DEFINE_POLYNOMREG(5)
DEFINE_POLYNOMREG(4)
DEFINE_POLYNOMREG(3)
DEFINE_POLYNOMREG(2)

template <eInstructionSet::domain IS, eInstructionSet::domain ISFma, typename T>
struct PolynomReg <IS, ISFma, 1, T>
{
    static IPSDK_FORCEINLINE
    typename RegType<IS, T>::Type
    act(const typename RegType<IS, T>::Type& in,
        const T& c0,
        const T& c1)
    {
        typename RegType<IS, T>::Type out;
        act(in, c0, c1, out);
        return out;
    }

    static IPSDK_FORCEINLINE
    void
    act(const typename RegType<IS, T>::Type& x,
        const T& c0,
        const T& c1,
        typename RegType<IS, T>::Type& out)
    {
        typedef typename RegType<IS, T>::Type Reg;
        Reg r0, r1;
        AssignReg<IS, T>::act(r0, c0);
        AssignReg<IS, T>::act(r1, c1);
        FmaddReg<ISFma, T>::act(r1, x, r0, out);
    }
};


} // end of namespace detail
} // end of namespace simd
} // end of namespace ipsdk

#endif // __IPSDKUTIL_INSTRUCTIONSET_ARITHMETIC_DETAIL_SSE2_FLOORREG_H__