xingri.gao | c18d447 | 2023-02-28 02:51:02 +0000 | [diff] [blame] | 1 | /* |
| 2 | ** FAAD2 - Freeware Advanced Audio (AAC) Decoder including SBR decoding |
| 3 | ** Copyright (C) 2003-2005 M. Bakker, Nero AG, http://www.nero.com |
| 4 | ** |
| 5 | ** This program is free software; you can redistribute it and/or modify |
| 6 | ** it under the terms of the GNU General Public License as published by |
| 7 | ** the Free Software Foundation; either version 2 of the License, or |
| 8 | ** (at your option) any later version. |
| 9 | ** |
| 10 | ** This program is distributed in the hope that it will be useful, |
| 11 | ** but WITHOUT ANY WARRANTY; without even the implied warranty of |
| 12 | ** MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the |
| 13 | ** GNU General Public License for more details. |
| 14 | ** |
| 15 | ** You should have received a copy of the GNU General Public License |
| 16 | ** along with this program; if not, write to the Free Software |
| 17 | ** Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA 02111-1307, USA. |
| 18 | ** |
| 19 | ** Any non-GPL usage of this software or parts of this software is strictly |
| 20 | ** forbidden. |
| 21 | ** |
| 22 | ** The "appropriate copyright message" mentioned in section 2c of the GPLv2 |
| 23 | ** must read: "Code from FAAD2 is copyright (c) Nero AG, www.nero.com" |
| 24 | ** |
| 25 | ** Commercial non-GPL licensing of this software is possible. |
| 26 | ** For more info contact Nero AG through Mpeg4AAClicense@nero.com. |
| 27 | ** |
| 28 | ** $Id: sbr_dec.c,v 1.44 2009/01/26 22:32:31 menno Exp $ |
| 29 | **/ |
| 30 | |
| 31 | #include <stdlib.h> |
| 32 | #include "common.h" |
| 33 | #include "structs.h" |
| 34 | |
| 35 | #ifdef SBR_DEC |
| 36 | |
| 37 | #include <string.h> |
| 38 | |
| 39 | #include "syntax.h" |
| 40 | #include "bits.h" |
| 41 | #include "sbr_syntax.h" |
| 42 | #include "sbr_qmf.h" |
| 43 | #include "sbr_hfgen.h" |
| 44 | #include "sbr_hfadj.h" |
| 45 | |
| 46 | #ifndef __unused |
| 47 | #define __unused __attribute__((unused)) |
| 48 | #endif |
| 49 | |
| 50 | /* static function declarations */ |
| 51 | static uint8_t sbr_save_prev_data(sbr_info *sbr, uint8_t ch); |
| 52 | static void sbr_save_matrix(sbr_info *sbr, uint8_t ch); |
| 53 | |
| 54 | |
| 55 | sbr_info *sbrDecodeInit(uint16_t framelength, uint8_t id_aac, |
| 56 | uint32_t sample_rate, uint8_t downSampledSBR |
| 57 | #ifdef DRM |
| 58 | , uint8_t IsDRM |
| 59 | #endif |
| 60 | ) |
| 61 | { |
| 62 | sbr_info *sbr = faad_malloc(sizeof(sbr_info)); |
| 63 | memset(sbr, 0, sizeof(sbr_info)); |
| 64 | |
| 65 | /* save id of the parent element */ |
| 66 | sbr->id_aac = id_aac; |
| 67 | sbr->sample_rate = sample_rate; |
| 68 | |
| 69 | sbr->bs_freq_scale = 2; |
| 70 | sbr->bs_alter_scale = 1; |
| 71 | sbr->bs_noise_bands = 2; |
| 72 | sbr->bs_limiter_bands = 2; |
| 73 | sbr->bs_limiter_gains = 2; |
| 74 | sbr->bs_interpol_freq = 1; |
| 75 | sbr->bs_smoothing_mode = 1; |
| 76 | sbr->bs_start_freq = 5; |
| 77 | sbr->bs_amp_res = 1; |
| 78 | sbr->bs_samplerate_mode = 1; |
| 79 | sbr->prevEnvIsShort[0] = -1; |
| 80 | sbr->prevEnvIsShort[1] = -1; |
| 81 | sbr->header_count = 0; |
| 82 | sbr->Reset = 1; |
| 83 | |
| 84 | #ifdef DRM |
| 85 | sbr->Is_DRM_SBR = IsDRM; |
| 86 | #endif |
| 87 | sbr->tHFGen = T_HFGEN; |
| 88 | sbr->tHFAdj = T_HFADJ; |
| 89 | |
| 90 | sbr->bsco = 0; |
| 91 | sbr->bsco_prev = 0; |
| 92 | sbr->M_prev = 0; |
| 93 | sbr->frame_len = framelength; |
| 94 | |
| 95 | /* force sbr reset */ |
| 96 | sbr->bs_start_freq_prev = -1; |
| 97 | |
| 98 | if (framelength == 960) { |
| 99 | sbr->numTimeSlotsRate = RATE * NO_TIME_SLOTS_960; |
| 100 | sbr->numTimeSlots = NO_TIME_SLOTS_960; |
| 101 | } else { |
| 102 | sbr->numTimeSlotsRate = RATE * NO_TIME_SLOTS; |
| 103 | sbr->numTimeSlots = NO_TIME_SLOTS; |
| 104 | } |
| 105 | |
| 106 | sbr->GQ_ringbuf_index[0] = 0; |
| 107 | sbr->GQ_ringbuf_index[1] = 0; |
| 108 | |
| 109 | if (id_aac == ID_CPE) { |
| 110 | /* stereo */ |
| 111 | uint8_t j; |
| 112 | sbr->qmfa[0] = qmfa_init(32); |
| 113 | sbr->qmfa[1] = qmfa_init(32); |
| 114 | sbr->qmfs[0] = qmfs_init((downSampledSBR) ? 32 : 64); |
| 115 | sbr->qmfs[1] = qmfs_init((downSampledSBR) ? 32 : 64); |
| 116 | |
| 117 | for (j = 0; j < 5; j++) { |
| 118 | sbr->G_temp_prev[0][j] = faad_malloc(64 * sizeof(real_t)); |
| 119 | sbr->G_temp_prev[1][j] = faad_malloc(64 * sizeof(real_t)); |
| 120 | sbr->Q_temp_prev[0][j] = faad_malloc(64 * sizeof(real_t)); |
| 121 | sbr->Q_temp_prev[1][j] = faad_malloc(64 * sizeof(real_t)); |
| 122 | } |
| 123 | |
| 124 | memset(sbr->Xsbr[0], 0, (sbr->numTimeSlotsRate + sbr->tHFGen) * 64 * sizeof(qmf_t)); |
| 125 | memset(sbr->Xsbr[1], 0, (sbr->numTimeSlotsRate + sbr->tHFGen) * 64 * sizeof(qmf_t)); |
| 126 | } else { |
| 127 | /* mono */ |
| 128 | uint8_t j; |
| 129 | sbr->qmfa[0] = qmfa_init(32); |
| 130 | sbr->qmfs[0] = qmfs_init((downSampledSBR) ? 32 : 64); |
| 131 | sbr->qmfs[1] = NULL; |
| 132 | |
| 133 | for (j = 0; j < 5; j++) { |
| 134 | sbr->G_temp_prev[0][j] = faad_malloc(64 * sizeof(real_t)); |
| 135 | sbr->Q_temp_prev[0][j] = faad_malloc(64 * sizeof(real_t)); |
| 136 | } |
| 137 | |
| 138 | memset(sbr->Xsbr[0], 0, (sbr->numTimeSlotsRate + sbr->tHFGen) * 64 * sizeof(qmf_t)); |
| 139 | } |
| 140 | |
yuliang.hu | 585b720 | 2024-09-06 17:28:13 +0800 | [diff] [blame^] | 141 | sbr->process_buf = calloc(38 * 64, sizeof(qmf_t)); |
| 142 | if (sbr->process_buf) { |
| 143 | sbr->process_buf_size = 38 * 64 * sizeof(qmf_t); |
| 144 | } else { |
| 145 | sbr->process_buf_size = 0; |
| 146 | } |
| 147 | sbr->process_buf2 = calloc(38 * 64, sizeof(qmf_t)); |
| 148 | if (sbr->process_buf2) { |
| 149 | sbr->process_buf_size2 = 38 * 64 * sizeof(qmf_t); |
| 150 | } else { |
| 151 | sbr->process_buf_size2 = 0; |
| 152 | } |
| 153 | |
xingri.gao | c18d447 | 2023-02-28 02:51:02 +0000 | [diff] [blame] | 154 | return sbr; |
| 155 | } |
| 156 | |
| 157 | void sbrDecodeEnd(sbr_info *sbr) |
| 158 | { |
| 159 | uint8_t j; |
| 160 | |
| 161 | if (sbr) { |
| 162 | qmfa_end(sbr->qmfa[0]); |
| 163 | qmfs_end(sbr->qmfs[0]); |
| 164 | if (sbr->qmfs[1] != NULL) { |
| 165 | qmfa_end(sbr->qmfa[1]); |
| 166 | qmfs_end(sbr->qmfs[1]); |
| 167 | } |
| 168 | |
| 169 | for (j = 0; j < 5; j++) { |
| 170 | if (sbr->G_temp_prev[0][j]) { |
| 171 | faad_free(sbr->G_temp_prev[0][j]); |
| 172 | } |
| 173 | if (sbr->Q_temp_prev[0][j]) { |
| 174 | faad_free(sbr->Q_temp_prev[0][j]); |
| 175 | } |
| 176 | if (sbr->G_temp_prev[1][j]) { |
| 177 | faad_free(sbr->G_temp_prev[1][j]); |
| 178 | } |
| 179 | if (sbr->Q_temp_prev[1][j]) { |
| 180 | faad_free(sbr->Q_temp_prev[1][j]); |
| 181 | } |
| 182 | } |
| 183 | |
| 184 | #ifdef PS_DEC |
| 185 | if (sbr->ps != NULL) { |
| 186 | ps_free(sbr->ps); |
| 187 | } |
| 188 | #endif |
| 189 | |
| 190 | #ifdef DRM_PS |
| 191 | if (sbr->drm_ps != NULL) { |
| 192 | drm_ps_free(sbr->drm_ps); |
| 193 | } |
| 194 | #endif |
| 195 | |
| 196 | faad_free(sbr); |
| 197 | } |
| 198 | } |
| 199 | |
| 200 | void sbrReset(sbr_info *sbr) |
| 201 | { |
| 202 | uint8_t j; |
| 203 | if (sbr->qmfa[0] != NULL) { |
| 204 | memset(sbr->qmfa[0]->x, 0, 2 * sbr->qmfa[0]->channels * 10 * sizeof(real_t)); |
| 205 | } |
| 206 | if (sbr->qmfa[1] != NULL) { |
| 207 | memset(sbr->qmfa[1]->x, 0, 2 * sbr->qmfa[1]->channels * 10 * sizeof(real_t)); |
| 208 | } |
| 209 | if (sbr->qmfs[0] != NULL) { |
| 210 | memset(sbr->qmfs[0]->v, 0, 2 * sbr->qmfs[0]->channels * 20 * sizeof(real_t)); |
| 211 | } |
| 212 | if (sbr->qmfs[1] != NULL) { |
| 213 | memset(sbr->qmfs[1]->v, 0, 2 * sbr->qmfs[1]->channels * 20 * sizeof(real_t)); |
| 214 | } |
| 215 | |
| 216 | for (j = 0; j < 5; j++) { |
| 217 | if (sbr->G_temp_prev[0][j] != NULL) { |
| 218 | memset(sbr->G_temp_prev[0][j], 0, 64 * sizeof(real_t)); |
| 219 | } |
| 220 | if (sbr->G_temp_prev[1][j] != NULL) { |
| 221 | memset(sbr->G_temp_prev[1][j], 0, 64 * sizeof(real_t)); |
| 222 | } |
| 223 | if (sbr->Q_temp_prev[0][j] != NULL) { |
| 224 | memset(sbr->Q_temp_prev[0][j], 0, 64 * sizeof(real_t)); |
| 225 | } |
| 226 | if (sbr->Q_temp_prev[1][j] != NULL) { |
| 227 | memset(sbr->Q_temp_prev[1][j], 0, 64 * sizeof(real_t)); |
| 228 | } |
| 229 | } |
| 230 | |
| 231 | memset(sbr->Xsbr[0], 0, (sbr->numTimeSlotsRate + sbr->tHFGen) * 64 * sizeof(qmf_t)); |
| 232 | memset(sbr->Xsbr[1], 0, (sbr->numTimeSlotsRate + sbr->tHFGen) * 64 * sizeof(qmf_t)); |
| 233 | |
| 234 | sbr->GQ_ringbuf_index[0] = 0; |
| 235 | sbr->GQ_ringbuf_index[1] = 0; |
| 236 | sbr->header_count = 0; |
| 237 | sbr->Reset = 1; |
| 238 | |
| 239 | sbr->L_E_prev[0] = 0; |
| 240 | sbr->L_E_prev[1] = 0; |
| 241 | sbr->bs_freq_scale = 2; |
| 242 | sbr->bs_alter_scale = 1; |
| 243 | sbr->bs_noise_bands = 2; |
| 244 | sbr->bs_limiter_bands = 2; |
| 245 | sbr->bs_limiter_gains = 2; |
| 246 | sbr->bs_interpol_freq = 1; |
| 247 | sbr->bs_smoothing_mode = 1; |
| 248 | sbr->bs_start_freq = 5; |
| 249 | sbr->bs_amp_res = 1; |
| 250 | sbr->bs_samplerate_mode = 1; |
| 251 | sbr->prevEnvIsShort[0] = -1; |
| 252 | sbr->prevEnvIsShort[1] = -1; |
| 253 | sbr->bsco = 0; |
| 254 | sbr->bsco_prev = 0; |
| 255 | sbr->M_prev = 0; |
| 256 | sbr->bs_start_freq_prev = -1; |
| 257 | |
| 258 | sbr->f_prev[0] = 0; |
| 259 | sbr->f_prev[1] = 0; |
| 260 | for (j = 0; j < MAX_M; j++) { |
| 261 | sbr->E_prev[0][j] = 0; |
| 262 | sbr->Q_prev[0][j] = 0; |
| 263 | sbr->E_prev[1][j] = 0; |
| 264 | sbr->Q_prev[1][j] = 0; |
| 265 | sbr->bs_add_harmonic_prev[0][j] = 0; |
| 266 | sbr->bs_add_harmonic_prev[1][j] = 0; |
| 267 | } |
| 268 | sbr->bs_add_harmonic_flag_prev[0] = 0; |
| 269 | sbr->bs_add_harmonic_flag_prev[1] = 0; |
| 270 | } |
| 271 | |
| 272 | static uint8_t sbr_save_prev_data(sbr_info *sbr, uint8_t ch) |
| 273 | { |
| 274 | uint8_t i; |
| 275 | |
| 276 | /* save data for next frame */ |
| 277 | sbr->kx_prev = sbr->kx; |
| 278 | sbr->M_prev = sbr->M; |
| 279 | sbr->bsco_prev = sbr->bsco; |
| 280 | |
| 281 | sbr->L_E_prev[ch] = sbr->L_E[ch]; |
| 282 | |
| 283 | /* sbr->L_E[ch] can become 0 on files with bit errors */ |
| 284 | if (sbr->L_E[ch] <= 0) { |
| 285 | return 19; |
| 286 | } |
| 287 | |
| 288 | sbr->f_prev[ch] = sbr->f[ch][sbr->L_E[ch] - 1]; |
| 289 | for (i = 0; i < MAX_M; i++) { |
| 290 | sbr->E_prev[ch][i] = sbr->E[ch][i][sbr->L_E[ch] - 1]; |
| 291 | sbr->Q_prev[ch][i] = sbr->Q[ch][i][sbr->L_Q[ch] - 1]; |
| 292 | } |
| 293 | |
| 294 | for (i = 0; i < MAX_M; i++) { |
| 295 | sbr->bs_add_harmonic_prev[ch][i] = sbr->bs_add_harmonic[ch][i]; |
| 296 | } |
| 297 | sbr->bs_add_harmonic_flag_prev[ch] = sbr->bs_add_harmonic_flag[ch]; |
| 298 | |
| 299 | if (sbr->l_A[ch] == sbr->L_E[ch]) { |
| 300 | sbr->prevEnvIsShort[ch] = 0; |
| 301 | } else { |
| 302 | sbr->prevEnvIsShort[ch] = -1; |
| 303 | } |
| 304 | |
| 305 | return 0; |
| 306 | } |
| 307 | |
| 308 | static void sbr_save_matrix(sbr_info *sbr, uint8_t ch) |
| 309 | { |
| 310 | uint8_t i; |
| 311 | |
| 312 | for (i = 0; i < sbr->tHFGen; i++) { |
| 313 | memmove(sbr->Xsbr[ch][i], sbr->Xsbr[ch][i + sbr->numTimeSlotsRate], 64 * sizeof(qmf_t)); |
| 314 | } |
| 315 | for (i = sbr->tHFGen; i < MAX_NTSRHFG; i++) { |
| 316 | memset(sbr->Xsbr[ch][i], 0, 64 * sizeof(qmf_t)); |
| 317 | } |
| 318 | } |
| 319 | |
| 320 | static uint8_t sbr_process_channel(sbr_info *sbr, real_t *channel_buf, qmf_t X[MAX_NTSR][64], |
| 321 | uint8_t ch, uint8_t dont_process, |
| 322 | const uint8_t downSampledSBR __unused) |
| 323 | { |
| 324 | int16_t k, l; |
| 325 | uint8_t ret = 0; |
| 326 | |
| 327 | #ifdef SBR_LOW_POWER |
| 328 | ALIGN real_t deg[64]; |
| 329 | #endif |
| 330 | |
| 331 | #ifdef DRM |
| 332 | if (sbr->Is_DRM_SBR) { |
| 333 | sbr->bsco = max((int32_t)sbr->maxAACLine * 32 / (int32_t)sbr->frame_len - (int32_t)sbr->kx, 0); |
| 334 | } else { |
| 335 | #endif |
| 336 | sbr->bsco = 0; |
| 337 | #ifdef DRM |
| 338 | } |
| 339 | #endif |
| 340 | |
| 341 | |
| 342 | //#define PRE_QMF_PRINT |
| 343 | #ifdef PRE_QMF_PRINT |
| 344 | { |
| 345 | int i; |
| 346 | for (i = 0; i < 1024; i++) { |
| 347 | printf("%d\n", channel_buf[i]); |
| 348 | } |
| 349 | } |
| 350 | #endif |
| 351 | |
| 352 | |
| 353 | /* subband analysis */ |
| 354 | if (dont_process) { |
| 355 | sbr_qmf_analysis_32(sbr, sbr->qmfa[ch], channel_buf, sbr->Xsbr[ch], sbr->tHFGen, 32); |
| 356 | } else { |
| 357 | sbr_qmf_analysis_32(sbr, sbr->qmfa[ch], channel_buf, sbr->Xsbr[ch], sbr->tHFGen, sbr->kx); |
| 358 | } |
| 359 | |
| 360 | if (!dont_process) { |
| 361 | #if 1 |
| 362 | /* insert high frequencies here */ |
| 363 | /* hf generation using patching */ |
| 364 | hf_generation(sbr, sbr->Xsbr[ch], sbr->Xsbr[ch] |
| 365 | #ifdef SBR_LOW_POWER |
| 366 | , deg |
| 367 | #endif |
| 368 | , ch); |
| 369 | #endif |
| 370 | |
| 371 | #if 0 //def SBR_LOW_POWER |
| 372 | for (l = sbr->t_E[ch][0]; l < sbr->t_E[ch][sbr->L_E[ch]]; l++) { |
| 373 | for (k = 0; k < sbr->kx; k++) { |
| 374 | QMF_RE(sbr->Xsbr[ch][sbr->tHFAdj + l][k]) = 0; |
| 375 | } |
| 376 | } |
| 377 | #endif |
| 378 | |
| 379 | #if 1 |
| 380 | /* hf adjustment */ |
| 381 | ret = hf_adjustment(sbr, sbr->Xsbr[ch] |
| 382 | #ifdef SBR_LOW_POWER |
| 383 | , deg |
| 384 | #endif |
| 385 | , ch); |
| 386 | #endif |
| 387 | if (ret > 0) { |
| 388 | dont_process = 1; |
| 389 | } |
| 390 | } |
| 391 | |
| 392 | if ((sbr->just_seeked != 0) || dont_process) { |
| 393 | for (l = 0; l < sbr->numTimeSlotsRate; l++) { |
| 394 | for (k = 0; k < 32; k++) { |
| 395 | QMF_RE(X[l][k]) = QMF_RE(sbr->Xsbr[ch][l + sbr->tHFAdj][k]); |
| 396 | #ifndef SBR_LOW_POWER |
| 397 | QMF_IM(X[l][k]) = QMF_IM(sbr->Xsbr[ch][l + sbr->tHFAdj][k]); |
| 398 | #endif |
| 399 | } |
| 400 | for (k = 32; k < 64; k++) { |
| 401 | QMF_RE(X[l][k]) = 0; |
| 402 | #ifndef SBR_LOW_POWER |
| 403 | QMF_IM(X[l][k]) = 0; |
| 404 | #endif |
| 405 | } |
| 406 | } |
| 407 | } else { |
| 408 | for (l = 0; l < sbr->numTimeSlotsRate; l++) { |
| 409 | uint8_t kx_band, M_band, bsco_band; |
| 410 | |
| 411 | if (l < sbr->t_E[ch][0]) { |
| 412 | kx_band = sbr->kx_prev; |
| 413 | M_band = sbr->M_prev; |
| 414 | bsco_band = sbr->bsco_prev; |
| 415 | } else { |
| 416 | kx_band = sbr->kx; |
| 417 | M_band = sbr->M; |
| 418 | bsco_band = sbr->bsco; |
| 419 | } |
| 420 | |
| 421 | #ifndef SBR_LOW_POWER |
| 422 | for (k = 0; k < kx_band + bsco_band; k++) { |
| 423 | QMF_RE(X[l][k]) = QMF_RE(sbr->Xsbr[ch][l + sbr->tHFAdj][k]); |
| 424 | QMF_IM(X[l][k]) = QMF_IM(sbr->Xsbr[ch][l + sbr->tHFAdj][k]); |
| 425 | } |
| 426 | for (k = kx_band + bsco_band; k < kx_band + M_band; k++) { |
| 427 | QMF_RE(X[l][k]) = QMF_RE(sbr->Xsbr[ch][l + sbr->tHFAdj][k]); |
| 428 | QMF_IM(X[l][k]) = QMF_IM(sbr->Xsbr[ch][l + sbr->tHFAdj][k]); |
| 429 | } |
| 430 | for (k = max(kx_band + bsco_band, kx_band + M_band); k < 64; k++) { |
| 431 | QMF_RE(X[l][k]) = 0; |
| 432 | QMF_IM(X[l][k]) = 0; |
| 433 | } |
| 434 | #else |
| 435 | for (k = 0; k < kx_band + bsco_band; k++) { |
| 436 | QMF_RE(X[l][k]) = QMF_RE(sbr->Xsbr[ch][l + sbr->tHFAdj][k]); |
| 437 | } |
| 438 | for (k = kx_band + bsco_band; k < min(kx_band + M_band, 63); k++) { |
| 439 | QMF_RE(X[l][k]) = QMF_RE(sbr->Xsbr[ch][l + sbr->tHFAdj][k]); |
| 440 | } |
| 441 | for (k = max(kx_band + bsco_band, kx_band + M_band); k < 64; k++) { |
| 442 | QMF_RE(X[l][k]) = 0; |
| 443 | } |
| 444 | QMF_RE(X[l][kx_band - 1 + bsco_band]) += |
| 445 | QMF_RE(sbr->Xsbr[ch][l + sbr->tHFAdj][kx_band - 1 + bsco_band]); |
| 446 | #endif |
| 447 | } |
| 448 | } |
| 449 | |
| 450 | return ret; |
| 451 | } |
| 452 | |
| 453 | uint8_t sbrDecodeCoupleFrame(sbr_info *sbr, real_t *left_chan, real_t *right_chan, |
| 454 | const uint8_t just_seeked, const uint8_t downSampledSBR) |
| 455 | { |
| 456 | uint8_t dont_process = 0; |
| 457 | uint8_t ret = 0; |
yuliang.hu | 585b720 | 2024-09-06 17:28:13 +0800 | [diff] [blame^] | 458 | ALIGN qmf_t (*X)[64] = NULL; |
xingri.gao | c18d447 | 2023-02-28 02:51:02 +0000 | [diff] [blame] | 459 | |
| 460 | if (sbr == NULL) { |
| 461 | return 20; |
| 462 | } |
yuliang.hu | 585b720 | 2024-09-06 17:28:13 +0800 | [diff] [blame^] | 463 | memset(sbr->process_buf, 0, sizeof(sbr->process_buf_size)); |
| 464 | X = (qmf_t (*)[64])sbr->process_buf; |
xingri.gao | c18d447 | 2023-02-28 02:51:02 +0000 | [diff] [blame] | 465 | |
| 466 | /* case can occur due to bit errors */ |
| 467 | if (sbr->id_aac != ID_CPE) { |
| 468 | return 21; |
| 469 | } |
| 470 | |
| 471 | if (sbr->ret || (sbr->header_count == 0)) { |
| 472 | /* don't process just upsample */ |
| 473 | dont_process = 1; |
| 474 | |
| 475 | /* Re-activate reset for next frame */ |
| 476 | if (sbr->ret && sbr->Reset) { |
| 477 | sbr->bs_start_freq_prev = -1; |
| 478 | } |
| 479 | } |
| 480 | |
| 481 | if (just_seeked) { |
| 482 | sbr->just_seeked = 1; |
| 483 | } else { |
| 484 | sbr->just_seeked = 0; |
| 485 | } |
| 486 | |
| 487 | sbr->ret += sbr_process_channel(sbr, left_chan, X, 0, dont_process, downSampledSBR); |
| 488 | /* subband synthesis */ |
| 489 | if (downSampledSBR) { |
Jiacai.Liu | ba4fdc0 | 2023-11-21 19:58:04 +0800 | [diff] [blame] | 490 | sbr_qmf_synthesis_32(sbr, sbr->qmfs[0], &X[0][0], left_chan); |
xingri.gao | c18d447 | 2023-02-28 02:51:02 +0000 | [diff] [blame] | 491 | } else { |
Jiacai.Liu | ba4fdc0 | 2023-11-21 19:58:04 +0800 | [diff] [blame] | 492 | sbr_qmf_synthesis_64(sbr, sbr->qmfs[0], &X[0][0], left_chan); |
xingri.gao | c18d447 | 2023-02-28 02:51:02 +0000 | [diff] [blame] | 493 | } |
| 494 | |
| 495 | sbr->ret += sbr_process_channel(sbr, right_chan, X, 1, dont_process, downSampledSBR); |
| 496 | /* subband synthesis */ |
| 497 | if (downSampledSBR) { |
Jiacai.Liu | ba4fdc0 | 2023-11-21 19:58:04 +0800 | [diff] [blame] | 498 | sbr_qmf_synthesis_32(sbr, sbr->qmfs[1], &X[0][0], right_chan); |
xingri.gao | c18d447 | 2023-02-28 02:51:02 +0000 | [diff] [blame] | 499 | } else { |
Jiacai.Liu | ba4fdc0 | 2023-11-21 19:58:04 +0800 | [diff] [blame] | 500 | sbr_qmf_synthesis_64(sbr, sbr->qmfs[1], &X[0][0], right_chan); |
xingri.gao | c18d447 | 2023-02-28 02:51:02 +0000 | [diff] [blame] | 501 | } |
| 502 | |
| 503 | if (sbr->bs_header_flag) { |
| 504 | sbr->just_seeked = 0; |
| 505 | } |
| 506 | |
| 507 | if (sbr->header_count != 0 && sbr->ret == 0) { |
| 508 | ret = sbr_save_prev_data(sbr, 0); |
| 509 | if (ret) { |
| 510 | return ret; |
| 511 | } |
| 512 | ret = sbr_save_prev_data(sbr, 1); |
| 513 | if (ret) { |
| 514 | return ret; |
| 515 | } |
| 516 | } |
| 517 | |
| 518 | sbr_save_matrix(sbr, 0); |
| 519 | sbr_save_matrix(sbr, 1); |
| 520 | |
| 521 | sbr->frame++; |
| 522 | |
| 523 | //#define POST_QMF_PRINT |
| 524 | #ifdef POST_QMF_PRINT |
| 525 | { |
| 526 | int i; |
| 527 | for (i = 0; i < 2048; i++) { |
| 528 | printf("%d\n", left_chan[i]); |
| 529 | } |
| 530 | for (i = 0; i < 2048; i++) { |
| 531 | printf("%d\n", right_chan[i]); |
| 532 | } |
| 533 | } |
| 534 | #endif |
| 535 | |
| 536 | return 0; |
| 537 | } |
| 538 | |
| 539 | uint8_t sbrDecodeSingleFrame(sbr_info *sbr, real_t *channel, |
| 540 | const uint8_t just_seeked, const uint8_t downSampledSBR) |
| 541 | { |
| 542 | uint8_t dont_process = 0; |
| 543 | uint8_t ret = 0; |
yuliang.hu | 585b720 | 2024-09-06 17:28:13 +0800 | [diff] [blame^] | 544 | ALIGN qmf_t (*X)[64] = NULL; |
xingri.gao | c18d447 | 2023-02-28 02:51:02 +0000 | [diff] [blame] | 545 | |
| 546 | if (sbr == NULL) { |
| 547 | return 20; |
| 548 | } |
yuliang.hu | 585b720 | 2024-09-06 17:28:13 +0800 | [diff] [blame^] | 549 | memset(sbr->process_buf, 0, sizeof(sbr->process_buf_size)); |
| 550 | X = (qmf_t (*)[64])sbr->process_buf; |
xingri.gao | c18d447 | 2023-02-28 02:51:02 +0000 | [diff] [blame] | 551 | |
| 552 | /* case can occur due to bit errors */ |
| 553 | if (sbr->id_aac != ID_SCE && sbr->id_aac != ID_LFE) { |
| 554 | return 21; |
| 555 | } |
| 556 | |
| 557 | if (sbr->ret || (sbr->header_count == 0)) { |
| 558 | /* don't process just upsample */ |
| 559 | dont_process = 1; |
| 560 | |
| 561 | /* Re-activate reset for next frame */ |
| 562 | if (sbr->ret && sbr->Reset) { |
| 563 | sbr->bs_start_freq_prev = -1; |
| 564 | } |
| 565 | } |
| 566 | |
| 567 | if (just_seeked) { |
| 568 | sbr->just_seeked = 1; |
| 569 | } else { |
| 570 | sbr->just_seeked = 0; |
| 571 | } |
| 572 | |
| 573 | sbr->ret += sbr_process_channel(sbr, channel, X, 0, dont_process, downSampledSBR); |
| 574 | /* subband synthesis */ |
| 575 | if (downSampledSBR) { |
Jiacai.Liu | ba4fdc0 | 2023-11-21 19:58:04 +0800 | [diff] [blame] | 576 | sbr_qmf_synthesis_32(sbr, sbr->qmfs[0], &X[0][0], channel); |
xingri.gao | c18d447 | 2023-02-28 02:51:02 +0000 | [diff] [blame] | 577 | } else { |
Jiacai.Liu | ba4fdc0 | 2023-11-21 19:58:04 +0800 | [diff] [blame] | 578 | sbr_qmf_synthesis_64(sbr, sbr->qmfs[0], &X[0][0], channel); |
xingri.gao | c18d447 | 2023-02-28 02:51:02 +0000 | [diff] [blame] | 579 | } |
| 580 | |
| 581 | if (sbr->bs_header_flag) { |
| 582 | sbr->just_seeked = 0; |
| 583 | } |
| 584 | |
| 585 | if (sbr->header_count != 0 && sbr->ret == 0) { |
| 586 | ret = sbr_save_prev_data(sbr, 0); |
| 587 | if (ret) { |
| 588 | return ret; |
| 589 | } |
| 590 | } |
| 591 | |
| 592 | sbr_save_matrix(sbr, 0); |
| 593 | |
| 594 | sbr->frame++; |
| 595 | |
| 596 | //#define POST_QMF_PRINT |
| 597 | #ifdef POST_QMF_PRINT |
| 598 | { |
| 599 | int i; |
| 600 | for (i = 0; i < 2048; i++) { |
| 601 | printf("%d\n", channel[i]); |
| 602 | } |
| 603 | } |
| 604 | #endif |
| 605 | |
| 606 | return 0; |
| 607 | } |
| 608 | |
| 609 | #if (defined(PS_DEC) || defined(DRM_PS)) |
| 610 | uint8_t sbrDecodeSingleFramePS(sbr_info *sbr, real_t *left_channel, real_t *right_channel, |
| 611 | const uint8_t just_seeked, const uint8_t downSampledSBR) |
| 612 | { |
| 613 | uint8_t l, k; |
| 614 | uint8_t dont_process = 0; |
| 615 | uint8_t ret = 0; |
yuliang.hu | 585b720 | 2024-09-06 17:28:13 +0800 | [diff] [blame^] | 616 | ALIGN qmf_t (*X_left)[64] = NULL; |
| 617 | ALIGN qmf_t (*X_right)[64] = NULL; |
xingri.gao | c18d447 | 2023-02-28 02:51:02 +0000 | [diff] [blame] | 618 | |
| 619 | if (sbr == NULL) { |
| 620 | return 20; |
| 621 | } |
yuliang.hu | 585b720 | 2024-09-06 17:28:13 +0800 | [diff] [blame^] | 622 | memset(sbr->process_buf, 0, sizeof(sbr->process_buf_size));/* must set this to 0 */ |
| 623 | memset(sbr->process_buf2, 0, sizeof(sbr->process_buf_size2)); |
| 624 | X_left = (qmf_t (*)[64])sbr->process_buf; |
| 625 | X_right = (qmf_t (*)[64])sbr->process_buf2; |
| 626 | |
xingri.gao | c18d447 | 2023-02-28 02:51:02 +0000 | [diff] [blame] | 627 | |
| 628 | /* case can occur due to bit errors */ |
| 629 | if (sbr->id_aac != ID_SCE && sbr->id_aac != ID_LFE) { |
| 630 | return 21; |
| 631 | } |
| 632 | |
| 633 | if (sbr->ret || (sbr->header_count == 0)) { |
| 634 | /* don't process just upsample */ |
| 635 | dont_process = 1; |
| 636 | |
| 637 | /* Re-activate reset for next frame */ |
| 638 | if (sbr->ret && sbr->Reset) { |
| 639 | sbr->bs_start_freq_prev = -1; |
| 640 | } |
| 641 | } |
| 642 | |
| 643 | if (just_seeked) { |
| 644 | sbr->just_seeked = 1; |
| 645 | } else { |
| 646 | sbr->just_seeked = 0; |
| 647 | } |
| 648 | |
| 649 | if (sbr->qmfs[1] == NULL) { |
| 650 | sbr->qmfs[1] = qmfs_init((downSampledSBR) ? 32 : 64); |
| 651 | } |
| 652 | |
| 653 | sbr->ret += sbr_process_channel(sbr, left_channel, X_left, 0, dont_process, downSampledSBR); |
| 654 | |
| 655 | /* copy some extra data for PS */ |
| 656 | for (l = sbr->numTimeSlotsRate; l < sbr->numTimeSlotsRate + 6; l++) { |
| 657 | for (k = 0; k < 5; k++) { |
| 658 | QMF_RE(X_left[l][k]) = QMF_RE(sbr->Xsbr[0][sbr->tHFAdj + l][k]); |
| 659 | QMF_IM(X_left[l][k]) = QMF_IM(sbr->Xsbr[0][sbr->tHFAdj + l][k]); |
| 660 | } |
| 661 | } |
| 662 | |
| 663 | /* perform parametric stereo */ |
| 664 | #ifdef DRM_PS |
| 665 | if (sbr->Is_DRM_SBR) { |
| 666 | drm_ps_decode(sbr->drm_ps, (sbr->ret > 0), X_left, X_right); |
| 667 | } else { |
| 668 | #endif |
| 669 | #ifdef PS_DEC |
| 670 | ps_decode(sbr->ps, X_left, X_right); |
| 671 | #endif |
| 672 | #ifdef DRM_PS |
| 673 | } |
| 674 | #endif |
| 675 | |
| 676 | /* subband synthesis */ |
| 677 | if (downSampledSBR) { |
Jiacai.Liu | ba4fdc0 | 2023-11-21 19:58:04 +0800 | [diff] [blame] | 678 | sbr_qmf_synthesis_32(sbr, sbr->qmfs[0], &X_left[0][0], left_channel); |
| 679 | sbr_qmf_synthesis_32(sbr, sbr->qmfs[1], &X_right[0][0], right_channel); |
xingri.gao | c18d447 | 2023-02-28 02:51:02 +0000 | [diff] [blame] | 680 | } else { |
Jiacai.Liu | ba4fdc0 | 2023-11-21 19:58:04 +0800 | [diff] [blame] | 681 | sbr_qmf_synthesis_64(sbr, sbr->qmfs[0], &X_left[0][0], left_channel); |
| 682 | sbr_qmf_synthesis_64(sbr, sbr->qmfs[1], &X_right[0][0], right_channel); |
xingri.gao | c18d447 | 2023-02-28 02:51:02 +0000 | [diff] [blame] | 683 | } |
| 684 | |
| 685 | if (sbr->bs_header_flag) { |
| 686 | sbr->just_seeked = 0; |
| 687 | } |
| 688 | |
| 689 | if (sbr->header_count != 0 && sbr->ret == 0) { |
| 690 | ret = sbr_save_prev_data(sbr, 0); |
| 691 | if (ret) { |
| 692 | return ret; |
| 693 | } |
| 694 | } |
| 695 | |
| 696 | sbr_save_matrix(sbr, 0); |
| 697 | |
| 698 | sbr->frame++; |
| 699 | |
| 700 | return 0; |
| 701 | } |
| 702 | #endif |
| 703 | |
| 704 | #endif |