InstructionSetArithmeticUtils.h

// InstructionSetUtils.h:
// ----------------------
//

#ifndef __IPSDKUTIL_INSTRUCTIONSETUTILS_H__
#define __IPSDKUTIL_INSTRUCTIONSETUTILS_H__

#include <IPSDKUtil/InstructionSet/Pack.h>
#include <IPSDKUtil/InstructionSet/Arithmetic/add.h>
#include <IPSDKUtil/InstructionSet/Arithmetic/sub.h>
#include <IPSDKUtil/InstructionSet/Arithmetic/mul.h>

namespace ipsdk {
namespace simd {


template <ipsdk::eInstructionSet::domain IS, typename T>
IPSDK_FORCEINLINE
void complexMultiply(const simd::Pack<IS, T>& packIn1Re,
                     const simd::Pack<IS, T>& packIn1Im,
                     const simd::Pack<IS, T>& packIn2Re,
                     const simd::Pack<IS, T>& packIn2Im,
                     simd::Pack<IS, T>& packOutRe,
                     simd::Pack<IS, T>& packOutIm)
{
    simd::Pack<IS, T> packOp1, packOp2;

    simd::mul<IS>(packIn1Re, packIn2Re, packOp1);
    simd::mul<IS>(packIn1Im, packIn2Im, packOp2);
    simd::sub<IS>(packOp1, packOp2, packOutRe);

    simd::mul<IS>(packIn1Re, packIn2Im, packOp1);
    simd::mul<IS>(packIn1Im, packIn2Re, packOp2);
    simd::add<IS>(packOp1, packOp2, packOutIm);
}


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

#endif // __IPSDKUTIL_INSTRUCTIONSETUTILS_H__