Chris@69: /***********************************************************************
Chris@69: Copyright (c) 2006-2011, Skype Limited. All rights reserved.
Chris@69: Redistribution and use in source and binary forms, with or without
Chris@69: modification, are permitted provided that the following conditions
Chris@69: are met:
Chris@69: - Redistributions of source code must retain the above copyright notice,
Chris@69: this list of conditions and the following disclaimer.
Chris@69: - Redistributions in binary form must reproduce the above copyright
Chris@69: notice, this list of conditions and the following disclaimer in the
Chris@69: documentation and/or other materials provided with the distribution.
Chris@69: - Neither the name of Internet Society, IETF or IETF Trust, nor the
Chris@69: names of specific contributors, may be used to endorse or promote
Chris@69: products derived from this software without specific prior written
Chris@69: permission.
Chris@69: THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS"
Chris@69: AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
Chris@69: IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
Chris@69: ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT OWNER OR CONTRIBUTORS BE
Chris@69: LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR
Chris@69: CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF
Chris@69: SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS
Chris@69: INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN
Chris@69: CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE)
Chris@69: ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
Chris@69: POSSIBILITY OF SUCH DAMAGE.
Chris@69: ***********************************************************************/
Chris@69: 
Chris@69: #ifndef SILK_MACROS_H
Chris@69: #define SILK_MACROS_H
Chris@69: 
Chris@69: #ifdef HAVE_CONFIG_H
Chris@69: #include "config.h"
Chris@69: #endif
Chris@69: 
Chris@69: #include "opus_types.h"
Chris@69: #include "opus_defines.h"
Chris@69: #include "arch.h"
Chris@69: 
Chris@69: /* This is an OPUS_INLINE header file for general platform. */
Chris@69: 
Chris@69: /* (a32 * (opus_int32)((opus_int16)(b32))) >> 16 output have to be 32bit int */
Chris@69: #if OPUS_FAST_INT64
Chris@69: #define silk_SMULWB(a32, b32)            ((opus_int32)(((a32) * (opus_int64)((opus_int16)(b32))) >> 16))
Chris@69: #else
Chris@69: #define silk_SMULWB(a32, b32)            ((((a32) >> 16) * (opus_int32)((opus_int16)(b32))) + ((((a32) & 0x0000FFFF) * (opus_int32)((opus_int16)(b32))) >> 16))
Chris@69: #endif
Chris@69: 
Chris@69: /* a32 + (b32 * (opus_int32)((opus_int16)(c32))) >> 16 output have to be 32bit int */
Chris@69: #if OPUS_FAST_INT64
Chris@69: #define silk_SMLAWB(a32, b32, c32)       ((opus_int32)((a32) + (((b32) * (opus_int64)((opus_int16)(c32))) >> 16)))
Chris@69: #else
Chris@69: #define silk_SMLAWB(a32, b32, c32)       ((a32) + ((((b32) >> 16) * (opus_int32)((opus_int16)(c32))) + ((((b32) & 0x0000FFFF) * (opus_int32)((opus_int16)(c32))) >> 16)))
Chris@69: #endif
Chris@69: 
Chris@69: /* (a32 * (b32 >> 16)) >> 16 */
Chris@69: #if OPUS_FAST_INT64
Chris@69: #define silk_SMULWT(a32, b32)            ((opus_int32)(((a32) * (opus_int64)((b32) >> 16)) >> 16))
Chris@69: #else
Chris@69: #define silk_SMULWT(a32, b32)            (((a32) >> 16) * ((b32) >> 16) + ((((a32) & 0x0000FFFF) * ((b32) >> 16)) >> 16))
Chris@69: #endif
Chris@69: 
Chris@69: /* a32 + (b32 * (c32 >> 16)) >> 16 */
Chris@69: #if OPUS_FAST_INT64
Chris@69: #define silk_SMLAWT(a32, b32, c32)       ((opus_int32)((a32) + (((b32) * ((opus_int64)(c32) >> 16)) >> 16)))
Chris@69: #else
Chris@69: #define silk_SMLAWT(a32, b32, c32)       ((a32) + (((b32) >> 16) * ((c32) >> 16)) + ((((b32) & 0x0000FFFF) * ((c32) >> 16)) >> 16))
Chris@69: #endif
Chris@69: 
Chris@69: /* (opus_int32)((opus_int16)(a3))) * (opus_int32)((opus_int16)(b32)) output have to be 32bit int */
Chris@69: #define silk_SMULBB(a32, b32)            ((opus_int32)((opus_int16)(a32)) * (opus_int32)((opus_int16)(b32)))
Chris@69: 
Chris@69: /* a32 + (opus_int32)((opus_int16)(b32)) * (opus_int32)((opus_int16)(c32)) output have to be 32bit int */
Chris@69: #define silk_SMLABB(a32, b32, c32)       ((a32) + ((opus_int32)((opus_int16)(b32))) * (opus_int32)((opus_int16)(c32)))
Chris@69: 
Chris@69: /* (opus_int32)((opus_int16)(a32)) * (b32 >> 16) */
Chris@69: #define silk_SMULBT(a32, b32)            ((opus_int32)((opus_int16)(a32)) * ((b32) >> 16))
Chris@69: 
Chris@69: /* a32 + (opus_int32)((opus_int16)(b32)) * (c32 >> 16) */
Chris@69: #define silk_SMLABT(a32, b32, c32)       ((a32) + ((opus_int32)((opus_int16)(b32))) * ((c32) >> 16))
Chris@69: 
Chris@69: /* a64 + (b32 * c32) */
Chris@69: #define silk_SMLAL(a64, b32, c32)        (silk_ADD64((a64), ((opus_int64)(b32) * (opus_int64)(c32))))
Chris@69: 
Chris@69: /* (a32 * b32) >> 16 */
Chris@69: #if OPUS_FAST_INT64
Chris@69: #define silk_SMULWW(a32, b32)            ((opus_int32)(((opus_int64)(a32) * (b32)) >> 16))
Chris@69: #else
Chris@69: #define silk_SMULWW(a32, b32)            silk_MLA(silk_SMULWB((a32), (b32)), (a32), silk_RSHIFT_ROUND((b32), 16))
Chris@69: #endif
Chris@69: 
Chris@69: /* a32 + ((b32 * c32) >> 16) */
Chris@69: #if OPUS_FAST_INT64
Chris@69: #define silk_SMLAWW(a32, b32, c32)       ((opus_int32)((a32) + (((opus_int64)(b32) * (c32)) >> 16)))
Chris@69: #else
Chris@69: #define silk_SMLAWW(a32, b32, c32)       silk_MLA(silk_SMLAWB((a32), (b32), (c32)), (b32), silk_RSHIFT_ROUND((c32), 16))
Chris@69: #endif
Chris@69: 
Chris@69: /* add/subtract with output saturated */
Chris@69: #define silk_ADD_SAT32(a, b)             ((((opus_uint32)(a) + (opus_uint32)(b)) & 0x80000000) == 0 ?                              \
Chris@69:                                         ((((a) & (b)) & 0x80000000) != 0 ? silk_int32_MIN : (a)+(b)) :   \
Chris@69:                                         ((((a) | (b)) & 0x80000000) == 0 ? silk_int32_MAX : (a)+(b)) )
Chris@69: 
Chris@69: #define silk_SUB_SAT32(a, b)             ((((opus_uint32)(a)-(opus_uint32)(b)) & 0x80000000) == 0 ?                                        \
Chris@69:                                         (( (a) & ((b)^0x80000000) & 0x80000000) ? silk_int32_MIN : (a)-(b)) :    \
Chris@69:                                         ((((a)^0x80000000) & (b)  & 0x80000000) ? silk_int32_MAX : (a)-(b)) )
Chris@69: 
Chris@69: #if defined(MIPSr1_ASM)
Chris@69: #include "mips/macros_mipsr1.h"
Chris@69: #endif
Chris@69: 
Chris@69: #include "ecintrin.h"
Chris@69: #ifndef OVERRIDE_silk_CLZ16
Chris@69: static OPUS_INLINE opus_int32 silk_CLZ16(opus_int16 in16)
Chris@69: {
Chris@69:     return 32 - EC_ILOG(in16<<16|0x8000);
Chris@69: }
Chris@69: #endif
Chris@69: 
Chris@69: #ifndef OVERRIDE_silk_CLZ32
Chris@69: static OPUS_INLINE opus_int32 silk_CLZ32(opus_int32 in32)
Chris@69: {
Chris@69:     return in32 ? 32 - EC_ILOG(in32) : 32;
Chris@69: }
Chris@69: #endif
Chris@69: 
Chris@69: /* Row based */
Chris@69: #define matrix_ptr(Matrix_base_adr, row, column, N) \
Chris@69:     (*((Matrix_base_adr) + ((row)*(N)+(column))))
Chris@69: #define matrix_adr(Matrix_base_adr, row, column, N) \
Chris@69:       ((Matrix_base_adr) + ((row)*(N)+(column)))
Chris@69: 
Chris@69: /* Column based */
Chris@69: #ifndef matrix_c_ptr
Chris@69: #   define matrix_c_ptr(Matrix_base_adr, row, column, M) \
Chris@69:     (*((Matrix_base_adr) + ((row)+(M)*(column))))
Chris@69: #endif
Chris@69: 
Chris@69: #ifdef OPUS_ARM_INLINE_ASM
Chris@69: #include "arm/macros_armv4.h"
Chris@69: #endif
Chris@69: 
Chris@69: #ifdef OPUS_ARM_INLINE_EDSP
Chris@69: #include "arm/macros_armv5e.h"
Chris@69: #endif
Chris@69: 
Chris@69: #ifdef OPUS_ARM_PRESUME_AARCH64_NEON_INTR
Chris@69: #include "arm/macros_arm64.h"
Chris@69: #endif
Chris@69: 
Chris@69: #endif /* SILK_MACROS_H */
Chris@69: