/src/aac/libFDK/include/cplx_mul.h
Line  | Count  | Source  | 
1  |  | /* -----------------------------------------------------------------------------  | 
2  |  | Software License for The Fraunhofer FDK AAC Codec Library for Android  | 
3  |  |  | 
4  |  | © Copyright  1995 - 2018 Fraunhofer-Gesellschaft zur Förderung der angewandten  | 
5  |  | Forschung e.V. All rights reserved.  | 
6  |  |  | 
7  |  |  1.    INTRODUCTION  | 
8  |  | The Fraunhofer FDK AAC Codec Library for Android ("FDK AAC Codec") is software | 
9  |  | that implements the MPEG Advanced Audio Coding ("AAC") encoding and decoding | 
10  |  | scheme for digital audio. This FDK AAC Codec software is intended to be used on  | 
11  |  | a wide variety of Android devices.  | 
12  |  |  | 
13  |  | AAC's HE-AAC and HE-AAC v2 versions are regarded as today's most efficient  | 
14  |  | general perceptual audio codecs. AAC-ELD is considered the best-performing  | 
15  |  | full-bandwidth communications codec by independent studies and is widely  | 
16  |  | deployed. AAC has been standardized by ISO and IEC as part of the MPEG  | 
17  |  | specifications.  | 
18  |  |  | 
19  |  | Patent licenses for necessary patent claims for the FDK AAC Codec (including  | 
20  |  | those of Fraunhofer) may be obtained through Via Licensing  | 
21  |  | (www.vialicensing.com) or through the respective patent owners individually for  | 
22  |  | the purpose of encoding or decoding bit streams in products that are compliant  | 
23  |  | with the ISO/IEC MPEG audio standards. Please note that most manufacturers of  | 
24  |  | Android devices already license these patent claims through Via Licensing or  | 
25  |  | directly from the patent owners, and therefore FDK AAC Codec software may  | 
26  |  | already be covered under those patent licenses when it is used for those  | 
27  |  | licensed purposes only.  | 
28  |  |  | 
29  |  | Commercially-licensed AAC software libraries, including floating-point versions  | 
30  |  | with enhanced sound quality, are also available from Fraunhofer. Users are  | 
31  |  | encouraged to check the Fraunhofer website for additional applications  | 
32  |  | information and documentation.  | 
33  |  |  | 
34  |  | 2.    COPYRIGHT LICENSE  | 
35  |  |  | 
36  |  | Redistribution and use in source and binary forms, with or without modification,  | 
37  |  | are permitted without payment of copyright license fees provided that you  | 
38  |  | satisfy the following conditions:  | 
39  |  |  | 
40  |  | You must retain the complete text of this software license in redistributions of  | 
41  |  | the FDK AAC Codec or your modifications thereto in source code form.  | 
42  |  |  | 
43  |  | You must retain the complete text of this software license in the documentation  | 
44  |  | and/or other materials provided with redistributions of the FDK AAC Codec or  | 
45  |  | your modifications thereto in binary form. You must make available free of  | 
46  |  | charge copies of the complete source code of the FDK AAC Codec and your  | 
47  |  | modifications thereto to recipients of copies in binary form.  | 
48  |  |  | 
49  |  | The name of Fraunhofer may not be used to endorse or promote products derived  | 
50  |  | from this library without prior written permission.  | 
51  |  |  | 
52  |  | You may not charge copyright license fees for anyone to use, copy or distribute  | 
53  |  | the FDK AAC Codec software or your modifications thereto.  | 
54  |  |  | 
55  |  | Your modified versions of the FDK AAC Codec must carry prominent notices stating  | 
56  |  | that you changed the software and the date of any change. For modified versions  | 
57  |  | of the FDK AAC Codec, the term "Fraunhofer FDK AAC Codec Library for Android"  | 
58  |  | must be replaced by the term "Third-Party Modified Version of the Fraunhofer FDK  | 
59  |  | AAC Codec Library for Android."  | 
60  |  |  | 
61  |  | 3.    NO PATENT LICENSE  | 
62  |  |  | 
63  |  | NO EXPRESS OR IMPLIED LICENSES TO ANY PATENT CLAIMS, including without  | 
64  |  | limitation the patents of Fraunhofer, ARE GRANTED BY THIS SOFTWARE LICENSE.  | 
65  |  | Fraunhofer provides no warranty of patent non-infringement with respect to this  | 
66  |  | software.  | 
67  |  |  | 
68  |  | You may use this FDK AAC Codec software or modifications thereto only for  | 
69  |  | purposes that are authorized by appropriate patent licenses.  | 
70  |  |  | 
71  |  | 4.    DISCLAIMER  | 
72  |  |  | 
73  |  | This FDK AAC Codec software is provided by Fraunhofer on behalf of the copyright  | 
74  |  | holders and contributors "AS IS" and WITHOUT ANY EXPRESS OR IMPLIED WARRANTIES,  | 
75  |  | including but not limited to the implied warranties of merchantability and  | 
76  |  | fitness for a particular purpose. IN NO EVENT SHALL THE COPYRIGHT HOLDER OR  | 
77  |  | CONTRIBUTORS BE LIABLE for any direct, indirect, incidental, special, exemplary,  | 
78  |  | or consequential damages, including but not limited to procurement of substitute  | 
79  |  | goods or services; loss of use, data, or profits, or business interruption,  | 
80  |  | however caused and on any theory of liability, whether in contract, strict  | 
81  |  | liability, or tort (including negligence), arising in any way out of the use of  | 
82  |  | this software, even if advised of the possibility of such damage.  | 
83  |  |  | 
84  |  | 5.    CONTACT INFORMATION  | 
85  |  |  | 
86  |  | Fraunhofer Institute for Integrated Circuits IIS  | 
87  |  | Attention: Audio and Multimedia Departments - FDK AAC LL  | 
88  |  | Am Wolfsmantel 33  | 
89  |  | 91058 Erlangen, Germany  | 
90  |  |  | 
91  |  | www.iis.fraunhofer.de/amm  | 
92  |  | amm-info@iis.fraunhofer.de  | 
93  |  | ----------------------------------------------------------------------------- */  | 
94  |  |  | 
95  |  | /******************* Library for basic calculation routines ********************  | 
96  |  |  | 
97  |  |    Author(s):  | 
98  |  |  | 
99  |  |    Description: fixed point intrinsics  | 
100  |  |  | 
101  |  | *******************************************************************************/  | 
102  |  |  | 
103  |  | #if !defined(CPLX_MUL_H)  | 
104  |  | #define CPLX_MUL_H  | 
105  |  |  | 
106  |  | #include "common_fix.h"  | 
107  |  |  | 
108  |  | #if defined(__arm__) || defined(_M_ARM)  | 
109  |  | #include "arm/cplx_mul_arm.h"  | 
110  |  |  | 
111  |  | #elif defined(__GNUC__) && defined(__mips__) && __mips_isa_rev < 6  | 
112  |  | #include "mips/cplx_mul_mips.h"  | 
113  |  |  | 
114  |  | #endif /* #if defined all cores: bfin, arm, etc. */  | 
115  |  |  | 
116  |  | /* #############################################################################  | 
117  |  |  */  | 
118  |  |  | 
119  |  | /* Fallback generic implementations */  | 
120  |  |  | 
121  |  | #if !defined(FUNCTION_cplxMultDiv2_32x16X2)  | 
122  |  | #define FUNCTION_cplxMultDiv2_32x16X2  | 
123  |  |  | 
124  |  | inline void cplxMultDiv2(FIXP_DBL *c_Re, FIXP_DBL *c_Im, const FIXP_DBL a_Re,  | 
125  |  |                          const FIXP_DBL a_Im, const FIXP_SGL b_Re,  | 
126  | 3.89G  |                          const FIXP_SGL b_Im) { | 
127  | 3.89G  |   *c_Re = fMultDiv2(a_Re, b_Re) - fMultDiv2(a_Im, b_Im);  | 
128  | 3.89G  |   *c_Im = fMultDiv2(a_Re, b_Im) + fMultDiv2(a_Im, b_Re);  | 
129  | 3.89G  | }  | 
130  |  | #endif  | 
131  |  |  | 
132  |  | #if !defined(FUNCTION_cplxMultDiv2_16x16X2)  | 
133  |  | #define FUNCTION_cplxMultDiv2_16x16X2  | 
134  |  |  | 
135  |  | inline void cplxMultDiv2(FIXP_DBL *c_Re, FIXP_DBL *c_Im, const FIXP_SGL a_Re,  | 
136  |  |                          const FIXP_SGL a_Im, const FIXP_SGL b_Re,  | 
137  | 0  |                          const FIXP_SGL b_Im) { | 
138  | 0  |   *c_Re = fMultDiv2(a_Re, b_Re) - fMultDiv2(a_Im, b_Im);  | 
139  | 0  |   *c_Im = fMultDiv2(a_Re, b_Im) + fMultDiv2(a_Im, b_Re);  | 
140  | 0  | }  | 
141  |  |  | 
142  |  | inline void cplxMultDiv2(FIXP_SGL *c_Re, FIXP_SGL *c_Im, const FIXP_SGL a_Re,  | 
143  |  |                          const FIXP_SGL a_Im, const FIXP_SGL b_Re,  | 
144  | 0  |                          const FIXP_SGL b_Im) { | 
145  | 0  |   *c_Re = FX_DBL2FX_SGL(fMultDiv2(a_Re, b_Re) - fMultDiv2(a_Im, b_Im));  | 
146  | 0  |   *c_Im = FX_DBL2FX_SGL(fMultDiv2(a_Re, b_Im) + fMultDiv2(a_Im, b_Re));  | 
147  | 0  | }  | 
148  |  | #endif  | 
149  |  |  | 
150  |  | #if !defined(FUNCTION_cplxMultDiv2_32x16)  | 
151  |  | #define FUNCTION_cplxMultDiv2_32x16  | 
152  |  |  | 
153  |  | inline void cplxMultDiv2(FIXP_DBL *c_Re, FIXP_DBL *c_Im, const FIXP_DBL a_Re,  | 
154  | 3.60G  |                          const FIXP_DBL a_Im, const FIXP_SPK w) { | 
155  | 3.60G  |   cplxMultDiv2(c_Re, c_Im, a_Re, a_Im, w.v.re, w.v.im);  | 
156  | 3.60G  | }  | 
157  |  | #endif  | 
158  |  |  | 
159  |  | #if !defined(FUNCTION_cplxMultDiv2_16x16)  | 
160  |  | #define FUNCTION_cplxMultDiv2_16x16  | 
161  |  |  | 
162  |  | inline void cplxMultDiv2(FIXP_DBL *c_Re, FIXP_DBL *c_Im, const FIXP_SGL a_Re,  | 
163  | 0  |                          const FIXP_SGL a_Im, const FIXP_SPK w) { | 
164  | 0  |   cplxMultDiv2(c_Re, c_Im, a_Re, a_Im, w.v.re, w.v.im);  | 
165  | 0  | }  | 
166  |  |  | 
167  |  | inline void cplxMultDiv2(FIXP_SGL *c_Re, FIXP_SGL *c_Im, const FIXP_SGL a_Re,  | 
168  | 0  |                          const FIXP_SGL a_Im, const FIXP_SPK w) { | 
169  | 0  |   cplxMultDiv2(c_Re, c_Im, a_Re, a_Im, w.v.re, w.v.im);  | 
170  | 0  | }  | 
171  |  | #endif  | 
172  |  |  | 
173  |  | #if !defined(FUNCTION_cplxMultSubDiv2_32x16X2)  | 
174  |  | #define FUNCTION_cplxMultSubDiv2_32x16X2  | 
175  |  |  | 
176  |  | inline void cplxMultSubDiv2(FIXP_DBL *c_Re, FIXP_DBL *c_Im, const FIXP_DBL a_Re,  | 
177  |  |                             const FIXP_DBL a_Im, const FIXP_SGL b_Re,  | 
178  | 0  |                             const FIXP_SGL b_Im) { | 
179  | 0  |   *c_Re -= fMultDiv2(a_Re, b_Re) - fMultDiv2(a_Im, b_Im);  | 
180  | 0  |   *c_Im -= fMultDiv2(a_Re, b_Im) + fMultDiv2(a_Im, b_Re);  | 
181  | 0  | }  | 
182  |  | #endif  | 
183  |  |  | 
184  |  | #if !defined(FUNCTION_cplxMultDiv2_32x32X2)  | 
185  |  | #define FUNCTION_cplxMultDiv2_32x32X2  | 
186  |  |  | 
187  |  | inline void cplxMultDiv2(FIXP_DBL *c_Re, FIXP_DBL *c_Im, const FIXP_DBL a_Re,  | 
188  |  |                          const FIXP_DBL a_Im, const FIXP_DBL b_Re,  | 
189  | 425M  |                          const FIXP_DBL b_Im) { | 
190  | 425M  |   *c_Re = fMultDiv2(a_Re, b_Re) - fMultDiv2(a_Im, b_Im);  | 
191  | 425M  |   *c_Im = fMultDiv2(a_Re, b_Im) + fMultDiv2(a_Im, b_Re);  | 
192  | 425M  | }  | 
193  |  | #endif  | 
194  |  |  | 
195  |  | #if !defined(FUNCTION_cplxMultDiv2_32x32)  | 
196  |  | #define FUNCTION_cplxMultDiv2_32x32  | 
197  |  |  | 
198  |  | inline void cplxMultDiv2(FIXP_DBL *c_Re, FIXP_DBL *c_Im, const FIXP_DBL a_Re,  | 
199  | 391k  |                          const FIXP_DBL a_Im, const FIXP_DPK w) { | 
200  | 391k  |   cplxMultDiv2(c_Re, c_Im, a_Re, a_Im, w.v.re, w.v.im);  | 
201  | 391k  | }  | 
202  |  | #endif  | 
203  |  |  | 
204  |  | #if !defined(FUNCTION_cplxMultSubDiv2_32x32X2)  | 
205  |  | #define FUNCTION_cplxMultSubDiv2_32x32X2  | 
206  |  |  | 
207  |  | inline void cplxMultSubDiv2(FIXP_DBL *c_Re, FIXP_DBL *c_Im, const FIXP_DBL a_Re,  | 
208  |  |                             const FIXP_DBL a_Im, const FIXP_DBL b_Re,  | 
209  | 0  |                             const FIXP_DBL b_Im) { | 
210  | 0  |   *c_Re -= fMultDiv2(a_Re, b_Re) - fMultDiv2(a_Im, b_Im);  | 
211  | 0  |   *c_Im -= fMultDiv2(a_Re, b_Im) + fMultDiv2(a_Im, b_Re);  | 
212  | 0  | }  | 
213  |  | #endif  | 
214  |  |  | 
215  |  |   /* #############################################################################  | 
216  |  |    */  | 
217  |  |  | 
218  |  | #if !defined(FUNCTION_cplxMult_32x16X2)  | 
219  |  | #define FUNCTION_cplxMult_32x16X2  | 
220  |  |  | 
221  |  | inline void cplxMult(FIXP_DBL *c_Re, FIXP_DBL *c_Im, const FIXP_DBL a_Re,  | 
222  |  |                      const FIXP_DBL a_Im, const FIXP_SGL b_Re,  | 
223  | 2.34G  |                      const FIXP_SGL b_Im) { | 
224  | 2.34G  |   *c_Re = fMult(a_Re, b_Re) - fMult(a_Im, b_Im);  | 
225  | 2.34G  |   *c_Im = fMult(a_Re, b_Im) + fMult(a_Im, b_Re);  | 
226  | 2.34G  | }  | 
227  |  | inline void cplxMult(FIXP_SGL *c_Re, FIXP_SGL *c_Im, const FIXP_SGL a_Re,  | 
228  |  |                      const FIXP_SGL a_Im, const FIXP_SGL b_Re,  | 
229  | 0  |                      const FIXP_SGL b_Im) { | 
230  | 0  |   *c_Re = FX_DBL2FX_SGL(fMult(a_Re, b_Re) - fMult(a_Im, b_Im));  | 
231  | 0  |   *c_Im = FX_DBL2FX_SGL(fMult(a_Re, b_Im) + fMult(a_Im, b_Re));  | 
232  | 0  | }  | 
233  |  | #endif  | 
234  |  |  | 
235  |  | #if !defined(FUNCTION_cplxMult_32x16)  | 
236  |  | #define FUNCTION_cplxMult_32x16  | 
237  |  |  | 
238  |  | inline void cplxMult(FIXP_DBL *c_Re, FIXP_DBL *c_Im, const FIXP_DBL a_Re,  | 
239  | 1.93G  |                      const FIXP_DBL a_Im, const FIXP_SPK w) { | 
240  | 1.93G  |   cplxMult(c_Re, c_Im, a_Re, a_Im, w.v.re, w.v.im);  | 
241  | 1.93G  | }  | 
242  |  | #endif  | 
243  |  |  | 
244  |  | #if !defined(FUNCTION_cplxMult_32x32X2)  | 
245  |  | #define FUNCTION_cplxMult_32x32X2  | 
246  |  |  | 
247  |  | inline void cplxMult(FIXP_DBL *c_Re, FIXP_DBL *c_Im, const FIXP_DBL a_Re,  | 
248  |  |                      const FIXP_DBL a_Im, const FIXP_DBL b_Re,  | 
249  | 0  |                      const FIXP_DBL b_Im) { | 
250  | 0  |   *c_Re = fMult(a_Re, b_Re) - fMult(a_Im, b_Im);  | 
251  | 0  |   *c_Im = fMult(a_Re, b_Im) + fMult(a_Im, b_Re);  | 
252  | 0  | }  | 
253  |  | #endif  | 
254  |  |  | 
255  |  | #if !defined(FUNCTION_cplxMult_32x32)  | 
256  |  | #define FUNCTION_cplxMult_32x32  | 
257  |  | inline void cplxMult(FIXP_DBL *c_Re, FIXP_DBL *c_Im, const FIXP_DBL a_Re,  | 
258  | 0  |                      const FIXP_DBL a_Im, const FIXP_DPK w) { | 
259  | 0  |   cplxMult(c_Re, c_Im, a_Re, a_Im, w.v.re, w.v.im);  | 
260  | 0  | }  | 
261  |  | #endif  | 
262  |  |  | 
263  |  |   /* #############################################################################  | 
264  |  |    */  | 
265  |  |  | 
266  |  | #endif /* CPLX_MUL_H */  |