• Home
  • Line#
  • Scopes#
  • Navigate#
  • Raw
  • Download
1 /******************************************************************************
2  *                                                                            *
3  * Copyright (C) 2018 The Android Open Source Project
4  *
5  * Licensed under the Apache License, Version 2.0 (the "License");
6  * you may not use this file except in compliance with the License.
7  * You may obtain a copy of the License at:
8  *
9  * http://www.apache.org/licenses/LICENSE-2.0
10  *
11  * Unless required by applicable law or agreed to in writing, software
12  * distributed under the License is distributed on an "AS IS" BASIS,
13  * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
14  * See the License for the specific language governing permissions and
15  * limitations under the License.
16  *
17  *****************************************************************************
18  * Originally developed and contributed by Ittiam Systems Pvt. Ltd, Bangalore
19  */
20 
21 #include <stdlib.h>
22 #include <stdio.h>
23 #include <string.h>
24 
25 #include <ixheaacd_type_def.h>
26 #include "ixheaacd_constants.h"
27 #include <ixheaacd_basic_ops32.h>
28 #include <ixheaacd_basic_ops16.h>
29 #include <ixheaacd_basic_ops40.h>
30 
31 #include "ixheaacd_defines.h"
32 #include <ixheaacd_aac_rom.h>
33 #include "ixheaacd_audioobjtypes.h"
34 
35 #include "ixheaacd_bitbuffer.h"
36 #include "ixheaacd_pulsedata.h"
37 #include "ixheaacd_pns.h"
38 #include "ixheaacd_lt_predict.h"
39 #include "ixheaacd_channelinfo.h"
40 #include "ixheaacd_cnst.h"
41 #include "ixheaacd_tns.h"
42 #include "ixheaacd_aac_imdct.h"
43 
44 static const WORD32 ixheaacd_codebook_Q30[8] = {
45     612922971,  747985734,  872956397,  978505219,
46     1057528322, 1146642451, 1282693056, 1470524861};
47 
48 #define SHIFT_VAL 8
49 #define SHIFT_VAL1 (15 - SHIFT_VAL)
50 
ixheaacd_lt_prediction(ia_aac_dec_channel_info_struct * ptr_aac_dec_channel_info,ltp_info * ltp,WORD32 * spec,ia_aac_dec_tables_struct * aac_tables_ptr,UWORD16 win_shape_prev,UWORD32 sr_index,UWORD32 object_type,UWORD32 frame_len,WORD32 * in_data,WORD32 * out_data)51 VOID ixheaacd_lt_prediction(
52     ia_aac_dec_channel_info_struct *ptr_aac_dec_channel_info, ltp_info *ltp,
53     WORD32 *spec, ia_aac_dec_tables_struct *aac_tables_ptr,
54     UWORD16 win_shape_prev, UWORD32 sr_index, UWORD32 object_type,
55     UWORD32 frame_len, WORD32 *in_data, WORD32 *out_data) {
56   ia_ics_info_struct *ptr_ics_info = &ptr_aac_dec_channel_info->str_ics_info;
57   WORD16 *lt_pred_stat = ptr_aac_dec_channel_info->ltp_buf;
58   UWORD16 win_shape = ptr_aac_dec_channel_info->str_ics_info.window_shape;
59   WORD16 sfb;
60   WORD16 bin, i, num_samples;
61   const WORD8 *swb_offset = aac_tables_ptr->scale_factor_bands_long[sr_index];
62   WORD32 *ptr_spec = &spec[0];
63   WORD32 *ptr_x_est = &out_data[0];
64 
65   if (512 == ptr_ics_info->frame_length) {
66     swb_offset = aac_tables_ptr->scale_fac_bands_512[sr_index];
67   } else if (480 == ptr_ics_info->frame_length) {
68     swb_offset = aac_tables_ptr->scale_fac_bands_480[sr_index];
69   }
70 
71   if (ptr_ics_info->window_sequence != EIGHT_SHORT_SEQUENCE) {
72     if (ltp->data_present) {
73       num_samples = frame_len << 1;
74 
75       for (i = 0; i < num_samples; i++) {
76         in_data[i] =
77             ixheaacd_shr32(ixheaacd_mult32x16in32_shl_sat(
78                                ixheaacd_codebook_Q30[ltp->coef],
79                                lt_pred_stat[num_samples + i - ltp->lag]),
80                            SHIFT_VAL);
81       }
82 
83       ixheaacd_filter_bank_ltp(aac_tables_ptr, ptr_ics_info->window_sequence,
84                                win_shape, win_shape_prev, in_data, out_data,
85                                object_type, frame_len);
86 
87       if (ptr_aac_dec_channel_info->str_tns_info.tns_data_present == 1)
88         ixheaacd_aac_tns_process(ptr_aac_dec_channel_info, 1, aac_tables_ptr,
89                                  object_type, 0, out_data);
90 
91       for (sfb = 0; sfb < ltp->last_band; sfb++) {
92         WORD8 sfb_width = swb_offset[sfb];
93         if (ltp->long_used[sfb]) {
94           for (bin = sfb_width - 1; bin >= 0; bin--) {
95             WORD32 temp = *ptr_spec;
96             temp = ixheaacd_add32_sat(temp,
97                                       ixheaacd_shr32(*ptr_x_est++, SHIFT_VAL1));
98             *ptr_spec++ = temp;
99           }
100         } else {
101           ptr_spec += sfb_width;
102           ptr_x_est += sfb_width;
103         }
104       }
105     }
106   }
107 }
108 
ixheaacd_filter_bank_ltp(ia_aac_dec_tables_struct * aac_tables_ptr,WORD16 window_sequence,WORD16 window_shape,WORD16 window_shape_prev,WORD32 * in_data,WORD32 * out_mdct,UWORD32 object_type,UWORD32 frame_len)109 VOID ixheaacd_filter_bank_ltp(ia_aac_dec_tables_struct *aac_tables_ptr,
110                               WORD16 window_sequence, WORD16 window_shape,
111                               WORD16 window_shape_prev, WORD32 *in_data,
112                               WORD32 *out_mdct, UWORD32 object_type,
113                               UWORD32 frame_len) {
114   WORD32 i;
115 
116   const WORD16 *window_long = NULL;
117   const WORD16 *window_long_prev = NULL;
118   const WORD16 *window_short = NULL;
119   const WORD16 *window_short_prev = NULL;
120 
121   UWORD16 nlong = frame_len;
122   UWORD16 nlong2 = frame_len << 1;
123   UWORD16 nshort = frame_len / 8;
124   UWORD16 nflat_ls = (nlong - nshort) / 2;
125   WORD32 imdct_scale = 0;
126   WORD32 expo = 0;
127 
128   if (object_type == AOT_ER_AAC_LD) {
129     if (!window_shape) {
130       if (512 == frame_len) {
131         window_long =
132             (WORD16 *)aac_tables_ptr->pstr_imdct_tables->window_sine_512;
133       } else {
134         window_long =
135             (WORD16 *)aac_tables_ptr->pstr_imdct_tables->window_sine_480;
136       }
137     } else {
138       if (512 == frame_len) {
139         window_long =
140             (WORD16 *)aac_tables_ptr->pstr_imdct_tables->low_overlap_win;
141       } else {
142         window_long =
143             (WORD16 *)aac_tables_ptr->pstr_imdct_tables->low_overlap_win_480;
144       }
145     }
146 
147     if (!window_shape_prev) {
148       if (512 == frame_len) {
149         window_long_prev =
150             (WORD16 *)aac_tables_ptr->pstr_imdct_tables->window_sine_512;
151       } else {
152         window_long_prev =
153             (WORD16 *)aac_tables_ptr->pstr_imdct_tables->window_sine_480;
154       }
155     } else {
156       if (512 == frame_len) {
157         window_long_prev =
158             (WORD16 *)aac_tables_ptr->pstr_imdct_tables->low_overlap_win;
159       } else {
160         window_long_prev =
161             (WORD16 *)aac_tables_ptr->pstr_imdct_tables->low_overlap_win_480;
162       }
163     }
164 
165     if (!window_shape)
166       window_short = aac_tables_ptr->pstr_imdct_tables->only_short_window_sine;
167     else
168       window_short = aac_tables_ptr->pstr_imdct_tables->only_short_window_kbd;
169     if (!window_shape_prev)
170       window_short_prev =
171           aac_tables_ptr->pstr_imdct_tables->only_short_window_sine;
172     else
173       window_short_prev =
174           aac_tables_ptr->pstr_imdct_tables->only_short_window_kbd;
175 
176   } else {
177     if (!window_shape)
178       window_long = aac_tables_ptr->pstr_imdct_tables->only_long_window_sine;
179     else
180       window_long = aac_tables_ptr->pstr_imdct_tables->only_long_window_kbd;
181     if (!window_shape_prev)
182       window_long_prev =
183           aac_tables_ptr->pstr_imdct_tables->only_long_window_sine;
184     else
185       window_long_prev =
186           aac_tables_ptr->pstr_imdct_tables->only_long_window_kbd;
187 
188     if (!window_shape)
189       window_short = aac_tables_ptr->pstr_imdct_tables->only_short_window_sine;
190     else
191       window_short = aac_tables_ptr->pstr_imdct_tables->only_short_window_kbd;
192     if (!window_shape_prev)
193       window_short_prev =
194           aac_tables_ptr->pstr_imdct_tables->only_short_window_sine;
195     else
196       window_short_prev =
197           aac_tables_ptr->pstr_imdct_tables->only_short_window_kbd;
198   }
199 
200   switch (window_sequence) {
201     case ONLY_LONG_SEQUENCE:
202 
203       if ((512 != nlong) && (480 != nlong)) {
204         for (i = 0; i<nlong>> 1; i++) {
205           in_data[i] =
206               ixheaacd_mult32x16in32_shl(in_data[i], window_long_prev[2 * i]);
207 
208           in_data[i + nlong] = ixheaacd_mult32x16in32_shl(
209               in_data[i + nlong], window_long[2 * i + 1]);
210         }
211         for (i = 0; i<nlong>> 1; i++) {
212           in_data[i + (nlong >> 1)] = ixheaacd_mult32x16in32_shl(
213               in_data[i + (nlong >> 1)], window_long_prev[nlong - 1 - 2 * i]);
214 
215           in_data[i + nlong + (nlong >> 1)] =
216               ixheaacd_mult32x16in32_shl(in_data[i + nlong + (nlong >> 1)],
217                                          window_long[nlong - 1 - 2 * i - 1]);
218         }
219 
220       } else {
221         WORD32 *win1, *win2, *win3;
222         WORD32 *ptr_in1, *ptr_in2;
223         win1 = (WORD32 *)window_long_prev;
224         win2 = (WORD32 *)window_long;
225         ptr_in1 = &in_data[0];
226         ptr_in2 = &in_data[nlong];
227         win3 = win2 + nlong - 1;
228 
229         for (i = nlong - 1; i >= 0; i--) {
230           WORD32 temp1 = ixheaacd_mult32_shl(*ptr_in1, *win1++);
231           WORD32 temp2 = ixheaacd_mult32_shl(*ptr_in2, win2[i]);
232 
233           *ptr_in1++ = temp1;
234           *ptr_in2++ = temp2;
235         }
236       }
237 
238       for (i = 0; i < nlong / 2; i++) {
239         out_mdct[nlong / 2 + i] =
240             ixheaacd_sub32(in_data[i], in_data[nlong - 1 - i]);
241         out_mdct[i] = (-ixheaacd_add32(in_data[nlong + i + nlong / 2],
242                                        in_data[nlong2 - nlong / 2 - 1 - i]));
243       }
244 
245       if (512 == nlong || (480 == nlong)) {
246         if (512 == nlong)
247           ixheaacd_inverse_transform_512(
248               out_mdct, in_data, &imdct_scale,
249               aac_tables_ptr->pstr_imdct_tables->cosine_array_1024,
250               aac_tables_ptr->pstr_imdct_tables, object_type);
251 
252         else
253           ixheaacd_mdct_480_ld(out_mdct, in_data, &imdct_scale, 1,
254                                aac_tables_ptr->pstr_imdct_tables, object_type);
255 
256         imdct_scale += 1;
257 
258         if (imdct_scale > 0) {
259           WORD32 *ptr_out_mdct = &out_mdct[0];
260 
261           for (i = nlong - 1; i >= 0; i -= 4) {
262             *ptr_out_mdct = ixheaacd_shl32(*ptr_out_mdct, imdct_scale);
263             ptr_out_mdct++;
264             *ptr_out_mdct = ixheaacd_shl32(*ptr_out_mdct, imdct_scale);
265             ptr_out_mdct++;
266             *ptr_out_mdct = ixheaacd_shl32(*ptr_out_mdct, imdct_scale);
267             ptr_out_mdct++;
268             *ptr_out_mdct = ixheaacd_shl32(*ptr_out_mdct, imdct_scale);
269             ptr_out_mdct++;
270           }
271         } else if (imdct_scale < 0) {
272           WORD32 *ptr_out_mdct = &out_mdct[0];
273           imdct_scale = -imdct_scale;
274           for (i = nlong - 1; i >= 0; i -= 4) {
275             *ptr_out_mdct = ixheaacd_shr32(*ptr_out_mdct, imdct_scale);
276             ptr_out_mdct++;
277             *ptr_out_mdct = ixheaacd_shr32(*ptr_out_mdct, imdct_scale);
278             ptr_out_mdct++;
279             *ptr_out_mdct = ixheaacd_shr32(*ptr_out_mdct, imdct_scale);
280             ptr_out_mdct++;
281             *ptr_out_mdct = ixheaacd_shr32(*ptr_out_mdct, imdct_scale);
282             ptr_out_mdct++;
283           }
284         }
285       }
286 
287       else if (1024 == nlong) {
288         expo = ixheaacd_calc_max_spectral_line_dec(out_mdct, 1024) - 1;
289 
290         expo = 8 - expo;
291 
292         imdct_scale = ixheaacd_inverse_transform(
293             out_mdct, in_data, aac_tables_ptr->pstr_imdct_tables, expo, 1024);
294 
295         ixheaacd_post_twiddle_dec(in_data, out_mdct,
296                                   aac_tables_ptr->pstr_imdct_tables, 1024);
297 
298         imdct_scale += 1;
299 
300         for (i = 0; i < nlong; i++) {
301           out_mdct[i] = ixheaacd_shl32_dir(in_data[i], imdct_scale);
302         }
303       }
304 
305       break;
306 
307     case LONG_START_SEQUENCE:
308 
309       for (i = 0; i<nlong>> 1; i++)
310         in_data[i] =
311             ixheaacd_mult32x16in32_shl(in_data[i], window_long_prev[2 * i]);
312 
313       for (i = 0; i<nlong>> 1; i++)
314         in_data[i + (nlong >> 1)] = ixheaacd_mult32x16in32_shl(
315             in_data[i + (nlong >> 1)], window_long_prev[nlong - 1 - 2 * i - 1]);
316 
317       for (i = 0; i<nshort>> 1; i++)
318         in_data[i + nlong + nflat_ls + (nshort >> 1)] =
319             ixheaacd_mult32x16in32_shl(
320                 in_data[i + nlong + nflat_ls + (nshort >> 1)],
321                 window_short[nshort - 1 - 2 * i - 1]);
322 
323       for (i = 0; i < nflat_ls; i++) in_data[i + nlong + nflat_ls + nshort] = 0;
324 
325       for (i = 0; i < nlong / 2; i++) {
326         out_mdct[nlong / 2 + i] =
327             ixheaacd_sub32(in_data[i], in_data[nlong - 1 - i]);
328         out_mdct[nlong / 2 - 1 - i] =
329             -ixheaacd_add32(in_data[nlong + i], in_data[nlong2 - 1 - i]);
330       }
331 
332       {
333         expo = ixheaacd_calc_max_spectral_line_dec(out_mdct, 1024) - 1;
334 
335         expo = 8 - expo;
336         imdct_scale = ixheaacd_inverse_transform(
337             out_mdct, in_data, aac_tables_ptr->pstr_imdct_tables, expo, 1024);
338 
339         ixheaacd_post_twiddle_dec(in_data, out_mdct,
340                                   aac_tables_ptr->pstr_imdct_tables, 1024);
341       }
342 
343       imdct_scale += 1;
344 
345       for (i = 0; i < nlong; i++) {
346         out_mdct[i] = ixheaacd_shl32_dir(in_data[i], imdct_scale);
347       }
348       break;
349 
350     case LONG_STOP_SEQUENCE:
351       for (i = 0; i < nflat_ls; i++) in_data[i] = 0;
352 
353       for (i = 0; i<nshort>> 1; i++)
354         in_data[i + nflat_ls] = ixheaacd_mult32x16in32_shl(
355             in_data[i + nflat_ls], window_short_prev[2 * i]);
356 
357       for (i = 0; i<nshort>> 1; i++)
358         in_data[i + nflat_ls + (nshort >> 1)] =
359             ixheaacd_mult32x16in32_shl(in_data[i + nflat_ls + (nshort >> 1)],
360                                        window_short_prev[127 - 2 * i]);
361 
362       for (i = 0; i<nlong>> 1; i++)
363         in_data[i + nlong] = ixheaacd_mult32x16in32_shl(in_data[i + nlong],
364                                                         window_long[2 * i + 1]);
365 
366       for (i = 0; i<nlong>> 1; i++)
367         in_data[i + nlong + (nlong >> 1)] =
368             ixheaacd_mult32x16in32_shl(in_data[i + nlong + (nlong >> 1)],
369                                        window_long[nlong - 1 - 2 * i - 1]);
370 
371       for (i = 0; i < nlong / 2; i++) {
372         out_mdct[nlong / 2 + i] =
373             ixheaacd_sub32(in_data[i], in_data[nlong - 1 - i]);
374         out_mdct[nlong / 2 - 1 - i] =
375             -ixheaacd_add32(in_data[nlong + i], in_data[nlong2 - 1 - i]);
376       }
377 
378       {
379         expo = ixheaacd_calc_max_spectral_line_dec(out_mdct, 1024) - 1;
380 
381         expo = 8 - expo;
382         imdct_scale = ixheaacd_inverse_transform(
383             out_mdct, in_data, aac_tables_ptr->pstr_imdct_tables, expo, 1024);
384 
385         ixheaacd_post_twiddle_dec(in_data, out_mdct,
386                                   aac_tables_ptr->pstr_imdct_tables, 1024);
387       }
388 
389       imdct_scale += 1;
390 
391       for (i = 0; i < nlong; i++) {
392         out_mdct[i] = ixheaacd_shl32_dir(in_data[i], imdct_scale);
393       }
394 
395       break;
396   }
397 }
398 
ixheaacd_lt_update_state(WORD16 * lt_pred_stat,WORD16 * time,WORD32 * overlap,WORD32 frame_len,WORD32 object_type,WORD32 stride,WORD16 window_sequence,WORD16 * p_window_next)399 VOID ixheaacd_lt_update_state(WORD16 *lt_pred_stat, WORD16 *time,
400                               WORD32 *overlap, WORD32 frame_len,
401                               WORD32 object_type, WORD32 stride,
402                               WORD16 window_sequence, WORD16 *p_window_next) {
403   WORD32 i;
404   if (object_type == AOT_ER_AAC_LD) {
405     WORD16 *ptr_ltp_state0 = &lt_pred_stat[0];
406     WORD16 *ptr_ltp_state_fl = &lt_pred_stat[frame_len + 0];
407     WORD16 *ptr_ltp_state_2fl = &lt_pred_stat[(frame_len * 2) + 0];
408     WORD16 *ptr_time_in = &time[0 * stride];
409 
410     for (i = 0; i < frame_len; i++) {
411       *ptr_ltp_state0++ = *ptr_ltp_state_fl;
412       *ptr_ltp_state_fl++ = *ptr_ltp_state_2fl;
413       *ptr_ltp_state_2fl++ = *ptr_time_in;
414       ptr_time_in += stride;
415     }
416 
417   } else {
418     WORD16 *ptr_ltp_state0 = &lt_pred_stat[0];
419     WORD16 *ptr_ltp_state_fl = &lt_pred_stat[frame_len + 0];
420     WORD16 *ptr_time_in = &time[0 * stride];
421 
422     for (i = 0; i < frame_len; i++) {
423       *ptr_ltp_state0++ = *ptr_ltp_state_fl;
424       *ptr_ltp_state_fl++ = *ptr_time_in;
425       ptr_time_in += stride;
426     }
427   }
428 
429   if ((window_sequence == ONLY_LONG_SEQUENCE) ||
430       (window_sequence == LONG_STOP_SEQUENCE)) {
431     if (512 == frame_len) {
432       WORD32 *window = (WORD32 *)p_window_next;
433 
434       for (i = 0; i < 256; i++) {
435         lt_pred_stat[(frame_len * 3) + i] =
436             ixheaacd_round16(ixheaacd_mult16x16in32_shl(
437                 (WORD16)ixheaacd_shl16(
438                     (WORD16)-ixheaacd_sat16(overlap[255 - i]), 1),
439                 (WORD16)ixheaacd_shr32(window[511 - i], 15)));
440 
441         lt_pred_stat[(frame_len * 3) + 256 + i] =
442             ixheaacd_round16(ixheaacd_mult16x16in32_shl(
443                 (WORD16)ixheaacd_shl16((WORD16)-ixheaacd_sat16(overlap[i]), 1),
444                 (WORD16)ixheaacd_shr32(window[255 - i], 15)));
445       }
446     } else if (480 == frame_len) {
447       WORD32 *window = (WORD32 *)p_window_next;
448 
449       for (i = 0; i < 240; i++) {
450         lt_pred_stat[(frame_len * 3) + i] =
451             ixheaacd_round16(ixheaacd_mult16x16in32_shl(
452                 (WORD16)ixheaacd_shl16(
453                     (WORD16)-ixheaacd_sat16(overlap[239 - i]), 1),
454                 (WORD16)ixheaacd_shr32(window[479 - i], 15)));
455 
456         lt_pred_stat[(frame_len * 3) + 240 + i] =
457             ixheaacd_round16(ixheaacd_mult16x16in32_shl(
458                 (WORD16)ixheaacd_shl16((WORD16)-ixheaacd_sat16(overlap[i]), 1),
459                 (WORD16)ixheaacd_shr32(window[239 - i], 15)));
460       }
461     } else {
462       for (i = 0; i < 512; i++) {
463         lt_pred_stat[(frame_len * 2) + i] = ixheaacd_round16(
464             ixheaacd_shl32_sat(ixheaacd_mult16x16in32_shl(
465                                    (WORD16)-ixheaacd_sat16(overlap[511 - i]),
466                                    p_window_next[2 * i + 1]),
467                                1));
468 
469         lt_pred_stat[(frame_len * 2) + 512 + i] =
470             ixheaacd_round16(ixheaacd_shl32_sat(
471                 ixheaacd_mult16x16in32_shl((WORD16)-ixheaacd_sat16(overlap[i]),
472                                            p_window_next[1023 - 2 * i - 1]),
473                 1));
474       }
475     }
476 
477   } else if (window_sequence == LONG_START_SEQUENCE) {
478     for (i = 0; i < 448; i++) {
479       lt_pred_stat[(frame_len * 2) + i] =
480           ixheaacd_shl16((WORD16)-ixheaacd_sat16(overlap[511 - i]), 1);
481     }
482     for (i = 0; i < 64; i++) {
483       lt_pred_stat[(frame_len * 2) + 448 + i] =
484           ixheaacd_round16(ixheaacd_shl32_sat(
485               ixheaacd_mult16x16in32_shl(
486                   (WORD16)-ixheaacd_sat16(overlap[511 - 448 - i]),
487                   p_window_next[2 * i + 1]),
488               1));
489     }
490     for (i = 0; i < 64; i++) {
491       lt_pred_stat[(frame_len * 2) + 512 + i] =
492           ixheaacd_round16(ixheaacd_shl32_sat(
493               ixheaacd_mult16x16in32_shl((WORD16)-ixheaacd_sat16(overlap[i]),
494                                          p_window_next[127 - 2 * i - 1]),
495               1));
496     }
497     for (i = 576; i < 1024; i++) {
498       lt_pred_stat[(frame_len * 2) + i] = 0;
499     }
500   } else {
501     for (i = 0; i < 448; i++) {
502       lt_pred_stat[(frame_len * 2) + i] =
503           ixheaacd_shl16(ixheaacd_sat16(overlap[i]), 1);
504     }
505     for (i = 0; i < 64; i++) {
506       lt_pred_stat[(frame_len * 2) + 448 + i] = ixheaacd_round16(
507           ixheaacd_shl32_sat(ixheaacd_mult16x16in32_shl(
508                                  (WORD16)-ixheaacd_sat16(overlap[511 - i]),
509                                  p_window_next[2 * i + 1]),
510                              1));
511     }
512     for (i = 0; i < 64; i++) {
513       lt_pred_stat[(frame_len * 2) + 512 + i] = ixheaacd_round16(
514           ixheaacd_shl32_sat(ixheaacd_mult16x16in32_shl(
515                                  (WORD16)-ixheaacd_sat16(overlap[448 + i]),
516                                  p_window_next[127 - 2 * i - 1]),
517                              1));
518     }
519     for (i = 576; i < 1024; i++) {
520       lt_pred_stat[(frame_len * 2) + i] = 0;
521     }
522   }
523 }
524