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
|
/*
* dca_internal.h
* Copyright (C) 2004 Gildas Bazin <gbazin@videolan.org>
* Copyright (C) 2000-2003 Michel Lespinasse <walken@zoy.org>
* Copyright (C) 1999-2000 Aaron Holtzman <aholtzma@ess.engr.uvic.ca>
*
* This file is part of libdca, a free DTS Coherent Acoustics stream decoder.
* See http://www.videolan.org/developers/libdca.html for updates.
*
* libdca is free software; you can redistribute it and/or modify
* it under the terms of the GNU General Public License as published by
* the Free Software Foundation; either version 2 of the License, or
* (at your option) any later version.
*
* libdca is distributed in the hope that it will be useful,
* but WITHOUT ANY WARRANTY; without even the implied warranty of
* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
* GNU General Public License for more details.
*
* You should have received a copy of the GNU General Public License
* along with this program; if not, write to the
* Free Software Foundation, Inc.,
* 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA.
*/
#define DCA_SUBFRAMES_MAX (16)
#define DCA_PRIM_CHANNELS_MAX (5)
#define DCA_SUBBANDS (32)
#define DCA_ABITS_MAX (32) /* Should be 28 */
#define DCA_SUBSUBFAMES_MAX (4)
#define DCA_LFE_MAX (3)
struct dca_state_s {
/* Frame header */
int frame_type; /* type of the current frame */
int samples_deficit; /* deficit sample count */
int crc_present; /* crc is present in the bitstream */
int sample_blocks; /* number of PCM sample blocks */
int frame_size; /* primary frame byte size */
int amode; /* audio channels arrangement */
int sample_rate; /* audio sampling rate */
int bit_rate; /* transmission bit rate */
int downmix; /* embedded downmix enabled */
int dynrange; /* embedded dynamic range flag */
int timestamp; /* embedded time stamp flag */
int aux_data; /* auxiliary data flag */
int hdcd; /* source material is mastered in HDCD */
int ext_descr; /* extension audio descriptor flag */
int ext_coding; /* extended coding flag */
int aspf; /* audio sync word insertion flag */
int lfe; /* low frequency effects flag */
int predictor_history; /* predictor history flag */
int header_crc; /* header crc check bytes */
int multirate_inter; /* multirate interpolator switch */
int version; /* encoder software revision */
int copy_history; /* copy history */
int source_pcm_res; /* source pcm resolution */
int front_sum; /* front sum/difference flag */
int surround_sum; /* surround sum/difference flag */
int dialog_norm; /* dialog normalisation parameter */
/* Primary audio coding header */
int subframes; /* number of subframes */
int prim_channels; /* number of primary audio channels */
/* subband activity count */
int subband_activity[DCA_PRIM_CHANNELS_MAX];
/* high frequency vq start subband */
int vq_start_subband[DCA_PRIM_CHANNELS_MAX];
/* joint intensity coding index */
int joint_intensity[DCA_PRIM_CHANNELS_MAX];
/* transient mode code book */
int transient_huffman[DCA_PRIM_CHANNELS_MAX];
/* scale factor code book */
int scalefactor_huffman[DCA_PRIM_CHANNELS_MAX];
/* bit allocation quantizer select */
int bitalloc_huffman[DCA_PRIM_CHANNELS_MAX];
/* quantization index codebook select */
int quant_index_huffman[DCA_PRIM_CHANNELS_MAX][DCA_ABITS_MAX];
/* scale factor adjustment */
float scalefactor_adj[DCA_PRIM_CHANNELS_MAX][DCA_ABITS_MAX];
/* Primary audio coding side information */
int subsubframes; /* number of subsubframes */
int partial_samples; /* partial subsubframe samples count */
/* prediction mode (ADPCM used or not) */
int prediction_mode[DCA_PRIM_CHANNELS_MAX][DCA_SUBBANDS];
/* prediction VQ coefs */
int prediction_vq[DCA_PRIM_CHANNELS_MAX][DCA_SUBBANDS];
/* bit allocation index */
int bitalloc[DCA_PRIM_CHANNELS_MAX][DCA_SUBBANDS];
/* transition mode (transients) */
int transition_mode[DCA_PRIM_CHANNELS_MAX][DCA_SUBBANDS];
/* scale factors (2 if transient)*/
int scale_factor[DCA_PRIM_CHANNELS_MAX][DCA_SUBBANDS][2];
/* joint subband scale factors codebook */
int joint_huff[DCA_PRIM_CHANNELS_MAX];
/* joint subband scale factors */
int joint_scale_factor[DCA_PRIM_CHANNELS_MAX][DCA_SUBBANDS];
/* stereo downmix coefficients */
int downmix_coef[DCA_PRIM_CHANNELS_MAX][2];
/* dynamic range coefficient */
int dynrange_coef;
/* VQ encoded high frequency subbands */
int high_freq_vq[DCA_PRIM_CHANNELS_MAX][DCA_SUBBANDS];
/* Low frequency effect data */
double lfe_data[2*DCA_SUBSUBFAMES_MAX*DCA_LFE_MAX * 2 /*history*/];
int lfe_scale_factor;
/* Subband samples history (for ADPCM) */
double subband_samples_hist[DCA_PRIM_CHANNELS_MAX][DCA_SUBBANDS][4];
double subband_fir_hist[DCA_PRIM_CHANNELS_MAX][512];
double subband_fir_noidea[DCA_PRIM_CHANNELS_MAX][64];
/* Audio output */
level_t clev; /* centre channel mix level */
level_t slev; /* surround channels mix level */
int output; /* type of output */
level_t level; /* output level */
sample_t bias; /* output bias */
sample_t * samples; /* pointer to the internal audio samples buffer */
int downmixed;
int dynrnge; /* apply dynamic range */
level_t dynrng; /* dynamic range */
void * dynrngdata; /* dynamic range callback funtion and data */
level_t (* dynrngcall) (level_t range, void * dynrngdata);
/* Bitstream handling */
uint32_t * buffer_start;
uint32_t bits_left;
uint32_t current_word;
int word_mode; /* 16/14 bits word format (1 -> 16, 0 -> 14) */
int bigendian_mode; /* endianness (1 -> be, 0 -> le) */
/* Current position in DTS frame */
int current_subframe;
int current_subsubframe;
/* Pre-calculated cosine modulation coefs for the QMF */
double cos_mod[544];
/* Debug flag */
int debug_flag;
};
#define LEVEL_PLUS6DB 2.0
#define LEVEL_PLUS3DB 1.4142135623730951
#define LEVEL_3DB 0.7071067811865476
#define LEVEL_45DB 0.5946035575013605
#define LEVEL_6DB 0.5
int dca_downmix_init (int input, int flags, level_t * level,
level_t clev, level_t slev);
int dca_downmix_coeff (level_t * coeff, int acmod, int output, level_t level,
level_t clev, level_t slev);
void dca_downmix (sample_t * samples, int acmod, int output, sample_t bias,
level_t clev, level_t slev);
void dca_upmix (sample_t * samples, int acmod, int output);
#define ROUND(x) ((int)((x) + ((x) > 0 ? 0.5 : -0.5)))
#ifndef LIBDCA_FIXED
typedef sample_t quantizer_t;
#define SAMPLE(x) (x)
#define LEVEL(x) (x)
#define MUL(a,b) ((a) * (b))
#define MUL_L(a,b) ((a) * (b))
#define MUL_C(a,b) ((a) * (b))
#define DIV(a,b) ((a) / (b))
#define BIAS(x) ((x) + bias)
#else /* LIBDCA_FIXED */
typedef int16_t quantizer_t;
#define SAMPLE(x) (sample_t)((x) * (1 << 30))
#define LEVEL(x) (level_t)((x) * (1 << 26))
#if 0
#define MUL(a,b) ((int)(((int64_t)(a) * (b) + (1 << 29)) >> 30))
#define MUL_L(a,b) ((int)(((int64_t)(a) * (b) + (1 << 25)) >> 26))
#elif 1
#define MUL(a,b) \
({ int32_t _ta=(a), _tb=(b), _tc; \
_tc=(_ta & 0xffff)*(_tb >> 16)+(_ta >> 16)*(_tb & 0xffff); (int32_t)(((_tc >> 14))+ (((_ta >> 16)*(_tb >> 16)) << 2 )); })
#define MUL_L(a,b) \
({ int32_t _ta=(a), _tb=(b), _tc; \
_tc=(_ta & 0xffff)*(_tb >> 16)+(_ta >> 16)*(_tb & 0xffff); (int32_t)((_tc >> 10) + (((_ta >> 16)*(_tb >> 16)) << 6)); })
#else
#define MUL(a,b) (((a) >> 15) * ((b) >> 15))
#define MUL_L(a,b) (((a) >> 13) * ((b) >> 13))
#endif
#define MUL_C(a,b) MUL_L (a, LEVEL (b))
#define DIV(a,b) ((((int64_t)LEVEL (a)) << 26) / (b))
#define BIAS(x) (x)
#endif
|
