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