xref: /aosp_15_r20/external/libxaac/encoder/ixheaace_mps_polyphase.c (revision 15dc779a375ca8b5125643b829a8aa4b70d7f451)

/******************************************************************************
 *                                                                            *
 * Copyright (C) 2023 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 <string.h>
#include <math.h>
#include "ixheaac_type_def.h"
#include "ixheaace_mps_common_fix.h"
#include "ixheaace_mps_defines.h"
#include "ixheaace_mps_common_define.h"
#include "ixheaace_bitbuffer.h"

#include "ixheaace_mps_struct_def.h"
#include "ixheaace_mps_sac_polyphase.h"
#include "ixheaace_mps_sac_hybfilter.h"
#include "ixheaace_mps_spatial_bitstream.h"
#include "ixheaace_mps_bitstream.h"
#include "ixheaace_mps_param_extract.h"
#include "ixheaace_mps_tree.h"
#include "ixheaace_mps_rom.h"

static VOID ixheaace_mps_bit_reversal(FLOAT32 *a, FLOAT32 *b) {
  FLOAT32 temp;
  temp = *a;
  *a = *b;
  *b = temp;
}

static VOID ixheaace_mps_bit_add_sub1(FLOAT32 *a, FLOAT32 *b) {
  FLOAT32 temp;
  temp = *a + *b;
  *b = *b - *a;
  *a = temp;
}

static VOID ixheaace_mps_bit_add_sub2(FLOAT32 *a, FLOAT32 *b) {
  FLOAT32 temp;
  temp = *a + *b;
  *b = *a - *b;
  *a = temp;
}

static VOID ixheaace_mps_fft32(FLOAT32 *subband) {
  FLOAT32 val_0, val_1, val_2, val_3, val_4, val_5, val_6, val_7, val_8, val_9, val_10, val_11,
      val_12, val_13, val_14, val_15;
  FLOAT32 val_1_0, val_1_1, val_1_2, val_1_3, val_1_4, val_1_5, val_1_6, val_1_7, val_1_8,
      val_1_9, val_1_10, val_1_11, val_1_12, val_1_13, val_1_14, val_1_15;
  FLOAT32 val_2_0, val_2_1, val_2_2, val_2_3, val_2_4, val_2_5, val_2_6, val_2_7, val_2_8,
      val_2_9, val_2_10, val_2_11, val_2_12, val_2_13, val_2_14, val_2_15;
  FLOAT32 val_3_0, val_3_1, val_3_2, val_3_3, val_3_4, val_3_5, val_3_6, val_3_7, val_3_8,
      val_3_9, val_3_10, val_3_11, val_3_12, val_3_13, val_3_14, val_3_15;
  FLOAT32 val_3_16, val_3_17, val_3_18, val_3_19, val_3_20, val_3_21, val_3_22, val_3_23,
      val_3_24, val_3_25, val_3_26, val_3_27, val_3_28, val_3_29, val_3_30, val_3_31;
  FLOAT32 val_4_0, val_4_1, val_4_2, val_4_3, val_4_4, val_4_5, val_4_6, val_4_7, val_4_8,
      val_4_9, val_4_10, val_4_11, val_4_12, val_4_13, val_4_14, val_4_15;

  val_2_0 = subband[2] - subband[34];
  val_2_1 = subband[3] - subband[35];
  val_3_0 = subband[0] + subband[32];
  val_3_1 = subband[1] + subband[33];
  val_3_2 = subband[2] + subband[34];
  val_3_3 = subband[3] + subband[35];

  val_2_2 = subband[6] - subband[38];
  val_2_3 = subband[7] - subband[39];
  val_3_4 = subband[4] + subband[36];
  val_3_5 = subband[5] + subband[37];
  val_3_6 = subband[6] + subband[38];
  val_3_7 = subband[7] + subband[39];

  val_2_4 = subband[10] - subband[42];
  val_2_5 = subband[11] - subband[43];
  val_3_8 = subband[8] + subband[40];
  val_3_9 = subband[9] + subband[41];
  val_3_10 = subband[10] + subband[42];
  val_3_11 = subband[11] + subband[43];

  val_2_6 = subband[14] - subband[46];
  val_2_7 = subband[15] - subband[47];
  val_3_12 = subband[12] + subband[44];
  val_3_13 = subband[13] + subband[45];
  val_3_14 = subband[14] + subband[46];
  val_3_15 = subband[15] + subband[47];

  val_2_8 = subband[18] - subband[50];
  val_2_9 = subband[19] - subband[51];
  val_3_16 = subband[16] + subband[48];
  val_3_17 = subband[17] + subband[49];
  val_3_18 = subband[18] + subband[50];
  val_3_19 = subband[19] + subband[51];

  val_2_10 = subband[22] - subband[54];
  val_2_11 = subband[23] - subband[55];
  val_3_20 = subband[20] + subband[52];
  val_3_21 = subband[21] + subband[53];
  val_3_22 = subband[22] + subband[54];
  val_3_23 = subband[23] + subband[55];

  val_2_12 = subband[26] - subband[58];
  val_2_13 = subband[27] - subband[59];
  val_3_24 = subband[24] + subband[56];
  val_3_25 = subband[25] + subband[57];
  val_3_26 = subband[26] + subband[58];
  val_3_27 = subband[27] + subband[59];

  val_2_14 = subband[30] - subband[62];
  val_2_15 = subband[31] - subband[63];
  val_3_28 = subband[28] + subband[60];
  val_3_29 = subband[29] + subband[61];
  val_3_30 = subband[30] + subband[62];
  val_3_31 = subband[31] + subband[63];

  val_4_1 = -(val_2_0 + val_2_14);
  val_4_2 = val_2_0 - val_2_14;
  val_4_0 = val_2_1 + val_2_15;
  val_4_3 = val_2_1 - val_2_15;
  val_4_5 = -(val_2_2 + val_2_12);
  val_4_6 = val_2_2 - val_2_12;
  val_4_4 = val_2_3 + val_2_13;
  val_4_7 = val_2_3 - val_2_13;
  val_4_9 = -(val_2_4 + val_2_10);
  val_4_10 = val_2_4 - val_2_10;
  val_4_8 = val_2_5 + val_2_11;
  val_4_11 = val_2_5 - val_2_11;
  val_4_13 = -(val_2_6 + val_2_8);
  val_4_14 = val_2_6 - val_2_8;
  val_4_12 = val_2_7 + val_2_9;
  val_4_15 = val_2_7 - val_2_9;

  val_2_0 = val_4_0 * fft_c[3] + val_4_4 * fft_c[2] + val_4_8 * fft_c[1] + val_4_12 * fft_c[0];

  val_2_4 = val_4_0 * fft_c[2] + val_4_4 * fft_c[0] + val_4_8 * fft_c[3] - val_4_12 * fft_c[1];
  val_2_8 = val_4_0 * fft_c[1] + val_4_4 * fft_c[3] - val_4_8 * fft_c[0] + val_4_12 * fft_c[2];

  val_2_12 = val_4_0 * fft_c[0] - val_4_4 * fft_c[1] + val_4_8 * fft_c[2] - val_4_12 * fft_c[3];

  val_2_1 = val_4_1 * fft_c[3] + val_4_5 * fft_c[2] + val_4_9 * fft_c[1] + val_4_13 * fft_c[0];

  val_2_5 = val_4_1 * fft_c[2] + val_4_5 * fft_c[0] + val_4_9 * fft_c[3] - val_4_13 * fft_c[1];
  val_2_9 = val_4_1 * fft_c[1] + val_4_5 * fft_c[3] - val_4_9 * fft_c[0] + val_4_13 * fft_c[2];

  val_2_13 = val_4_1 * fft_c[0] - val_4_5 * fft_c[1] + val_4_9 * fft_c[2] - val_4_13 * fft_c[3];

  val_2_2 = val_4_2 * fft_c[0] + val_4_6 * fft_c[1] + val_4_10 * fft_c[2] + val_4_14 * fft_c[3];

  val_2_6 = val_4_2 * fft_c[1] - val_4_6 * fft_c[3] - val_4_10 * fft_c[0] - val_4_14 * fft_c[2];

  val_2_10 = val_4_2 * fft_c[2] - val_4_6 * fft_c[0] + val_4_10 * fft_c[3] + val_4_14 * fft_c[1];

  val_2_14 = val_4_2 * fft_c[3] - val_4_6 * fft_c[2] + val_4_10 * fft_c[1] - val_4_14 * fft_c[0];

  val_2_3 = val_4_3 * fft_c[0] + val_4_7 * fft_c[1] + val_4_11 * fft_c[2] + val_4_15 * fft_c[3];

  val_2_7 = val_4_3 * fft_c[1] - val_4_7 * fft_c[3] - val_4_11 * fft_c[0] - val_4_15 * fft_c[2];

  val_2_11 = val_4_3 * fft_c[2] - val_4_7 * fft_c[0] + val_4_11 * fft_c[3] + val_4_15 * fft_c[1];

  val_2_15 = val_4_3 * fft_c[3] - val_4_7 * fft_c[2] + val_4_11 * fft_c[1] - val_4_15 * fft_c[0];

  val_4_0 = val_2_0 + val_2_2;
  val_4_14 = val_2_0 - val_2_2;
  val_4_1 = val_2_1 + val_2_3;
  val_4_15 = val_2_1 - val_2_3;
  val_4_2 = val_2_4 + val_2_6;
  val_4_12 = val_2_4 - val_2_6;
  val_4_3 = val_2_5 + val_2_7;
  val_4_13 = val_2_5 - val_2_7;
  val_4_4 = val_2_8 + val_2_10;
  val_4_10 = val_2_8 - val_2_10;
  val_4_5 = val_2_9 + val_2_11;
  val_4_11 = val_2_9 - val_2_11;
  val_4_6 = val_2_12 + val_2_14;
  val_4_8 = val_2_12 - val_2_14;
  val_4_7 = val_2_13 + val_2_15;
  val_4_9 = val_2_13 - val_2_15;

  val_1_0 = val_3_0 + val_3_16;
  val_1_1 = val_3_1 + val_3_17;
  val_1_2 = val_3_2 + val_3_18;
  val_1_3 = val_3_3 + val_3_19;
  val_1_4 = val_3_4 + val_3_20;
  val_1_5 = val_3_5 + val_3_21;
  val_1_6 = val_3_6 + val_3_22;
  val_1_7 = val_3_7 + val_3_23;
  val_1_8 = val_3_8 + val_3_24;
  val_1_9 = val_3_9 + val_3_25;
  val_1_10 = val_3_10 + val_3_26;
  val_1_11 = val_3_11 + val_3_27;
  val_1_12 = val_3_12 + val_3_28;
  val_1_13 = val_3_13 + val_3_29;
  val_1_14 = val_3_14 + val_3_30;
  val_1_15 = val_3_15 + val_3_31;

  val_0 = val_1_0 + val_1_8;
  val_2 = val_1_0 - val_1_8;
  val_1 = val_1_1 + val_1_9;
  val_3 = val_1_1 - val_1_9;
  val_4 = val_1_2 + val_1_10;
  val_6 = val_1_2 - val_1_10;
  val_5 = val_1_3 + val_1_11;
  val_7 = val_1_3 - val_1_11;
  val_8 = val_1_4 + val_1_12;
  val_10 = val_1_4 - val_1_12;
  val_9 = val_1_5 + val_1_13;
  val_11 = val_1_5 - val_1_13;
  val_12 = val_1_6 + val_1_14;
  val_14 = val_1_6 - val_1_14;
  val_13 = val_1_7 + val_1_15;
  val_15 = val_1_7 - val_1_15;

  val_2_0 = val_0 + val_8;
  val_2_4 = val_0 - val_8;
  val_2_1 = val_1 + val_9;
  val_2_5 = val_1 - val_9;
  val_2_8 = val_2 - val_11;
  val_2_10 = val_2 + val_11;
  val_2_9 = val_3 + val_10;
  val_2_11 = val_3 - val_10;
  val_2_2 = val_4 + val_12;
  val_2_7 = val_4 - val_12;
  val_2_3 = val_5 + val_13;
  val_2_6 = val_13 - val_5;

  val_1 = val_6 + val_14;
  val_2 = val_6 - val_14;
  val_0 = val_7 + val_15;
  val_3 = val_7 - val_15;

  val_2_12 = (val_0 + val_2) * MPS_INV_SQRT2;
  val_2_14 = (val_0 - val_2) * MPS_INV_SQRT2;
  val_2_13 = (val_3 - val_1) * MPS_INV_SQRT2;
  val_2_15 = (val_1 + val_3) * -MPS_INV_SQRT2;

  val_1_0 = val_3_0 - val_3_16;
  val_1_1 = val_3_1 - val_3_17;
  val_1_2 = val_3_2 - val_3_18;
  val_1_3 = val_3_3 - val_3_19;
  val_1_4 = val_3_4 - val_3_20;
  val_1_5 = val_3_5 - val_3_21;
  val_1_6 = val_3_6 - val_3_22;
  val_1_7 = val_3_7 - val_3_23;
  val_1_8 = val_3_8 - val_3_24;
  val_1_9 = val_3_9 - val_3_25;
  val_1_10 = val_3_10 - val_3_26;
  val_1_11 = val_3_11 - val_3_27;
  val_1_12 = val_3_12 - val_3_28;
  val_1_13 = val_3_13 - val_3_29;
  val_1_14 = val_3_14 - val_3_30;
  val_1_15 = val_3_15 - val_3_31;

  val_3_0 = val_2_0 + val_2_2;
  val_3_16 = val_2_0 - val_2_2;
  val_3_1 = val_2_1 + val_2_3;
  val_3_17 = val_2_1 - val_2_3;
  val_3_8 = val_2_4 - val_2_6;
  val_3_24 = val_2_4 + val_2_6;
  val_3_9 = val_2_5 - val_2_7;
  val_3_25 = val_2_5 + val_2_7;
  val_3_12 = val_2_8 + val_2_14;
  val_3_28 = val_2_8 - val_2_14;
  val_3_13 = val_2_9 + val_2_15;
  val_3_29 = val_2_9 - val_2_15;
  val_3_4 = val_2_10 + val_2_12;
  val_3_20 = val_2_10 - val_2_12;
  val_3_5 = val_2_11 + val_2_13;
  val_3_21 = val_2_11 - val_2_13;

  val_9 = (val_1_2 + val_1_14) * -MPS_COS_3PI_DIV8;
  val_10 = (val_1_2 - val_1_14) * MPS_COS_PI_DIV8;
  val_8 = (val_1_3 + val_1_15) * MPS_COS_3PI_DIV8;
  val_11 = (val_1_3 - val_1_15) * MPS_COS_PI_DIV8;
  val_5 = (val_1_4 + val_1_12) * -MPS_INV_SQRT2;
  val_6 = (val_1_4 - val_1_12) * MPS_INV_SQRT2;
  val_4 = (val_1_5 + val_1_13) * MPS_INV_SQRT2;
  val_7 = (val_1_5 - val_1_13) * MPS_INV_SQRT2;
  val_13 = (val_1_6 + val_1_10) * -MPS_COS_PI_DIV8;
  val_14 = (val_1_6 - val_1_10) * MPS_COS_3PI_DIV8;
  val_12 = (val_1_7 + val_1_11) * MPS_COS_PI_DIV8;
  val_15 = (val_1_7 - val_1_11) * MPS_COS_3PI_DIV8;

  val_1_2 = val_8 * MPS_SQRT2PLUS1 - val_12 * MPS_SQRT2MINUS1;
  val_1_3 = val_9 * MPS_SQRT2PLUS1 - val_13 * MPS_SQRT2MINUS1;
  val_1_4 = val_10 * MPS_SQRT2MINUS1 - val_14 * MPS_SQRT2PLUS1;
  val_1_5 = val_11 * MPS_SQRT2MINUS1 - val_15 * MPS_SQRT2PLUS1;

  val_8 += val_12;
  val_9 += val_13;
  val_10 += val_14;
  val_11 += val_15;
  val_1_6 = val_1_0 + val_4;
  val_1_10 = val_1_0 - val_4;
  val_1_7 = val_1_1 + val_5;
  val_1_11 = val_1_1 - val_5;

  val_1_12 = val_6 - val_1_9;
  val_1_14 = val_6 + val_1_9;
  val_1_13 = val_1_8 + val_7;
  val_1_15 = val_1_8 - val_7;

  val_0 = val_1_6 - val_1_14;
  val_2 = val_1_6 + val_1_14;
  val_1 = val_1_7 + val_1_15;
  val_3 = val_1_7 - val_1_15;
  val_4 = val_1_10 + val_1_12;
  val_6 = val_1_10 - val_1_12;
  val_5 = val_1_11 + val_1_13;
  val_7 = val_1_11 - val_1_13;

  val_1_10 = val_8 + val_10;
  val_10 = val_8 - val_10;
  val_1_11 = val_9 + val_11;
  val_11 = val_9 - val_11;

  val_12 = val_1_2 + val_1_4;
  val_14 = val_1_2 - val_1_4;
  val_13 = val_1_3 + val_1_5;
  val_15 = val_1_3 - val_1_5;

  val_3_2 = val_2 + val_1_10;
  val_3_18 = val_2 - val_1_10;
  val_3_3 = val_3 + val_1_11;
  val_3_19 = val_3 - val_1_11;
  val_3_6 = val_0 + val_12;
  val_3_22 = val_0 - val_12;
  val_3_7 = val_1 + val_13;
  val_3_23 = val_1 - val_13;
  val_3_14 = val_4 + val_10;
  val_3_30 = val_4 - val_10;
  val_3_15 = val_5 + val_11;
  val_3_31 = val_5 - val_11;
  val_3_10 = val_6 + val_14;
  val_3_26 = val_6 - val_14;
  val_3_11 = val_7 + val_15;
  val_3_27 = val_7 - val_15;

  val_1_0 = subband[0] - subband[32];
  val_1_1 = subband[1] - subband[33];
  val_1_2 = subband[4] - subband[36];
  val_1_3 = subband[5] - subband[37];
  val_1_4 = subband[8] - subband[40];
  val_1_5 = subband[9] - subband[41];
  val_1_6 = subband[12] - subband[44];
  val_1_7 = subband[13] - subband[45];
  val_1_8 = subband[16] - subband[48];
  val_1_9 = subband[17] - subband[49];
  val_1_10 = subband[20] - subband[52];
  val_1_11 = subband[21] - subband[53];
  val_1_12 = subband[24] - subband[56];
  val_1_13 = subband[25] - subband[57];
  val_1_14 = subband[28] - subband[60];
  val_1_15 = subband[29] - subband[61];

  val_9 = (val_1_2 + val_1_14) * -MPS_COS_3PI_DIV8;
  val_10 = (val_1_2 - val_1_14) * MPS_COS_PI_DIV8;
  val_8 = (val_1_3 + val_1_15) * MPS_COS_3PI_DIV8;
  val_11 = (val_1_3 - val_1_15) * MPS_COS_PI_DIV8;
  val_5 = (val_1_4 + val_1_12) * -MPS_INV_SQRT2;
  val_6 = (val_1_4 - val_1_12) * MPS_INV_SQRT2;
  val_4 = (val_1_5 + val_1_13) * MPS_INV_SQRT2;
  val_7 = (val_1_5 - val_1_13) * MPS_INV_SQRT2;
  val_13 = (val_1_6 + val_1_10) * -MPS_COS_PI_DIV8;
  val_14 = (val_1_6 - val_1_10) * MPS_COS_3PI_DIV8;
  val_12 = (val_1_7 + val_1_11) * MPS_COS_PI_DIV8;
  val_15 = (val_1_7 - val_1_11) * MPS_COS_3PI_DIV8;

  val_1_2 = val_8 * MPS_SQRT2PLUS1 - val_12 * MPS_SQRT2MINUS1;
  val_1_3 = val_9 * MPS_SQRT2PLUS1 - val_13 * MPS_SQRT2MINUS1;
  val_1_4 = val_10 * MPS_SQRT2MINUS1 - val_14 * MPS_SQRT2PLUS1;
  val_1_5 = val_11 * MPS_SQRT2MINUS1 - val_15 * MPS_SQRT2PLUS1;

  val_8 += val_12;
  val_9 += val_13;
  val_10 += val_14;
  val_11 += val_15;
  val_1_6 = val_1_0 + val_4;
  val_1_10 = val_1_0 - val_4;
  val_1_7 = val_1_1 + val_5;
  val_1_11 = val_1_1 - val_5;

  val_1_12 = val_6 - val_1_9;
  val_1_14 = val_6 + val_1_9;
  val_1_13 = val_1_8 + val_7;
  val_1_15 = val_1_8 - val_7;

  val_0 = val_1_6 - val_1_14;
  val_2 = val_1_6 + val_1_14;
  val_1 = val_1_7 + val_1_15;
  val_3 = val_1_7 - val_1_15;
  val_4 = val_1_10 + val_1_12;
  val_6 = val_1_10 - val_1_12;
  val_5 = val_1_11 + val_1_13;
  val_7 = val_1_11 - val_1_13;

  val_1_10 = val_8 + val_10;
  val_10 = val_8 - val_10;
  val_1_11 = val_9 + val_11;
  val_11 = val_9 - val_11;

  val_12 = val_1_2 + val_1_4;
  val_14 = val_1_2 - val_1_4;
  val_13 = val_1_3 + val_1_5;
  val_15 = val_1_3 - val_1_5;

  val_1_0 = val_2 + val_1_10;
  val_1_8 = val_2 - val_1_10;
  val_1_1 = val_3 + val_1_11;
  val_1_9 = val_3 - val_1_11;
  val_1_2 = val_0 + val_12;
  val_1_10 = val_0 - val_12;
  val_1_3 = val_1 + val_13;
  val_1_11 = val_1 - val_13;
  val_1_6 = val_4 + val_10;
  val_1_14 = val_4 - val_10;
  val_1_7 = val_5 + val_11;
  val_1_15 = val_5 - val_11;
  val_1_4 = val_6 + val_14;
  val_1_12 = val_6 - val_14;
  val_1_5 = val_7 + val_15;
  val_1_13 = val_7 - val_15;

  *subband++ = val_3_0;
  *subband++ = val_3_1;
  *subband++ = val_1_0 + val_4_0;
  *subband++ = val_1_1 + val_4_1;
  *subband++ = val_3_2;
  *subband++ = val_3_3;
  *subband++ = val_1_2 + val_4_2;
  *subband++ = val_1_3 + val_4_3;
  *subband++ = val_3_4;
  *subband++ = val_3_5;
  *subband++ = val_1_4 + val_4_4;
  *subband++ = val_1_5 + val_4_5;
  *subband++ = val_3_6;
  *subband++ = val_3_7;
  *subband++ = val_1_6 + val_4_6;
  *subband++ = val_1_7 + val_4_7;
  *subband++ = val_3_8;
  *subband++ = val_3_9;
  *subband++ = val_1_8 + val_4_8;
  *subband++ = val_1_9 + val_4_9;
  *subband++ = val_3_10;
  *subband++ = val_3_11;
  *subband++ = val_1_10 + val_4_10;
  *subband++ = val_1_11 + val_4_11;
  *subband++ = val_3_12;
  *subband++ = val_3_13;
  *subband++ = val_1_12 + val_4_12;
  *subband++ = val_1_13 + val_4_13;
  *subband++ = val_3_14;
  *subband++ = val_3_15;
  *subband++ = val_1_14 + val_4_14;
  *subband++ = val_1_15 + val_4_15;
  *subband++ = val_3_16;
  *subband++ = val_3_17;
  *subband++ = val_1_0 - val_4_0;
  *subband++ = val_1_1 - val_4_1;
  *subband++ = val_3_18;
  *subband++ = val_3_19;
  *subband++ = val_1_2 - val_4_2;
  *subband++ = val_1_3 - val_4_3;
  *subband++ = val_3_20;
  *subband++ = val_3_21;
  *subband++ = val_1_4 - val_4_4;
  *subband++ = val_1_5 - val_4_5;
  *subband++ = val_3_22;
  *subband++ = val_3_23;
  *subband++ = val_1_6 - val_4_6;
  *subband++ = val_1_7 - val_4_7;
  *subband++ = val_3_24;
  *subband++ = val_3_25;
  *subband++ = val_1_8 - val_4_8;
  *subband++ = val_1_9 - val_4_9;
  *subband++ = val_3_26;
  *subband++ = val_3_27;
  *subband++ = val_1_10 - val_4_10;
  *subband++ = val_1_11 - val_4_11;
  *subband++ = val_3_28;
  *subband++ = val_3_29;
  *subband++ = val_1_12 - val_4_12;
  *subband++ = val_1_13 - val_4_13;
  *subband++ = val_3_30;
  *subband++ = val_3_31;
  *subband++ = val_1_14 - val_4_14;
  *subband++ = val_1_15 - val_4_15;
}

static VOID ixheaace_mps_cos_mod(
    FLOAT32 *subband, ixheaace_mps_sac_pstr_qmf_synth_filter_bank pstr_qmf_synth_filter_bank) {
  WORD32 idx, length;
  FLOAT32 weight_imag, weight_real;
  FLOAT32 real_1, real_2;
  FLOAT32 imag_1, imag_2;
  FLOAT32 accu1, accu2;

  length = 32;

  for (idx = 0; idx < (length / 2); idx++) {
    real_1 = subband[2 * idx];
    imag_2 = subband[2 * idx + 1];
    real_2 = subband[2 * length - 2 - 2 * idx];
    imag_1 = subband[2 * length - 1 - 2 * idx];

    weight_imag = pstr_qmf_synth_filter_bank->sin_twiddle[idx];
    weight_real = pstr_qmf_synth_filter_bank->cos_twiddle[idx];

    accu1 = imag_1 * weight_imag + real_1 * weight_real;
    accu2 = imag_1 * weight_real - real_1 * weight_imag;

    subband[2 * idx] = accu1;
    subband[2 * idx + 1] = accu2;

    weight_imag = pstr_qmf_synth_filter_bank->sin_twiddle[length - 1 - idx];
    weight_real = pstr_qmf_synth_filter_bank->cos_twiddle[length - 1 - idx];

    accu1 = imag_2 * weight_imag + real_2 * weight_real;
    accu2 = imag_2 * weight_real - real_2 * weight_imag;

    subband[2 * length - 2 - 2 * idx] = accu1;
    subband[2 * length - 1 - 2 * idx] = accu2;
  }

  ixheaace_mps_fft32(subband);

  weight_imag = pstr_qmf_synth_filter_bank->alt_sin_twiddle[0];
  weight_real = pstr_qmf_synth_filter_bank->alt_sin_twiddle[length];

  for (idx = 0; idx < length / 2; idx++) {
    real_1 = subband[2 * idx];
    imag_1 = subband[2 * idx + 1];
    real_2 = subband[2 * length - 2 - 2 * idx];
    imag_2 = subband[2 * length - 1 - 2 * idx];

    accu1 = real_1 * weight_real + imag_1 * weight_imag;
    accu2 = real_1 * weight_imag - imag_1 * weight_real;

    subband[2 * idx] = accu1;
    subband[2 * length - 1 - 2 * idx] = accu2;

    weight_imag = pstr_qmf_synth_filter_bank->alt_sin_twiddle[idx + 1];
    weight_real = pstr_qmf_synth_filter_bank->alt_sin_twiddle[length - 1 - idx];

    accu1 = real_2 * weight_imag + imag_2 * weight_real;
    accu2 = real_2 * weight_real - imag_2 * weight_imag;

    subband[2 * length - 2 - 2 * idx] = accu1;
    subband[2 * idx + 1] = accu2;
  }
}

static VOID ixheaace_mps_sin_mod(
    FLOAT32 *subband, ixheaace_mps_sac_pstr_qmf_synth_filter_bank pstr_qmf_synth_filter_bank) {
  WORD32 idx, length;
  FLOAT32 weight_real, weight_imag;
  FLOAT32 real_1, imag_1, real_2, imag_2;
  FLOAT32 accu1, accu2;

  length = 32;

  for (idx = 0; idx < length / 2; idx++) {
    real_1 = subband[2 * idx];
    imag_2 = subband[2 * idx + 1];
    real_2 = subband[2 * length - 2 - 2 * idx];
    imag_1 = subband[2 * length - 1 - 2 * idx];

    weight_real = pstr_qmf_synth_filter_bank->sin_twiddle[idx];
    weight_imag = pstr_qmf_synth_filter_bank->cos_twiddle[idx];

    accu1 = imag_1 * weight_imag + real_1 * weight_real;
    accu2 = imag_1 * weight_real - real_1 * weight_imag;

    subband[2 * idx + 1] = accu1;
    subband[2 * idx] = accu2;

    weight_real = pstr_qmf_synth_filter_bank->sin_twiddle[length - 1 - idx];
    weight_imag = pstr_qmf_synth_filter_bank->cos_twiddle[length - 1 - idx];

    accu1 = imag_2 * weight_imag + real_2 * weight_real;
    accu2 = imag_2 * weight_real - real_2 * weight_imag;

    subband[2 * length - 1 - 2 * idx] = accu1;
    subband[2 * length - 2 - 2 * idx] = accu2;
  }

  ixheaace_mps_fft32(subband);

  weight_imag = pstr_qmf_synth_filter_bank->alt_sin_twiddle[0];
  weight_real = pstr_qmf_synth_filter_bank->alt_sin_twiddle[length];

  for (idx = 0; idx < length / 2; idx++) {
    real_1 = subband[2 * idx];
    imag_1 = subband[2 * idx + 1];
    real_2 = subband[2 * length - 2 - 2 * idx];
    imag_2 = subband[2 * length - 1 - 2 * idx];

    accu1 = -(real_1 * weight_real + imag_1 * weight_imag);
    accu2 = -(real_1 * weight_imag - imag_1 * weight_real);

    subband[2 * length - 1 - 2 * idx] = accu1;
    subband[2 * idx] = accu2;

    weight_imag = pstr_qmf_synth_filter_bank->alt_sin_twiddle[idx + 1];
    weight_real = pstr_qmf_synth_filter_bank->alt_sin_twiddle[length - 1 - idx];

    accu1 = -(real_2 * weight_imag + imag_2 * weight_real);
    accu2 = -(real_2 * weight_real - imag_2 * weight_imag);

    subband[2 * idx + 1] = accu1;
    subband[2 * length - 2 - 2 * idx] = accu2;
  }
}

static VOID ixheaace_mps_inverse_modulation(
    FLOAT32 *qmf_real, FLOAT32 *qmf_imag,
    ixheaace_mps_sac_pstr_qmf_synth_filter_bank pstr_qmf_synth_filter_bank) {
  WORD32 idx, no_synthesis_channels, length;

  FLOAT32 r1, i1, r2, i2;

  no_synthesis_channels = NUM_QMF_BANDS;

  length = no_synthesis_channels / 2;

  ixheaace_mps_cos_mod(qmf_real, pstr_qmf_synth_filter_bank);
  ixheaace_mps_sin_mod(qmf_imag, pstr_qmf_synth_filter_bank);

  for (idx = 0; idx < length; idx++) {
    r1 = qmf_real[idx];
    i2 = qmf_imag[no_synthesis_channels - 1 - idx];
    r2 = qmf_real[no_synthesis_channels - idx - 1];
    i1 = qmf_imag[idx];

    qmf_real[idx] = (r1 - i1);
    qmf_imag[no_synthesis_channels - 1 - idx] = -(r1 + i1);
    qmf_real[no_synthesis_channels - idx - 1] = (r2 - i2);
    qmf_imag[idx] = -(r2 + i2);
  }
}

static VOID ixheaace_mps_fct3_4(FLOAT32 *subband) {
  FLOAT32 val_00, val_01, val_10, val_11;

  subband[1] *= MPS_INV_SQRT2;

  val_00 = subband[0] + subband[1];
  val_01 = subband[0] - subband[1];

  val_10 = subband[2] * MPS_COS_6_PI_BY_16 + subband[3] * MPS_SIN_6_PI_BY_16;
  val_11 = subband[2] * MPS_SIN_6_PI_BY_16 - subband[3] * MPS_COS_6_PI_BY_16;

  subband[0] = val_00 + val_10;
  subband[3] = val_00 - val_10;
  subband[1] = val_01 + val_11;
  subband[2] = val_01 - val_11;
}

static VOID ixheaace_mps_fst3_4r(FLOAT32 *subband) {
  FLOAT32 val_00, val_01, val_10, val_11;

  subband[2] *= MPS_INV_SQRT2;

  val_00 = subband[3] + subband[2];
  val_01 = subband[3] - subband[2];

  val_10 = subband[1] * MPS_COS_6_PI_BY_16 + subband[0] * MPS_SIN_6_PI_BY_16;
  val_11 = subband[0] * MPS_COS_6_PI_BY_16 - subband[1] * MPS_SIN_6_PI_BY_16;

  subband[3] = val_00 + val_10;
  subband[0] = val_10 - val_00;
  subband[2] = val_11 - val_01;
  subband[1] = val_11 + val_01;
}

static VOID ixheaace_mps_fct4_4r(FLOAT32 *subband) {
  FLOAT32 val_00, val_01, val_10, val_11;
  subband[1] *= MPS_INV_SQRT2;

  val_00 = subband[0] + subband[1];
  val_01 = subband[0] - subband[1];

  subband[2] *= MPS_INV_SQRT2;

  val_11 = subband[3] - subband[2];
  val_10 = subband[3] + subband[2];

  subband[3] = val_00 * MPS_COS_PI_BY_16 + val_10 * MPS_SIN_PI_BY_16;
  subband[0] = val_00 * MPS_SIN_PI_BY_16 - val_10 * MPS_COS_PI_BY_16;

  subband[2] = val_01 * MPS_COS_3_PI_BY_16 - val_11 * MPS_SIN_3_PI_BY_16;
  subband[1] = val_01 * MPS_SIN_3_PI_BY_16 + val_11 * MPS_COS_3_PI_BY_16;
}

static VOID ixheaace_mps_fst4_4(FLOAT32 *subband) {
  FLOAT32 val_00, val_01, val_10, val_11;

  subband[1] *= MPS_INV_SQRT2;

  val_10 = subband[0] + subband[1];
  val_11 = subband[0] - subband[1];

  subband[2] *= MPS_INV_SQRT2;

  val_01 = subband[3] - subband[2];
  val_00 = subband[3] + subband[2];

  subband[0] = val_00 * MPS_COS_PI_BY_16 + val_10 * MPS_SIN_PI_BY_16;
  subband[3] = val_10 * MPS_COS_PI_BY_16 - val_00 * MPS_SIN_PI_BY_16;

  subband[1] = val_11 * MPS_SIN_3_PI_BY_16 - val_01 * MPS_COS_3_PI_BY_16;
  subband[2] = val_01 * MPS_SIN_3_PI_BY_16 + val_11 * MPS_COS_3_PI_BY_16;
}

static VOID ixheaace_mps_fct3_64(FLOAT32 *subband) {
  WORD32 idx;
  const FLOAT32 *t_ptr;
  FLOAT32 cos_val, sine_val;
  FLOAT32 xp;

  ixheaace_mps_bit_reversal(&subband[1], &subband[32]);
  ixheaace_mps_bit_reversal(&subband[2], &subband[16]);
  ixheaace_mps_bit_reversal(&subband[3], &subband[48]);
  ixheaace_mps_bit_reversal(&subband[4], &subband[8]);
  ixheaace_mps_bit_reversal(&subband[5], &subband[40]);
  ixheaace_mps_bit_reversal(&subband[6], &subband[24]);
  ixheaace_mps_bit_reversal(&subband[7], &subband[56]);
  ixheaace_mps_bit_reversal(&subband[9], &subband[36]);
  ixheaace_mps_bit_reversal(&subband[10], &subband[20]);
  ixheaace_mps_bit_reversal(&subband[11], &subband[52]);
  ixheaace_mps_bit_reversal(&subband[13], &subband[44]);
  ixheaace_mps_bit_reversal(&subband[14], &subband[28]);
  ixheaace_mps_bit_reversal(&subband[15], &subband[60]);
  ixheaace_mps_bit_reversal(&subband[17], &subband[34]);
  ixheaace_mps_bit_reversal(&subband[19], &subband[50]);
  ixheaace_mps_bit_reversal(&subband[21], &subband[42]);
  ixheaace_mps_bit_reversal(&subband[22], &subband[26]);
  ixheaace_mps_bit_reversal(&subband[23], &subband[58]);
  ixheaace_mps_bit_reversal(&subband[25], &subband[38]);
  ixheaace_mps_bit_reversal(&subband[27], &subband[54]);
  ixheaace_mps_bit_reversal(&subband[29], &subband[46]);
  ixheaace_mps_bit_reversal(&subband[31], &subband[62]);
  ixheaace_mps_bit_reversal(&subband[35], &subband[49]);
  ixheaace_mps_bit_reversal(&subband[37], &subband[41]);
  ixheaace_mps_bit_reversal(&subband[39], &subband[57]);
  ixheaace_mps_bit_reversal(&subband[43], &subband[53]);
  ixheaace_mps_bit_reversal(&subband[47], &subband[61]);
  ixheaace_mps_bit_reversal(&subband[55], &subband[59]);

  ixheaace_mps_bit_add_sub1(&subband[33], &subband[62]);
  ixheaace_mps_bit_add_sub1(&subband[34], &subband[60]);
  ixheaace_mps_bit_add_sub1(&subband[35], &subband[61]);
  ixheaace_mps_bit_add_sub1(&subband[36], &subband[56]);
  ixheaace_mps_bit_add_sub1(&subband[37], &subband[57]);
  ixheaace_mps_bit_add_sub1(&subband[38], &subband[58]);
  ixheaace_mps_bit_add_sub1(&subband[39], &subband[59]);
  ixheaace_mps_bit_add_sub1(&subband[40], &subband[48]);
  ixheaace_mps_bit_add_sub1(&subband[41], &subband[49]);
  ixheaace_mps_bit_add_sub1(&subband[42], &subband[50]);
  ixheaace_mps_bit_add_sub1(&subband[43], &subband[51]);
  ixheaace_mps_bit_add_sub1(&subband[44], &subband[52]);
  ixheaace_mps_bit_add_sub1(&subband[45], &subband[53]);
  ixheaace_mps_bit_add_sub1(&subband[46], &subband[54]);
  ixheaace_mps_bit_add_sub1(&subband[47], &subband[55]);
  ixheaace_mps_bit_add_sub1(&subband[17], &subband[30]);
  ixheaace_mps_bit_add_sub1(&subband[18], &subband[28]);
  ixheaace_mps_bit_add_sub1(&subband[19], &subband[29]);
  ixheaace_mps_bit_add_sub1(&subband[20], &subband[24]);
  ixheaace_mps_bit_add_sub1(&subband[21], &subband[25]);
  ixheaace_mps_bit_add_sub1(&subband[22], &subband[26]);
  ixheaace_mps_bit_add_sub1(&subband[23], &subband[27]);
  ixheaace_mps_bit_add_sub1(&subband[9], &subband[14]);
  ixheaace_mps_bit_add_sub1(&subband[10], &subband[12]);
  ixheaace_mps_bit_add_sub1(&subband[11], &subband[13]);
  ixheaace_mps_bit_add_sub1(&subband[41], &subband[46]);
  ixheaace_mps_bit_add_sub1(&subband[42], &subband[44]);
  ixheaace_mps_bit_add_sub1(&subband[43], &subband[45]);
  ixheaace_mps_bit_add_sub1(&subband[54], &subband[49]);
  ixheaace_mps_bit_add_sub1(&subband[52], &subband[50]);
  ixheaace_mps_bit_add_sub1(&subband[53], &subband[51]);
  ixheaace_mps_bit_add_sub1(&subband[5], &subband[6]);
  ixheaace_mps_bit_add_sub1(&subband[21], &subband[22]);
  ixheaace_mps_bit_add_sub1(&subband[26], &subband[25]);
  ixheaace_mps_bit_add_sub1(&subband[37], &subband[38]);
  ixheaace_mps_bit_add_sub1(&subband[58], &subband[57]);

  ixheaace_mps_fct3_4(subband);
  ixheaace_mps_fct4_4r(subband + 4);
  ixheaace_mps_fct3_4(subband + 8);
  ixheaace_mps_fst3_4r(subband + 12);
  ixheaace_mps_fct3_4(subband + 16);
  ixheaace_mps_fct4_4r(subband + 20);
  ixheaace_mps_fst4_4(subband + 24);
  ixheaace_mps_fst3_4r(subband + 28);
  ixheaace_mps_fct3_4(subband + 32);
  ixheaace_mps_fct4_4r(subband + 36);
  ixheaace_mps_fct3_4(subband + 40);
  ixheaace_mps_fst3_4r(subband + 44);
  ixheaace_mps_fct3_4(subband + 48);
  ixheaace_mps_fst3_4r(subband + 52);
  ixheaace_mps_fst4_4(subband + 56);
  ixheaace_mps_fst3_4r(subband + 60);

  for (idx = 0; idx < 4; idx++) {
    ixheaace_mps_bit_add_sub2(&subband[idx], &subband[7 - idx]);
    ixheaace_mps_bit_add_sub2(&subband[16 + idx], &subband[23 - idx]);
    ixheaace_mps_bit_add_sub1(&subband[31 - idx], &subband[24 + idx]);
    ixheaace_mps_bit_add_sub2(&subband[32 + idx], &subband[39 - idx]);
    ixheaace_mps_bit_add_sub1(&subband[63 - idx], &subband[56 + idx]);
  }
  t_ptr = trig_data_fct4_8;

  for (idx = 0; idx < 4; idx++) {
    cos_val = *t_ptr++;
    sine_val = *t_ptr++;

    xp = subband[8 + idx] * cos_val + subband[15 - idx] * sine_val;
    subband[8 + idx] = subband[8 + idx] * sine_val - subband[15 - idx] * cos_val;
    subband[15 - idx] = xp;
    xp = subband[40 + idx] * cos_val + subband[47 - idx] * sine_val;
    subband[40 + idx] = subband[40 + idx] * sine_val - subband[47 - idx] * cos_val;
    subband[47 - idx] = xp;
    xp = subband[48 + idx] * sine_val + subband[55 - idx] * cos_val;
    subband[55 - idx] = subband[48 + idx] * cos_val - subband[55 - idx] * sine_val;
    subband[48 + idx] = xp;
  }

  for (idx = 0; idx < 8; idx++) {
    ixheaace_mps_bit_add_sub2(&subband[idx], &subband[15 - idx]);
    ixheaace_mps_bit_add_sub2(&subband[32 + idx], &subband[47 - idx]);
    ixheaace_mps_bit_add_sub1(&subband[63 - idx], &subband[48 + idx]);
  }

  t_ptr = trig_data_fct4_16;

  for (idx = 0; idx < 8; idx++) {
    cos_val = *t_ptr++;
    sine_val = *t_ptr++;

    xp = subband[16 + idx] * cos_val + subband[31 - idx] * sine_val;
    subband[16 + idx] = subband[16 + idx] * sine_val - subband[31 - idx] * cos_val;
    subband[31 - idx] = xp;
  }

  for (idx = 0; idx < 16; idx++) {
    ixheaace_mps_bit_add_sub2(&subband[idx], &subband[31 - idx]);
  }

  t_ptr = trig_data_fct4_32;

  for (idx = 0; idx < 16; idx++) {
    cos_val = *t_ptr++;
    sine_val = *t_ptr++;

    xp = subband[32 + idx] * cos_val + subband[63 - idx] * sine_val;
    subband[32 + idx] = subband[32 + idx] * sine_val - subband[63 - idx] * cos_val;
    subband[63 - idx] = xp;
  }

  for (idx = 0; idx < 32; idx++) {
    ixheaace_mps_bit_add_sub2(&subband[idx], &subband[63 - idx]);
  }
}

static VOID ixheaace_mps_fst3_64(FLOAT32 *subband) {
  WORD32 k;

  for (k = 0; k < 32; k++) {
    ixheaace_mps_bit_reversal(&subband[k], &subband[63 - k]);
  }

  ixheaace_mps_fct3_64(subband);

  for (k = 1; k < NUM_QMF_BANDS; k += 2) {
    subband[k] = -subband[k];
  }
}

static VOID ixheaace_mps_forward_modulation(const FLOAT32 *time_in, FLOAT32 *real_subband,
                                            FLOAT32 *imag_subband) {
  WORD32 idx;
  real_subband[0] = time_in[0];

  for (idx = 1; idx < NUM_QMF_BANDS; idx++) {
    real_subband[idx] = time_in[idx] - time_in[(NUM_QMF_BANDS * 2) - idx];
    imag_subband[idx - 1] = time_in[idx] + time_in[(NUM_QMF_BANDS * 2) - idx];
  }

  imag_subband[63] = time_in[NUM_QMF_BANDS];

  ixheaace_mps_fct3_64(real_subband);

  ixheaace_mps_fst3_64(imag_subband);
}

VOID ixheaace_mps_515_calculate_sbr_syn_filterbank(
    FLOAT32 *real_subband, FLOAT32 *imag_subband, FLOAT32 *time_sig, WORD32 channel,
    ixheaace_mps_sac_pstr_qmf_synth_filter_bank pstr_qmf_synth_filter_bank, WORD32 slots,
    FLOAT32 *sbr_qmf_states_synthesis_per) {
  WORD32 idx, subband, slot;
  const FLOAT32 *ptr_filter;
  FLOAT32 accumlate;
  FLOAT32 *synth_buf;

  synth_buf = &(sbr_qmf_states_synthesis_per[channel * QMF_FILTER_STATE_SYN_SIZE]);

  for (slot = 0; slot < slots; slot++) {
    ptr_filter = pstr_qmf_synth_filter_bank->p_filter;
    ixheaace_mps_inverse_modulation(real_subband + NUM_QMF_BANDS * slot,
                                    imag_subband + NUM_QMF_BANDS * slot,
                                    pstr_qmf_synth_filter_bank);

    for (subband = 0; subband < NUM_QMF_BANDS; subband++) {
      real_subband[NUM_QMF_BANDS * slot + subband] =
          real_subband[NUM_QMF_BANDS * slot + subband] * (-0.015625f);
      imag_subband[NUM_QMF_BANDS * slot + subband] =
          imag_subband[NUM_QMF_BANDS * slot + subband] * (-0.015625f);
    }

    for (subband = 0; subband < NUM_QMF_BANDS; subband++) {
      FLOAT32 new_sample;
      new_sample = imag_subband[NUM_QMF_BANDS * slot + 63 - subband];

      for (idx = 0; idx < 5; idx++) {
        accumlate = synth_buf[2 * idx * NUM_QMF_BANDS + subband] + (*ptr_filter++) * new_sample;
        synth_buf[2 * idx * NUM_QMF_BANDS + subband] = accumlate;
      }
    }

    for (idx = 0; idx < 5; idx++) {
      accumlate = synth_buf[2 * idx * NUM_QMF_BANDS + NUM_QMF_BANDS + (NUM_QMF_BANDS - 1)] +
                  (*ptr_filter++) * real_subband[NUM_QMF_BANDS * slot];
      synth_buf[2 * idx * NUM_QMF_BANDS + NUM_QMF_BANDS + (NUM_QMF_BANDS - 1)] = accumlate;
    }
    time_sig[0] = accumlate;

    ptr_filter -= 10;

    for (subband = 0; subband < 63; subband++) {
      FLOAT32 new_sample;

      new_sample = real_subband[NUM_QMF_BANDS * slot + 63 - subband];

      for (idx = 0; idx < 5; idx++) {
        accumlate = synth_buf[2 * idx * NUM_QMF_BANDS + NUM_QMF_BANDS + subband] +
                    (*--ptr_filter) * new_sample;
        synth_buf[2 * idx * NUM_QMF_BANDS + NUM_QMF_BANDS + subband] = accumlate;
      }

      time_sig[NUM_QMF_BANDS - 1 - subband] = accumlate;
    }

    time_sig += NUM_QMF_BANDS;

    memmove((synth_buf + NUM_QMF_BANDS), synth_buf, (640 - NUM_QMF_BANDS) * sizeof(FLOAT32));

    memset((synth_buf), 0, NUM_QMF_BANDS * sizeof(FLOAT32));
  }
}

VOID ixheaace_mps_515_calculate_ana_filterbank(
    ixheaace_mps_sac_sbr_encoder_ana_filter_bank *pstr_sbr_enc_ana_filter_bank, FLOAT32 *time_sig,
    FLOAT32 *real_subband, FLOAT32 *imag_subband,
    ixheaace_mps_sac_pstr_qmf_ana_filter_bank pstr_qmf_synth_filter_bank) {
  WORD32 idx;
  FLOAT32 temp_array[128];
  const FLOAT32 *ptr_filter = pstr_qmf_synth_filter_bank->p_filter;
  FLOAT32 accumlate;

  for (idx = 0; idx < 576; idx++) {
    pstr_sbr_enc_ana_filter_bank->x[idx] = pstr_sbr_enc_ana_filter_bank->x[idx + NUM_QMF_BANDS];
  }
  for (idx = 0; idx < NUM_QMF_BANDS; idx++) {
    pstr_sbr_enc_ana_filter_bank->x[idx + 576] = time_sig[idx];
  }

  for (idx = 0; idx < NUM_QMF_BANDS; idx++) {
    accumlate = 0.0f;
    accumlate += (*ptr_filter++ * pstr_sbr_enc_ana_filter_bank->x[idx]);
    accumlate += (*ptr_filter++ * pstr_sbr_enc_ana_filter_bank->x[idx + (NUM_QMF_BANDS * 2)]);
    accumlate += (*ptr_filter++ * pstr_sbr_enc_ana_filter_bank->x[idx + (NUM_QMF_BANDS * 4)]);
    accumlate += (*ptr_filter++ * pstr_sbr_enc_ana_filter_bank->x[idx + (NUM_QMF_BANDS * 6)]);
    accumlate += (*ptr_filter++ * pstr_sbr_enc_ana_filter_bank->x[idx + (NUM_QMF_BANDS * 8)]);
    temp_array[127 - idx] = accumlate;
  }

  accumlate = 0.0f;
  for (idx = 0; idx < 5; idx++) {
    accumlate +=
        (*ptr_filter++ * pstr_sbr_enc_ana_filter_bank->x[127 + idx * (NUM_QMF_BANDS * 2)]);
  }
  temp_array[0] = accumlate;

  ptr_filter -= 10;
  for (idx = 0; idx < 63; idx++) {
    accumlate = 0.0f;
    accumlate += (*--ptr_filter * pstr_sbr_enc_ana_filter_bank->x[NUM_QMF_BANDS + idx]);
    accumlate += (*--ptr_filter *
                  pstr_sbr_enc_ana_filter_bank->x[NUM_QMF_BANDS + idx + (NUM_QMF_BANDS * 2)]);
    accumlate += (*--ptr_filter *
                  pstr_sbr_enc_ana_filter_bank->x[NUM_QMF_BANDS + idx + (NUM_QMF_BANDS * 4)]);
    accumlate += (*--ptr_filter *
                  pstr_sbr_enc_ana_filter_bank->x[NUM_QMF_BANDS + idx + (NUM_QMF_BANDS * 6)]);
    accumlate += (*--ptr_filter *
                  pstr_sbr_enc_ana_filter_bank->x[NUM_QMF_BANDS + idx + (NUM_QMF_BANDS * 8)]);
    temp_array[63 - idx] = accumlate;
  }

  ixheaace_mps_forward_modulation(temp_array, real_subband, imag_subband);
}