annotate src/opus-1.3/silk/macros.h @ 69:7aeed7906520

Add Opus sources and macOS builds
author Chris Cannam
date Wed, 23 Jan 2019 13:48:08 +0000
parents
children
rev   line source
Chris@69 1 /***********************************************************************
Chris@69 2 Copyright (c) 2006-2011, Skype Limited. All rights reserved.
Chris@69 3 Redistribution and use in source and binary forms, with or without
Chris@69 4 modification, are permitted provided that the following conditions
Chris@69 5 are met:
Chris@69 6 - Redistributions of source code must retain the above copyright notice,
Chris@69 7 this list of conditions and the following disclaimer.
Chris@69 8 - Redistributions in binary form must reproduce the above copyright
Chris@69 9 notice, this list of conditions and the following disclaimer in the
Chris@69 10 documentation and/or other materials provided with the distribution.
Chris@69 11 - Neither the name of Internet Society, IETF or IETF Trust, nor the
Chris@69 12 names of specific contributors, may be used to endorse or promote
Chris@69 13 products derived from this software without specific prior written
Chris@69 14 permission.
Chris@69 15 THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS"
Chris@69 16 AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
Chris@69 17 IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
Chris@69 18 ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT OWNER OR CONTRIBUTORS BE
Chris@69 19 LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR
Chris@69 20 CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF
Chris@69 21 SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS
Chris@69 22 INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN
Chris@69 23 CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE)
Chris@69 24 ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
Chris@69 25 POSSIBILITY OF SUCH DAMAGE.
Chris@69 26 ***********************************************************************/
Chris@69 27
Chris@69 28 #ifndef SILK_MACROS_H
Chris@69 29 #define SILK_MACROS_H
Chris@69 30
Chris@69 31 #ifdef HAVE_CONFIG_H
Chris@69 32 #include "config.h"
Chris@69 33 #endif
Chris@69 34
Chris@69 35 #include "opus_types.h"
Chris@69 36 #include "opus_defines.h"
Chris@69 37 #include "arch.h"
Chris@69 38
Chris@69 39 /* This is an OPUS_INLINE header file for general platform. */
Chris@69 40
Chris@69 41 /* (a32 * (opus_int32)((opus_int16)(b32))) >> 16 output have to be 32bit int */
Chris@69 42 #if OPUS_FAST_INT64
Chris@69 43 #define silk_SMULWB(a32, b32) ((opus_int32)(((a32) * (opus_int64)((opus_int16)(b32))) >> 16))
Chris@69 44 #else
Chris@69 45 #define silk_SMULWB(a32, b32) ((((a32) >> 16) * (opus_int32)((opus_int16)(b32))) + ((((a32) & 0x0000FFFF) * (opus_int32)((opus_int16)(b32))) >> 16))
Chris@69 46 #endif
Chris@69 47
Chris@69 48 /* a32 + (b32 * (opus_int32)((opus_int16)(c32))) >> 16 output have to be 32bit int */
Chris@69 49 #if OPUS_FAST_INT64
Chris@69 50 #define silk_SMLAWB(a32, b32, c32) ((opus_int32)((a32) + (((b32) * (opus_int64)((opus_int16)(c32))) >> 16)))
Chris@69 51 #else
Chris@69 52 #define silk_SMLAWB(a32, b32, c32) ((a32) + ((((b32) >> 16) * (opus_int32)((opus_int16)(c32))) + ((((b32) & 0x0000FFFF) * (opus_int32)((opus_int16)(c32))) >> 16)))
Chris@69 53 #endif
Chris@69 54
Chris@69 55 /* (a32 * (b32 >> 16)) >> 16 */
Chris@69 56 #if OPUS_FAST_INT64
Chris@69 57 #define silk_SMULWT(a32, b32) ((opus_int32)(((a32) * (opus_int64)((b32) >> 16)) >> 16))
Chris@69 58 #else
Chris@69 59 #define silk_SMULWT(a32, b32) (((a32) >> 16) * ((b32) >> 16) + ((((a32) & 0x0000FFFF) * ((b32) >> 16)) >> 16))
Chris@69 60 #endif
Chris@69 61
Chris@69 62 /* a32 + (b32 * (c32 >> 16)) >> 16 */
Chris@69 63 #if OPUS_FAST_INT64
Chris@69 64 #define silk_SMLAWT(a32, b32, c32) ((opus_int32)((a32) + (((b32) * ((opus_int64)(c32) >> 16)) >> 16)))
Chris@69 65 #else
Chris@69 66 #define silk_SMLAWT(a32, b32, c32) ((a32) + (((b32) >> 16) * ((c32) >> 16)) + ((((b32) & 0x0000FFFF) * ((c32) >> 16)) >> 16))
Chris@69 67 #endif
Chris@69 68
Chris@69 69 /* (opus_int32)((opus_int16)(a3))) * (opus_int32)((opus_int16)(b32)) output have to be 32bit int */
Chris@69 70 #define silk_SMULBB(a32, b32) ((opus_int32)((opus_int16)(a32)) * (opus_int32)((opus_int16)(b32)))
Chris@69 71
Chris@69 72 /* a32 + (opus_int32)((opus_int16)(b32)) * (opus_int32)((opus_int16)(c32)) output have to be 32bit int */
Chris@69 73 #define silk_SMLABB(a32, b32, c32) ((a32) + ((opus_int32)((opus_int16)(b32))) * (opus_int32)((opus_int16)(c32)))
Chris@69 74
Chris@69 75 /* (opus_int32)((opus_int16)(a32)) * (b32 >> 16) */
Chris@69 76 #define silk_SMULBT(a32, b32) ((opus_int32)((opus_int16)(a32)) * ((b32) >> 16))
Chris@69 77
Chris@69 78 /* a32 + (opus_int32)((opus_int16)(b32)) * (c32 >> 16) */
Chris@69 79 #define silk_SMLABT(a32, b32, c32) ((a32) + ((opus_int32)((opus_int16)(b32))) * ((c32) >> 16))
Chris@69 80
Chris@69 81 /* a64 + (b32 * c32) */
Chris@69 82 #define silk_SMLAL(a64, b32, c32) (silk_ADD64((a64), ((opus_int64)(b32) * (opus_int64)(c32))))
Chris@69 83
Chris@69 84 /* (a32 * b32) >> 16 */
Chris@69 85 #if OPUS_FAST_INT64
Chris@69 86 #define silk_SMULWW(a32, b32) ((opus_int32)(((opus_int64)(a32) * (b32)) >> 16))
Chris@69 87 #else
Chris@69 88 #define silk_SMULWW(a32, b32) silk_MLA(silk_SMULWB((a32), (b32)), (a32), silk_RSHIFT_ROUND((b32), 16))
Chris@69 89 #endif
Chris@69 90
Chris@69 91 /* a32 + ((b32 * c32) >> 16) */
Chris@69 92 #if OPUS_FAST_INT64
Chris@69 93 #define silk_SMLAWW(a32, b32, c32) ((opus_int32)((a32) + (((opus_int64)(b32) * (c32)) >> 16)))
Chris@69 94 #else
Chris@69 95 #define silk_SMLAWW(a32, b32, c32) silk_MLA(silk_SMLAWB((a32), (b32), (c32)), (b32), silk_RSHIFT_ROUND((c32), 16))
Chris@69 96 #endif
Chris@69 97
Chris@69 98 /* add/subtract with output saturated */
Chris@69 99 #define silk_ADD_SAT32(a, b) ((((opus_uint32)(a) + (opus_uint32)(b)) & 0x80000000) == 0 ? \
Chris@69 100 ((((a) & (b)) & 0x80000000) != 0 ? silk_int32_MIN : (a)+(b)) : \
Chris@69 101 ((((a) | (b)) & 0x80000000) == 0 ? silk_int32_MAX : (a)+(b)) )
Chris@69 102
Chris@69 103 #define silk_SUB_SAT32(a, b) ((((opus_uint32)(a)-(opus_uint32)(b)) & 0x80000000) == 0 ? \
Chris@69 104 (( (a) & ((b)^0x80000000) & 0x80000000) ? silk_int32_MIN : (a)-(b)) : \
Chris@69 105 ((((a)^0x80000000) & (b) & 0x80000000) ? silk_int32_MAX : (a)-(b)) )
Chris@69 106
Chris@69 107 #if defined(MIPSr1_ASM)
Chris@69 108 #include "mips/macros_mipsr1.h"
Chris@69 109 #endif
Chris@69 110
Chris@69 111 #include "ecintrin.h"
Chris@69 112 #ifndef OVERRIDE_silk_CLZ16
Chris@69 113 static OPUS_INLINE opus_int32 silk_CLZ16(opus_int16 in16)
Chris@69 114 {
Chris@69 115 return 32 - EC_ILOG(in16<<16|0x8000);
Chris@69 116 }
Chris@69 117 #endif
Chris@69 118
Chris@69 119 #ifndef OVERRIDE_silk_CLZ32
Chris@69 120 static OPUS_INLINE opus_int32 silk_CLZ32(opus_int32 in32)
Chris@69 121 {
Chris@69 122 return in32 ? 32 - EC_ILOG(in32) : 32;
Chris@69 123 }
Chris@69 124 #endif
Chris@69 125
Chris@69 126 /* Row based */
Chris@69 127 #define matrix_ptr(Matrix_base_adr, row, column, N) \
Chris@69 128 (*((Matrix_base_adr) + ((row)*(N)+(column))))
Chris@69 129 #define matrix_adr(Matrix_base_adr, row, column, N) \
Chris@69 130 ((Matrix_base_adr) + ((row)*(N)+(column)))
Chris@69 131
Chris@69 132 /* Column based */
Chris@69 133 #ifndef matrix_c_ptr
Chris@69 134 # define matrix_c_ptr(Matrix_base_adr, row, column, M) \
Chris@69 135 (*((Matrix_base_adr) + ((row)+(M)*(column))))
Chris@69 136 #endif
Chris@69 137
Chris@69 138 #ifdef OPUS_ARM_INLINE_ASM
Chris@69 139 #include "arm/macros_armv4.h"
Chris@69 140 #endif
Chris@69 141
Chris@69 142 #ifdef OPUS_ARM_INLINE_EDSP
Chris@69 143 #include "arm/macros_armv5e.h"
Chris@69 144 #endif
Chris@69 145
Chris@69 146 #ifdef OPUS_ARM_PRESUME_AARCH64_NEON_INTR
Chris@69 147 #include "arm/macros_arm64.h"
Chris@69 148 #endif
Chris@69 149
Chris@69 150 #endif /* SILK_MACROS_H */
Chris@69 151