cannam@154: /* Copyright (C) 2013 Xiph.Org Foundation and contributors */ cannam@154: /* cannam@154: Redistribution and use in source and binary forms, with or without cannam@154: modification, are permitted provided that the following conditions cannam@154: are met: cannam@154: cannam@154: - Redistributions of source code must retain the above copyright cannam@154: notice, this list of conditions and the following disclaimer. cannam@154: cannam@154: - Redistributions in binary form must reproduce the above copyright cannam@154: notice, this list of conditions and the following disclaimer in the cannam@154: documentation and/or other materials provided with the distribution. cannam@154: cannam@154: THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS cannam@154: ``AS IS'' AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT cannam@154: LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR cannam@154: A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT OWNER cannam@154: OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, cannam@154: EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, cannam@154: PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR cannam@154: PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF cannam@154: LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING cannam@154: NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS cannam@154: SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE. cannam@154: */ cannam@154: cannam@154: #ifndef FIXED_ARMv4_H cannam@154: #define FIXED_ARMv4_H cannam@154: cannam@154: /** 16x32 multiplication, followed by a 16-bit shift right. Results fits in 32 bits */ cannam@154: #undef MULT16_32_Q16 cannam@154: static OPUS_INLINE opus_val32 MULT16_32_Q16_armv4(opus_val16 a, opus_val32 b) cannam@154: { cannam@154: unsigned rd_lo; cannam@154: int rd_hi; cannam@154: __asm__( cannam@154: "#MULT16_32_Q16\n\t" cannam@154: "smull %0, %1, %2, %3\n\t" cannam@154: : "=&r"(rd_lo), "=&r"(rd_hi) cannam@154: : "%r"(b),"r"(SHL32(a,16)) cannam@154: ); cannam@154: return rd_hi; cannam@154: } cannam@154: #define MULT16_32_Q16(a, b) (MULT16_32_Q16_armv4(a, b)) cannam@154: cannam@154: cannam@154: /** 16x32 multiplication, followed by a 15-bit shift right. Results fits in 32 bits */ cannam@154: #undef MULT16_32_Q15 cannam@154: static OPUS_INLINE opus_val32 MULT16_32_Q15_armv4(opus_val16 a, opus_val32 b) cannam@154: { cannam@154: unsigned rd_lo; cannam@154: int rd_hi; cannam@154: __asm__( cannam@154: "#MULT16_32_Q15\n\t" cannam@154: "smull %0, %1, %2, %3\n\t" cannam@154: : "=&r"(rd_lo), "=&r"(rd_hi) cannam@154: : "%r"(b), "r"(SHL32(a,16)) cannam@154: ); cannam@154: /*We intentionally don't OR in the high bit of rd_lo for speed.*/ cannam@154: return SHL32(rd_hi,1); cannam@154: } cannam@154: #define MULT16_32_Q15(a, b) (MULT16_32_Q15_armv4(a, b)) cannam@154: cannam@154: cannam@154: /** 16x32 multiply, followed by a 15-bit shift right and 32-bit add. cannam@154: b must fit in 31 bits. cannam@154: Result fits in 32 bits. */ cannam@154: #undef MAC16_32_Q15 cannam@154: #define MAC16_32_Q15(c, a, b) ADD32(c, MULT16_32_Q15(a, b)) cannam@154: cannam@154: /** 16x32 multiply, followed by a 16-bit shift right and 32-bit add. cannam@154: Result fits in 32 bits. */ cannam@154: #undef MAC16_32_Q16 cannam@154: #define MAC16_32_Q16(c, a, b) ADD32(c, MULT16_32_Q16(a, b)) cannam@154: cannam@154: /** 32x32 multiplication, followed by a 31-bit shift right. Results fits in 32 bits */ cannam@154: #undef MULT32_32_Q31 cannam@154: #define MULT32_32_Q31(a,b) (opus_val32)((((opus_int64)(a)) * ((opus_int64)(b)))>>31) cannam@154: cannam@154: #endif