ModuloReg.h

Go to the documentation of this file.

// ModuloReg.h:
// -------------------
//

#ifndef __IPSDKUTIL_INSTRUCTIONSET_ARITHMETIC_DETAIL_COMMON_MODULOREG_H__
#define __IPSDKUTIL_INSTRUCTIONSET_ARITHMETIC_DETAIL_COMMON_MODULOREG_H__

#include <IPSDKUtil/IPSDKUtilExports.h>
#include <IPSDKUtil/InstructionSet/RegType.h>
#include <IPSDKUtil/InstructionSet/RegMaskType.h>
#include <IPSDKUtil/InstructionSet/Arithmetic/detail/DivReg.h>
#include <IPSDKUtil/InstructionSet/Arithmetic/detail/FloorReg.h>
#include <IPSDKUtil/InstructionSet/Arithmetic/detail/MaxReg.h>
#include <IPSDKUtil/InstructionSet/Arithmetic/detail/MinReg.h>
#include <IPSDKUtil/InstructionSet/Arithmetic/detail/MulReg.h>
#include <IPSDKUtil/InstructionSet/Arithmetic/detail/SubReg.h>
#include <IPSDKUtil/InstructionSet/Logical/detail/IfElseReg.h>
#include <IPSDKUtil/InstructionSet/Comparison/detail/IsGreaterEqualReg.h>
#include <IPSDKUtil/InstructionSet/detail/AssignReg.h>
#include <IPSDKUtil/Tools/ForceInline.h>

namespace ipsdk {
namespace simd {
namespace detail {


template <eInstructionSet::domain IS, typename T>
struct ModuloReg<
    IS,
    T,
    typename boost::enable_if<typename boost::is_integral<T>::type>::type
>
{
    static IPSDK_FORCEINLINE
    typename RegType<IS, T>::Type
    act(const typename RegType<IS, T>::Type& in1,
        const typename RegType<IS, T>::Type& in2)
    {
        typename RegType<IS, T>::Type out;
        act(in1, in2, out);
        return out;
    }

    static IPSDK_FORCEINLINE
    void
    act(const typename RegType<IS, T>::Type& in,
        const typename RegType<IS, T>::Type& m,
        typename RegType<IS, T>::Type& out)
    {
        typename RegType<IS, T>::Type cste0, cste1;
        AssignReg<IS, T>::act(cste0, 0);
        AssignReg<IS, T>::act(cste1, 1);

        // process integer division
        typename RegType<IS, T>::Type n, nM1;
        DivReg<IS, T>::act(in, m, n);
        SubReg<IS, T>::act(n, cste1, nM1);

        // compute two solutions for modulo operation :
        // -> in - n*mod if in >= 0
        // -> in - (n-1)*mod if in < 0
        typename RegMaskType<IS, T>::Type test;
        IsGreaterEqualReg<IS, T>::act(in, cste0, test);
        typename RegType<IS, T>::Type res;
        IfElseReg<IS, T>::act(test, n, nM1, res);
        MulReg<IS, T>::act(res, m, res);
        SubReg<IS, T>::act(in, res, out);
    }
};

template <eInstructionSet::domain IS, typename T>
struct ModuloReg<
    IS,
    T,
    typename boost::enable_if<
        typename boost::mpl::not_<
            typename boost::is_integral<T>::type
        >::type
    >::type
>
{
    static IPSDK_FORCEINLINE
    typename RegType<IS, T>::Type
    act(const typename RegType<IS, T>::Type& in1,
        const typename RegType<IS, T>::Type& in2)
    {
        typename RegType<IS, T>::Type out;
        act(in1, in2, out);
        return out;
    }

    static IPSDK_FORCEINLINE
    void
    act(const typename RegType<IS, T>::Type& in,
        const typename RegType<IS, T>::Type& m,
        typename RegType<IS, T>::Type& out)
    {
        // process integer division
        typename RegType<IS, T>::Type n;
        DivReg<IS, T>::act(in, m, n);
        FloorReg<IS, T>::act(n, n);

        typename RegType<IS, T>::Type res;
        MulReg<IS, T>::act(n, m, res);
        SubReg<IS, T>::act(in, res, res);

        // bound result to range [0; m[
        typename RegType<IS, T>::Type zero;
        AssignReg<IS, T>::act(zero, 0);
        MaxReg<IS, T>::act(res, zero, res);

        typename RegMaskType<IS, T>::Type test;
        IsGreaterEqualReg<IS, T>::act(res, m, test);
        IfElseReg<IS, T>::act(test, m, res, out);
    }
};


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

#endif // __IPSDKUTIL_INSTRUCTIONSET_ARITHMETIC_DETAIL_COMMON_MODULOREG_H__