|
Chris@69
|
1 /* Copyright (c) 2003-2008 Timothy B. Terriberry
|
|
Chris@69
|
2 Copyright (c) 2008 Xiph.Org Foundation */
|
|
Chris@69
|
3 /*
|
|
Chris@69
|
4 Redistribution and use in source and binary forms, with or without
|
|
Chris@69
|
5 modification, are permitted provided that the following conditions
|
|
Chris@69
|
6 are met:
|
|
Chris@69
|
7
|
|
Chris@69
|
8 - Redistributions of source code must retain the above copyright
|
|
Chris@69
|
9 notice, this list of conditions and the following disclaimer.
|
|
Chris@69
|
10
|
|
Chris@69
|
11 - Redistributions in binary form must reproduce the above copyright
|
|
Chris@69
|
12 notice, this list of conditions and the following disclaimer in the
|
|
Chris@69
|
13 documentation and/or other materials provided with the distribution.
|
|
Chris@69
|
14
|
|
Chris@69
|
15 THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS
|
|
Chris@69
|
16 ``AS IS'' AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT
|
|
Chris@69
|
17 LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR
|
|
Chris@69
|
18 A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT OWNER
|
|
Chris@69
|
19 OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL,
|
|
Chris@69
|
20 EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO,
|
|
Chris@69
|
21 PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR
|
|
Chris@69
|
22 PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF
|
|
Chris@69
|
23 LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING
|
|
Chris@69
|
24 NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS
|
|
Chris@69
|
25 SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
|
|
Chris@69
|
26 */
|
|
Chris@69
|
27
|
|
Chris@69
|
28 /*Some common macros for potential platform-specific optimization.*/
|
|
Chris@69
|
29 #include "opus_types.h"
|
|
Chris@69
|
30 #include <math.h>
|
|
Chris@69
|
31 #include <limits.h>
|
|
Chris@69
|
32 #include "arch.h"
|
|
Chris@69
|
33 #if !defined(_ecintrin_H)
|
|
Chris@69
|
34 # define _ecintrin_H (1)
|
|
Chris@69
|
35
|
|
Chris@69
|
36 /*Some specific platforms may have optimized intrinsic or OPUS_INLINE assembly
|
|
Chris@69
|
37 versions of these functions which can substantially improve performance.
|
|
Chris@69
|
38 We define macros for them to allow easy incorporation of these non-ANSI
|
|
Chris@69
|
39 features.*/
|
|
Chris@69
|
40
|
|
Chris@69
|
41 /*Modern gcc (4.x) can compile the naive versions of min and max with cmov if
|
|
Chris@69
|
42 given an appropriate architecture, but the branchless bit-twiddling versions
|
|
Chris@69
|
43 are just as fast, and do not require any special target architecture.
|
|
Chris@69
|
44 Earlier gcc versions (3.x) compiled both code to the same assembly
|
|
Chris@69
|
45 instructions, because of the way they represented ((_b)>(_a)) internally.*/
|
|
Chris@69
|
46 # define EC_MINI(_a,_b) ((_a)+(((_b)-(_a))&-((_b)<(_a))))
|
|
Chris@69
|
47
|
|
Chris@69
|
48 /*Count leading zeros.
|
|
Chris@69
|
49 This macro should only be used for implementing ec_ilog(), if it is defined.
|
|
Chris@69
|
50 All other code should use EC_ILOG() instead.*/
|
|
Chris@69
|
51 #if defined(_MSC_VER) && (_MSC_VER >= 1400)
|
|
Chris@69
|
52 # include <intrin.h>
|
|
Chris@69
|
53 /*In _DEBUG mode this is not an intrinsic by default.*/
|
|
Chris@69
|
54 # pragma intrinsic(_BitScanReverse)
|
|
Chris@69
|
55
|
|
Chris@69
|
56 static __inline int ec_bsr(unsigned long _x){
|
|
Chris@69
|
57 unsigned long ret;
|
|
Chris@69
|
58 _BitScanReverse(&ret,_x);
|
|
Chris@69
|
59 return (int)ret;
|
|
Chris@69
|
60 }
|
|
Chris@69
|
61 # define EC_CLZ0 (1)
|
|
Chris@69
|
62 # define EC_CLZ(_x) (-ec_bsr(_x))
|
|
Chris@69
|
63 #elif defined(ENABLE_TI_DSPLIB)
|
|
Chris@69
|
64 # include "dsplib.h"
|
|
Chris@69
|
65 # define EC_CLZ0 (31)
|
|
Chris@69
|
66 # define EC_CLZ(_x) (_lnorm(_x))
|
|
Chris@69
|
67 #elif __GNUC_PREREQ(3,4)
|
|
Chris@69
|
68 # if INT_MAX>=2147483647
|
|
Chris@69
|
69 # define EC_CLZ0 ((int)sizeof(unsigned)*CHAR_BIT)
|
|
Chris@69
|
70 # define EC_CLZ(_x) (__builtin_clz(_x))
|
|
Chris@69
|
71 # elif LONG_MAX>=2147483647L
|
|
Chris@69
|
72 # define EC_CLZ0 ((int)sizeof(unsigned long)*CHAR_BIT)
|
|
Chris@69
|
73 # define EC_CLZ(_x) (__builtin_clzl(_x))
|
|
Chris@69
|
74 # endif
|
|
Chris@69
|
75 #endif
|
|
Chris@69
|
76
|
|
Chris@69
|
77 #if defined(EC_CLZ)
|
|
Chris@69
|
78 /*Note that __builtin_clz is not defined when _x==0, according to the gcc
|
|
Chris@69
|
79 documentation (and that of the BSR instruction that implements it on x86).
|
|
Chris@69
|
80 The majority of the time we can never pass it zero.
|
|
Chris@69
|
81 When we need to, it can be special cased.*/
|
|
Chris@69
|
82 # define EC_ILOG(_x) (EC_CLZ0-EC_CLZ(_x))
|
|
Chris@69
|
83 #else
|
|
Chris@69
|
84 int ec_ilog(opus_uint32 _v);
|
|
Chris@69
|
85 # define EC_ILOG(_x) (ec_ilog(_x))
|
|
Chris@69
|
86 #endif
|
|
Chris@69
|
87 #endif
|
