wrapper/libxaac
1
0
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity Actions
libxaac/encoder/ixheaace_mps_polyphase.c
bmdivya100655 088122e748
Combined Workspace for Encoder and Decoder (#33) 
* Combined Workspace for Encoder and Decoder

* Addressed review comments and some minor edits

---------

Co-authored-by: Divya B M <100655@ittiam.com>
2023-05-22 20:11:52 +05:30

1023 lines
35 KiB
C
Raw Permalink Blame History

/******************************************************************************
* *
* 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);
}
Reference in a new issue View git blame Copy permalink