blob: ffc84052d3dc7e0e659b576c4f32c7fda4f9d511 [file] [log] [blame]
xingri.gaoc18d4472023-02-28 02:51:02 +00001/*
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 */
51static uint8_t sbr_save_prev_data(sbr_info *sbr, uint8_t ch);
52static void sbr_save_matrix(sbr_info *sbr, uint8_t ch);
53
54
55sbr_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.hu585b7202024-09-06 17:28:13 +0800141 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.gaoc18d4472023-02-28 02:51:02 +0000154 return sbr;
155}
156
157void 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
200void 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
272static 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
308static 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
320static 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
453uint8_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.hu585b7202024-09-06 17:28:13 +0800458 ALIGN qmf_t (*X)[64] = NULL;
xingri.gaoc18d4472023-02-28 02:51:02 +0000459
460 if (sbr == NULL) {
461 return 20;
462 }
yuliang.hu585b7202024-09-06 17:28:13 +0800463 memset(sbr->process_buf, 0, sizeof(sbr->process_buf_size));
464 X = (qmf_t (*)[64])sbr->process_buf;
xingri.gaoc18d4472023-02-28 02:51:02 +0000465
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.Liuba4fdc02023-11-21 19:58:04 +0800490 sbr_qmf_synthesis_32(sbr, sbr->qmfs[0], &X[0][0], left_chan);
xingri.gaoc18d4472023-02-28 02:51:02 +0000491 } else {
Jiacai.Liuba4fdc02023-11-21 19:58:04 +0800492 sbr_qmf_synthesis_64(sbr, sbr->qmfs[0], &X[0][0], left_chan);
xingri.gaoc18d4472023-02-28 02:51:02 +0000493 }
494
495 sbr->ret += sbr_process_channel(sbr, right_chan, X, 1, dont_process, downSampledSBR);
496 /* subband synthesis */
497 if (downSampledSBR) {
Jiacai.Liuba4fdc02023-11-21 19:58:04 +0800498 sbr_qmf_synthesis_32(sbr, sbr->qmfs[1], &X[0][0], right_chan);
xingri.gaoc18d4472023-02-28 02:51:02 +0000499 } else {
Jiacai.Liuba4fdc02023-11-21 19:58:04 +0800500 sbr_qmf_synthesis_64(sbr, sbr->qmfs[1], &X[0][0], right_chan);
xingri.gaoc18d4472023-02-28 02:51:02 +0000501 }
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
539uint8_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.hu585b7202024-09-06 17:28:13 +0800544 ALIGN qmf_t (*X)[64] = NULL;
xingri.gaoc18d4472023-02-28 02:51:02 +0000545
546 if (sbr == NULL) {
547 return 20;
548 }
yuliang.hu585b7202024-09-06 17:28:13 +0800549 memset(sbr->process_buf, 0, sizeof(sbr->process_buf_size));
550 X = (qmf_t (*)[64])sbr->process_buf;
xingri.gaoc18d4472023-02-28 02:51:02 +0000551
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.Liuba4fdc02023-11-21 19:58:04 +0800576 sbr_qmf_synthesis_32(sbr, sbr->qmfs[0], &X[0][0], channel);
xingri.gaoc18d4472023-02-28 02:51:02 +0000577 } else {
Jiacai.Liuba4fdc02023-11-21 19:58:04 +0800578 sbr_qmf_synthesis_64(sbr, sbr->qmfs[0], &X[0][0], channel);
xingri.gaoc18d4472023-02-28 02:51:02 +0000579 }
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))
610uint8_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.hu585b7202024-09-06 17:28:13 +0800616 ALIGN qmf_t (*X_left)[64] = NULL;
617 ALIGN qmf_t (*X_right)[64] = NULL;
xingri.gaoc18d4472023-02-28 02:51:02 +0000618
619 if (sbr == NULL) {
620 return 20;
621 }
yuliang.hu585b7202024-09-06 17:28:13 +0800622 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.gaoc18d4472023-02-28 02:51:02 +0000627
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.Liuba4fdc02023-11-21 19:58:04 +0800678 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.gaoc18d4472023-02-28 02:51:02 +0000680 } else {
Jiacai.Liuba4fdc02023-11-21 19:58:04 +0800681 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.gaoc18d4472023-02-28 02:51:02 +0000683 }
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