/****************************************************************************** * * * Copyright (C) 2018 The Android Open Source Project * * Licensed under the Apache License, Version 2.0 (the "License"); * you may not use this file except in compliance with the License. * You may obtain a copy of the License at: * * http://www.apache.org/licenses/LICENSE-2.0 * * Unless required by applicable law or agreed to in writing, software * distributed under the License is distributed on an "AS IS" BASIS, * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. * See the License for the specific language governing permissions and * limitations under the License. * ***************************************************************************** * Originally developed and contributed by Ittiam Systems Pvt. Ltd, Bangalore */ #include "ixheaacd_sbr_common.h" #include "ixheaac_type_def.h" #include "ixheaac_constants.h" #include "ixheaac_basic_ops32.h" #include "ixheaac_basic_ops16.h" #include "ixheaac_basic_ops40.h" #include "ixheaac_basic_ops.h" #include "ixheaacd_defines.h" #include "ixheaacd_common_rom.h" #include "ixheaacd_basic_funcs.h" #include "ixheaacd_aac_rom.h" #include "ixheaacd_aac_imdct.h" #include "ixheaacd_intrinsics.h" #include "ixheaac_basic_op.h" #include "ixheaacd_function_selector.h" #include "ixheaacd_audioobjtypes.h" #include "ixheaacd_tns.h" #define DIG_REV(i, m, j) \ do { \ unsigned _ = (i); \ _ = ((_ & 0x33333333) << 2) | ((_ & ~0x33333333) >> 2); \ _ = ((_ & 0x0F0F0F0F) << 4) | ((_ & ~0x0F0F0F0F) >> 4); \ _ = ((_ & 0x00FF00FF) << 8) | ((_ & ~0x00FF00FF) >> 8); \ _ = ((_ & 0x0000FFFF) << 16) | ((_ & ~0x0000FFFF) >> 16); \ (j) = _ >> (m); \ } while (0) #define MPYHIRC(x, y) \ \ (((WORD32)((WORD16)(x >> 16) * (UWORD16)(y & 0x0000FFFF) + 0x4000) >> \ 15) + \ ((WORD32)((WORD16)(x >> 16) * (WORD16)((y) >> 16)) << 1)) #define MPYLUHS(x, y) \ ((WORD32)((UWORD16)(x & 0x0000FFFF) * (WORD16)(y >> 16))) #define MPYLIRC(x, y) \ (((WORD32)((WORD16)(x) * (UWORD16)(y & 0x0000FFFF) + 0x4000) >> 15)+ \ ((WORD32)((WORD16)(x) * (WORD16)((y) >> 16)) << 1)) WORD32 rev_dig[] = { 0, 8, 2, 10 }; #define MDCT_LEN 480 #define FFT15X2 30 #define MDCT_LEN_BY2 240 #define FFT5 5 #define FFT16 16 #define FFT4 4 #define FFT3 3 #define FFT15 15 #define FFT16X2 32 #define MDCT_LEN_960 960 WORD32 ixheaacd_fft5out[FFT15X2]; static PLATFORM_INLINE WORD32 ixheaacd_mult32x16lin32(WORD32 a, WORD32 b) { WORD32 result; WORD64 temp_result; temp_result = (WORD64)a * (WORD64)(((b & 0xFFFF) << 16) >> 16); result = (WORD32)(temp_result >> 16); return (result); } static PLATFORM_INLINE WORD32 ixheaacd_mac32x16lin32(WORD32 a, WORD32 b, WORD32 c) { WORD32 result; result = a + ixheaacd_mult32x16lin32(b, c); return (result); } static PLATFORM_INLINE WORD32 ixheaacd_mult32x16lin32_sat(WORD32 a, WORD32 b) { WORD32 result; WORD64 temp_result; temp_result = (WORD64)a * (WORD64)(((b & 0xFFFF) << 16) >> 16); if (temp_result < (WORD64)MIN_32) result = MIN_32; else if (temp_result > (WORD64)MAX_32) result = MAX_32; else result = (WORD32)(temp_result); return (result); } WORD16 ixheaacd_neg_expo_inc_dec(WORD16 neg_expo) { return (neg_expo + 2); } WORD16 ixheaacd_neg_expo_inc_arm(WORD16 neg_expo) { return (neg_expo + 3); } VOID ixheaacd_pretwiddle_compute_960_dec( WORD32 *spec_data1, WORD32 *spec_data2, WORD32 *out_ptr, ia_aac_dec_imdct_tables_struct *ptr_imdct_tables, WORD npoints4, WORD32 neg_expo) { WORD32 i; WORD32 tempr, tempi; WORD16 c, c1, s, s1; WORD32 *out_ptr1 = out_ptr + ((npoints4 << 2) - 1); const WORD16 *cos_sin_ptr = ptr_imdct_tables->cosine_array_240; for (i = 0; i < npoints4; i++) { c = *cos_sin_ptr++; s = *cos_sin_ptr++; tempr = *spec_data1++; tempi = *spec_data2--; *out_ptr = ixheaac_mac32x16in32(ixheaac_mult32x16in32(tempr, c), tempi, s); *out_ptr = ixheaac_shl32(*out_ptr, neg_expo); out_ptr++; *out_ptr = ixheaac_sub32(ixheaac_mult32x16in32(tempi, c), ixheaac_mult32x16in32(tempr, s)); *out_ptr = ixheaac_shl32(*out_ptr, neg_expo); out_ptr++; c1 = *cos_sin_ptr++; s1 = *cos_sin_ptr++; tempi = *spec_data1++; tempr = *spec_data2--; *out_ptr1 = ixheaac_sub32(ixheaac_mult32x16in32(tempi, c1), ixheaac_mult32x16in32(tempr, s1)); *out_ptr1 = ixheaac_shl32(*out_ptr1, neg_expo); out_ptr1--; *out_ptr1 = ixheaac_mac32x16in32(ixheaac_mult32x16in32(tempr, c1), tempi, s1); *out_ptr1 = ixheaac_shl32(*out_ptr1, neg_expo); out_ptr1--; } } VOID ixheaacd_pretwiddle_compute_dec( WORD32 *spec_data1, WORD32 *spec_data2, WORD32 *out_ptr, ia_aac_dec_imdct_tables_struct *ptr_imdct_tables, WORD npoints4, WORD32 neg_expo) { WORD32 i; WORD32 tempr, tempi; WORD32 tempr1, tempi1; WORD32 npoints2 = npoints4 * 2; WORD32 *out_ptr1 = out_ptr + (npoints2 << 1) - 1; const WORD16 *cos_sin_ptr = ptr_imdct_tables->cosine_array_2048_256; WORD16 cos = 0, cos1 = 0, sin = 0, sin1 = 0; if (neg_expo < 0) { neg_expo = -neg_expo; if (npoints4 == 256) { cos = *cos_sin_ptr++; sin = *cos_sin_ptr++; } else if (npoints4 == 32) { cos = *cos_sin_ptr++; sin = *cos_sin_ptr; cos_sin_ptr += 15; } tempr = *spec_data1++; tempi = *spec_data2--; *out_ptr = ixheaac_mac32x16in32(ixheaac_mult32x16in32(tempr, cos), tempi, sin); *out_ptr = ixheaac_shl32(*out_ptr, neg_expo); out_ptr++; *out_ptr = ixheaac_sub32(ixheaac_mult32x16in32(tempi, cos), ixheaac_mult32x16in32(tempr, sin)); *out_ptr = ixheaac_shl32(*out_ptr, neg_expo); out_ptr++; for (i = 0; i < npoints4 - 1; i++) { if (npoints4 == 256) { sin = *cos_sin_ptr++; cos = *cos_sin_ptr++; } else if (npoints4 == 32) { sin = *cos_sin_ptr++; cos = *cos_sin_ptr; cos_sin_ptr += 15; } tempi1 = *spec_data1++; tempr = *spec_data1++; tempr1 = *spec_data2--; tempi = *spec_data2--; *out_ptr1 = ixheaac_sub32(ixheaac_mult32x16in32(tempi1, cos), ixheaac_mult32x16in32(tempr1, sin)); *out_ptr1 = ixheaac_shl32(*out_ptr1, neg_expo); out_ptr1--; *out_ptr1 = ixheaac_mac32x16in32(ixheaac_mult32x16in32(tempr1, cos), tempi1, sin); *out_ptr1 = ixheaac_shl32(*out_ptr1, neg_expo); out_ptr1--; *out_ptr = ixheaac_mac32x16in32(ixheaac_mult32x16in32(tempr, sin), tempi, cos); *out_ptr = ixheaac_shl32(*out_ptr, neg_expo); out_ptr++; *out_ptr = ixheaac_sub32(ixheaac_mult32x16in32(tempi, sin), ixheaac_mult32x16in32(tempr, cos)); *out_ptr = ixheaac_shl32(*out_ptr, neg_expo); out_ptr++; } cos1 = *cos_sin_ptr++; sin1 = *cos_sin_ptr; tempr1 = *spec_data2; tempi1 = *spec_data1; *out_ptr1 = ixheaac_sub32(ixheaac_mult32x16in32(tempi1, cos1), ixheaac_mult32x16in32(tempr1, sin1)); *out_ptr1 = ixheaac_shl32(*out_ptr1, neg_expo); out_ptr1--; *out_ptr1 = ixheaac_mac32x16in32(ixheaac_mult32x16in32(tempr1, cos1), tempi1, sin1); *out_ptr1 = ixheaac_shl32(*out_ptr1, neg_expo); out_ptr1--; } else { if (npoints4 == 256) { cos = *cos_sin_ptr++; sin = *cos_sin_ptr++; } else if (npoints4 == 32) { cos = *cos_sin_ptr++; sin = *cos_sin_ptr; cos_sin_ptr += 15; } tempr = *spec_data1++; tempi = *spec_data2--; *out_ptr = ixheaac_mac32x16in32(ixheaac_mult32x16in32(tempr, cos), tempi, sin); *out_ptr = ixheaac_shr32(*out_ptr, neg_expo); out_ptr++; *out_ptr = ixheaac_sub32(ixheaac_mult32x16in32(tempi, cos), ixheaac_mult32x16in32(tempr, sin)); *out_ptr = ixheaac_shr32(*out_ptr, neg_expo); out_ptr++; for (i = 0; i < npoints4 - 1; i++) { if (npoints4 == 256) { sin = *cos_sin_ptr++; cos = *cos_sin_ptr++; } else if (npoints4 == 32) { sin = *cos_sin_ptr++; cos = *cos_sin_ptr; cos_sin_ptr += 15; } tempi1 = *spec_data1++; tempr = *spec_data1++; tempr1 = *spec_data2--; tempi = *spec_data2--; *out_ptr1 = ixheaac_sub32(ixheaac_mult32x16in32(tempi1, cos), ixheaac_mult32x16in32(tempr1, sin)); *out_ptr1 = ixheaac_shr32(*out_ptr1, neg_expo); out_ptr1--; *out_ptr1 = ixheaac_mac32x16in32(ixheaac_mult32x16in32(tempr1, cos), tempi1, sin); *out_ptr1 = ixheaac_shr32(*out_ptr1, neg_expo); out_ptr1--; *out_ptr = ixheaac_mac32x16in32(ixheaac_mult32x16in32(tempr, sin), tempi, cos); *out_ptr = ixheaac_shr32(*out_ptr, neg_expo); out_ptr++; *out_ptr = ixheaac_sub32(ixheaac_mult32x16in32(tempi, sin), ixheaac_mult32x16in32(tempr, cos)); *out_ptr = ixheaac_shr32(*out_ptr, neg_expo); out_ptr++; } cos1 = *cos_sin_ptr++; sin1 = *cos_sin_ptr; tempr1 = *spec_data2; tempi1 = *spec_data1; *out_ptr1 = ixheaac_sub32(ixheaac_mult32x16in32(tempi1, cos1), ixheaac_mult32x16in32(tempr1, sin1)); *out_ptr1 = ixheaac_shr32(*out_ptr1, neg_expo); out_ptr1--; *out_ptr1 = ixheaac_mac32x16in32(ixheaac_mult32x16in32(tempr1, cos1), tempi1, sin1); *out_ptr1 = ixheaac_shr32(*out_ptr1, neg_expo); out_ptr1--; } } VOID ixheaacd_post_twiddle_dec(WORD32 out_ptr[], WORD32 spec_data[], ia_aac_dec_imdct_tables_struct *ptr_imdct_tables, WORD npoints) { WORD i; WORD16 cos, cos1, sin, sin1; WORD32 *spec_data1 = spec_data + npoints - 1; WORD32 *out_ptr1 = out_ptr + npoints - 1; WORD16 adjust = 50, adjust1 = -50; const WORD16 *cos_sin_ptr = ptr_imdct_tables->cosine_array_2048_256; if (npoints == 1024) { WORD32 tempr, tempi, outi, outr, temp1, temp2; tempr = *spec_data++; tempi = *spec_data++; cos = *cos_sin_ptr; cos_sin_ptr++; sin = *cos_sin_ptr; cos_sin_ptr++; outi = ixheaac_sub32(ixheaac_mult32x16in32(tempr, sin), ixheaac_mult32x16in32(tempi, cos)); outr = ixheaac_mac32x16in32(ixheaac_mult32x16in32(tempr, cos), tempi, sin); temp1 = ixheaac_mult32x16in32(outi, adjust1); temp2 = ixheaac_mult32x16in32(outr, adjust); outr = outr + temp1; outi = outi + temp2; *out_ptr1-- = outi; *out_ptr++ = outr; for (i = 0; i < (npoints / 2 - 2); i++) { sin = *cos_sin_ptr++; cos = *cos_sin_ptr++; tempi = *spec_data1--; tempr = *spec_data1--; outi = ixheaac_sub32(ixheaac_mult32x16in32(tempr, sin), ixheaac_mult32x16in32(tempi, cos)); outr = ixheaac_mac32x16in32(ixheaac_mult32x16in32(tempr, cos), tempi, sin); temp1 = ixheaac_mult32x16in32(outi, adjust1); temp2 = ixheaac_mult32x16in32(outr, adjust); outr = outr + temp1; outi = outi + temp2; *out_ptr++ = outi; *out_ptr1-- = outr; i++; tempr = *spec_data++; tempi = *spec_data++; outi = ixheaac_sub32(ixheaac_mult32x16in32(tempr, cos), ixheaac_mult32x16in32(tempi, sin)); outr = ixheaac_mac32x16in32(ixheaac_mult32x16in32(tempr, sin), tempi, cos); temp1 = ixheaac_mult32x16in32(outi, adjust1); temp2 = ixheaac_mult32x16in32(outr, adjust); outr = outr + temp1; outi = outi + temp2; *out_ptr1-- = outi; *out_ptr++ = outr; } cos1 = *cos_sin_ptr++; sin1 = *cos_sin_ptr; tempi = *spec_data1--; tempr = *spec_data1--; outi = ixheaac_sub32(ixheaac_mult32x16in32(tempr, sin1), ixheaac_mult32x16in32(tempi, cos1)); outr = ixheaac_mac32x16in32(ixheaac_mult32x16in32(tempr, cos1), tempi, sin1); temp1 = ixheaac_mult32x16in32(outi, adjust1); temp2 = ixheaac_mult32x16in32(outr, adjust); outr = outr + temp1; outi = outi + temp2; *out_ptr++ = outi; *out_ptr1-- = outr; } else if (npoints == 128) { WORD32 tempr, tempi, outi, outr, temp1, temp2; tempr = *spec_data++; tempi = *spec_data++; cos = *cos_sin_ptr++; sin = *cos_sin_ptr; cos_sin_ptr += 15; outi = ixheaac_sub32(ixheaac_mult32x16in32(tempr, sin), ixheaac_mult32x16in32(tempi, cos)); outr = ixheaac_mac32x16in32(ixheaac_mult32x16in32(tempr, cos), tempi, sin); temp1 = ixheaac_mult32x16in32(outi, -(201 << 1)); temp2 = ixheaac_mult32x16in32(outr, 201 << 1); outr = outr + temp1; outi = outi + temp2; *out_ptr1-- = outi; *out_ptr++ = outr; for (i = 0; i < (npoints / 2 - 2); i++) { sin = *cos_sin_ptr++; cos = *cos_sin_ptr; cos_sin_ptr += 15; tempi = *spec_data1--; tempr = *spec_data1--; outi = ixheaac_sub32(ixheaac_mult32x16in32(tempr, sin), ixheaac_mult32x16in32(tempi, cos)); outr = ixheaac_mac32x16in32(ixheaac_mult32x16in32(tempr, cos), tempi, sin); temp1 = ixheaac_mult32x16in32(outi, -(201 << 1)); temp2 = ixheaac_mult32x16in32(outr, 201 << 1); outr = outr + temp1; outi = outi + temp2; *out_ptr++ = outi; *out_ptr1-- = outr; i++; tempr = *spec_data++; tempi = *spec_data++; outi = ixheaac_sub32(ixheaac_mult32x16in32(tempr, cos), ixheaac_mult32x16in32(tempi, sin)); outr = ixheaac_mac32x16in32(ixheaac_mult32x16in32(tempr, sin), tempi, cos); temp1 = ixheaac_mult32x16in32(outi, -(201 << 1)); temp2 = ixheaac_mult32x16in32(outr, 201 << 1); outr = outr + temp1; outi = outi + temp2; *out_ptr1-- = outi; *out_ptr++ = outr; } cos1 = *cos_sin_ptr++; sin1 = *cos_sin_ptr; tempi = *spec_data1--; tempr = *spec_data1--; outi = ixheaac_sub32(ixheaac_mult32x16in32(tempr, sin1), ixheaac_mult32x16in32(tempi, cos1)); outr = ixheaac_mac32x16in32(ixheaac_mult32x16in32(tempr, cos1), tempi, sin1); temp1 = ixheaac_mult32x16in32(outi, -(201 << 1)); temp2 = ixheaac_mult32x16in32(outr, 201 << 1); outr = outr + temp1; outi = outi + temp2; *out_ptr++ = outi; *out_ptr1-- = outr; } } VOID ixheaacd_post_twid_overlap_add_dec( WORD32 pcm_out[], WORD32 spec_data[], ia_aac_dec_imdct_tables_struct *ptr_imdct_tables, WORD npoints, WORD32 *ptr_overlap_buf, WORD16 q_shift, const WORD16 *window, WORD16 ch_fac) { WORD i; WORD16 cos, cos1, sin, sin1; WORD32 size = npoints / 2; WORD32 *pcmout1 = pcm_out + (ch_fac * size); const WORD16 *cos_sin_ptr = ptr_imdct_tables->cosine_array_2048_256; pcm_out = pcmout1 - ch_fac; spec_data += size; if (q_shift > 0) { WORD32 tempr, tempi, outr, outi, win1, accu, temp1, temp2; WORD16 adjust, adjust1; WORD32 overlap_data; tempr = *(spec_data - size); tempi = *(spec_data - size + 1); adjust = 50; adjust1 = -50; cos = *cos_sin_ptr++; sin = *cos_sin_ptr++; outi = ixheaac_sub32(ixheaac_mult32x16in32(tempr, sin), ixheaac_mult32x16in32(tempi, cos)); outr = ixheaac_add32(ixheaac_mult32x16in32(tempr, cos), ixheaac_mult32x16in32(tempi, sin)); overlap_data = *ptr_overlap_buf; temp1 = ixheaac_mult32x16in32(outi, adjust1); temp2 = ixheaac_mult32x16in32(outr, adjust); outr = outr + temp1; outi = outi + temp2; *ptr_overlap_buf++ = ixheaac_shr32_sat(outr, 16 - q_shift); win1 = *((WORD32 *)window + size - 1); accu = ixheaac_sub32_sat( ixheaac_shl32_sat(ixheaacd_mult32x16lin32(outi, win1), q_shift), ixheaacd_mult32x16lin32_sat(overlap_data, (WORD16)(win1 >> 16))); *pcm_out = accu; pcm_out -= ch_fac; accu = ixheaac_sub32_sat( ixheaac_shl32_sat( ixheaac_mult32x16hin32(ixheaac_negate32_sat(outi), win1), q_shift), ixheaacd_mult32x16lin32_sat(overlap_data, (WORD16)(win1))); *pcmout1 = accu; pcmout1 += ch_fac; for (i = size - 2; i != 0;) { sin = *cos_sin_ptr++; cos = *cos_sin_ptr++; tempr = *(spec_data + i); tempi = *(spec_data + i + 1); outr = ixheaac_add32(ixheaac_mult32x16in32(tempr, cos), ixheaac_mult32x16in32(tempi, sin)); outi = ixheaac_sub32(ixheaac_mult32x16in32(tempr, sin), ixheaac_mult32x16in32(tempi, cos)); temp1 = ixheaac_mult32x16in32(outi, adjust1); temp2 = ixheaac_mult32x16in32(outr, adjust); outr = outr + temp1; outi = outi + temp2; overlap_data = *ptr_overlap_buf; *ptr_overlap_buf++ = ixheaac_shr32_sat(outi, 16 - q_shift); win1 = *((WORD32 *)window + i); accu = ixheaac_sub32_sat( ixheaac_shl32_sat(ixheaacd_mult32x16lin32(outr, win1), q_shift), ixheaacd_mult32x16lin32_sat(overlap_data, (WORD16)(win1 >> 16))); *pcm_out = accu; pcm_out -= ch_fac; accu = ixheaac_sub32_sat( ixheaac_shl32_sat( ixheaac_mult32x16hin32(ixheaac_negate32_sat(outr), win1), q_shift), ixheaacd_mult32x16lin32_sat(overlap_data, (WORD16)win1)); *pcmout1 = accu; pcmout1 += ch_fac; tempr = *(spec_data - i); tempi = *(spec_data - i + 1); i -= 2; outi = ixheaac_sub32(ixheaac_mult32x16in32(tempr, cos), ixheaac_mult32x16in32(tempi, sin)); outr = ixheaac_add32(ixheaac_mult32x16in32(tempr, sin), ixheaac_mult32x16in32(tempi, cos)); overlap_data = *ptr_overlap_buf; temp1 = ixheaac_mult32x16in32(outi, adjust1); temp2 = ixheaac_mult32x16in32(outr, adjust); outr = outr + temp1; outi = outi + temp2; *ptr_overlap_buf++ = ixheaac_shr32_sat(outr, 16 - q_shift); win1 = *((WORD32 *)window + i + 1); accu = ixheaac_sub32_sat( ixheaac_shl32_sat(ixheaacd_mult32x16lin32(outi, win1), q_shift), ixheaacd_mult32x16lin32_sat(overlap_data, (WORD16)(win1 >> 16))); *pcm_out = accu; pcm_out -= ch_fac; accu = ixheaac_sub32_sat( ixheaac_shl32_sat( ixheaac_mult32x16hin32(ixheaac_negate32_sat(outi), win1), q_shift), ixheaacd_mult32x16lin32_sat(overlap_data, (WORD16)(win1))); *pcmout1 = accu; pcmout1 += ch_fac; } cos1 = *cos_sin_ptr++; sin1 = *cos_sin_ptr; tempr = *(spec_data + i); tempi = *(spec_data + i + 1); outr = ixheaac_add32(ixheaac_mult32x16in32(tempr, cos1), ixheaac_mult32x16in32(tempi, sin1)); outi = ixheaac_sub32(ixheaac_mult32x16in32(tempr, sin1), ixheaac_mult32x16in32(tempi, cos1)); temp1 = ixheaac_mult32x16in32(outi, adjust1); temp2 = ixheaac_mult32x16in32(outr, adjust); outr = outr + temp1; outi = outi + temp2; overlap_data = *ptr_overlap_buf; *ptr_overlap_buf++ = ixheaac_shr32_sat(outi, 16 - q_shift); win1 = *((WORD32 *)window + i); accu = ixheaac_sub32_sat( ixheaac_shl32_sat(ixheaacd_mult32x16lin32(outr, win1), q_shift), ixheaacd_mult32x16lin32_sat(overlap_data, (WORD16)(win1 >> 16))); *pcm_out = accu; pcm_out -= ch_fac; accu = ixheaac_sub32_sat( ixheaac_shl32_sat( ixheaac_mult32x16hin32(ixheaac_negate32_sat(outr), win1), q_shift), ixheaacd_mult32x16lin32_sat(overlap_data, (WORD16)win1)); *pcmout1 = accu; pcmout1 += ch_fac; } else { q_shift = -q_shift; { WORD32 tempr, tempi, temp1, temp2, outr, outi, win1, accu; WORD16 adjust, adjust1; WORD16 overlap_data; tempr = *(spec_data - size); tempi = *(spec_data - size + 1); adjust = 50; adjust1 = -50; cos = *cos_sin_ptr++; sin = *cos_sin_ptr++; outi = ixheaac_sub32(ixheaac_mult32x16in32(tempr, sin), ixheaac_mult32x16in32(tempi, cos)); outr = ixheaac_add32(ixheaac_mult32x16in32(tempr, cos), ixheaac_mult32x16in32(tempi, sin)); overlap_data = *ptr_overlap_buf; temp1 = ixheaac_mult32x16in32(outi, adjust1); temp2 = ixheaac_mult32x16in32(outr, adjust); outr = outr + temp1; outi = outi + temp2; *ptr_overlap_buf++ = ixheaac_shr32_sat(outr, 16 + q_shift); win1 = *((WORD32 *)window + size - 1); accu = ixheaac_sub32_sat( ixheaac_shr32(ixheaacd_mult32x16lin32(outi, win1), q_shift), ixheaacd_mult32x16lin32_sat(overlap_data, (WORD16)(win1 >> 16))); *pcm_out = accu; pcm_out -= ch_fac; accu = ixheaac_sub32_sat( ixheaac_shr32( ixheaac_mult32x16hin32(ixheaac_negate32_sat(outi), win1), q_shift), ixheaacd_mult32x16lin32_sat(overlap_data, (WORD16)(win1))); *pcmout1 = accu; pcmout1 += ch_fac; for (i = size - 2; i != 0;) { sin = *cos_sin_ptr++; cos = *cos_sin_ptr++; tempr = *(spec_data + i); tempi = *(spec_data + i + 1); outr = ixheaac_add32(ixheaac_mult32x16in32(tempr, cos), ixheaac_mult32x16in32(tempi, sin)); outi = ixheaac_sub32(ixheaac_mult32x16in32(tempr, sin), ixheaac_mult32x16in32(tempi, cos)); overlap_data = *ptr_overlap_buf; temp1 = ixheaac_mult32x16in32(outi, adjust1); temp2 = ixheaac_mult32x16in32(outr, adjust); outr = outr + temp1; outi = outi + temp2; *ptr_overlap_buf++ = ixheaac_shr32_sat(outi, 16 + q_shift); win1 = *((WORD32 *)window + i); accu = ixheaac_sub32_sat( ixheaac_shr32(ixheaacd_mult32x16lin32(outr, win1), q_shift), ixheaacd_mult32x16lin32_sat(overlap_data, (WORD16)(win1 >> 16))); *pcm_out = accu; pcm_out -= ch_fac; accu = ixheaac_sub32_sat( ixheaac_shr32( ixheaac_mult32x16hin32(ixheaac_negate32_sat(outr), win1), q_shift), ixheaacd_mult32x16lin32_sat(overlap_data, (WORD16)win1)); *pcmout1 = accu; pcmout1 += ch_fac; tempr = *(spec_data - i); tempi = *(spec_data - i + 1); i -= 2; outi = ixheaac_sub32(ixheaac_mult32x16in32(tempr, cos), ixheaac_mult32x16in32(tempi, sin)); outr = ixheaac_add32(ixheaac_mult32x16in32(tempr, sin), ixheaac_mult32x16in32(tempi, cos)); overlap_data = *ptr_overlap_buf; temp1 = ixheaac_mult32x16in32(outi, adjust1); temp2 = ixheaac_mult32x16in32(outr, adjust); outr = outr + temp1; outi = outi + temp2; *ptr_overlap_buf++ = ixheaac_shr32_sat(outr, 16 + q_shift); win1 = *((WORD32 *)window + i + 1); accu = ixheaac_sub32_sat( ixheaac_shr32(ixheaacd_mult32x16lin32(outi, win1), q_shift), ixheaacd_mult32x16lin32_sat(overlap_data, (WORD16)(win1 >> 16))); *pcm_out = accu; pcm_out -= ch_fac; accu = ixheaac_sub32_sat( ixheaac_shr32( ixheaac_mult32x16hin32(ixheaac_negate32_sat(outi), win1), q_shift), ixheaacd_mult32x16lin32_sat(overlap_data, (WORD16)(win1))); *pcmout1 = accu; pcmout1 += ch_fac; } cos1 = *cos_sin_ptr++; sin1 = *cos_sin_ptr++; tempr = *(spec_data + i); tempi = *(spec_data + i + 1); outr = ixheaac_add32(ixheaac_mult32x16in32(tempr, cos1), ixheaac_mult32x16in32(tempi, sin1)); outi = ixheaac_sub32(ixheaac_mult32x16in32(tempr, sin1), ixheaac_mult32x16in32(tempi, cos1)); overlap_data = *ptr_overlap_buf; temp1 = ixheaac_mult32x16in32(outi, adjust1); temp2 = ixheaac_mult32x16in32(outr, adjust); outr = outr + temp1; outi = outi + temp2; *ptr_overlap_buf++ = ixheaac_shr32_sat(outi, 16 + q_shift); win1 = *((WORD32 *)window + i); accu = ixheaac_sub32_sat( ixheaac_shr32(ixheaacd_mult32x16lin32(outr, win1), q_shift), ixheaacd_mult32x16lin32_sat(overlap_data, (WORD16)(win1 >> 16))); *pcm_out = accu; pcm_out -= ch_fac; accu = ixheaac_sub32_sat( ixheaac_shr32( ixheaac_mult32x16hin32(ixheaac_negate32_sat(outr), win1), q_shift), ixheaacd_mult32x16lin32_sat(overlap_data, (WORD16)win1)); *pcmout1 = accu; pcmout1 += ch_fac; } } } VOID ixheaacd_imdct_using_fft_dec( ia_aac_dec_imdct_tables_struct *ptr_imdct_tables, WORD32 npoints, WORD32 *ptr_x, WORD32 *ptr_y) { WORD32 i, j, k, k1, n_stages; WORD32 x0r, x0i, x1r, x1i, x2r, x2i, x3r, x3i, x4r, x4i, x5r, x5i, x6r, x6i, x7r, x7i; WORD32 del, nodespacing, in_loop_cnt, tmp, twiddle_val, *ptr_tmp; const WORD32 *ptr_twiddle; WORD8 *ptr_dig_rev_table; n_stages = ixheaac_norm32(npoints); n_stages = (30 - n_stages) / 3; ptr_tmp = ptr_y; ptr_twiddle = ptr_imdct_tables->fft_twiddle; ptr_dig_rev_table = ((npoints << 1) == 1024) ? ptr_imdct_tables->dig_rev_table8_long : ptr_imdct_tables->dig_rev_table8_short; for (i = npoints; i != 0; i -= 8) { WORD32 *data = ptr_x; data = data + (*ptr_dig_rev_table++ << 1); x0r = *data; x0i = *(data + 1); data += (npoints >> 1); x2r = *data; x2i = *(data + 1); data += (npoints >> 1); x4r = *data; x4i = *(data + 1); data += (npoints >> 1); x6r = *data; x6i = *(data + 1); data -= 5 * (npoints >> 2); x0r = x0r + x4r; x0i = x0i + x4i; x4r = x0r - (x4r << 1); x4i = x0i - (x4i << 1); x2r = x2r + x6r; x2i = x2i + x6i; x6r = x2r - (x6r << 1); x6i = x2i - (x6i << 1); x0r = x0r + x2r; x0i = x0i + x2i; x2r = x0r - (x2r << 1); x2i = x0i - (x2i << 1); x4r = x4r + x6i; x4i = x4i - x6r; tmp = x6r; x6r = x4r - (x6i << 1); x6i = x4i + (tmp << 1); x1r = *data; x1i = *(data + 1); data += (npoints >> 1); x3r = *data; x3i = *(data + 1); data += (npoints >> 1); x5r = *data; x5i = *(data + 1); data += (npoints >> 1); x7r = *data; x7i = *(data + 1); data -= 7 * (npoints >> 2); x1r = x1r + x5r; x1i = x1i + x5i; x5r = x1r - (x5r << 1); x5i = x1i - (x5i << 1); x3r = x3r + x7r; x3i = x3i + x7i; x7r = x3r - (x7r << 1); x7i = x3i - (x7i << 1); x1r = x1r + x3r; x1i = x1i + x3i; x3r = x1r - (x3r << 1); x3i = x1i - (x3i << 1); x5r = x5r + x5i; x5i = x5r - (x5i << 1); x7r = x7r + x7i; x7i = x7r - (x7i << 1); x7i = x5r - x7i; x5r = x7i - (x5r << 1); x5i = x7r - x5i; x7r = x5i - (x7r << 1); x7i = x7i << 1; x5r = x5r << 1; x5i = x5i << 1; x7r = x7r << 1; x0r = x0r + x1r; x0i = x0i + x1i; x1r = x0r - (x1r << 1); x1i = x0i - (x1i << 1); x2r = x2r + x3i; tmp = x2r - (x3i << 1); x2i = x2i - x3r; x3i = x2i + (x3r << 1); *ptr_tmp = x0r; *(ptr_tmp + 1) = x0i; ptr_tmp += 4; *ptr_tmp = x2r; *(ptr_tmp + 1) = x2i; ptr_tmp += 4; *ptr_tmp = x1r; *(ptr_tmp + 1) = x1i; ptr_tmp += 4; *ptr_tmp = tmp; *(ptr_tmp + 1) = x3i; ptr_tmp -= 10; tmp = 0x5A82; x7i = x4r + (ixheaacd_mult32x16lin32(x7i, tmp)); x4r = x7i - (x4r << 1); x7r = x4i + (ixheaacd_mult32x16lin32(x7r, tmp)); x4i = x7r - (x4i << 1); x5i = x6r + (ixheaacd_mult32x16lin32(x5i, tmp)); x6r = x5i - (x6r << 1); x5r = x6i + (ixheaacd_mult32x16lin32(x5r, tmp)); x6i = x5r - (x6i << 1); *ptr_tmp = x7i; *(ptr_tmp + 1) = x7r; ptr_tmp += 4; *ptr_tmp = x5i; *(ptr_tmp + 1) = x5r; ptr_tmp += 4; *ptr_tmp = -x4r; *(ptr_tmp + 1) = -x4i; ptr_tmp += 4; *ptr_tmp = -x6r; *(ptr_tmp + 1) = -x6i; ptr_tmp += 2; } del = 8; nodespacing = 64; in_loop_cnt = npoints >> 6; for (k1 = n_stages - 2; k1 > 0; k1--) { WORD32 *data = ptr_y; const WORD32 *twiddles; for (i = 0; i != npoints; i += 8 * del) { data = ptr_y + (i << 1); x0r = *data; x0i = *(data + 1); data += (del << 2); x2r = *data; x2i = *(data + 1); data += (del << 2); x4r = *data; x4i = *(data + 1); data += (del << 2); x6r = *data; x6i = *(data + 1); data -= 5 * (del << 1); x0r = x0r + x4r; x0i = x0i + x4i; x4r = x0r - (x4r << 1); x4i = x0i - (x4i << 1); x2r = x2r + x6r; x2i = x2i + x6i; x6r = x2r - (x6r << 1); x6i = x2i - (x6i << 1); x0r = x0r + x2r; x0i = x0i + x2i; x2r = x0r - (x2r << 1); x2i = x0i - (x2i << 1); x4r = x4r + x6i; x4i = x4i - x6r; tmp = x6r; x6r = x4r - (x6i << 1); x6i = x4i + (tmp << 1); x1r = *data; x1i = *(data + 1); data += (del << 2); x3r = *data; x3i = *(data + 1); data += (del << 2); x5r = *data; x5i = *(data + 1); data += (del << 2); x7r = *data; x7i = *(data + 1); data -= 7 * (del << 1); x1r = x1r + x5r; x1i = x1i + x5i; x5r = x1r - (x5r << 1); x5i = x1i - (x5i << 1); x3r = x3r + x7r; x3i = x3i + x7i; x7r = x3r - (x7r << 1); x7i = x3i - (x7i << 1); x1r = x1r + x3r; x1i = x1i + x3i; x3r = x1r - (x3r << 1); x3i = x1i - (x3i << 1); x5r = x5r + x5i; x5i = x5r - (x5i << 1); x7r = x7r + x7i; x7i = x7r - (x7i << 1); x7i = x5r - x7i; x5r = x7i - (x5r << 1); x5i = x7r - x5i; x7r = x5i - (x7r << 1); x7i = x7i << 1; x5r = x5r << 1; x5i = x5i << 1; x7r = x7r << 1; x0r = x0r + x1r; x0i = x0i + x1i; x1r = x0r - (x1r << 1); x1i = x0i - (x1i << 1); x2r = x2r + x3i; tmp = x2r - (x3i << 1); x2i = x2i - x3r; x3i = x2i + (x3r << 1); *data = x0r; *(data + 1) = x0i; data += (del << 2); *data = x2r; *(data + 1) = x2i; data += (del << 2); *data = x1r; *(data + 1) = x1i; data += (del << 2); *data = tmp; *(data + 1) = x3i; data -= 5 * (del << 1); tmp = 0x5A82; x7i = x4r + (ixheaacd_mult32x16lin32(x7i, tmp)); x4r = x7i - (x4r << 1); x7r = x4i + (ixheaacd_mult32x16lin32(x7r, tmp)); x4i = x7r - (x4i << 1); x5i = x6r + (ixheaacd_mult32x16lin32(x5i, tmp)); x6r = x5i - (x6r << 1); x5r = x6i + (ixheaacd_mult32x16lin32(x5r, tmp)); x6i = x5r - (x6i << 1); *data = x7i; *(data + 1) = x7r; data += (del << 2); *data = x5i; *(data + 1) = x5r; data += (del << 2); *data = -x4r; *(data + 1) = -x4i; data += (del << 2); *data = -x6r; *(data + 1) = -x6i; data -= 7 * (del << 1); } twiddles = ptr_twiddle; data = ptr_y; for (j = nodespacing; j < nodespacing * del; j += nodespacing) { data = data + 2; for (k = in_loop_cnt; k != 0; k--) { data += (del << 2); x2r = *data; x2i = *(data + 1); data += (del << 2); x4r = *data; x4i = *(data + 1); data += (del << 2); x6r = *data; x6i = *(data + 1); data -= 6 * (del << 1); twiddles += (j >> 2); twiddle_val = *(twiddles); tmp = (ixheaacd_mult32x16lin32(x2r, twiddle_val) - ixheaac_mult32x16hin32(x2i, twiddle_val)); x2i = (ixheaacd_mac32x16lin32( ixheaac_mult32x16hin32(x2r, twiddle_val), x2i, twiddle_val)) << 1; x2r = tmp << 1; twiddles += (j >> 2); twiddle_val = *(twiddles); tmp = (ixheaacd_mult32x16lin32(x4r, twiddle_val) - ixheaac_mult32x16hin32(x4i, twiddle_val)); x4i = (ixheaacd_mac32x16lin32( ixheaac_mult32x16hin32(x4r, twiddle_val), x4i, twiddle_val)) << 1; x4r = tmp << 1; twiddles += (j >> 2); twiddle_val = *(twiddles); tmp = (ixheaacd_mult32x16lin32(x6r, twiddle_val) - ixheaac_mult32x16hin32(x6i, twiddle_val)); x6i = (ixheaacd_mac32x16lin32( ixheaac_mult32x16hin32(x6r, twiddle_val), x6i, twiddle_val)) << 1; x6r = tmp << 1; x0r = *data; x0i = *(data + 1); data += (del << 1); x0r = x0r + x4r; x0i = x0i + x4i; x4r = x0r - (x4r << 1); x4i = x0i - (x4i << 1); x2r = x2r + x6r; x2i = x2i + x6i; x6r = x2r - (x6r << 1); x6i = x2i - (x6i << 1); x0r = x0r + x2r; x0i = x0i + x2i; x2r = x0r - (x2r << 1); x2i = x0i - (x2i << 1); x4r = x4r + x6i; x4i = x4i - x6r; tmp = x6r; x6r = x4r - (x6i << 1); x6i = x4i + (tmp << 1); x1r = *data; x1i = *(data + 1); data += (del << 2); twiddles -= 5 * (j >> 3); twiddle_val = *(twiddles); tmp = (ixheaacd_mult32x16lin32(x1r, twiddle_val) - ixheaac_mult32x16hin32(x1i, twiddle_val)); x1i = (ixheaacd_mac32x16lin32( ixheaac_mult32x16hin32(x1r, twiddle_val), x1i, twiddle_val)) << 1; x1r = tmp << 1; x3r = *data; x3i = *(data + 1); data += (del << 2); twiddles += (j >> 2); twiddle_val = *(twiddles); tmp = (ixheaacd_mult32x16lin32(x3r, twiddle_val) - ixheaac_mult32x16hin32(x3i, twiddle_val)); x3i = (ixheaacd_mac32x16lin32( ixheaac_mult32x16hin32(x3r, twiddle_val), x3i, twiddle_val)); x3r = tmp; x5r = *data; x5i = *(data + 1); data += (del << 2); twiddles += (j >> 2); twiddle_val = *(twiddles); tmp = (ixheaacd_mult32x16lin32(x5r, twiddle_val) - ixheaac_mult32x16hin32(x5i, twiddle_val)); x5i = (ixheaacd_mac32x16lin32( ixheaac_mult32x16hin32(x5r, twiddle_val), x5i, twiddle_val)); x5r = tmp; x7r = *data; x7i = *(data + 1); data -= 7 * (del << 1); twiddles += (j >> 2); twiddle_val = *(twiddles); twiddles -= 7 * (j >> 3); tmp = (ixheaacd_mult32x16lin32(x7r, twiddle_val) - ixheaac_mult32x16hin32(x7i, twiddle_val)); x7i = (ixheaacd_mac32x16lin32( ixheaac_mult32x16hin32(x7r, twiddle_val), x7i, twiddle_val)); x7r = tmp; x1r = x1r + (x5r << 1); x1i = x1i + (x5i << 1); x5r = x1r - (x5r << 2); x5i = x1i - (x5i << 2); x3r = x3r + x7r; x3i = x3i + x7i; x7r = x3r - (x7r << 1); x7i = x3i - (x7i << 1); x1r = x1r + (x3r << 1); x1i = x1i + (x3i << 1); x3r = x1r - (x3r << 2); x3i = x1i - (x3i << 2); x5r = x5r + x5i; x5i = x5r - (x5i << 1); x7r = x7r + x7i; x7i = x7r - (x7i << 1); x7i = x5r - (x7i << 1); x5r = x7i - (x5r << 1); x5i = (x7r << 1) - x5i; x7r = x5i - (x7r << 2); x7i = x7i << 1; x5r = x5r << 1; x5i = x5i << 1; x7r = x7r << 1; x0r = x0r + x1r; x0i = x0i + x1i; x1r = x0r - (x1r << 1); x1i = x0i - (x1i << 1); x2r = x2r + x3i; tmp = x2r - (x3i << 1); x2i = x2i - x3r; x3i = x2i + (x3r << 1); *data = x0r; *(data + 1) = x0i; data += (del << 2); *data = x2r; *(data + 1) = x2i; data += (del << 2); *data = x1r; *(data + 1) = x1i; data += (del << 2); *data = tmp; *(data + 1) = x3i; data -= 5 * (del << 1); tmp = 0x5A82; x7i = x4r + (ixheaacd_mult32x16lin32(x7i, tmp)); x4r = x7i - (x4r << 1); x7r = x4i + (ixheaacd_mult32x16lin32(x7r, tmp)); x4i = x7r - (x4i << 1); x5i = x6r + (ixheaacd_mult32x16lin32(x5i, tmp)); x6r = x5i - (x6r << 1); x5r = x6i + (ixheaacd_mult32x16lin32(x5r, tmp)); x6i = x5r - (x6i << 1); *data = x7i; *(data + 1) = x7r; data += (del << 2); *data = x5i; *(data + 1) = x5r; data += (del << 2); *data = -x4r; *(data + 1) = -x4i; data += (del << 2); *data = -x6r; *(data + 1) = -x6i; data -= 7 * (del << 1); data += (del << 4); } data -= npoints << 1; } nodespacing >>= 3; del <<= 3; in_loop_cnt >>= 3; } { WORD32 *data = ptr_y; const WORD32 *twiddles; twiddles = ptr_twiddle; data = ptr_y; data = data - 2; for (j = 0; j < nodespacing * del; j += nodespacing) { data = data + 2; { data += (del << 2); x2r = *data; x2i = *(data + 1); data += (del << 2); x4r = *data; x4i = *(data + 1); data += (del << 2); x6r = *data; x6i = *(data + 1); data -= 6 * (del << 1); twiddles += (j >> 2); twiddle_val = *(twiddles); tmp = (ixheaacd_mult32x16lin32(x2r, twiddle_val) - ixheaac_mult32x16hin32(x2i, twiddle_val)); x2i = (ixheaacd_mac32x16lin32( ixheaac_mult32x16hin32(x2r, twiddle_val), x2i, twiddle_val)) << 1; x2r = tmp << 1; twiddles += (j >> 2); twiddle_val = *(twiddles); tmp = (ixheaacd_mult32x16lin32(x4r, twiddle_val) - ixheaac_mult32x16hin32(x4i, twiddle_val)); x4i = (ixheaacd_mac32x16lin32( ixheaac_mult32x16hin32(x4r, twiddle_val), x4i, twiddle_val)) << 1; x4r = tmp << 1; twiddles += (j >> 2); twiddle_val = *(twiddles); tmp = (ixheaacd_mult32x16lin32(x6r, twiddle_val) - ixheaac_mult32x16hin32(x6i, twiddle_val)); x6i = (ixheaacd_mac32x16lin32( ixheaac_mult32x16hin32(x6r, twiddle_val), x6i, twiddle_val)) << 1; x6r = tmp << 1; x0r = *data; x0i = *(data + 1); data += (del << 1); x0r = x0r + x4r; x0i = x0i + x4i; x4r = x0r - (x4r << 1); x4i = x0i - (x4i << 1); x2r = x2r + x6r; x2i = x2i + x6i; x6r = x2r - (x6r << 1); x6i = x2i - (x6i << 1); x0r = x0r + x2r; x0i = x0i + x2i; x2r = x0r - (x2r << 1); x2i = x0i - (x2i << 1); x4r = x4r + x6i; x4i = x4i - x6r; tmp = x6r; x6r = x4r - (x6i << 1); x6i = x4i + (tmp << 1); x1r = *data; x1i = *(data + 1); data += (del << 2); twiddles -= 5 * (j >> 3); twiddle_val = *(twiddles); tmp = (ixheaacd_mult32x16lin32(x1r, twiddle_val) - ixheaac_mult32x16hin32(x1i, twiddle_val)); x1i = (ixheaacd_mac32x16lin32( ixheaac_mult32x16hin32(x1r, twiddle_val), x1i, twiddle_val)) << 1; x1r = tmp << 1; x3r = *data; x3i = *(data + 1); data += (del << 2); twiddles += (j >> 2); twiddle_val = *(twiddles); tmp = (ixheaacd_mult32x16lin32(x3r, twiddle_val) - ixheaac_mult32x16hin32(x3i, twiddle_val)); x3i = (ixheaacd_mac32x16lin32( ixheaac_mult32x16hin32(x3r, twiddle_val), x3i, twiddle_val)); x3r = tmp; x5r = *data; x5i = *(data + 1); data += (del << 2); twiddles += (j >> 2); twiddle_val = *(twiddles); tmp = (ixheaacd_mult32x16lin32(x5r, twiddle_val) - ixheaac_mult32x16hin32(x5i, twiddle_val)); x5i = (ixheaacd_mac32x16lin32( ixheaac_mult32x16hin32(x5r, twiddle_val), x5i, twiddle_val)); x5r = tmp; x7r = *data; x7i = *(data + 1); data -= 7 * (del << 1); twiddles += (j >> 2); twiddle_val = *(twiddles); twiddles -= 7 * (j >> 3); tmp = (ixheaacd_mult32x16lin32(x7r, twiddle_val) - ixheaac_mult32x16hin32(x7i, twiddle_val)); x7i = (ixheaacd_mac32x16lin32( ixheaac_mult32x16hin32(x7r, twiddle_val), x7i, twiddle_val)); x7r = tmp; x1r = x1r + (x5r << 1); x1i = x1i + (x5i << 1); x5r = x1r - (x5r << 2); x5i = x1i - (x5i << 2); x3r = x3r + x7r; x3i = x3i + x7i; x7r = x3r - (x7r << 1); x7i = x3i - (x7i << 1); x1r = x1r + (x3r << 1); x1i = x1i + (x3i << 1); x3r = x1r - (x3r << 2); x3i = x1i - (x3i << 2); x5r = x5r + x5i; x5i = x5r - (x5i << 1); x7r = x7r + x7i; x7i = x7r - (x7i << 1); x7i = x5r - (x7i << 1); x5r = x7i - (x5r << 1); x5i = (x7r << 1) - x5i; x7r = x5i - (x7r << 2); x7i = x7i << 1; x5r = x5r << 1; x5i = x5i << 1; x7r = x7r << 1; x0r = x0r + x1r; x0i = x0i + x1i; x1r = x0r - (x1r << 1); x1i = x0i - (x1i << 1); x2r = x2r + x3i; tmp = x2r - (x3i << 1); x2i = x2i - x3r; x3i = x2i + (x3r << 1); *data = x0r; *(data + 1) = x0i; data += (del << 2); *data = x2r; *(data + 1) = x2i; data += (del << 2); *data = x1r; *(data + 1) = x1i; data += (del << 2); *data = tmp; *(data + 1) = x3i; data -= 5 * (del << 1); tmp = 0x5A82; x7i = x4r + (ixheaacd_mult32x16lin32(x7i, tmp)); x4r = x7i - (x4r << 1); x7r = x4i + (ixheaacd_mult32x16lin32(x7r, tmp)); x4i = x7r - (x4i << 1); x5i = x6r + (ixheaacd_mult32x16lin32(x5i, tmp)); x6r = x5i - (x6r << 1); x5r = x6i + (ixheaacd_mult32x16lin32(x5r, tmp)); x6i = x5r - (x6i << 1); *data = x7i; *(data + 1) = x7r; data += (del << 2); *data = x5i; *(data + 1) = x5r; data += (del << 2); *data = -x4r; *(data + 1) = -x4i; data += (del << 2); *data = -x6r; *(data + 1) = -x6i; data -= 7 * (del << 1); data += (del << 4); } data -= npoints << 1; } nodespacing >>= 3; del <<= 3; in_loop_cnt >>= 3; } } VOID ixheaacd_inverse_transform_960( WORD32 spec_data[], WORD32 scratch[], ia_aac_dec_imdct_tables_struct *ptr_imdct_tables, WORD32 expo, WORD32 *imdct_scale) { WORD32 n; WORD32 Nd2; WORD16 const_mltfac; WORD32 neg_expo; WORD32 i; n = 120; Nd2 = n >> 1; neg_expo = 4; ixheaacd_pre_twiddle_120(spec_data, scratch, n, ptr_imdct_tables->cosine_array_240, neg_expo - expo); ixheaacd_fft_120(ptr_imdct_tables, Nd2, spec_data, scratch); neg_expo += 2; *imdct_scale = neg_expo + 1; ixheaacd_post_twiddle_120(spec_data, scratch, ptr_imdct_tables->cosine_array_240, n); const_mltfac = 17476; for (i = 0; i < 120; i++) { spec_data[i] = ixheaac_mult32x16in32_shl(spec_data[i], const_mltfac); } } WORD32 ixheaacd_inverse_transform( WORD32 spec_data[], WORD32 scratch[], ia_aac_dec_imdct_tables_struct *ptr_imdct_tables, WORD32 expo, WORD32 npoints) { (*ixheaacd_pretwiddle_compute)(spec_data, spec_data + npoints - 1, scratch, ptr_imdct_tables, (npoints >> 2), expo); (*ixheaacd_imdct_using_fft)(ptr_imdct_tables, npoints >> 1, scratch, spec_data); expo += 2; return expo; } VOID ixheaacd_mdct_960(WORD32 *inp, WORD32 *scratch, WORD32 *mdct_scale, WORD32 mdct_flag, ia_aac_dec_imdct_tables_struct *imdct_tables_ptr) { WORD32 expo, neg_expo = 0, k; WORD16 const_mltfac = 17476; expo = (*ixheaacd_calc_max_spectral_line)(inp, MDCT_LEN_960) - 1; ; memcpy(scratch, inp, sizeof(WORD32) * MDCT_LEN_960); neg_expo = 7 - expo; ixheaacd_pre_twiddle_960(inp, scratch, MDCT_LEN_960, imdct_tables_ptr->cosine_array_1920, neg_expo); ixheaacd_fft_960(inp, scratch, imdct_tables_ptr); ixheaacd_post_twiddle_960(inp, scratch, imdct_tables_ptr->cosine_array_1920, MDCT_LEN_960); if (0 == mdct_flag) { WORD32 *data = inp; for (k = MDCT_LEN_960 - 1; k >= 0; k -= 2) { *data = ixheaac_mult32x16in32_shl(*data, const_mltfac); data++; *data = ixheaac_mult32x16in32_shl(*data, const_mltfac); data++; } } *mdct_scale = neg_expo + 1 + 1 + 1; } VOID ixheaacd_mdct_480_ld(WORD32 *inp, WORD32 *scratch, WORD32 *mdct_scale, WORD32 mdct_flag, ia_aac_dec_imdct_tables_struct *imdct_tables_ptr, WORD32 object_type) { WORD32 expo, neg_expo = 0, k; WORD32 const_mltfac = 1145324612; expo = (*ixheaacd_calc_max_spectral_line)(inp, MDCT_LEN) - 1; ; memcpy(scratch, inp, sizeof(WORD32) * MDCT_LEN); neg_expo = 7 - expo; ixheaacd_pre_twiddle(inp, scratch, 480, imdct_tables_ptr->cosine_array_960, neg_expo); ixheaacd_fft_480_ld(inp, scratch, imdct_tables_ptr); if (object_type == AOT_ER_AAC_LD) { ixheaacd_post_twiddle_ld(inp, scratch, imdct_tables_ptr->cosine_array_960, 480); } else if (object_type == AOT_ER_AAC_ELD) { ixheaacd_post_twiddle_eld(inp + (480), scratch, imdct_tables_ptr->cosine_array_960, 480); } if (0 == mdct_flag) { WORD32 *data = inp; if (object_type != AOT_ER_AAC_ELD) { for (k = MDCT_LEN - 1; k >= 0; k -= 2) { *data = ixheaac_mult32_shl(*data, const_mltfac); data++; *data = ixheaac_mult32_shl(*data, const_mltfac); data++; } neg_expo += 1; } else { data = inp + 480; for (k = (MDCT_LEN << 1) - 1; k >= 0; k -= 2) { *data = ixheaac_mult32_shl(*data, const_mltfac); data++; *data = ixheaac_mult32_shl(*data, const_mltfac); data++; } neg_expo += 1; } } *mdct_scale = neg_expo + 3; } VOID ixheaacd_inverse_transform_512( WORD32 data[], WORD32 temp[], WORD32 *imdct_scale, WORD32 *cos_sin_ptr, ia_aac_dec_imdct_tables_struct *imdct_tables_ptr, WORD32 object_type) { WORD32 n; WORD32 npoints_2; WORD16 expo, neg_expo; n = 512; npoints_2 = n >> 1; expo = (*ixheaacd_calc_max_spectral_line)(data, n) - 1; memcpy(temp, data, sizeof(WORD32) * n); neg_expo = 7 - expo; ixheaacd_pre_twiddle(data, temp, n, cos_sin_ptr, neg_expo); (*ixheaacd_fft32x32_ld)(imdct_tables_ptr, npoints_2, data, temp); neg_expo = (*ixheaacd_neg_expo_inc)(neg_expo); *imdct_scale = neg_expo + 1; if (object_type == AOT_ER_AAC_ELD) ixheaacd_post_twiddle_eld((data + n), temp, cos_sin_ptr, n); else ixheaacd_post_twiddle_ld((data), temp, cos_sin_ptr, n); } VOID ixheaacd_fft_960(WORD32 *inp, WORD32 *op, ia_aac_dec_imdct_tables_struct *imdct_tables_ptr) { WORD32 i; WORD32 *buf1, *buf2; WORD16 *re_arr_tab_sml_480_ptr; (*ixheaacd_aac_ld_dec_rearrange_960)(inp, op, 480, imdct_tables_ptr->re_arr_tab_32); buf1 = op; buf2 = inp; for (i = 0; i < FFT15; i++) { ixheaacd_fft_32_points(imdct_tables_ptr->w_32, 32, buf1, buf2); buf1 += (FFT16X2 * 2); buf2 += (FFT16X2 * 2); } re_arr_tab_sml_480_ptr = imdct_tables_ptr->re_arr_tab_sml_480; buf1 = inp; for (i = 0; i < FFT16 * 2; i++) { ixheaacd_ld_dec_fft_15_opt(buf1, op, ixheaacd_fft5out, re_arr_tab_sml_480_ptr); buf1 += 2; re_arr_tab_sml_480_ptr += FFT15; } } VOID ixheaacd_fft_32_points(WORD16 *ptr_w, WORD32 npoints, WORD32* ptr_x, WORD32* ptr_y) { WORD32 i, j, l1, l2, h2, predj, tw_offset, stride, fft_jmp; WORD32 xt0_0, yt0_0, xt1_0, yt1_0, xt2_0, yt2_0; WORD32 xh0_0, xh1_0, xh20_0, xh21_0, xl0_0, xl1_0, xl20_0, xl21_0; WORD32 x_0, x_1, x_l1_0, x_l1_1, x_l2_0, x_l2_1; WORD32 x_h2_0, x_h2_1; WORD16 si10, si20, si30, co10, co20, co30; WORD16 *w; WORD32 *x, *x2, *x0; WORD32 *y0, *y1, *y2, *y3; WORD32 n0, j0; WORD32 radix; WORD32 norm; radix = 2; norm = 25; stride = 32; tw_offset = 0; fft_jmp = 192; while (stride > radix) { j = 0; fft_jmp >>= 2; h2 = stride >> 1; l1 = stride; l2 = stride + (stride >> 1); x = ptr_x; w = ptr_w + tw_offset; tw_offset += fft_jmp; for (i = 0; i < npoints; i += 4) { co10 = w[j + 1]; si10 = w[j + 0]; co20 = w[j + 3]; si20 = w[j + 2]; co30 = w[j + 5]; si30 = w[j + 4]; x_0 = x[0]; x_1 = x[1]; x_l1_0 = x[l1]; x_l1_1 = x[l1 + 1]; x_l2_0 = x[l2]; x_l2_1 = x[l2 + 1]; x_h2_0 = x[h2]; x_h2_1 = x[h2 + 1]; xh0_0 = ixheaac_add32_sat(x_0, x_l1_0); xh1_0 = ixheaac_add32_sat(x_1, x_l1_1); xl0_0 = ixheaac_sub32_sat(x_0, x_l1_0); xl1_0 = ixheaac_sub32_sat(x_1, x_l1_1); xh20_0 = ixheaac_add32_sat(x_h2_0, x_l2_0); xh21_0 = ixheaac_add32_sat(x_h2_1, x_l2_1); xl20_0 = ixheaac_sub32_sat(x_h2_0, x_l2_0); xl21_0 = ixheaac_sub32_sat(x_h2_1, x_l2_1); x0 = x; x2 = x0; j += 6; x += 2; predj = (j - fft_jmp); if (!predj) x += fft_jmp; if (!predj) j = 0; x0[0] = ixheaac_add32_sat(xh0_0, xh20_0); x0[1] = ixheaac_add32_sat(xh1_0, xh21_0); xt0_0 = ixheaac_sub32_sat(xh0_0, xh20_0); yt0_0 = ixheaac_sub32_sat(xh1_0, xh21_0); xt1_0 = ixheaac_add32_sat(xl0_0, xl21_0); yt2_0 = ixheaac_add32_sat(xl1_0, xl20_0); xt2_0 = ixheaac_sub32_sat(xl0_0, xl21_0); yt1_0 = ixheaac_sub32_sat(xl1_0, xl20_0); x2[h2] = ixheaac_add32_sat(MPYLIRC(si10, yt1_0), MPYLIRC(co10, xt1_0)); x2[h2 + 1] = ixheaac_sub32_sat(MPYLIRC(co10, yt1_0), MPYLIRC(si10, xt1_0)); x2[l1] = ixheaac_add32_sat(MPYLIRC(si20, yt0_0), MPYLIRC(co20, xt0_0)); x2[l1 + 1] = ixheaac_sub32_sat(MPYLIRC(co20, yt0_0), MPYLIRC(si20, xt0_0)); yt0_0 = MPYLIRC(si20, yt0_0); x2[l2] = ixheaac_add32_sat(MPYLIRC(si30, yt2_0), MPYLIRC(co30, xt2_0)); x2[l2 + 1] = ixheaac_sub32_sat(MPYLIRC(co30, yt2_0), MPYLIRC(si30, xt2_0)); yt2_0 = MPYLIRC(si30, yt2_0); } stride >>= 2; } y0 = ptr_y; y2 = ptr_y + (int)npoints; x0 = ptr_x; x2 = ptr_x + (int)(npoints >> 1); y1 = y0 + (int)(npoints >> 2); y3 = y2 + (int)(npoints >> 2); l1 = norm + 1; j0 = 8; n0 = npoints >> 1; j = 0; for (i = 0; i < 4; i++) { int t1, t2; h2 = rev_dig[i]; t1 = h2 << 1; t2 = t1 + 1; y0[t1] = ixheaac_add32_sat(x0[0], x0[2]); y2[t1] = ixheaac_sub32_sat(x0[0], x0[2]); y0[t2] = ixheaac_add32_sat(x0[1], x0[3]); y2[t2] = ixheaac_sub32_sat(x0[1], x0[3]); y1[t1] = ixheaac_add32_sat(x0[4], x0[6]); y3[t1] = ixheaac_sub32_sat(x0[4], x0[6]); y1[t2] = ixheaac_add32_sat(x0[5], x0[7]); y3[t2] = ixheaac_sub32_sat(x0[5], x0[7]); x0 += 8; t1 += 2; t2 += 2; y0[t1] = ixheaac_add32_sat(x2[0], x2[2]); y2[t1] = ixheaac_sub32_sat(x2[0], x2[2]); y0[t2] = ixheaac_add32_sat(x2[1], x2[3]); y2[t2] = ixheaac_sub32_sat(x2[1], x2[3]); y1[t1] = ixheaac_add32_sat(x2[4], x2[6]); y3[t1] = ixheaac_sub32_sat(x2[4], x2[6]); y1[t2] = ixheaac_add32_sat(x2[5], x2[7]); y3[t2] = ixheaac_sub32_sat(x2[5], x2[7]); x2 += 8; j += j0; if (j == n0) { j += n0; x0 += (int)npoints >> 1; x2 += (int)npoints >> 1; } } } VOID ixheaacd_dec_rearrange_short(WORD32 *ip, WORD32 *op, WORD32 mdct_len_2, WORD16 *re_arr_tab) { WORD32 n, i = 0; for (n = 0; n < mdct_len_2; n++) { WORD32 idx = re_arr_tab[n] << 1; op[i++] = ip[idx]; op[i++] = ip[idx + 1]; } } VOID ixheaacd_ld_dec_fft_15_opt(WORD32 *inp, WORD32 *op, WORD32 *fft3out, WORD16 *ptr_re_arr_tab_sml_240) { WORD32 i, n, idx; WORD32 *buf1, *buf2, *buf1a; WORD32 add_r, sub_r; WORD32 add_i, sub_i; WORD32 x_01_r, x_01_i, temp; WORD32 p1, p2, p3, p4; WORD32 sinmu = 1859775393; WORD32 c_51 = 2042378317; WORD32 c_52 = -1652318768; WORD32 c_53 = -780119100; WORD32 c_54 = 1200479854; WORD32 c_55 = -1342177280; WORD32 r1, r2, r3, r4; WORD32 s1, s2, s3, s4, t, temp1, temp2; WORD32 *fft3outptr = fft3out; WORD32 xr_0, xr_1, xr_2; WORD32 xi_0, xi_1, xi_2; buf2 = fft3out; buf1 = buf1a = fft3out; n = 0; { *buf1++ = inp[0]; *buf1++ = inp[1]; *buf1++ = inp[192]; *buf1++ = inp[193]; *buf1++ = inp[384]; *buf1++ = inp[385]; *buf1++ = inp[576]; *buf1++ = inp[577]; *buf1++ = inp[768]; *buf1++ = inp[769]; r1 = ixheaac_add32_sat(buf1a[2], buf1a[8]); r4 = ixheaac_sub32_sat(buf1a[2], buf1a[8]); r3 = ixheaac_add32_sat(buf1a[4], buf1a[6]); r2 = ixheaac_sub32_sat(buf1a[4], buf1a[6]); t = ixheaac_mult32_shl(ixheaac_sub32_sat(r1, r3), c_54); r1 = ixheaac_add32_sat(r1, r3); temp1 = ixheaac_add32_sat(buf1a[0], r1); r1 = ixheaac_add32_sat(temp1, (ixheaac_mult32_shl(r1, c_55) << 1)); r3 = ixheaac_sub32_sat(r1, t); r1 = ixheaac_add32_sat(r1, t); t = ixheaac_mult32_shl(ixheaac_add32_sat(r4, r2), c_51); r4 = ixheaac_add32_sat(t, (ixheaac_mult32_shl(r4, c_52) << 1)); r2 = ixheaac_add32_sat(t, ixheaac_mult32_shl(r2, c_53)); s1 = ixheaac_add32_sat(buf1a[3], buf1a[9]); s4 = ixheaac_sub32_sat(buf1a[3], buf1a[9]); s3 = ixheaac_add32_sat(buf1a[5], buf1a[7]); s2 = ixheaac_sub32_sat(buf1a[5], buf1a[7]); t = ixheaac_mult32_shl(ixheaac_sub32_sat(s1, s3), c_54); s1 = ixheaac_add32_sat(s1, s3); temp2 = ixheaac_add32_sat(buf1a[1], s1); s1 = ixheaac_add32_sat(temp2, (ixheaac_mult32_shl(s1, c_55) << 1)); s3 = ixheaac_sub32_sat(s1, t); s1 = ixheaac_add32_sat(s1, t); t = ixheaac_mult32_shl(ixheaac_add32_sat(s4, s2), c_51); s4 = ixheaac_add32_sat(t, (ixheaac_mult32_shl(s4, c_52) << 1)); s2 = ixheaac_add32_sat(t, (ixheaac_mult32_shl(s2, c_53))); *buf2++ = temp1; *buf2++ = temp2; *buf2++ = ixheaac_add32_sat(r1, s2); *buf2++ = ixheaac_sub32_sat(s1, r2); *buf2++ = ixheaac_sub32_sat(r3, s4); *buf2++ = ixheaac_add32_sat(s3, r4); *buf2++ = ixheaac_add32_sat(r3, s4); *buf2++ = ixheaac_sub32_sat(s3, r4); *buf2++ = ixheaac_sub32_sat(r1, s2); *buf2++ = ixheaac_add32_sat(s1, r2); buf1a = buf1; *buf1++ = inp[320]; *buf1++ = inp[321]; *buf1++ = inp[512]; *buf1++ = inp[513]; *buf1++ = inp[704]; *buf1++ = inp[705]; *buf1++ = inp[896]; *buf1++ = inp[897]; *buf1++ = inp[128]; *buf1++ = inp[129]; r1 = ixheaac_add32_sat(buf1a[2], buf1a[8]); r4 = ixheaac_sub32_sat(buf1a[2], buf1a[8]); r3 = ixheaac_add32_sat(buf1a[4], buf1a[6]); r2 = ixheaac_sub32_sat(buf1a[4], buf1a[6]); t = ixheaac_mult32_shl(ixheaac_sub32_sat(r1, r3), c_54); r1 = ixheaac_add32_sat(r1, r3); temp1 = ixheaac_add32_sat(buf1a[0], r1); r1 = ixheaac_add32_sat(temp1, (ixheaac_mult32_shl(r1, c_55) << 1)); r3 = ixheaac_sub32_sat(r1, t); r1 = ixheaac_add32_sat(r1, t); t = ixheaac_mult32_shl(ixheaac_add32_sat(r4, r2), c_51); r4 = ixheaac_add32_sat(t, (ixheaac_mult32_shl(r4, c_52) << 1)); r2 = ixheaac_add32_sat(t, ixheaac_mult32_shl(r2, c_53)); s1 = ixheaac_add32_sat(buf1a[3], buf1a[9]); s4 = ixheaac_sub32_sat(buf1a[3], buf1a[9]); s3 = ixheaac_add32_sat(buf1a[5], buf1a[7]); s2 = ixheaac_sub32_sat(buf1a[5], buf1a[7]); t = ixheaac_mult32_shl(ixheaac_sub32_sat(s1, s3), c_54); s1 = ixheaac_add32_sat(s1, s3); temp2 = ixheaac_add32_sat(buf1a[1], s1); s1 = ixheaac_add32_sat(temp2, (ixheaac_mult32_shl(s1, c_55) << 1)); s3 = ixheaac_sub32_sat(s1, t); s1 = ixheaac_add32_sat(s1, t); t = ixheaac_mult32_shl(ixheaac_add32_sat(s4, s2), c_51); s4 = ixheaac_add32_sat(t, (ixheaac_mult32_shl(s4, c_52) << 1)); s2 = ixheaac_add32_sat(t, ixheaac_mult32_shl(s2, c_53)); *buf2++ = temp1; *buf2++ = temp2; *buf2++ = ixheaac_add32_sat(r1, s2); *buf2++ = ixheaac_sub32_sat(s1, r2); *buf2++ = ixheaac_sub32_sat(r3, s4); *buf2++ = ixheaac_add32_sat(s3, r4); *buf2++ = ixheaac_add32_sat(r3, s4); *buf2++ = ixheaac_sub32_sat(s3, r4); *buf2++ = ixheaac_sub32_sat(r1, s2); *buf2++ = ixheaac_add32_sat(s1, r2); buf1a = buf1; *buf1++ = inp[640]; *buf1++ = inp[641]; *buf1++ = inp[832]; *buf1++ = inp[833]; *buf1++ = inp[64]; *buf1++ = inp[65]; *buf1++ = inp[256]; *buf1++ = inp[257]; *buf1++ = inp[448]; *buf1++ = inp[449]; r1 = ixheaac_add32_sat(buf1a[2], buf1a[8]); r4 = ixheaac_sub32_sat(buf1a[2], buf1a[8]); r3 = ixheaac_add32_sat(buf1a[4], buf1a[6]); r2 = ixheaac_sub32_sat(buf1a[4], buf1a[6]); t = ixheaac_mult32_shl(ixheaac_sub32_sat(r1, r3), c_54); r1 = ixheaac_add32_sat(r1, r3); temp1 = ixheaac_add32_sat(buf1a[0], r1); r1 = ixheaac_add32_sat(temp1, (ixheaac_mult32_shl(r1, c_55) << 1)); r3 = ixheaac_sub32_sat(r1, t); r1 = ixheaac_add32_sat(r1, t); t = ixheaac_mult32_shl(ixheaac_add32_sat(r4, r2), c_51); r4 = ixheaac_add32_sat(t, ixheaac_mult32_shl(r4, c_52) << 1); r2 = ixheaac_add32_sat(t, ixheaac_mult32_shl(r2, c_53)); s1 = ixheaac_add32_sat(buf1a[3], buf1a[9]); s4 = ixheaac_sub32_sat(buf1a[3], buf1a[9]); s3 = ixheaac_add32_sat(buf1a[5], buf1a[7]); s2 = ixheaac_sub32_sat(buf1a[5], buf1a[7]); t = ixheaac_mult32_shl(ixheaac_sub32_sat(s1, s3), c_54); s1 = ixheaac_add32_sat(s1, s3); temp2 = ixheaac_add32_sat(buf1a[1], s1); s1 = ixheaac_add32_sat(temp2, (ixheaac_mult32_shl(s1, c_55) << 1)); s3 = ixheaac_sub32_sat(s1, t); s1 = ixheaac_add32_sat(s1, t); t = ixheaac_mult32_shl(ixheaac_add32_sat(s4, s2), c_51); s4 = ixheaac_add32_sat(t, (ixheaac_mult32_shl(s4, c_52) << 1)); s2 = ixheaac_add32_sat(t, ixheaac_mult32_shl(s2, c_53)); *buf2++ = temp1; *buf2++ = temp2; *buf2++ = ixheaac_add32_sat(r1, s2); *buf2++ = ixheaac_sub32_sat(s1, r2); *buf2++ = ixheaac_sub32_sat(r3, s4); *buf2++ = ixheaac_add32_sat(s3, r4); *buf2++ = ixheaac_add32_sat(r3, s4); *buf2++ = ixheaac_sub32_sat(s3, r4); *buf2++ = ixheaac_sub32_sat(r1, s2); *buf2++ = ixheaac_add32_sat(s1, r2); buf1a = buf1; } n = 0; for (i = 0; i < FFT5; i++) { xr_0 = fft3outptr[0]; xi_0 = fft3outptr[1]; xr_1 = fft3outptr[10]; xi_1 = fft3outptr[11]; xr_2 = fft3outptr[20]; xi_2 = fft3outptr[21]; x_01_r = ixheaac_add32_sat(xr_0, xr_1); x_01_i = ixheaac_add32_sat(xi_0, xi_1); add_r = ixheaac_add32_sat(xr_1, xr_2); add_i = ixheaac_add32_sat(xi_1, xi_2); sub_r = ixheaac_sub32_sat(xr_1, xr_2); sub_i = ixheaac_sub32_sat(xi_1, xi_2); p1 = add_r >> 1; p2 = ixheaac_mult32_shl(sub_i, sinmu); p3 = ixheaac_mult32_shl(sub_r, sinmu); p4 = add_i >> 1; temp = ixheaac_sub32_sat(xr_0, p1); temp1 = ixheaac_add32_sat(xi_0, p3); temp2 = ixheaac_sub32_sat(xi_0, p3); idx = ptr_re_arr_tab_sml_240[n++] << 1; op[idx] = ixheaac_add32_sat(x_01_r, xr_2); op[idx + 1] = ixheaac_add32_sat(x_01_i, xi_2); idx = ptr_re_arr_tab_sml_240[n++] << 1; op[idx] = ixheaac_add32_sat(temp, p2); op[idx + 1] = ixheaac_sub32_sat(temp2, p4); idx = ptr_re_arr_tab_sml_240[n++] << 1; op[idx] = ixheaac_sub32_sat(temp, p2); op[idx + 1] = ixheaac_sub32_sat(temp1, p4); fft3outptr += 2; } return; } VOID ixheaacd_fft_120(ia_aac_dec_imdct_tables_struct *imdct_tables_ptr, WORD32 npoints, WORD32* ptr_x, WORD32* ptr_y) { WORD32 i; WORD32 *buf1, *buf2; WORD32 *inp, *op; inp = ptr_x; op = ptr_y; ixheaacd_dec_rearrange_short(inp, op, 60, imdct_tables_ptr->re_arr_tab_4); buf1 = op; buf2 = inp; for (i = 0; i < FFT15; i++) { { WORD32 x_0, x_1, x_2, x_3, x_4, x_5, x_6, x_7; WORD32 *y0, *y1, *y2, *y3; WORD32 *x0; WORD32 xh0_0, xh1_0, xh0_1, xh1_1, xl0_0, xl1_0, xl0_1, xl1_1; WORD32 h2; WORD32 n00, n01, n10, n11, n20, n21, n30, n31; ptr_x = buf1; ptr_y = buf2; npoints = 4; h2 = 0; y0 = ptr_y; y2 = ptr_y + (WORD32)npoints; x0 = ptr_x; y1 = y0 + (WORD32)(npoints >> 1); y3 = y2 + (WORD32)(npoints >> 1); x_0 = x0[0]; x_1 = x0[1]; x_2 = x0[2]; x_3 = x0[3]; x_4 = x0[4]; x_5 = x0[5]; x_6 = x0[6]; x_7 = x0[7]; x0 += 8; xh0_0 = ixheaac_add32_sat(x_0, x_4); xh1_0 = ixheaac_add32_sat(x_1, x_5); xl0_0 = ixheaac_sub32_sat(x_0, x_4); xl1_0 = ixheaac_sub32_sat(x_1, x_5); xh0_1 = ixheaac_add32_sat(x_2, x_6); xh1_1 = ixheaac_add32_sat(x_3, x_7); xl0_1 = ixheaac_sub32_sat(x_2, x_6); xl1_1 = ixheaac_sub32_sat(x_3, x_7); n00 = ixheaac_add32_sat(xh0_0, xh0_1); n01 = ixheaac_add32_sat(xh1_0, xh1_1); n10 = ixheaac_add32_sat(xl0_0, xl1_1); n11 = ixheaac_sub32_sat(xl1_0, xl0_1); n20 = ixheaac_sub32_sat(xh0_0, xh0_1); n21 = ixheaac_sub32_sat(xh1_0, xh1_1); n30 = ixheaac_sub32_sat(xl0_0, xl1_1); n31 = ixheaac_add32_sat(xl1_0, xl0_1); y0[2 * h2] = n00; y0[2 * h2 + 1] = n01; y1[2 * h2] = n10; y1[2 * h2 + 1] = n11; y2[2 * h2] = n20; y2[2 * h2 + 1] = n21; y3[2 * h2] = n30; y3[2 * h2 + 1] = n31; } buf1 += (FFT4 * 2); buf2 += (FFT4 * 2); } ixheaacd_dec_rearrange_short(inp, op, 60, imdct_tables_ptr->re_arr_tab_15_4); buf1 = op; buf2 = inp; for (i = 0; i < FFT4; i++) { ixheaacd_fft_960_15(buf1, buf2, imdct_tables_ptr); buf1 += (FFT15 * 2); buf2 += (FFT15 * 2); } ixheaacd_dec_rearrange_short(inp, op, 60, imdct_tables_ptr->re_arr_tab_120); } VOID ixheaacd_fft_960_15(WORD32 *inp, WORD32 *op, ia_aac_dec_imdct_tables_struct *imdct_tables_ptr) { WORD32 i; WORD32 *buf1, *buf2; ixheaacd_dec_rearrange_short(inp, op, FFT15, imdct_tables_ptr->re_arr_tab_5); buf1 = op; buf2 = inp; for (i = 0; i < FFT3; i++) { ixheaacd_fft_5(buf1, buf2); buf1 += (FFT5 * 2); buf2 += (FFT5 * 2); } ixheaacd_dec_rearrange_short(inp, op, FFT15, imdct_tables_ptr->re_arr_tab_3); buf1 = op; buf2 = inp; for (i = 0; i < FFT5; i++) { ixheaacd_fft_3(buf1, buf2); buf1 += (FFT3 * 2); buf2 += (FFT3 * 2); } ixheaacd_dec_rearrange_short(inp, op, FFT15, imdct_tables_ptr->re_arr_tab_sml); } VOID ixheaacd_fft_3(WORD32 *inp, WORD32 *op) { WORD32 add_r, sub_r; WORD32 add_i, sub_i; WORD32 x_01_r, x_01_i, temp; WORD32 p1, p2, p3, p4; WORD32 sinmu = 1859775393; x_01_r = ixheaac_add32_sat(inp[0], inp[2]); x_01_i = ixheaac_add32_sat(inp[1], inp[3]); add_r = ixheaac_add32_sat(inp[2], inp[4]); add_i = ixheaac_add32_sat(inp[3], inp[5]); sub_r = ixheaac_sub32_sat(inp[2], inp[4]); sub_i = ixheaac_sub32_sat(inp[3], inp[5]); p1 = add_r >> 1; p2 = ixheaac_mult32_shl(sub_i, sinmu); p3 = ixheaac_mult32_shl(sub_r, sinmu); p4 = add_i >> 1; temp = ixheaac_sub32_sat(inp[0], p1); op[0] = ixheaac_add32_sat(x_01_r, inp[4]); op[1] = ixheaac_add32_sat(x_01_i, inp[5]); op[2] = ixheaac_add32_sat(temp, p2); op[3] = ixheaac_sub32_sat(ixheaac_sub32_sat(inp[1], p3), p4); op[4] = ixheaac_sub32_sat(temp, p2); op[5] = ixheaac_sub32_sat(ixheaac_add32_sat(inp[1], p3), p4); } VOID ixheaacd_fft_5(WORD32 *inp, WORD32 *op) { WORD32 c_51 = 2042378317; WORD32 c_52 = -1652318768; WORD32 c_53 = -780119100; WORD32 c_54 = 1200479854; WORD32 c_55 = -1342177280; WORD32 r1, r2, r3, r4; WORD32 s1, s2, s3, s4, t, temp1, temp2; r1 = ixheaac_add32_sat(inp[2], inp[8]); r4 = ixheaac_sub32_sat(inp[2], inp[8]); r3 = ixheaac_add32_sat(inp[4], inp[6]); r2 = ixheaac_sub32_sat(inp[4], inp[6]); t = ixheaac_mult32_shl(ixheaac_sub32_sat(r1, r3), c_54); r1 = ixheaac_add32_sat(r1, r3); temp1 = ixheaac_add32_sat(inp[0], r1); r1 = ixheaac_add32_sat( temp1, ixheaac_shl32_sat((ixheaac_mult32_shl(r1, c_55)), 1)); r3 = ixheaac_sub32_sat(r1, t); r1 = ixheaac_add32_sat(r1, t); t = ixheaac_mult32_shl(ixheaac_add32_sat(r4, r2), c_51); r4 = ixheaac_add32_sat( t, ixheaac_shl32_sat(ixheaac_mult32_shl(r4, c_52), 1)); r2 = ixheaac_add32_sat(t, ixheaac_mult32_shl(r2, c_53)); s1 = ixheaac_add32_sat(inp[3], inp[9]); s4 = ixheaac_sub32_sat(inp[3], inp[9]); s3 = ixheaac_add32_sat(inp[5], inp[7]); s2 = ixheaac_sub32_sat(inp[5], inp[7]); t = ixheaac_mult32_shl(ixheaac_sub32_sat(s1, s3), c_54); s1 = ixheaac_add32_sat(s1, s3); temp2 = ixheaac_add32_sat(inp[1], s1); s1 = ixheaac_add32_sat( temp2, ixheaac_shl32_sat((ixheaac_mult32_shl(s1, c_55)), 1)); s3 = ixheaac_sub32_sat(s1, t); s1 = ixheaac_add32_sat(s1, t); t = ixheaac_mult32_shl(ixheaac_add32_sat(s4, s2), c_51); s4 = ixheaac_add32_sat( t, ixheaac_shl32_sat((ixheaac_mult32_shl(s4, c_52)), 1)); s2 = ixheaac_add32_sat(t, (ixheaac_mult32_shl(s2, c_53))); op[0] = temp1; op[1] = temp2; op[2] = ixheaac_add32_sat(r1, s2); op[3] = ixheaac_sub32_sat(s1, r2); op[4] = ixheaac_sub32_sat(r3, s4); op[5] = ixheaac_add32_sat(s3, r4); op[6] = ixheaac_add32_sat(r3, s4); op[7] = ixheaac_sub32_sat(s3, r4); op[8] = ixheaac_sub32_sat(r1, s2); op[9] = ixheaac_add32_sat(s1, r2); } VOID ixheaacd_fft_480_ld(WORD32 *inp, WORD32 *op, ia_aac_dec_imdct_tables_struct *imdct_tables_ptr) { WORD32 i; WORD32 *buf1, *buf2; UWORD8 *re_arr_tab_sml_240_ptr; (*ixheaacd_aac_ld_dec_rearrange)(inp, op, MDCT_LEN_BY2, imdct_tables_ptr->re_arr_tab_16); buf1 = op; buf2 = inp; for (i = 0; i < FFT15; i++) { (*ixheaacd_fft32x32_ld2)(imdct_tables_ptr, 16, buf1, buf2); buf1 += (FFT16X2); buf2 += (FFT16X2); } re_arr_tab_sml_240_ptr = imdct_tables_ptr->re_arr_tab_sml_240; buf1 = inp; for (i = 0; i < FFT16; i++) { (*ixheaacd_fft_15_ld)(buf1, op, ixheaacd_fft5out, re_arr_tab_sml_240_ptr); re_arr_tab_sml_240_ptr += FFT15; buf1 += 2; } } VOID ixheaacd_pre_twiddle_960(WORD32 *xptr, WORD32 *data, WORD32 n, WORD32 *cos_sin_ptr, WORD32 neg_expo) { WORD npoints_4, i; WORD32 tempr, tempi, temp; WORD32 c, c1, s, s1; WORD32 *in_ptr1, *in_ptr2; WORD32 *xprt1 = xptr + (n - 1); npoints_4 = n >> 2; in_ptr1 = data; in_ptr2 = data + n - 1; for (i = 0; i < npoints_4; i++) { c = *cos_sin_ptr++; s = *cos_sin_ptr++; tempr = *in_ptr1++; tempi = *in_ptr2--; temp = -ixheaac_add32(ixheaac_mult32x32in32(tempr, c), ixheaac_mult32x32in32(tempi, s)); *xptr++ = ixheaac_shr32_dir_sat(temp, neg_expo); temp = -ixheaac_sub32(ixheaac_mult32x32in32(tempi, c), ixheaac_mult32x32in32(tempr, s)); *xptr++ = ixheaac_shr32_dir_sat(temp, neg_expo); c1 = *cos_sin_ptr++; s1 = *cos_sin_ptr++; tempi = *in_ptr1++; tempr = *in_ptr2--; temp = -ixheaac_sub32(ixheaac_mult32x32in32(tempi, c1), ixheaac_mult32x32in32(tempr, s1)); *xprt1-- = ixheaac_shr32_dir_sat(temp, neg_expo); temp = -ixheaac_add32(ixheaac_mult32x32in32(tempr, c1), ixheaac_mult32x32in32(tempi, s1)); *xprt1-- = ixheaac_shr32_dir_sat(temp, neg_expo); } } VOID ixheaacd_pre_twiddle_120(WORD32 *xptr, WORD32 *data, WORD32 n, WORD16 *cos_sin_ptr, WORD32 neg_expo) { WORD npoints_4, i; WORD32 tempr, tempi, temp; WORD16 c, c1, s, s1; WORD32 *in_ptr1, *in_ptr2; WORD32 *xprt1 = xptr + (n - 1); npoints_4 = n >> 2; in_ptr1 = data; in_ptr2 = data + n - 1; for (i = 0; i < npoints_4; i++) { c = *cos_sin_ptr++; s = *cos_sin_ptr++; tempr = *in_ptr1++; tempi = *in_ptr2--; temp = -ixheaac_add32(ixheaac_mult32x16in32(tempr, c), ixheaac_mult32x16in32(tempi, s)); *xptr++ = ixheaac_shr32_dir_sat(temp, neg_expo); temp = -ixheaac_sub32(ixheaac_mult32x16in32(tempi, c), ixheaac_mult32x16in32(tempr, s)); *xptr++ = ixheaac_shr32_dir_sat(temp, neg_expo); c1 = *cos_sin_ptr++; s1 = *cos_sin_ptr++; tempi = *in_ptr1++; tempr = *in_ptr2--; temp = -ixheaac_sub32(ixheaac_mult32x16in32(tempi, c1), ixheaac_mult32x16in32(tempr, s1)); *xprt1-- = ixheaac_shr32_dir_sat(temp, neg_expo); temp = -ixheaac_add32(ixheaac_mult32x16in32(tempr, c1), ixheaac_mult32x16in32(tempi, s1)); *xprt1-- = ixheaac_shr32_dir_sat(temp, neg_expo); } } VOID ixheaacd_pre_twiddle(WORD32 *xptr, WORD32 *data, WORD32 n, WORD32 *cos_sin_ptr, WORD32 neg_expo) { WORD npoints_4, i; WORD32 tempr, tempi, temp; WORD32 c, c1, s, s1; WORD32 *in_ptr1, *in_ptr2; npoints_4 = n >> 2; in_ptr1 = data; in_ptr2 = data + n - 1; if (neg_expo >= 0) { for (i = npoints_4 - 1; i >= 0; i--) { c = *cos_sin_ptr++; c1 = *cos_sin_ptr++; s = *cos_sin_ptr++; s1 = *cos_sin_ptr++; tempr = *in_ptr1; tempi = *in_ptr2; in_ptr1 += 2; in_ptr2 -= 2; temp = -ixheaac_add32(ixheaac_mult32(tempr, c), ixheaac_mult32(tempi, s)); *xptr++ = ixheaac_shr32(temp, neg_expo); temp = ixheaac_sub32(ixheaac_mult32(tempr, s), ixheaac_mult32(tempi, c)); *xptr++ = ixheaac_shr32(temp, neg_expo); tempr = *in_ptr1; tempi = *in_ptr2; in_ptr1 += 2; in_ptr2 -= 2; temp = -ixheaac_add32(ixheaac_mult32(tempr, c1), ixheaac_mult32(tempi, s1)); *xptr++ = ixheaac_shr32(temp, neg_expo); temp = ixheaac_sub32(ixheaac_mult32(tempr, s1), ixheaac_mult32(tempi, c1)); *xptr++ = ixheaac_shr32(temp, neg_expo); } } else { neg_expo = -neg_expo; for (i = npoints_4 - 1; i >= 0; i--) { c = *cos_sin_ptr++; c1 = *cos_sin_ptr++; s = *cos_sin_ptr++; s1 = *cos_sin_ptr++; tempr = *in_ptr1; tempi = *in_ptr2; in_ptr1 += 2; in_ptr2 -= 2; temp = -ixheaac_add32(ixheaac_mult32(tempr, c), ixheaac_mult32(tempi, s)); *xptr++ = ixheaac_shl32(temp, neg_expo); temp = ixheaac_sub32(ixheaac_mult32(tempr, s), ixheaac_mult32(tempi, c)); *xptr++ = ixheaac_shl32(temp, neg_expo); tempr = *in_ptr1; tempi = *in_ptr2; in_ptr1 += 2; in_ptr2 -= 2; temp = -ixheaac_add32(ixheaac_mult32(tempr, c1), ixheaac_mult32(tempi, s1)); *xptr++ = ixheaac_shl32(temp, neg_expo); temp = ixheaac_sub32(ixheaac_mult32(tempr, s1), ixheaac_mult32(tempi, c1)); *xptr++ = ixheaac_shl32(temp, neg_expo); } } } VOID ixheaacd_post_twiddle_120(WORD32 out[], WORD32 x[], const WORD16 *cos_sin_ptr, WORD m) { WORD i; WORD16 c, c1, s, s1; WORD32 tempr, tempi, temp; WORD32 *in_ptr2 = x + (m - 1); WORD32 *in_ptr1 = x; WORD32 *xptr = out; WORD32 *xptr1 = out + (m - 1); for (i = 0; i < m; i += 4) { c = *cos_sin_ptr++; s = *cos_sin_ptr++; c1 = *cos_sin_ptr++; s1 = *cos_sin_ptr++; tempr = *in_ptr1++; tempi = *in_ptr1++; temp = -ixheaac_sub32_sat(ixheaac_mult32x16in32(tempr, s), ixheaac_mult32x16in32(tempi, c)); *xptr1-- = temp; temp = -ixheaac_add32_sat(ixheaac_mult32x16in32(tempr, c), ixheaac_mult32x16in32(tempi, s)); *xptr++ = temp; tempi = *in_ptr2--; tempr = *in_ptr2--; temp = -ixheaac_sub32_sat(ixheaac_mult32x16in32(tempr, s1), ixheaac_mult32x16in32(tempi, c1)); *xptr++ = temp; temp = -ixheaac_add32_sat(ixheaac_mult32x16in32(tempr, c1), ixheaac_mult32x16in32(tempi, s1)); *xptr1-- = temp; } } VOID ixheaacd_post_twiddle_960(WORD32 out[], WORD32 x[], const WORD32 *cos_sin_ptr, WORD m) { WORD i; WORD32 c, c1, s, s1; WORD32 tempr, tempi, temp; WORD32 *in_ptr2 = x + (m - 1); WORD32 *in_ptr1 = x; WORD32 *xptr = out; WORD32 *xptr1 = out + (m - 1); for (i = 0; i < m; i += 4) { c = *cos_sin_ptr++; s = *cos_sin_ptr++; c1 = *cos_sin_ptr++; s1 = *cos_sin_ptr++; tempr = *in_ptr1++; tempi = *in_ptr1++; temp = -ixheaac_sub32_sat(ixheaac_mult32x32in32(tempr, s), ixheaac_mult32x32in32(tempi, c)); *xptr1-- = temp; temp = -ixheaac_add32_sat(ixheaac_mult32x32in32(tempr, c), ixheaac_mult32x32in32(tempi, s)); *xptr++ = temp; tempi = *in_ptr2--; tempr = *in_ptr2--; temp = -ixheaac_sub32_sat(ixheaac_mult32x32in32(tempr, s1), ixheaac_mult32x32in32(tempi, c1)); *xptr++ = temp; temp = -ixheaac_add32_sat(ixheaac_mult32x32in32(tempr, c1), ixheaac_mult32x32in32(tempi, s1)); *xptr1-- = temp; } } VOID ixheaacd_post_twiddle_ld(WORD32 out[], WORD32 x[], const WORD32 *cos_sin_ptr, WORD m) { WORD i; WORD32 *ptr_x = &x[0]; WORD32 *ptr_out, *ptr_out1; ptr_out = &out[0]; ptr_out1 = &out[m - 1]; for (i = (m >> 2) - 1; i >= 0; i--) { WORD32 c, c1, s, s1; WORD32 re, im; c = *cos_sin_ptr++; c1 = *cos_sin_ptr++; s = *cos_sin_ptr++; s1 = *cos_sin_ptr++; re = *ptr_x++; im = *ptr_x++; *ptr_out1 = ixheaac_sub32(ixheaac_mult32(im, c), ixheaac_mult32(re, s)); *ptr_out = -ixheaac_add32(ixheaac_mult32(re, c), ixheaac_mult32(im, s)); ptr_out += 2; ptr_out1 -= 2; re = *ptr_x++; im = *ptr_x++; *ptr_out1 = ixheaac_sub32(ixheaac_mult32(im, c1), ixheaac_mult32(re, s1)); *ptr_out = -ixheaac_add32(ixheaac_mult32(re, c1), ixheaac_mult32(im, s1)); ptr_out += 2; ptr_out1 -= 2; } } VOID ixheaacd_post_twiddle_eld(WORD32 out[], WORD32 x[], const WORD32 *cos_sin_ptr, WORD m) { WORD i = 0; WORD32 *ptr_x = &x[0]; WORD32 *ptr_out_767, *ptr_out_256; WORD32 *ptr_out_768, *ptr_out_255; WORD32 *ptr_out_0, *ptr_out_1279; WORD32 tempr, tempi; ptr_out_767 = &out[m + (m >> 1) - 1 - 2 * i]; ptr_out_256 = &out[(m >> 1) + 2 * i]; ptr_out_768 = &out[m + (m >> 1) + 2 * i]; ptr_out_255 = &out[(m >> 1) - 1 - 2 * i]; for (i = 0; i < (m >> 3); i++) { WORD32 c, c1, s, s1; WORD32 re, im; c = *cos_sin_ptr++; c1 = *cos_sin_ptr++; s = *cos_sin_ptr++; s1 = *cos_sin_ptr++; re = *ptr_x++; im = *ptr_x++; tempi = ixheaac_sub32(ixheaac_mult32(im, c), ixheaac_mult32(re, s)); tempr = -ixheaac_add32(ixheaac_mult32(re, c), ixheaac_mult32(im, s)); *ptr_out_767 = tempr; *ptr_out_256 = tempi; *ptr_out_768 = *ptr_out_767; *ptr_out_255 = -*ptr_out_256; ptr_out_256 += 2; ptr_out_767 -= 2; ptr_out_768 += 2; ptr_out_255 -= 2; re = *ptr_x++; im = *ptr_x++; tempi = ixheaac_sub32(ixheaac_mult32(im, c1), ixheaac_mult32(re, s1)); tempr = -ixheaac_add32(ixheaac_mult32(re, c1), ixheaac_mult32(im, s1)); *ptr_out_767 = tempr; *ptr_out_256 = tempi; *ptr_out_768 = *ptr_out_767; *ptr_out_255 = -*ptr_out_256; ptr_out_256 += 2; ptr_out_767 -= 2; ptr_out_768 += 2; ptr_out_255 -= 2; } ptr_out_0 = &out[2 * 2 * i - (m >> 1)]; ptr_out_1279 = &out[m + m + (m >> 1) - 1 - 2 * 2 * i]; for (; i < (m >> 2); i++) { WORD32 c, c1, s, s1; WORD32 re, im; c = *cos_sin_ptr++; c1 = *cos_sin_ptr++; s = *cos_sin_ptr++; s1 = *cos_sin_ptr++; re = *ptr_x++; im = *ptr_x++; tempi = ixheaac_sub32(ixheaac_mult32(im, c), ixheaac_mult32(re, s)); tempr = -ixheaac_add32(ixheaac_mult32(re, c), ixheaac_mult32(im, s)); *ptr_out_767 = tempr; *ptr_out_256 = tempi; *ptr_out_0 = -*ptr_out_767; *ptr_out_1279 = *ptr_out_256; ptr_out_256 += 2; ptr_out_767 -= 2; ptr_out_0 += 2; ptr_out_1279 -= 2; re = *ptr_x++; im = *ptr_x++; tempi = ixheaac_sub32(ixheaac_mult32(im, c1), ixheaac_mult32(re, s1)); tempr = -ixheaac_add32(ixheaac_mult32(re, c1), ixheaac_mult32(im, s1)); *ptr_out_767 = tempr; *ptr_out_256 = tempi; *ptr_out_0 = -*ptr_out_767; *ptr_out_1279 = *ptr_out_256; ptr_out_256 += 2; ptr_out_767 -= 2; ptr_out_0 += 2; ptr_out_1279 -= 2; } } VOID ixheaacd_fft32x32_ld_dec(ia_aac_dec_imdct_tables_struct *imdct_tables_ptr, WORD32 npoints, WORD32 *ptr_x, WORD32 *ptr_y) { WORD32 i, j, l1, l2, h2, predj, tw_offset, stride, fft_jmp, k; WORD32 xt0_0, yt0_0, xt1_0, yt1_0, xt2_0, yt2_0; WORD32 xh0_0, xh1_0, xh20_0, xh21_0, xl0_0, xl1_0, xl20_0, xl21_0; WORD32 xh0_1, xh1_1, xl0_1, xl1_1; WORD32 x_0, x_1, x_2, x_3, x_l1_0, x_l1_1, x_l2_0, x_l2_1; WORD32 xh0_2, xh1_2, xl0_2, xl1_2, xh0_3, xh1_3, xl0_3, xl1_3; WORD32 x_4, x_5, x_6, x_7, x_h2_0, x_h2_1; WORD32 x_8, x_9, x_a, x_b, x_c, x_d, x_e, x_f; WORD32 si10, si20, si30, co10, co20, co30; WORD32 *w; WORD32 *x, *x2, *x0; WORD32 *y0, *y1, *y2, *y3; WORD32 n00, n10, n20, n30, n01, n11, n21, n31; WORD32 n02, n12, n22, n32, n03, n13, n23, n33; WORD32 n0, j0; WORD32 radix; WORD32 norm; WORD32 m; WORD32 *ptr_w; if (npoints == 256) ptr_w = imdct_tables_ptr->w_256; else if (npoints == 32) ptr_w = (WORD32*)imdct_tables_ptr->w_32; else ptr_w = imdct_tables_ptr->w_16; for (i = 31, m = 1; (npoints & (1 << i)) == 0; i--, m++) ; radix = m & 1 ? 2 : 4; norm = m - 2; stride = npoints; tw_offset = 0; fft_jmp = 6 * stride; while (stride > radix) { j = 0; fft_jmp >>= 2; h2 = stride >> 1; l1 = stride; l2 = stride + (stride >> 1); x = ptr_x; w = ptr_w + tw_offset; tw_offset += fft_jmp; stride >>= 2; for (i = 0; i < npoints; i += 4) { co10 = w[j + 1]; si10 = w[j + 0]; co20 = w[j + 3]; si20 = w[j + 2]; co30 = w[j + 5]; si30 = w[j + 4]; x_0 = x[0]; x_1 = x[1]; x_l1_0 = x[l1]; x_l1_1 = x[l1 + 1]; x_l2_0 = x[l2]; x_l2_1 = x[l2 + 1]; x_h2_0 = x[h2]; x_h2_1 = x[h2 + 1]; xh0_0 = ixheaac_add32_sat(x_0, x_l1_0); xh1_0 = ixheaac_add32_sat(x_1, x_l1_1); xl0_0 = ixheaac_sub32_sat(x_0, x_l1_0); xl1_0 = ixheaac_sub32_sat(x_1, x_l1_1); xh20_0 = ixheaac_add32_sat(x_h2_0, x_l2_0); xh21_0 = ixheaac_add32_sat(x_h2_1, x_l2_1); xl20_0 = ixheaac_sub32_sat(x_h2_0, x_l2_0); xl21_0 = ixheaac_sub32_sat(x_h2_1, x_l2_1); x0 = x; x2 = x0; j += 6; x += 2; predj = (j - fft_jmp); if (!predj) x += fft_jmp; if (!predj) j = 0; x0[0] = ixheaac_add32_sat(xh0_0, xh20_0); x0[1] = ixheaac_add32_sat(xh1_0, xh21_0); xt0_0 = ixheaac_sub32_sat(xh0_0, xh20_0); yt0_0 = ixheaac_sub32_sat(xh1_0, xh21_0); xt1_0 = ixheaac_add32_sat(xl0_0, xl21_0); yt2_0 = ixheaac_add32_sat(xl1_0, xl20_0); xt2_0 = ixheaac_sub32_sat(xl0_0, xl21_0); yt1_0 = ixheaac_sub32_sat(xl1_0, xl20_0); x2[h2] = MPYHIRC(si10, yt1_0) + MPYHIRC(co10, xt1_0) + (((MPYLUHS(si10, yt1_0) + MPYLUHS(co10, xt1_0) + 0x8000) >> 16) << 1); x2[h2 + 1] = MPYHIRC(co10, yt1_0) - MPYHIRC(si10, xt1_0) + (((MPYLUHS(co10, yt1_0) - MPYLUHS(si10, xt1_0) + 0x8000) >> 16) << 1); x2[l1] = MPYHIRC(si20, yt0_0) + MPYHIRC(co20, xt0_0) + (((MPYLUHS(si20, yt0_0) + MPYLUHS(co20, xt0_0) + 0x8000) >> 16) << 1); x2[l1 + 1] = MPYHIRC(co20, yt0_0) - MPYHIRC(si20, xt0_0) + (((MPYLUHS(co20, yt0_0) - MPYLUHS(si20, xt0_0) + 0x8000) >> 16) << 1); x2[l2] = MPYHIRC(si30, yt2_0) + MPYHIRC(co30, xt2_0) + (((MPYLUHS(si30, yt2_0) + MPYLUHS(co30, xt2_0) + 0x8000) >> 16) << 1); x2[l2 + 1] = MPYHIRC(co30, yt2_0) - MPYHIRC(si30, xt2_0) + (((MPYLUHS(co30, yt2_0) - MPYLUHS(si30, xt2_0) + 0x8000) >> 16) << 1); } } y0 = ptr_y; y2 = ptr_y + (WORD32)npoints; x0 = ptr_x; x2 = ptr_x + (WORD32)(npoints >> 1); if (radix == 2) { y1 = y0 + (WORD32)(npoints >> 2); y3 = y2 + (WORD32)(npoints >> 2); l1 = norm + 1; j0 = 8; n0 = npoints >> 1; } else { y1 = y0 + (WORD32)(npoints >> 1); y3 = y2 + (WORD32)(npoints >> 1); l1 = norm + 2; j0 = 4; n0 = npoints >> 2; } j = 0; k = 0; for (i = 0; i < npoints; i += 8) { if (npoints == 32) h2 = rev_dig[k++]; else DIG_REV(j, l1, h2); x_0 = x0[0]; x_1 = x0[1]; x_2 = x0[2]; x_3 = x0[3]; x_4 = x0[4]; x_5 = x0[5]; x_6 = x0[6]; x_7 = x0[7]; x0 += 8; xh0_0 = ixheaac_add32_sat(x_0, x_4); xh1_0 = ixheaac_add32_sat(x_1, x_5); xl0_0 = ixheaac_sub32_sat(x_0, x_4); xl1_0 = ixheaac_sub32_sat(x_1, x_5); xh0_1 = ixheaac_add32_sat(x_2, x_6); xh1_1 = ixheaac_add32_sat(x_3, x_7); xl0_1 = ixheaac_sub32_sat(x_2, x_6); xl1_1 = ixheaac_sub32_sat(x_3, x_7); n00 = ixheaac_add32_sat(xh0_0, xh0_1); n01 = ixheaac_add32_sat(xh1_0, xh1_1); n10 = ixheaac_add32_sat(xl0_0, xl1_1); n11 = ixheaac_sub32_sat(xl1_0, xl0_1); n20 = ixheaac_sub32_sat(xh0_0, xh0_1); n21 = ixheaac_sub32_sat(xh1_0, xh1_1); n30 = ixheaac_sub32_sat(xl0_0, xl1_1); n31 = ixheaac_add32_sat(xl1_0, xl0_1); if (radix == 2) { n00 = ixheaac_add32_sat(x_0, x_2); n01 = ixheaac_add32_sat(x_1, x_3); n20 = ixheaac_sub32_sat(x_0, x_2); n21 = ixheaac_sub32_sat(x_1, x_3); n10 = ixheaac_add32_sat(x_4, x_6); n11 = ixheaac_add32_sat(x_5, x_7); n30 = ixheaac_sub32_sat(x_4, x_6); n31 = ixheaac_sub32_sat(x_5, x_7); } y0[2 * h2] = n00; y0[2 * h2 + 1] = n01; y1[2 * h2] = n10; y1[2 * h2 + 1] = n11; y2[2 * h2] = n20; y2[2 * h2 + 1] = n21; y3[2 * h2] = n30; y3[2 * h2 + 1] = n31; x_8 = x2[0]; x_9 = x2[1]; x_a = x2[2]; x_b = x2[3]; x_c = x2[4]; x_d = x2[5]; x_e = x2[6]; x_f = x2[7]; x2 += 8; xh0_2 = ixheaac_add32_sat(x_8, x_c); xh1_2 = ixheaac_add32_sat(x_9, x_d); xl0_2 = ixheaac_sub32_sat(x_8, x_c); xl1_2 = ixheaac_sub32_sat(x_9, x_d); xh0_3 = ixheaac_add32_sat(x_a, x_e); xh1_3 = ixheaac_add32_sat(x_b, x_f); xl0_3 = ixheaac_sub32_sat(x_a, x_e); xl1_3 = ixheaac_sub32_sat(x_b, x_f); n02 = ixheaac_add32_sat(xh0_2, xh0_3); n03 = ixheaac_add32_sat(xh1_2, xh1_3); n12 = ixheaac_add32_sat(xl0_2, xl1_3); n13 = ixheaac_sub32_sat(xl1_2, xl0_3); n22 = ixheaac_sub32_sat(xh0_2, xh0_3); n23 = ixheaac_sub32_sat(xh1_2, xh1_3); n32 = ixheaac_sub32_sat(xl0_2, xl1_3); n33 = ixheaac_add32_sat(xl1_2, xl0_3); if (radix == 2) { n02 = ixheaac_add32_sat(x_8, x_a); n03 = ixheaac_add32_sat(x_9, x_b); n22 = ixheaac_sub32_sat(x_8, x_a); n23 = ixheaac_sub32_sat(x_9, x_b); n12 = ixheaac_add32_sat(x_c, x_e); n13 = ixheaac_add32_sat(x_d, x_f); n32 = ixheaac_sub32_sat(x_c, x_e); n33 = ixheaac_sub32_sat(x_d, x_f); } y0[2 * h2 + 2] = n02; y0[2 * h2 + 3] = n03; y1[2 * h2 + 2] = n12; y1[2 * h2 + 3] = n13; y2[2 * h2 + 2] = n22; y2[2 * h2 + 3] = n23; y3[2 * h2 + 2] = n32; y3[2 * h2 + 3] = n33; j += j0; if (j == n0) { j += n0; x0 += (WORD32)npoints >> 1; x2 += (WORD32)npoints >> 1; } } } VOID ixheaacd_rearrange_dec(WORD32 *ip, WORD32 *op, WORD32 mdct_len_2, UWORD8 *re_arr_tab) { WORD32 n, i = 0; for (n = 0; n < mdct_len_2; n++) { WORD32 idx = re_arr_tab[n] << 1; op[i++] = ip[idx]; op[i++] = ip[idx + 1]; } } VOID ixheaacd_fft_15_ld_dec(WORD32 *inp, WORD32 *op, WORD32 *fft3out, UWORD8 *re_arr_tab_sml_240_ptr) { WORD32 i, n, idx; WORD32 *buf1, *buf2, *buf1a; WORD32 add_r, sub_r; WORD32 add_i, sub_i; WORD32 x01_real, x_01_imag, temp; WORD32 p1, p2, p3, p4; WORD32 sinmu = 1859775393; WORD32 cos_51 = 2042378317; WORD32 cos_52 = -1652318768; WORD32 cos_53 = -780119100; WORD32 cos_54 = 1200479854; WORD32 cos_55 = -1342177280; WORD32 r1, r2, r3, r4; WORD32 s1, s2, s3, s4, t, temp1, temp2; WORD32 *fft3outptr = fft3out; WORD32 xr_0, xr_1, xr_2; WORD32 xi_0, xi_1, xi_2; buf2 = fft3out; buf1 = buf1a = fft3out; n = 0; { *buf1++ = inp[0]; *buf1++ = inp[1]; *buf1++ = inp[96]; *buf1++ = inp[97]; *buf1++ = inp[192]; *buf1++ = inp[193]; *buf1++ = inp[288]; *buf1++ = inp[289]; *buf1++ = inp[384]; *buf1++ = inp[385]; r1 = ixheaac_add32_sat(buf1a[2], buf1a[8]); r4 = ixheaac_sub32_sat(buf1a[2], buf1a[8]); r3 = ixheaac_add32_sat(buf1a[4], buf1a[6]); r2 = ixheaac_sub32_sat(buf1a[4], buf1a[6]); t = ixheaac_mult32_shl(ixheaac_sub32_sat(r1, r3), cos_54); r1 = ixheaac_add32_sat(r1, r3); temp1 = ixheaac_add32_sat(buf1a[0], r1); r1 = ixheaac_add32_sat( temp1, ixheaac_shl32_sat((ixheaac_mult32_shl(r1, cos_55)), 1)); r3 = ixheaac_sub32_sat(r1, t); r1 = ixheaac_add32_sat(r1, t); t = ixheaac_mult32_shl(ixheaac_add32_sat(r4, r2), cos_51); r4 = ixheaac_add32_sat( t, ixheaac_shl32_sat(ixheaac_mult32_shl(r4, cos_52), 1)); r2 = ixheaac_add32_sat(t, ixheaac_mult32_shl(r2, cos_53)); s1 = ixheaac_add32_sat(buf1a[3], buf1a[9]); s4 = ixheaac_sub32_sat(buf1a[3], buf1a[9]); s3 = ixheaac_add32_sat(buf1a[5], buf1a[7]); s2 = ixheaac_sub32_sat(buf1a[5], buf1a[7]); t = ixheaac_mult32_shl(ixheaac_sub32_sat(s1, s3), cos_54); s1 = ixheaac_add32_sat(s1, s3); temp2 = ixheaac_add32_sat(buf1a[1], s1); s1 = ixheaac_add32_sat( temp2, ixheaac_shl32_sat((ixheaac_mult32_shl(s1, cos_55)), 1)); s3 = ixheaac_sub32_sat(s1, t); s1 = ixheaac_add32_sat(s1, t); t = ixheaac_mult32_shl(ixheaac_add32_sat(s4, s2), cos_51); s4 = ixheaac_add32_sat( t, ixheaac_shl32_sat((ixheaac_mult32_shl(s4, cos_52)), 1)); s2 = ixheaac_add32_sat(t, (ixheaac_mult32_shl(s2, cos_53))); *buf2++ = temp1; *buf2++ = temp2; *buf2++ = ixheaac_add32_sat(r1, s2); *buf2++ = ixheaac_sub32_sat(s1, r2); *buf2++ = ixheaac_sub32_sat(r3, s4); *buf2++ = ixheaac_add32_sat(s3, r4); *buf2++ = ixheaac_add32_sat(r3, s4); *buf2++ = ixheaac_sub32_sat(s3, r4); *buf2++ = ixheaac_sub32_sat(r1, s2); *buf2++ = ixheaac_add32_sat(s1, r2); buf1a = buf1; *buf1++ = inp[160]; *buf1++ = inp[161]; *buf1++ = inp[256]; *buf1++ = inp[257]; *buf1++ = inp[352]; *buf1++ = inp[353]; *buf1++ = inp[448]; *buf1++ = inp[449]; *buf1++ = inp[64]; *buf1++ = inp[65]; r1 = ixheaac_add32_sat(buf1a[2], buf1a[8]); r4 = ixheaac_sub32_sat(buf1a[2], buf1a[8]); r3 = ixheaac_add32_sat(buf1a[4], buf1a[6]); r2 = ixheaac_sub32_sat(buf1a[4], buf1a[6]); t = ixheaac_mult32_shl(ixheaac_sub32_sat(r1, r3), cos_54); r1 = ixheaac_add32_sat(r1, r3); temp1 = ixheaac_add32_sat(buf1a[0], r1); r1 = ixheaac_add32_sat( temp1, ixheaac_shl32_sat((ixheaac_mult32_shl(r1, cos_55)), 1)); r3 = ixheaac_sub32_sat(r1, t); r1 = ixheaac_add32_sat(r1, t); t = ixheaac_mult32_shl(ixheaac_add32_sat(r4, r2), cos_51); r4 = ixheaac_add32_sat( t, ixheaac_shl32_sat(ixheaac_mult32_shl(r4, cos_52), 1)); r2 = ixheaac_add32_sat(t, ixheaac_mult32_shl(r2, cos_53)); s1 = ixheaac_add32_sat(buf1a[3], buf1a[9]); s4 = ixheaac_sub32_sat(buf1a[3], buf1a[9]); s3 = ixheaac_add32_sat(buf1a[5], buf1a[7]); s2 = ixheaac_sub32_sat(buf1a[5], buf1a[7]); t = ixheaac_mult32_shl(ixheaac_sub32_sat(s1, s3), cos_54); s1 = ixheaac_add32_sat(s1, s3); temp2 = ixheaac_add32_sat(buf1a[1], s1); s1 = ixheaac_add32_sat( temp2, ixheaac_shl32_sat((ixheaac_mult32_shl(s1, cos_55)), 1)); s3 = ixheaac_sub32_sat(s1, t); s1 = ixheaac_add32_sat(s1, t); t = ixheaac_mult32_shl(ixheaac_add32_sat(s4, s2), cos_51); s4 = ixheaac_add32_sat( t, ixheaac_shl32_sat((ixheaac_mult32_shl(s4, cos_52)), 1)); s2 = ixheaac_add32_sat(t, (ixheaac_mult32_shl(s2, cos_53))); *buf2++ = temp1; *buf2++ = temp2; *buf2++ = ixheaac_add32_sat(r1, s2); *buf2++ = ixheaac_sub32_sat(s1, r2); *buf2++ = ixheaac_sub32_sat(r3, s4); *buf2++ = ixheaac_add32_sat(s3, r4); *buf2++ = ixheaac_add32_sat(r3, s4); *buf2++ = ixheaac_sub32_sat(s3, r4); *buf2++ = ixheaac_sub32_sat(r1, s2); *buf2++ = ixheaac_add32_sat(s1, r2); buf1a = buf1; ; *buf1++ = inp[320]; *buf1++ = inp[321]; *buf1++ = inp[416]; *buf1++ = inp[417]; *buf1++ = inp[32]; *buf1++ = inp[33]; *buf1++ = inp[128]; *buf1++ = inp[129]; *buf1++ = inp[224]; *buf1++ = inp[225]; r1 = ixheaac_add32_sat(buf1a[2], buf1a[8]); r4 = ixheaac_sub32_sat(buf1a[2], buf1a[8]); r3 = ixheaac_add32_sat(buf1a[4], buf1a[6]); r2 = ixheaac_sub32_sat(buf1a[4], buf1a[6]); t = ixheaac_mult32_shl(ixheaac_sub32_sat(r1, r3), cos_54); r1 = ixheaac_add32_sat(r1, r3); temp1 = ixheaac_add32_sat(buf1a[0], r1); r1 = ixheaac_add32_sat( temp1, ixheaac_shl32_sat((ixheaac_mult32_shl(r1, cos_55)), 1)); r3 = ixheaac_sub32_sat(r1, t); r1 = ixheaac_add32_sat(r1, t); t = ixheaac_mult32_shl(ixheaac_add32_sat(r4, r2), cos_51); r4 = ixheaac_add32_sat( t, ixheaac_shl32_sat(ixheaac_mult32_shl(r4, cos_52), 1)); r2 = ixheaac_add32_sat(t, ixheaac_mult32_shl(r2, cos_53)); s1 = ixheaac_add32_sat(buf1a[3], buf1a[9]); s4 = ixheaac_sub32_sat(buf1a[3], buf1a[9]); s3 = ixheaac_add32_sat(buf1a[5], buf1a[7]); s2 = ixheaac_sub32_sat(buf1a[5], buf1a[7]); t = ixheaac_mult32_shl(ixheaac_sub32_sat(s1, s3), cos_54); s1 = ixheaac_add32_sat(s1, s3); temp2 = ixheaac_add32_sat(buf1a[1], s1); s1 = ixheaac_add32_sat( temp2, ixheaac_shl32_sat((ixheaac_mult32_shl(s1, cos_55)), 1)); s3 = ixheaac_sub32_sat(s1, t); s1 = ixheaac_add32_sat(s1, t); t = ixheaac_mult32_shl(ixheaac_add32_sat(s4, s2), cos_51); s4 = ixheaac_add32_sat( t, ixheaac_shl32_sat((ixheaac_mult32_shl(s4, cos_52)), 1)); s2 = ixheaac_add32_sat(t, (ixheaac_mult32_shl(s2, cos_53))); *buf2++ = temp1; *buf2++ = temp2; *buf2++ = ixheaac_add32_sat(r1, s2); *buf2++ = ixheaac_sub32_sat(s1, r2); *buf2++ = ixheaac_sub32_sat(r3, s4); *buf2++ = ixheaac_add32_sat(s3, r4); *buf2++ = ixheaac_add32_sat(r3, s4); *buf2++ = ixheaac_sub32_sat(s3, r4); *buf2++ = ixheaac_sub32_sat(r1, s2); *buf2++ = ixheaac_add32_sat(s1, r2); buf1a = buf1; ; } n = 0; for (i = 0; i < FFT5; i++) { xr_0 = fft3outptr[0]; xi_0 = fft3outptr[1]; xr_1 = fft3outptr[10]; xi_1 = fft3outptr[11]; xr_2 = fft3outptr[20]; xi_2 = fft3outptr[21]; x01_real = ixheaac_add32_sat(xr_0, xr_1); x_01_imag = ixheaac_add32_sat(xi_0, xi_1); add_r = ixheaac_add32_sat(xr_1, xr_2); add_i = ixheaac_add32_sat(xi_1, xi_2); sub_r = ixheaac_sub32_sat(xr_1, xr_2); sub_i = ixheaac_sub32_sat(xi_1, xi_2); p1 = add_r >> 1; p2 = ixheaac_mult32_shl(sub_i, sinmu); p3 = ixheaac_mult32_shl(sub_r, sinmu); p4 = add_i >> 1; temp = ixheaac_sub32_sat(xr_0, p1); temp1 = ixheaac_add32_sat(xi_0, p3); temp2 = ixheaac_sub32_sat(xi_0, p3); idx = re_arr_tab_sml_240_ptr[n++] << 1; op[idx] = ixheaac_add32_sat(x01_real, xr_2); op[idx + 1] = ixheaac_add32_sat(x_01_imag, xi_2); idx = re_arr_tab_sml_240_ptr[n++] << 1; op[idx] = ixheaac_add32_sat(temp, p2); op[idx + 1] = ixheaac_sub32_sat(temp2, p4); idx = re_arr_tab_sml_240_ptr[n++] << 1; op[idx] = ixheaac_sub32_sat(temp, p2); op[idx + 1] = ixheaac_sub32_sat(temp1, p4); fft3outptr += 2; } }