/******************************************************************************
* *
* 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 <string.h>
#include "ixheaacd_sbr_common.h"
#include <ixheaacd_type_def.h>
#include "ixheaacd_constants.h"
#include "ixheaacd_basic_ops32.h"
#include "ixheaacd_basic_ops16.h"
#include "ixheaacd_basic_ops40.h"
#include "ixheaacd_basic_ops.h"
#include "ixheaacd_intrinsics.h"
#include "ixheaacd_common_rom.h"
#include "ixheaacd_bitbuffer.h"
#include "ixheaacd_sbrdecsettings.h"
#include "ixheaacd_sbr_scale.h"
#include "ixheaacd_lpp_tran.h"
#include "ixheaacd_env_extr_part.h"
#include "ixheaacd_sbr_rom.h"
#include "ixheaacd_hybrid.h"
#include "ixheaacd_ps_dec.h"
#include "ixheaacd_env_extr.h"
#include "ixheaacd_qmf_dec.h"
#include <ixheaacd_basic_op.h>
#include "ixheaacd_env_calc.h"
#include "ixheaacd_interface.h"
#include "ixheaacd_function_selector.h"
#include "ixheaacd_audioobjtypes.h"
#define mult16x16_16(a, b) ixheaacd_mult16((a), (b))
#define mac16x16(a, b, c) ixheaacd_mac16x16in32((a), (b), (c))
#define mpy_32x16(a, b) fixmuldiv2_32x16b((a), (b))
#define mpy_16x16(a, b) ixheaacd_mult16x16in32((a), (b))
#define mpy_32x32(a, b) ixheaacd_mult32((a), (b))
#define mpy_32x16H_n(a, b) ixheaacd_mult32x16hin32((a), (b))
#define msu16x16(a, b, c) msu16x16in32((a), (b), (c))
#define DCT3_LEN (32)
#define DCT2_LEN (64)
#define LP_SHIFT_VAL 7
#define HQ_SHIFT_64 4
#define RADIXSHIFT 1
#define ROUNDING_SPECTRA 1
#define HQ_SHIFT_VAL 4
VOID ixheaacd_dct2_64(WORD32 *x, WORD32 *X,
ia_qmf_dec_tables_struct *qmf_dec_tables_ptr,
WORD16 *filter_states) {
ixheaacd_pretwdct2(x, X);
ixheaacd_sbr_imdct_using_fft(qmf_dec_tables_ptr->w1024, 32, X, x,
qmf_dec_tables_ptr->dig_rev_table2_128,
qmf_dec_tables_ptr->dig_rev_table2_128,
qmf_dec_tables_ptr->dig_rev_table2_128,
qmf_dec_tables_ptr->dig_rev_table2_128);
ixheaacd_fftposttw(x, qmf_dec_tables_ptr);
ixheaacd_posttwdct2(x, filter_states, qmf_dec_tables_ptr);
return;
}
VOID ixheaacd_cplx_anal_qmffilt(const WORD16 *time_sample_buf,
ia_sbr_scale_fact_struct *sbr_scale_factor,
WORD32 **qmf_real, WORD32 **qmf_imag,
ia_sbr_qmf_filter_bank_struct *qmf_bank,
ia_qmf_dec_tables_struct *qmf_dec_tables_ptr,
WORD32 ch_fac, WORD32 low_pow_flag,
WORD audio_object_type) {
WORD32 i, k;
WORD32 num_time_slots = qmf_bank->num_time_slots;
WORD32 analysis_buffer[4 * NO_ANALYSIS_CHANNELS];
WORD16 *filter_states = qmf_bank->core_samples_buffer;
WORD16 *fp1, *fp2, *tmp;
WORD16 *filter_1;
WORD16 *filter_2;
WORD16 *filt_ptr;
if (audio_object_type != AOT_ER_AAC_ELD &&
audio_object_type != AOT_ER_AAC_LD) {
qmf_bank->filter_pos +=
(qmf_dec_tables_ptr->qmf_c - qmf_bank->analy_win_coeff);
qmf_bank->analy_win_coeff = qmf_dec_tables_ptr->qmf_c;
} else {
qmf_bank->filter_pos +=
(qmf_dec_tables_ptr->qmf_c_eld3 - qmf_bank->analy_win_coeff);
qmf_bank->analy_win_coeff = qmf_dec_tables_ptr->qmf_c_eld3;
}
filter_1 = qmf_bank->filter_pos;
if (audio_object_type != AOT_ER_AAC_ELD &&
audio_object_type != AOT_ER_AAC_LD) {
filter_2 = filter_1 + 64;
} else {
filter_2 = filter_1 + 32;
}
sbr_scale_factor->st_lb_scale = 0;
sbr_scale_factor->lb_scale = -10;
if (!low_pow_flag) {
if (audio_object_type != AOT_ER_AAC_ELD &&
audio_object_type != AOT_ER_AAC_LD) {
sbr_scale_factor->lb_scale = -8;
} else {
sbr_scale_factor->lb_scale = -9;
}
qmf_bank->cos_twiddle =
(WORD16 *)qmf_dec_tables_ptr->sbr_sin_cos_twiddle_l32;
qmf_bank->alt_sin_twiddle =
(WORD16 *)qmf_dec_tables_ptr->sbr_alt_sin_twiddle_l32;
if (audio_object_type != AOT_ER_AAC_ELD &&
audio_object_type != AOT_ER_AAC_LD) {
qmf_bank->t_cos = (WORD16 *)qmf_dec_tables_ptr->sbr_t_cos_sin_l32;
} else {
qmf_bank->t_cos =
(WORD16 *)qmf_dec_tables_ptr->ixheaacd_sbr_t_cos_sin_l32_eld;
}
}
fp1 = qmf_bank->anal_filter_states;
fp2 = qmf_bank->anal_filter_states + NO_ANALYSIS_CHANNELS;
if (audio_object_type == AOT_ER_AAC_ELD ||
audio_object_type == AOT_ER_AAC_LD) {
filter_2 = qmf_bank->filter_2;
fp1 = qmf_bank->fp1_anal;
fp2 = qmf_bank->fp2_anal;
}
for (i = 0; i < num_time_slots; i++) {
for (k = 0; k < NO_ANALYSIS_CHANNELS; k++)
filter_states[NO_ANALYSIS_CHANNELS - 1 - k] = time_sample_buf[ch_fac * k];
if (audio_object_type != AOT_ER_AAC_ELD &&
audio_object_type != AOT_ER_AAC_LD) {
ixheaacd_sbr_qmfanal32_winadds(fp1, fp2, filter_1, filter_2,
analysis_buffer, filter_states,
time_sample_buf, ch_fac);
} else {
ixheaacd_sbr_qmfanal32_winadds_eld(fp1, fp2, filter_1, filter_2,
analysis_buffer, filter_states,
time_sample_buf, ch_fac);
}
time_sample_buf += NO_ANALYSIS_CHANNELS * ch_fac;
filter_states -= NO_ANALYSIS_CHANNELS;
if (filter_states < qmf_bank->anal_filter_states) {
filter_states = qmf_bank->anal_filter_states + 288;
}
tmp = fp1;
fp1 = fp2;
fp2 = tmp;
if (audio_object_type != AOT_ER_AAC_ELD &&
audio_object_type != AOT_ER_AAC_LD) {
filter_1 += 64;
filter_2 += 64;
} else {
filter_1 += 32;
filter_2 += 32;
}
filt_ptr = filter_1;
filter_1 = filter_2;
filter_2 = filt_ptr;
if (audio_object_type != AOT_ER_AAC_ELD &&
audio_object_type != AOT_ER_AAC_LD) {
if (filter_2 > (qmf_bank->analy_win_coeff + 640)) {
filter_1 = (WORD16 *)qmf_bank->analy_win_coeff;
filter_2 = (WORD16 *)qmf_bank->analy_win_coeff + 64;
}
} else {
if (filter_2 > (qmf_bank->analy_win_coeff + 320)) {
filter_1 = (WORD16 *)qmf_bank->analy_win_coeff;
filter_2 = (WORD16 *)qmf_bank->analy_win_coeff + 32;
}
}
if (!low_pow_flag) {
ixheaacd_fwd_modulation(analysis_buffer, qmf_real[i], qmf_imag[i],
qmf_bank, qmf_dec_tables_ptr);
} else {
ixheaacd_dct3_32(
(WORD32 *)analysis_buffer, qmf_real[i], qmf_dec_tables_ptr->dct23_tw,
qmf_dec_tables_ptr->post_fft_tbl, qmf_dec_tables_ptr->w_16,
qmf_dec_tables_ptr->dig_rev_table4_16);
}
}
qmf_bank->filter_pos = filter_1;
qmf_bank->core_samples_buffer = filter_states;
if (audio_object_type == AOT_ER_AAC_ELD || audio_object_type == AOT_ER_AAC_LD)
{
qmf_bank->fp1_anal = fp1;
qmf_bank->fp2_anal = fp2;
qmf_bank->filter_2 = filter_2;
}
}
VOID ixheaacd_inv_modulation_lp(WORD32 *qmf_real, WORD16 *filter_states,
ia_sbr_qmf_filter_bank_struct *syn_qmf,
ia_qmf_dec_tables_struct *qmf_dec_tables_ptr) {
WORD32 L = syn_qmf->no_channels;
const WORD32 M = (L >> 1);
WORD32 *dct_in = qmf_real;
WORD32 time_out[2 * NO_SYNTHESIS_CHANNELS];
WORD32 ui_rem = ((WORD32)(&time_out[0]) % 8);
WORD32 *ptime_out = (pVOID)((WORD8 *)&time_out[0] + 8 - ui_rem);
if (L == 64)
ixheaacd_dec_DCT2_64_asm(dct_in, ptime_out, qmf_dec_tables_ptr->w1024,
qmf_dec_tables_ptr->dig_rev_table2_128,
qmf_dec_tables_ptr->post_fft_tbl,
qmf_dec_tables_ptr->dct23_tw, filter_states + M);
else
ixheaacd_dct2_32(dct_in, time_out, qmf_dec_tables_ptr, filter_states);
filter_states[3 * M] = 0;
}
VOID ixheaacd_inv_emodulation(WORD32 *qmf_real,
ia_sbr_qmf_filter_bank_struct *syn_qmf,
ia_qmf_dec_tables_struct *qmf_dec_tables_ptr) {
ixheaacd_cos_sin_mod(qmf_real, syn_qmf, (WORD16 *)qmf_dec_tables_ptr->w1024,
(WORD32 *)qmf_dec_tables_ptr->dig_rev_table2_128);
}
VOID ixheaacd_esbr_cos_sin_mod(WORD32 *subband,
ia_sbr_qmf_filter_bank_struct *qmf_bank,
WORD32 *p_twiddle, WORD32 *p_dig_rev_tbl) {
WORD32 z;
WORD32 temp[128];
WORD32 scaleshift = 0;
WORD32 M_2;
WORD32 M = ixheaacd_shr32(qmf_bank->no_channels, 1);
const WORD32 *p_sin;
const WORD32 *p_sin_cos;
WORD32 subband_tmp[128];
p_sin_cos = qmf_bank->esbr_cos_twiddle;
ixheaacd_esbr_cos_sin_mod_loop1(subband, M, p_sin_cos, subband_tmp);
M_2 = ixheaacd_shr32(M, 1);
if (M == 32) {
ixheaacd_esbr_radix4bfly(p_twiddle, subband_tmp, 1, 8);
ixheaacd_esbr_radix4bfly(p_twiddle + 48, subband_tmp, 4, 2);
ixheaacd_postradixcompute2(subband, subband_tmp, p_dig_rev_tbl, 32);
ixheaacd_esbr_radix4bfly(p_twiddle, &subband_tmp[64], 1, 8);
ixheaacd_esbr_radix4bfly(p_twiddle + 48, &subband_tmp[64], 4, 2);
ixheaacd_postradixcompute2(&subband[64], &subband_tmp[64], p_dig_rev_tbl,
32);
}
else if (M == 16) {
ixheaacd_esbr_radix4bfly(p_twiddle, subband_tmp, 1, 4);
ixheaacd_postradixcompute4(subband, subband_tmp, p_dig_rev_tbl, 16);
ixheaacd_esbr_radix4bfly(p_twiddle, &subband_tmp[64], 1, 4);
ixheaacd_postradixcompute4(&subband[64], &subband_tmp[64], p_dig_rev_tbl,
16);
}
else if (M == 12) {
for (z = 0; z < (qmf_bank->no_channels >> 1); z++) {
temp[z] = subband_tmp[2 * z];
temp[12 + z] = subband_tmp[2 * z + 1];
}
ixheaacd_complex_fft_p3(temp, &temp[12], 12, -1, &scaleshift);
for (z = 0; z < (qmf_bank->no_channels >> 1); z++) {
subband[2 * z] = temp[z];
subband[2 * z + 1] = temp[z + 12];
}
scaleshift = 0;
for (z = 0; z < (qmf_bank->no_channels >> 1); z++) {
temp[z] = subband_tmp[64 + 2 * z];
temp[12 + z] = subband_tmp[64 + 2 * z + 1];
}
ixheaacd_complex_fft_p3(temp, &temp[12], 12, -1, &scaleshift);
for (z = 0; z < (qmf_bank->no_channels >> 1); z++) {
subband[64 + 2 * z] = temp[z];
subband[64 + 2 * z + 1] = temp[z + 12];
}
}
else {
for (z = 0; z < (qmf_bank->no_channels >> 1); z++) {
temp[z] = subband_tmp[2 * z];
temp[8 + z] = subband_tmp[2 * z + 1];
}
(*ixheaacd_complex_fft_p2)(temp, &temp[8], 8, -1, &scaleshift);
for (z = 0; z < (qmf_bank->no_channels >> 1); z++) {
subband[2 * z] = temp[z] << scaleshift;
subband[2 * z + 1] = temp[z + 8] << scaleshift;
}
scaleshift = 0;
for (z = 0; z < (qmf_bank->no_channels >> 1); z++) {
temp[z] = subband_tmp[64 + 2 * z];
temp[8 + z] = subband_tmp[64 + 2 * z + 1];
}
(*ixheaacd_complex_fft_p2)(temp, &temp[8], 8, -1, &scaleshift);
for (z = 0; z < (qmf_bank->no_channels >> 1); z++) {
subband[64 + 2 * z] = temp[z] << scaleshift;
subband[64 + 2 * z + 1] = temp[8 + z] << scaleshift;
}
}
p_sin = qmf_bank->esbr_alt_sin_twiddle;
ixheaacd_esbr_cos_sin_mod_loop2(subband, p_sin, M);
}