Line data Source code
1 : /******************************************************************************************************
2 :
3 : (C) 2022-2025 IVAS codec Public Collaboration with portions copyright Dolby International AB, Ericsson AB,
4 : Fraunhofer-Gesellschaft zur Foerderung der angewandten Forschung e.V., Huawei Technologies Co. LTD.,
5 : Koninklijke Philips N.V., Nippon Telegraph and Telephone Corporation, Nokia Technologies Oy, Orange,
6 : Panasonic Holdings Corporation, Qualcomm Technologies, Inc., VoiceAge Corporation, and other
7 : contributors to this repository. All Rights Reserved.
8 :
9 : This software is protected by copyright law and by international treaties.
10 : The IVAS codec Public Collaboration consisting of Dolby International AB, Ericsson AB,
11 : Fraunhofer-Gesellschaft zur Foerderung der angewandten Forschung e.V., Huawei Technologies Co. LTD.,
12 : Koninklijke Philips N.V., Nippon Telegraph and Telephone Corporation, Nokia Technologies Oy, Orange,
13 : Panasonic Holdings Corporation, Qualcomm Technologies, Inc., VoiceAge Corporation, and other
14 : contributors to this repository retain full ownership rights in their respective contributions in
15 : the software. This notice grants no license of any kind, including but not limited to patent
16 : license, nor is any license granted by implication, estoppel or otherwise.
17 :
18 : Contributors are required to enter into the IVAS codec Public Collaboration agreement before making
19 : contributions.
20 :
21 : This software is provided "AS IS", without any express or implied warranties. The software is in the
22 : development stage. It is intended exclusively for experts who have experience with such software and
23 : solely for the purpose of inspection. All implied warranties of non-infringement, merchantability
24 : and fitness for a particular purpose are hereby disclaimed and excluded.
25 :
26 : Any dispute, controversy or claim arising under or in relation to providing this software shall be
27 : submitted to and settled by the final, binding jurisdiction of the courts of Munich, Germany in
28 : accordance with the laws of the Federal Republic of Germany excluding its conflict of law rules and
29 : the United Nations Convention on Contracts on the International Sales of Goods.
30 :
31 : *******************************************************************************************************/
32 :
33 : /*====================================================================================
34 : EVS Codec 3GPP TS26.443 Nov 04, 2021. Version 12.14.0 / 13.10.0 / 14.6.0 / 15.4.0 / 16.3.0
35 : ====================================================================================*/
36 :
37 : #include <assert.h>
38 : #include <stdint.h>
39 : #include "options.h"
40 : #include <math.h>
41 : #include "cnst.h"
42 : #include "prot_fx.h"
43 : #include "rom_com.h"
44 : #include "wmc_auto.h"
45 :
46 :
47 : const PWord16 *getSineWindowTable( Word16 length );
48 :
49 :
50 681141 : void mdct_window_sine_IVAS_updated(
51 : PWord16 *window, /* Qx */
52 : const Word32 Fs, /* Q0 */
53 : const Word16 n, /* Q0 */
54 : const Word16 window_type, /* Q0 */
55 : const Word16 element_mode /* Q0 */
56 : )
57 : {
58 681141 : IF( element_mode == EVS_MONO )
59 : {
60 : const PWord16 *table;
61 0 : table = getSineWindowTable( n );
62 0 : FOR( Word32 i = 0; i < n / 2; i++ )
63 : {
64 0 : window[i].v.re = table[i].v.re; /* Qx */
65 0 : move16();
66 0 : window[i].v.im = table[i].v.im; /* Qx */
67 0 : move16();
68 : }
69 : // PMT("getSineWindowTable needs to be updated for IVAS")
70 : }
71 : ELSE
72 : {
73 681141 : const Word16 *window_table = 0;
74 681141 : Word16 buf_in_size = 0;
75 681141 : move16();
76 : Word16 temp[420];
77 681141 : set16_fx( temp, 0, 420 );
78 681141 : SWITCH( window_type )
79 : {
80 436653 : case FULL_OVERLAP:
81 436653 : window_table = tcx_mdct_window_48_fx; /* Q15 */
82 436653 : buf_in_size = 420;
83 436653 : move16();
84 436653 : BREAK;
85 122244 : case HALF_OVERLAP:
86 122244 : window_table = tcx_mdct_window_half_48_fx; /* Q15 */
87 122244 : buf_in_size = 180;
88 122244 : BREAK;
89 122244 : case TRANSITION_OVERLAP:
90 : case MIN_OVERLAP:
91 122244 : window_table = tcx_mdct_window_trans_48_fx; /* Q15 */
92 122244 : buf_in_size = 60;
93 122244 : move16();
94 122244 : BREAK;
95 :
96 0 : default:
97 0 : assert( 0 && "Unsupported window type" );
98 : BREAK;
99 : }
100 :
101 681141 : IF( EQ_32( Fs, 48000 ) )
102 : {
103 148378 : Copy( window_table, temp, n ); /* Q15 */
104 : }
105 : ELSE
106 : {
107 532763 : lerp( window_table, temp, n, buf_in_size );
108 : }
109 :
110 45907437 : FOR( Word32 i = 0; i < n / 2; i++ )
111 : {
112 45226296 : window[i].v.re = temp[n - 1 - i]; /* Qx */
113 45226296 : move16();
114 45226296 : window[i].v.im = temp[i]; /* Qx */
115 45226296 : move16();
116 : }
117 : }
118 681141 : }
119 :
120 : ////Use mdct_window_sine_IVAS_updated for IVAS (EVS path covered as well)
121 808 : void mdct_window_sine(
122 : PWord16 *window, /* Qx */
123 : Word16 n /* Q0 */
124 : )
125 : {
126 : {
127 : const PWord16 *table;
128 808 : table = getSineWindowTable( n );
129 43264 : FOR( Word32 i = 0; i < n / 2; i++ )
130 : {
131 42456 : window[i].v.re = table[i].v.re; /* Qx */
132 42456 : move16();
133 42456 : window[i].v.im = table[i].v.im; /* Qx */
134 42456 : move16();
135 : }
136 : // PMT("getSineWindowTable needs to be updated for IVAS")
137 : }
138 808 : }
139 :
140 :
141 126360 : void mdct_window_aldo(
142 : Word16 *window1, /* Q15 */
143 : PWord16 *window1_trunc, /* Q15 */
144 : PWord16 *window2, /* Q15 */
145 : Word16 n /* Q0 */
146 : )
147 : {
148 : Word16 i, n0, n1, n2, d, tmp;
149 : const Word16 *p1, *p2;
150 :
151 : /* set table pointers and decimation factor */
152 126360 : SWITCH( n )
153 : {
154 0 : case 320 / 2:
155 0 : p1 = window_48kHz_fx + 2; /* Q15 */
156 0 : p2 = window_48kHz_fx + 1110 - 3; /* Q15 */
157 0 : d = 6;
158 0 : move16();
159 0 : BREAK;
160 12372 : case 512 / 2:
161 12372 : p1 = window_256kHz; /* Q15 */
162 12372 : p2 = window_256kHz + 592 - 1; /* Q15 */
163 12372 : d = 2;
164 12372 : move16();
165 12372 : BREAK;
166 23765 : case 640 / 2:
167 23765 : p1 = window_48kHz_fx + 1; /* Q15 */
168 23765 : p2 = window_48kHz_fx + 1110 - 2; /* Q15 */
169 23765 : d = 3;
170 23765 : move16();
171 23765 : BREAK;
172 16280 : case 1024 / 2:
173 16280 : p1 = window_256kHz; /* Q15 */
174 16280 : p2 = window_256kHz + 592 - 1; /* Q15 */
175 16280 : d = 1;
176 16280 : move16();
177 16280 : BREAK;
178 24655 : case 1280 / 2:
179 24655 : p1 = window_48kHz_fx + 1; /* Q15 */
180 24655 : p2 = window_48kHz_fx + 1110 - 2; /* Q15 */
181 24655 : d = 3;
182 24655 : move16();
183 24655 : BREAK;
184 49288 : case 1920 / 2:
185 49288 : p1 = window_48kHz_fx; /* Q15 */
186 49288 : p2 = window_48kHz_fx + 1110 - 1; /* Q15 */
187 49288 : d = 1;
188 49288 : move16();
189 49288 : BREAK;
190 0 : default:
191 0 : assert( 0 );
192 : return;
193 : }
194 :
195 : /* set lengths */
196 126360 : n0 = shr( imult1616( n, 9 ), 5 );
197 126360 : n1 = shr( imult1616( n, 23 ), 5 ); /* left slope length */
198 126360 : n2 = shr( imult1616( n, 14 ), 5 ); /* right slope length */
199 :
200 : /* first part (long slope) */
201 126360 : IF( NE_16( n, 1280 / 2 ) )
202 : {
203 18783419 : FOR( i = 0; i < n0; i++ )
204 : {
205 18681714 : *window1 = *p1; /* Q15 */
206 18681714 : move16();
207 18681714 : window1++;
208 18681714 : p1 += d;
209 : }
210 :
211 101705 : tmp = shr( n, 1 );
212 14631927 : FOR( ; i < tmp; i++ )
213 : {
214 14530222 : window1_trunc->v.im = *p1; /* Q15 */
215 14530222 : move16();
216 14530222 : window1_trunc++;
217 14530222 : p1 += d;
218 : }
219 :
220 101705 : test();
221 101705 : if ( EQ_16( n, 512 / 2 ) || EQ_16( n, 320 / 2 ) )
222 12372 : p1++;
223 :
224 14631927 : FOR( ; i < n1; i++ )
225 : {
226 14530222 : window1_trunc--;
227 14530222 : window1_trunc->v.re = *p1; /* Q15 */
228 14530222 : move16();
229 14530222 : p1 += d;
230 : }
231 : }
232 : ELSE
233 : {
234 24655 : const Word16 *pi = window_8_16_32kHz_fx;
235 :
236 2243605 : FOR( i = 0; i < n0; i += 2 )
237 : {
238 2218950 : *window1 = *p1; /* Q15 */
239 2218950 : move16();
240 2218950 : window1++;
241 2218950 : p1 += d;
242 :
243 2218950 : *window1 = *pi; /* Q15 */
244 2218950 : move16();
245 2218950 : window1++;
246 2218950 : pi++;
247 : }
248 :
249 24655 : tmp = shr( n, 1 );
250 1750505 : FOR( ; i < tmp; i += 2 )
251 : {
252 1725850 : window1_trunc->v.im = *p1; /* Q15 */
253 1725850 : move16();
254 1725850 : window1_trunc++;
255 1725850 : p1 += d;
256 :
257 1725850 : window1_trunc->v.im = *pi; /* Q15 */
258 1725850 : move16();
259 1725850 : window1_trunc++;
260 1725850 : pi++;
261 : }
262 :
263 1750505 : FOR( ; i < n1; i += 2 )
264 : {
265 1725850 : window1_trunc--;
266 1725850 : window1_trunc->v.re = *pi; /* Q15 */
267 1725850 : move16();
268 1725850 : pi++;
269 :
270 1725850 : window1_trunc--;
271 1725850 : window1_trunc->v.re = *p1; /* Q15 */
272 1725850 : move16();
273 1725850 : p1 += d;
274 : }
275 : }
276 :
277 : /* second part (short slope) */
278 126360 : IF( NE_16( n, 1280 / 2 ) )
279 : {
280 101705 : tmp = shr( n2, 1 );
281 14631927 : FOR( i = 0; i < tmp; i++ )
282 : {
283 14530222 : window2->v.im = *p2; /* Q15 */
284 14530222 : move16();
285 14530222 : window2++;
286 14530222 : p2 -= d;
287 : }
288 :
289 101705 : test();
290 101705 : if ( EQ_16( n, 512 / 2 ) || EQ_16( n, 320 / 2 ) )
291 12372 : p2--;
292 :
293 14631927 : FOR( ; i < n2; i++ )
294 : {
295 14530222 : window2--;
296 14530222 : window2->v.re = *p2; /* Q15 */
297 14530222 : move16();
298 14530222 : p2 -= d;
299 : }
300 : }
301 : ELSE
302 : {
303 24655 : const Word16 *pi = window_8_16_32kHz_fx + 370 - 1;
304 :
305 24655 : tmp = shr( n2, 1 );
306 1750505 : FOR( i = 0; i < tmp; i += 2 )
307 : {
308 1725850 : window2->v.im = *p2; /* Q15 */
309 1725850 : move16();
310 1725850 : window2++;
311 1725850 : p2 -= d;
312 :
313 1725850 : window2->v.im = *pi; /* Q15 */
314 1725850 : move16();
315 1725850 : window2++;
316 1725850 : pi--;
317 : }
318 :
319 1750505 : FOR( ; i < n2; i += 2 )
320 : {
321 1725850 : window2--;
322 1725850 : window2->v.re = *pi; /* Q15 */
323 1725850 : move16();
324 1725850 : pi--;
325 :
326 1725850 : window2--;
327 1725850 : window2->v.re = *p2; /* Q15 */
328 1725850 : move16();
329 1725850 : p2 -= d;
330 : }
331 : }
332 126360 : }
|
