1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
341
342
343
344
345
346
347
348
349
350
351
352
353
354
355
356
357
358
359
360
361
362
363
364
365
366
367
368
369
370
371
372
373
374
375
376
377
378
379
380
381
382
383
384
385
386
387
388
389
390
391
392
393
394
395
396
397
398
399
400
401
402
403
404
405
406
407
408
409
410
411
412
413
414
415
416
417
418
419
420
421
422
423
424
425
426
427
428
429
430
431
432
433
434
435
436
437
438
439
440
441
442
443
444
445
446
447
448
449
450
451
452
453
454
455
456
457
|
/* ------------------------------------------------------------------
* Copyright (C) 1998-2009 PacketVideo
*
* 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.
* -------------------------------------------------------------------
*/
/*
Pathname: intensity_right.c
------------------------------------------------------------------------------
REVISION HISTORY
Description: Modified per review comments.
Description: Noticed that the code could be more efficient by
using some other method for storing the sign. The code was changed to
use a signed Int to store the table, and an adjustment of the q-format to
reflect the difference between the data being shifted by 16 and the table
being stored in q-15 format.
Description: Updated pseudocode
Description: When the multiplication of two 16-bits variables is stored in
an 32-bits variable, the result should be typecasted explicitly
to Int32 before it is stored.
*(pCoefRight++) = (Int32) tempInt2 * multiplier;
Description:
------------------------------------------------------------------------------
INPUT AND OUTPUT DEFINITIONS
Inputs:
scalefactor = Multiplier used to scale the data extracted from the left
channel for use on the right.
[const Int]
coef_per_win = Number of coefficients per window.
(128 for short, 1024 for long)
[const Int]
sfb_per_win = Number of scalefactor bands per window. This should be
a number divisible by four.
[const Int]
wins_in_group = The number of windows in the group being decoded.
This number falls within the range 1-8.
[const Int]
band_length = The length of the scalefactor band being decoded.
This value is divisible by 4.
[const Int]
codebook = Value that denotes which Huffman codebook was used for
the encoding of this grouped scalefactor band.
[const Int]
ms_used = Flag that denotes whether M/S is active for this band.
[const Bool]
q_formatLeft = The Q-format for the left channel's fixed-point spectral
coefficients, on a per-scalefactor band, non-grouped basis.
[const Int *, length MAXBANDS]
q_formatRight = The Q-format for the right channel's fixed-point spectral
coefficients, on a per-scalefactor band, non-grouped basis.
[Int *, length MAXBANDS]
coefLeft = Array containing the fixed-point spectral coefficients
for the left channel.
[const Int32 *, length 1024]
coefRight = Array containing the fixed-point spectral coefficients
for the right channel.
[Int32 *, length 1024]
Local Stores/Buffers/Pointers Needed:
intensity_factor = Table which stores the values of
0.5^(0), 0.5^(1/4), 0.5^(2/4) and 0.5^(3/4)
[UInt, length 4]
Global Stores/Buffers/Pointers Needed:
None
Outputs:
None
Pointers and Buffers Modified:
coefRight[] Contains the new spectral information
q_formatRight[] Q-format may be updated with changed fixed-point
data in coefRight.
Local Stores Modified:
None
Global Stores Modified:
None
------------------------------------------------------------------------------
FUNCTION DESCRIPTION
This function applies Intensity Stereo, generating data on the right channel
that is derived from data on the Left. A scalefactor is applied using the
following formula...
RightCh = LeftCh*0.5^(scalefactor/4)
This function works for one scalefactor band, which may belong to a group.
(i.e. the same scalefactor band repeated across multiple windows belonging
to one group.)
------------------------------------------------------------------------------
REQUIREMENTS
codebook must be either INTENSITY_HCB or INTENSITY_HCB2 when this function
is called.
ms_used must be 1 when TRUE, 0 when FALSE.
wins_in_group falls within the range [1-8]
------------------------------------------------------------------------------
REFERENCES
(1) ISO/IEC 14496-3:1999(E)
Part 3
Subpart 4.6.7.2.3 Decoding Process (Intensity Stereo)
(2) MPEG-2 NBC Audio Decoder
"This software module was originally developed by AT&T, Dolby
Laboratories, Fraunhofer Gesellschaft IIS in the course of development
of the MPEG-2 NBC/MPEG-4 Audio standard ISO/IEC 13818-7, 14496-1,2 and
3. This software module is an implementation of a part of one or more
MPEG-2 NBC/MPEG-4 Audio tools as specified by the MPEG-2 NBC/MPEG-4
Audio standard. ISO/IEC gives users of the MPEG-2 NBC/MPEG-4 Audio
standards free license to this software module or modifications thereof
for use in hardware or software products claiming conformance to the
MPEG-2 NBC/MPEG-4 Audio standards. Those intending to use this software
module in hardware or software products are advised that this use may
infringe existing patents. The original developer of this software
module and his/her company, the subsequent editors and their companies,
and ISO/IEC have no liability for use of this software module or
modifications thereof in an implementation. Copyright is not released
for non MPEG-2 NBC/MPEG-4 Audio conforming products.The original
developer retains full right to use the code for his/her own purpose,
assign or donate the code to a third party and to inhibit third party
from using the code for non MPEG-2 NBC/MPEG-4 Audio conforming products.
This copyright notice must be included in all copies or derivative
works."
Copyright(c)1996.
------------------------------------------------------------------------------
PSEUDO-CODE
multiplier = codebook AND 0x1;
multiplier = multiplier XOR ms_used;
multiplier = multiplier << 1;
multiplier = multiplier - 1;
multiplier = multiplier * intensity_factor[scalefactor & 0x3];
scf_div_4 = (scalefactor >> 2);
nextWinPtrUpdate = (coef_per_win - band_length);
FOR (win_indx = wins_in_group; win_indx > 0; win_indx--)
*(pQformatRight) = scf_div_4 + *(pQformatLeft) - 1;
FOR (tempInt = band_length; tempInt > 0; tempInt--)
tempInt2 = (Int)(*(pCoefLeft) >> 16);
*(pCoefRight) = tempInt2 * multiplier;
pCoefRight = pCoefRight + 1;
pCoefLeft = pCoefLeft + 1;
ENDFOR
pCoefRight = pCoefRight + nextWinPtrUpdate;
pCoefLeft = pCoefLeft + nextWinPtrUpdate;
pQformatRight = pQformatRight + sfb_per_win;
pQformatLeft = pQformatLeft + sfb_per_win;
ENDFOR
------------------------------------------------------------------------------
RESOURCES USED
When the code is written for a specific target processor the
resources used should be documented below.
STACK USAGE: [stack count for this module] + [variable to represent
stack usage for each subroutine called]
where: [stack usage variable] = stack usage for [subroutine
name] (see [filename].ext)
DATA MEMORY USED: x words
PROGRAM MEMORY USED: x words
CLOCK CYCLES: [cycle count equation for this module] + [variable
used to represent cycle count for each subroutine
called]
where: [cycle count variable] = cycle count for [subroutine
name] (see [filename].ext)
------------------------------------------------------------------------------
*/
/*----------------------------------------------------------------------------
; INCLUDES
----------------------------------------------------------------------------*/
#include "pv_audio_type_defs.h"
#include "intensity_right.h"
#include "e_huffmanconst.h"
#include "fxp_mul32.h"
#include "aac_mem_funcs.h"
/*----------------------------------------------------------------------------
; MACROS
; Define module specific macros here
----------------------------------------------------------------------------*/
/*----------------------------------------------------------------------------
; DEFINES
; Include all pre-processor statements here. Include conditional
; compile variables also.
----------------------------------------------------------------------------*/
/*----------------------------------------------------------------------------
; LOCAL FUNCTION DEFINITIONS
; Function Prototype declaration
----------------------------------------------------------------------------*/
/*----------------------------------------------------------------------------
; LOCAL STORE/BUFFER/POINTER DEFINITIONS
; Variable declaration - defined here and used outside this module
----------------------------------------------------------------------------*/
const Int16 intensity_factor[4] =
{
32767, /* (0.5^0.00)*2^15 - 1 (minus 1 for storage as type Int) */
27554, /* (0.5^0.25)*2^15 */
23170, /* (0.5^0.50)*2^15 */
19484
}; /* (0.5^0.75)*2^15 */
/*----------------------------------------------------------------------------
; EXTERNAL FUNCTION REFERENCES
; Declare functions defined elsewhere and referenced in this module
----------------------------------------------------------------------------*/
/*----------------------------------------------------------------------------
; EXTERNAL GLOBAL STORE/BUFFER/POINTER REFERENCES
; Declare variables used in this module but defined elsewhere
----------------------------------------------------------------------------*/
/*----------------------------------------------------------------------------
; FUNCTION CODE
----------------------------------------------------------------------------*/
void intensity_right(
const Int scalefactor,
const Int coef_per_win,
const Int sfb_per_win,
const Int wins_in_group,
const Int band_length,
const Int codebook,
const Bool ms_used,
const Int q_formatLeft[],
Int q_formatRight[],
const Int32 coefLeft[],
Int32 coefRight[])
{
const Int32 *pCoefLeft = coefLeft;
Int32 *pCoefRight = coefRight;
const Int *pQformatLeft = q_formatLeft;
Int *pQformatRight = q_formatRight;
Int multiplier;
Int scf_div_4;
Int nextWinPtrUpdate;
/*
* The sign of the intensity multiplier obeys the following table...
*
* codebook | ms_used | multiplier
* --------------------------------------
* INTENSITY_HCB | TRUE | -1
* INTENSITY_HCB | FALSE | +1
* INTENSITY_HCB2 | TRUE | +1
* INTENSITY_HCB2 | FALSE | -1
*
* In binary, the above table is represented as...
*
* codebook | ms_used | multiplier
* --------------------------------------
* 1111b | 1 | -1
* 1111b | 0 | +1
* 1110b | 1 | +1
* 1110b | 0 | -1
*
*/
/*
* Deriving the correct value for "multiplier" is illustrated
* below for all 4 possible combinations of codebook and ms_used
*/
/*
* 1111b AND 0x1 = 1b
* 1111b AND 0x1 = 1b
* 1110b AND 0x1 = 0b
* 1110b AND 0x1 = 0b
*/
multiplier = (codebook & 0x1);
/*
* 1b XOR 1 = 0b
* 1b XOR 0 = 1b
* 0b XOR 1 = 1b
* 0b XOR 0 = 0b
*/
multiplier ^= ms_used;
/*
* 0b << 1 = 0
* 1b << 1 = 2
* 1b << 1 = 2
* 0b << 1 = 0
*/
multiplier <<= 1;
/*
* 0 - 1 = -1
* 2 - 1 = +1
* 2 - 1 = +1
* 0 - 1 = -1
*/
multiplier--;
multiplier *= intensity_factor[scalefactor & 0x3];
scf_div_4 = (scalefactor >> 2);
/*
* Step through all the windows in this group, replacing this
* band in each window's spectrum with
* left-channel correlated data
*/
nextWinPtrUpdate = (coef_per_win - band_length);
for (Int win_indx = wins_in_group; win_indx > 0; win_indx--)
{
/*
* Calculate q_formatRight
*
* q_formatLeft must be included, since the values
* on the right-channel are derived from the values
* on the left-channel.
*
* scalefactor/4 is included, since the intensity
* formula is RightCh = LeftCh*0.5^(scalefactor/4)
*
* powers of 0.5 increase the q_format by 1.
* (Another way to multiply something by 0.5^(x)
* is to increase its q-format by x.)
*
* Finally the q-format must be decreased by 1.
* The reason for this is because the table is stored
* in q-15 format, but we are shifting by 16 to do
* a 16 x 16 multiply.
*/
*(pQformatRight) = scf_div_4 + *(pQformatLeft);
/*
* reconstruct right intensity values
*
* to make things faster, this for loop
* can be partially unrolled, since band_length is a multiple
* of four.
*/
if (multiplier == 32767)
{
Int32 tempInt2 = *(pCoefLeft++);
Int32 tempInt22 = *(pCoefLeft++);
for (Int tempInt = band_length >> 1; tempInt > 0; tempInt--)
{
*(pCoefRight++) = tempInt2;
*(pCoefRight++) = tempInt22;
tempInt2 = *(pCoefLeft++);
tempInt22 = *(pCoefLeft++);
}
}
else
{
Int32 tempInt2 = *(pCoefLeft++);
Int32 tempInt22 = *(pCoefLeft++);
for (Int tempInt = band_length >> 1; tempInt > 0; tempInt--)
{
*(pCoefRight++) = fxp_mul32_by_16(tempInt2, multiplier) << 1;
*(pCoefRight++) = fxp_mul32_by_16(tempInt22, multiplier) << 1;
tempInt2 = *(pCoefLeft++);
tempInt22 = *(pCoefLeft++);
}
}
/*
* Set pCoefRight and pCoefLeft to the beginning of
* this sfb in the next window in the group.
*/
pCoefRight += nextWinPtrUpdate;
pCoefLeft += (nextWinPtrUpdate - 2);
/*
* Update pQformatRight and pQformatLeft to this sfb in
* in the next window in the group.
*/
pQformatRight += sfb_per_win;
pQformatLeft += sfb_per_win;
} /* for (win_indx) */
} /* void intensity_right */
|