Convert tabs to spaces in the opus and celt code.
[opus.git] / celt / _kiss_fft_guts.h
1 /*Copyright (c) 2003-2004, Mark Borgerding
2
3   All rights reserved.
4
5   Redistribution and use in source and binary forms, with or without
6    modification, are permitted provided that the following conditions are met:
7
8     * Redistributions of source code must retain the above copyright notice,
9        this list of conditions and the following disclaimer.
10     * Redistributions in binary form must reproduce the above copyright notice,
11        this list of conditions and the following disclaimer in the
12        documentation and/or other materials provided with the distribution.
13
14   THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS"
15   AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
16   IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
17   ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT OWNER OR CONTRIBUTORS BE
18   LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR
19   CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF
20   SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS
21   INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN
22   CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE)
23   ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
24   POSSIBILITY OF SUCH DAMAGE.*/
25
26 #ifndef KISS_FFT_GUTS_H
27 #define KISS_FFT_GUTS_H
28
29 #define MIN(a,b) ((a)<(b) ? (a):(b))
30 #define MAX(a,b) ((a)>(b) ? (a):(b))
31
32 /* kiss_fft.h
33    defines kiss_fft_scalar as either short or a float type
34    and defines
35    typedef struct { kiss_fft_scalar r; kiss_fft_scalar i; }kiss_fft_cpx; */
36 #include "kiss_fft.h"
37
38 /*
39   Explanation of macros dealing with complex math:
40
41    C_MUL(m,a,b)         : m = a*b
42    C_FIXDIV( c , div )  : if a fixed point impl., c /= div. noop otherwise
43    C_SUB( res, a,b)     : res = a - b
44    C_SUBFROM( res , a)  : res -= a
45    C_ADDTO( res , a)    : res += a
46  * */
47 #ifdef FIXED_POINT
48 #include "arch.h"
49
50
51 # define SAMPPROD long long
52 #define SAMP_MAX 2147483647
53 #define TWID_MAX 32767
54 #define TRIG_UPSCALE 1
55
56 #define SAMP_MIN -SAMP_MAX
57
58
59 #   define S_MUL(a,b) MULT16_32_Q15(b, a)
60
61 #   define C_MUL(m,a,b) \
62       do{ (m).r = SUB32(S_MUL((a).r,(b).r) , S_MUL((a).i,(b).i)); \
63           (m).i = ADD32(S_MUL((a).r,(b).i) , S_MUL((a).i,(b).r)); }while(0)
64
65 #   define C_MULC(m,a,b) \
66       do{ (m).r = ADD32(S_MUL((a).r,(b).r) , S_MUL((a).i,(b).i)); \
67           (m).i = SUB32(S_MUL((a).i,(b).r) , S_MUL((a).r,(b).i)); }while(0)
68
69 #   define C_MUL4(m,a,b) \
70       do{ (m).r = SHR(SUB32(S_MUL((a).r,(b).r) , S_MUL((a).i,(b).i)),2); \
71           (m).i = SHR(ADD32(S_MUL((a).r,(b).i) , S_MUL((a).i,(b).r)),2); }while(0)
72
73 #   define C_MULBYSCALAR( c, s ) \
74       do{ (c).r =  S_MUL( (c).r , s ) ;\
75           (c).i =  S_MUL( (c).i , s ) ; }while(0)
76
77 #   define DIVSCALAR(x,k) \
78         (x) = S_MUL(  x, (TWID_MAX-((k)>>1))/(k)+1 )
79
80 #   define C_FIXDIV(c,div) \
81         do {    DIVSCALAR( (c).r , div);  \
82                 DIVSCALAR( (c).i  , div); }while (0)
83
84 #define  C_ADD( res, a,b)\
85     do {(res).r=ADD32((a).r,(b).r);  (res).i=ADD32((a).i,(b).i); \
86     }while(0)
87 #define  C_SUB( res, a,b)\
88     do {(res).r=SUB32((a).r,(b).r);  (res).i=SUB32((a).i,(b).i); \
89     }while(0)
90 #define C_ADDTO( res , a)\
91     do {(res).r = ADD32((res).r, (a).r);  (res).i = ADD32((res).i,(a).i);\
92     }while(0)
93
94 #define C_SUBFROM( res , a)\
95     do {(res).r = ADD32((res).r,(a).r);  (res).i = SUB32((res).i,(a).i); \
96     }while(0)
97
98 #else  /* not FIXED_POINT*/
99
100 #   define S_MUL(a,b) ( (a)*(b) )
101 #define C_MUL(m,a,b) \
102     do{ (m).r = (a).r*(b).r - (a).i*(b).i;\
103         (m).i = (a).r*(b).i + (a).i*(b).r; }while(0)
104 #define C_MULC(m,a,b) \
105     do{ (m).r = (a).r*(b).r + (a).i*(b).i;\
106         (m).i = (a).i*(b).r - (a).r*(b).i; }while(0)
107
108 #define C_MUL4(m,a,b) C_MUL(m,a,b)
109
110 #   define C_FIXDIV(c,div) /* NOOP */
111 #   define C_MULBYSCALAR( c, s ) \
112     do{ (c).r *= (s);\
113         (c).i *= (s); }while(0)
114 #endif
115
116 #ifndef CHECK_OVERFLOW_OP
117 #  define CHECK_OVERFLOW_OP(a,op,b) /* noop */
118 #endif
119
120 #ifndef C_ADD
121 #define  C_ADD( res, a,b)\
122     do { \
123             CHECK_OVERFLOW_OP((a).r,+,(b).r)\
124             CHECK_OVERFLOW_OP((a).i,+,(b).i)\
125             (res).r=(a).r+(b).r;  (res).i=(a).i+(b).i; \
126     }while(0)
127 #define  C_SUB( res, a,b)\
128     do { \
129             CHECK_OVERFLOW_OP((a).r,-,(b).r)\
130             CHECK_OVERFLOW_OP((a).i,-,(b).i)\
131             (res).r=(a).r-(b).r;  (res).i=(a).i-(b).i; \
132     }while(0)
133 #define C_ADDTO( res , a)\
134     do { \
135             CHECK_OVERFLOW_OP((res).r,+,(a).r)\
136             CHECK_OVERFLOW_OP((res).i,+,(a).i)\
137             (res).r += (a).r;  (res).i += (a).i;\
138     }while(0)
139
140 #define C_SUBFROM( res , a)\
141     do {\
142             CHECK_OVERFLOW_OP((res).r,-,(a).r)\
143             CHECK_OVERFLOW_OP((res).i,-,(a).i)\
144             (res).r -= (a).r;  (res).i -= (a).i; \
145     }while(0)
146 #endif /* C_ADD defined */
147
148 #ifdef FIXED_POINT
149 /*#  define KISS_FFT_COS(phase)  TRIG_UPSCALE*floor(MIN(32767,MAX(-32767,.5+32768 * cos (phase))))
150 #  define KISS_FFT_SIN(phase)  TRIG_UPSCALE*floor(MIN(32767,MAX(-32767,.5+32768 * sin (phase))))*/
151 #  define KISS_FFT_COS(phase)  floor(.5+TWID_MAX*cos (phase))
152 #  define KISS_FFT_SIN(phase)  floor(.5+TWID_MAX*sin (phase))
153 #  define HALF_OF(x) ((x)>>1)
154 #elif defined(USE_SIMD)
155 #  define KISS_FFT_COS(phase) _mm_set1_ps( cos(phase) )
156 #  define KISS_FFT_SIN(phase) _mm_set1_ps( sin(phase) )
157 #  define HALF_OF(x) ((x)*_mm_set1_ps(.5f))
158 #else
159 #  define KISS_FFT_COS(phase) (kiss_fft_scalar) cos(phase)
160 #  define KISS_FFT_SIN(phase) (kiss_fft_scalar) sin(phase)
161 #  define HALF_OF(x) ((x)*.5f)
162 #endif
163
164 #define  kf_cexp(x,phase) \
165         do{ \
166                 (x)->r = KISS_FFT_COS(phase);\
167                 (x)->i = KISS_FFT_SIN(phase);\
168         }while(0)
169
170 #define  kf_cexp2(x,phase) \
171    do{ \
172       (x)->r = TRIG_UPSCALE*celt_cos_norm((phase));\
173       (x)->i = TRIG_UPSCALE*celt_cos_norm((phase)-32768);\
174 }while(0)
175
176 #endif /* KISS_FFT_GUTS_H */