/src/libxaac/decoder/ixheaacd_aac_imdct.c

Source
/******************************************************************************
 *                                                                            *
 * 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;
  }
}

Coverage Report

Created: 2025-10-13 07:04