FFmpeg
 All Data Structures Namespaces Files Functions Variables Typedefs Enumerations Enumerator Macros Groups Pages
ttaenc.c
Go to the documentation of this file.
1 /*
2  * TTA (The Lossless True Audio) encoder
3  *
4  * This file is part of FFmpeg.
5  *
6  * FFmpeg is free software; you can redistribute it and/or
7  * modify it under the terms of the GNU Lesser General Public
8  * License as published by the Free Software Foundation; either
9  * version 2.1 of the License, or (at your option) any later version.
10  *
11  * FFmpeg is distributed in the hope that it will be useful,
12  * but WITHOUT ANY WARRANTY; without even the implied warranty of
13  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
14  * Lesser General Public License for more details.
15  *
16  * You should have received a copy of the GNU Lesser General Public
17  * License along with FFmpeg; if not, write to the Free Software
18  * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
19  */
20 
21 #define BITSTREAM_WRITER_LE
22 #include "ttadata.h"
23 #include "avcodec.h"
24 #include "put_bits.h"
25 #include "internal.h"
26 #include "libavutil/crc.h"
27 
28 typedef struct TTAEncContext {
29  const AVCRC *crc_table;
30  int bps;
33 
35 {
36  TTAEncContext *s = avctx->priv_data;
37 
39 
40  switch (avctx->sample_fmt) {
41  case AV_SAMPLE_FMT_U8:
42  avctx->bits_per_raw_sample = 8;
43  break;
44  case AV_SAMPLE_FMT_S16:
45  avctx->bits_per_raw_sample = 16;
46  break;
47  case AV_SAMPLE_FMT_S32:
48  if (avctx->bits_per_raw_sample > 24)
49  av_log(avctx, AV_LOG_WARNING, "encoding as 24 bits-per-sample\n");
50  avctx->bits_per_raw_sample = 24;
51  }
52 
53  s->bps = avctx->bits_per_raw_sample >> 3;
54  avctx->frame_size = 256 * avctx->sample_rate / 245;
55 
56  s->ch_ctx = av_malloc_array(avctx->channels, sizeof(*s->ch_ctx));
57  if (!s->ch_ctx)
58  return AVERROR(ENOMEM);
59 
60  return 0;
61 }
62 
63 static inline void ttafilter_process(TTAFilter *c, int32_t *in)
64 {
65  register int32_t *dl = c->dl, *qm = c->qm, *dx = c->dx, sum = c->round;
66 
67  if (c->error < 0) {
68  qm[0] -= dx[0]; qm[1] -= dx[1]; qm[2] -= dx[2]; qm[3] -= dx[3];
69  qm[4] -= dx[4]; qm[5] -= dx[5]; qm[6] -= dx[6]; qm[7] -= dx[7];
70  } else if (c->error > 0) {
71  qm[0] += dx[0]; qm[1] += dx[1]; qm[2] += dx[2]; qm[3] += dx[3];
72  qm[4] += dx[4]; qm[5] += dx[5]; qm[6] += dx[6]; qm[7] += dx[7];
73  }
74 
75  sum += dl[0] * qm[0] + dl[1] * qm[1] + dl[2] * qm[2] + dl[3] * qm[3] +
76  dl[4] * qm[4] + dl[5] * qm[5] + dl[6] * qm[6] + dl[7] * qm[7];
77 
78  dx[0] = dx[1]; dx[1] = dx[2]; dx[2] = dx[3]; dx[3] = dx[4];
79  dl[0] = dl[1]; dl[1] = dl[2]; dl[2] = dl[3]; dl[3] = dl[4];
80 
81  dx[4] = ((dl[4] >> 30) | 1);
82  dx[5] = ((dl[5] >> 30) | 2) & ~1;
83  dx[6] = ((dl[6] >> 30) | 2) & ~1;
84  dx[7] = ((dl[7] >> 30) | 4) & ~3;
85 
86  dl[4] = -dl[5]; dl[5] = -dl[6];
87  dl[6] = *in - dl[7]; dl[7] = *in;
88  dl[5] += dl[6]; dl[4] += dl[5];
89 
90  *in -= (sum >> c->shift);
91  c->error = *in;
92 }
93 
94 static int32_t get_sample(const AVFrame *frame, int sample,
95  enum AVSampleFormat format)
96 {
97  int32_t ret;
98 
99  if (format == AV_SAMPLE_FMT_U8) {
100  ret = frame->data[0][sample] - 0x80;
101  } else if (format == AV_SAMPLE_FMT_S16) {
102  const int16_t *ptr = (const int16_t *)frame->data[0];
103  ret = ptr[sample];
104  } else {
105  const int32_t *ptr = (const int32_t *)frame->data[0];
106  ret = ptr[sample] >> 8;
107  }
108 
109  return ret;
110 }
111 
112 static int tta_encode_frame(AVCodecContext *avctx, AVPacket *avpkt,
113  const AVFrame *frame, int *got_packet_ptr)
114 {
115  TTAEncContext *s = avctx->priv_data;
116  PutBitContext pb;
117  int ret, i, out_bytes, cur_chan = 0, res = 0, samples = 0;
118 
119  if ((ret = ff_alloc_packet2(avctx, avpkt, frame->nb_samples * 2 * avctx->channels * s->bps, 0)) < 0)
120  return ret;
121  init_put_bits(&pb, avpkt->data, avpkt->size);
122 
123  // init per channel states
124  for (i = 0; i < avctx->channels; i++) {
125  s->ch_ctx[i].predictor = 0;
127  ff_tta_rice_init(&s->ch_ctx[i].rice, 10, 10);
128  }
129 
130  for (i = 0; i < frame->nb_samples * avctx->channels; i++) {
131  TTAChannel *c = &s->ch_ctx[cur_chan];
132  TTAFilter *filter = &c->filter;
133  TTARice *rice = &c->rice;
134  uint32_t k, unary, outval;
135  int32_t value, temp;
136 
137  value = get_sample(frame, samples++, avctx->sample_fmt);
138 
139  if (avctx->channels > 1) {
140  if (cur_chan < avctx->channels - 1)
141  value = res = get_sample(frame, samples, avctx->sample_fmt) - value;
142  else
143  value -= res / 2;
144  }
145 
146  temp = value;
147 #define PRED(x, k) (int32_t)((((uint64_t)(x) << (k)) - (x)) >> (k))
148  switch (s->bps) {
149  case 1: value -= PRED(c->predictor, 4); break;
150  case 2:
151  case 3: value -= PRED(c->predictor, 5); break;
152  }
153  c->predictor = temp;
154 
155  ttafilter_process(filter, &value);
156  outval = (value > 0) ? (value << 1) - 1: -value << 1;
157 
158  k = rice->k0;
159 
160  rice->sum0 += outval - (rice->sum0 >> 4);
161  if (rice->k0 > 0 && rice->sum0 < ff_tta_shift_16[rice->k0])
162  rice->k0--;
163  else if (rice->sum0 > ff_tta_shift_16[rice->k0 + 1])
164  rice->k0++;
165 
166  if (outval >= ff_tta_shift_1[k]) {
167  outval -= ff_tta_shift_1[k];
168  k = rice->k1;
169 
170  rice->sum1 += outval - (rice->sum1 >> 4);
171  if (rice->k1 > 0 && rice->sum1 < ff_tta_shift_16[rice->k1])
172  rice->k1--;
173  else if (rice->sum1 > ff_tta_shift_16[rice->k1 + 1])
174  rice->k1++;
175 
176  unary = 1 + (outval >> k);
177  do {
178  if (unary > 31) {
179  put_bits(&pb, 31, 0x7FFFFFFF);
180  unary -= 31;
181  } else {
182  put_bits(&pb, unary, (1 << unary) - 1);
183  unary = 0;
184  }
185  } while (unary);
186  }
187 
188  put_bits(&pb, 1, 0);
189 
190  if (k)
191  put_bits(&pb, k, outval & (ff_tta_shift_1[k] - 1));
192 
193  if (cur_chan < avctx->channels - 1)
194  cur_chan++;
195  else
196  cur_chan = 0;
197  }
198 
199  flush_put_bits(&pb);
200  out_bytes = put_bits_count(&pb) >> 3;
201  put_bits32(&pb, av_crc(s->crc_table, UINT32_MAX, avpkt->data, out_bytes) ^ UINT32_MAX);
202  flush_put_bits(&pb);
203 
204  avpkt->pts = frame->pts;
205  avpkt->size = out_bytes + 4;
206  avpkt->duration = ff_samples_to_time_base(avctx, frame->nb_samples);
207  *got_packet_ptr = 1;
208  return 0;
209 }
210 
212 {
213  TTAEncContext *s = avctx->priv_data;
214  av_freep(&s->ch_ctx);
215  return 0;
216 }
217 
219  .name = "tta",
220  .long_name = NULL_IF_CONFIG_SMALL("TTA (True Audio)"),
221  .type = AVMEDIA_TYPE_AUDIO,
222  .id = AV_CODEC_ID_TTA,
223  .priv_data_size = sizeof(TTAEncContext),
225  .close = tta_encode_close,
226  .encode2 = tta_encode_frame,
228  .sample_fmts = (const enum AVSampleFormat[]){ AV_SAMPLE_FMT_U8,
232 };
static void av_unused put_bits32(PutBitContext *s, uint32_t value)
Write exactly 32 bits into a bitstream.
Definition: put_bits.h:202
const char * s
Definition: avisynth_c.h:631
This structure describes decoded (raw) audio or video data.
Definition: frame.h:171
static av_cold int tta_encode_init(AVCodecContext *avctx)
Definition: ttaenc.c:34
static void put_bits(Jpeg2000EncoderContext *s, int val, int n)
put n times val bit
Definition: j2kenc.c:167
#define AV_LOG_WARNING
Something somehow does not look correct.
Definition: log.h:182
else temp
Definition: vf_mcdeint.c:257
AVCodec ff_tta_encoder
Definition: ttaenc.c:218
static av_cold int init(AVCodecContext *avctx)
Definition: avrndec.c:35
void ff_tta_rice_init(TTARice *c, uint32_t k0, uint32_t k1)
Definition: ttadata.c:40
int size
Definition: avcodec.h:1423
int32_t predictor
Definition: ttadata.h:39
int bits_per_raw_sample
Bits per sample/pixel of internal libavcodec pixel/sample format.
Definition: avcodec.h:3004
#define sample
AVCodec.
Definition: avcodec.h:3473
static av_cold int tta_encode_close(AVCodecContext *avctx)
Definition: ttaenc.c:211
int32_t error
Definition: ttadata.h:28
int32_t round
Definition: ttadata.h:28
enum AVSampleFormat sample_fmt
audio sample format
Definition: avcodec.h:2269
#define av_cold
Definition: attributes.h:74
AV_SAMPLE_FMT_U8
int64_t pts
Presentation timestamp in time_base units (time when frame should be shown to user).
Definition: frame.h:257
static AVFrame * frame
uint8_t * data
Definition: avcodec.h:1422
int32_t qm[MAX_ORDER]
Definition: ttadata.h:29
int32_t shift
Definition: ttadata.h:28
signed 32 bits
Definition: samplefmt.h:63
int duration
Duration of this packet in AVStream->time_base units, 0 if unknown.
Definition: avcodec.h:1440
#define av_log(a,...)
const AVCRC * crc_table
Definition: ttaenc.c:29
#define AVERROR(e)
Definition: error.h:43
#define PRED(x, k)
#define NULL_IF_CONFIG_SMALL(x)
Return NULL if CONFIG_SMALL is true, otherwise the argument without modification. ...
Definition: internal.h:175
static int32_t get_sample(const AVFrame *frame, int sample, enum AVSampleFormat format)
Definition: ttaenc.c:94
const char * name
Name of the codec implementation.
Definition: avcodec.h:3480
Libavcodec external API header.
uint32_t sum1
Definition: ttadata.h:35
static int put_bits_count(PutBitContext *s)
Definition: put_bits.h:85
const uint32_t ff_tta_shift_1[]
Definition: ttadata.c:23
#define AV_CODEC_CAP_SMALL_LAST_FRAME
Codec can be fed a final frame with a smaller size.
Definition: avcodec.h:886
ret
Definition: avfilter.c:974
GLsizei GLboolean const GLfloat * value
Definition: opengl_enc.c:109
const uint32_t *const ff_tta_shift_16
Definition: ttadata.c:36
int32_t
uint32_t av_crc(const AVCRC *ctx, uint32_t crc, const uint8_t *buffer, size_t length)
Calculate the CRC of a block.
Definition: crc.c:356
uint32_t k1
Definition: ttadata.h:35
int frame_size
Number of samples per channel in an audio frame.
Definition: avcodec.h:2281
TTAChannel * ch_ctx
Definition: ttaenc.c:31
AVSampleFormat
Audio sample formats.
Definition: samplefmt.h:59
int sample_rate
samples per second
Definition: avcodec.h:2261
const uint8_t ff_tta_filter_configs[]
Definition: ttadata.c:38
main external API structure.
Definition: avcodec.h:1501
int32_t dx[MAX_ORDER]
Definition: ttadata.h:30
uint8_t pi<< 24) CONV_FUNC_GROUP(AV_SAMPLE_FMT_FLT, float, AV_SAMPLE_FMT_U8, uint8_t,(*(constuint8_t *) pi-0x80)*(1.0f/(1<< 7))) CONV_FUNC_GROUP(AV_SAMPLE_FMT_DBL, double, AV_SAMPLE_FMT_U8, uint8_t,(*(constuint8_t *) pi-0x80)*(1.0/(1<< 7))) CONV_FUNC_GROUP(AV_SAMPLE_FMT_U8, uint8_t, AV_SAMPLE_FMT_S16, int16_t,(*(constint16_t *) pi >>8)+0x80) CONV_FUNC_GROUP(AV_SAMPLE_FMT_FLT, float, AV_SAMPLE_FMT_S16, int16_t,*(constint16_t *) pi *(1.0f/(1<< 15))) CONV_FUNC_GROUP(AV_SAMPLE_FMT_DBL, double, AV_SAMPLE_FMT_S16, int16_t,*(constint16_t *) pi *(1.0/(1<< 15))) CONV_FUNC_GROUP(AV_SAMPLE_FMT_U8, uint8_t, AV_SAMPLE_FMT_S32, int32_t,(*(constint32_t *) pi >>24)+0x80) CONV_FUNC_GROUP(AV_SAMPLE_FMT_FLT, float, AV_SAMPLE_FMT_S32, int32_t,*(constint32_t *) pi *(1.0f/(1U<< 31))) CONV_FUNC_GROUP(AV_SAMPLE_FMT_DBL, double, AV_SAMPLE_FMT_S32, int32_t,*(constint32_t *) pi *(1.0/(1U<< 31))) CONV_FUNC_GROUP(AV_SAMPLE_FMT_U8, uint8_t, AV_SAMPLE_FMT_FLT, float, av_clip_uint8(lrintf(*(constfloat *) pi *(1<< 7))+0x80)) CONV_FUNC_GROUP(AV_SAMPLE_FMT_S16, int16_t, AV_SAMPLE_FMT_FLT, float, av_clip_int16(lrintf(*(constfloat *) pi *(1<< 15)))) CONV_FUNC_GROUP(AV_SAMPLE_FMT_S32, int32_t, AV_SAMPLE_FMT_FLT, float, av_clipl_int32(llrintf(*(constfloat *) pi *(1U<< 31)))) CONV_FUNC_GROUP(AV_SAMPLE_FMT_U8, uint8_t, AV_SAMPLE_FMT_DBL, double, av_clip_uint8(lrint(*(constdouble *) pi *(1<< 7))+0x80)) CONV_FUNC_GROUP(AV_SAMPLE_FMT_S16, int16_t, AV_SAMPLE_FMT_DBL, double, av_clip_int16(lrint(*(constdouble *) pi *(1<< 15)))) CONV_FUNC_GROUP(AV_SAMPLE_FMT_S32, int32_t, AV_SAMPLE_FMT_DBL, double, av_clipl_int32(llrint(*(constdouble *) pi *(1U<< 31))))#defineSET_CONV_FUNC_GROUP(ofmt, ifmt) staticvoidset_generic_function(AudioConvert *ac){}voidff_audio_convert_free(AudioConvert **ac){if(!*ac) return;ff_dither_free(&(*ac) ->dc);av_freep(ac);}AudioConvert *ff_audio_convert_alloc(AVAudioResampleContext *avr, enumAVSampleFormatout_fmt, enumAVSampleFormatin_fmt, intchannels, intsample_rate, intapply_map){AudioConvert *ac;intin_planar, out_planar;ac=av_mallocz(sizeof(*ac));if(!ac) returnNULL;ac->avr=avr;ac->out_fmt=out_fmt;ac->in_fmt=in_fmt;ac->channels=channels;ac->apply_map=apply_map;if(avr->dither_method!=AV_RESAMPLE_DITHER_NONE &&av_get_packed_sample_fmt(out_fmt)==AV_SAMPLE_FMT_S16 &&av_get_bytes_per_sample(in_fmt)>2){ac->dc=ff_dither_alloc(avr, out_fmt, in_fmt, channels, sample_rate, apply_map);if(!ac->dc){av_free(ac);returnNULL;}returnac;}in_planar=ff_sample_fmt_is_planar(in_fmt, channels);out_planar=ff_sample_fmt_is_planar(out_fmt, channels);if(in_planar==out_planar){ac->func_type=CONV_FUNC_TYPE_FLAT;ac->planes=in_planar?ac->channels:1;}elseif(in_planar) ac->func_type=CONV_FUNC_TYPE_INTERLEAVE;elseac->func_type=CONV_FUNC_TYPE_DEINTERLEAVE;set_generic_function(ac);if(ARCH_AARCH64) ff_audio_convert_init_aarch64(ac);if(ARCH_ARM) ff_audio_convert_init_arm(ac);if(ARCH_X86) ff_audio_convert_init_x86(ac);returnac;}intff_audio_convert(AudioConvert *ac, AudioData *out, AudioData *in){intuse_generic=1;intlen=in->nb_samples;intp;if(ac->dc){av_log(ac->avr, AV_LOG_TRACE,"%dsamples-audio_convert:%sto%s(dithered)\n", len, av_get_sample_fmt_name(ac->in_fmt), av_get_sample_fmt_name(ac->out_fmt));returnff_convert_dither(ac-> in
static void ttafilter_process(TTAFilter *c, int32_t *in)
Definition: ttaenc.c:63
int ff_alloc_packet2(AVCodecContext *avctx, AVPacket *avpkt, int64_t size, int64_t min_size)
Check AVPacket size and/or allocate data.
Definition: utils.c:1780
uint8_t * data[AV_NUM_DATA_POINTERS]
pointer to the picture/channel planes.
Definition: frame.h:182
const AVCRC * av_crc_get_table(AVCRCId crc_id)
Get an initialized standard CRC table.
Definition: crc.c:342
uint32_t k0
Definition: ttadata.h:35
static void filter(MpegAudioContext *s, int ch, const short *samples, int incr)
common internal api header.
static void flush_put_bits(PutBitContext *s)
Pad the end of the output stream with zeros.
Definition: put_bits.h:101
signed 16 bits
Definition: samplefmt.h:62
static double c[64]
#define AV_CODEC_CAP_LOSSLESS
Codec is lossless.
Definition: avcodec.h:943
static void init_put_bits(PutBitContext *s, uint8_t *buffer, int buffer_size)
Initialize the PutBitContext s.
Definition: put_bits.h:48
void * priv_data
Definition: avcodec.h:1543
TTAFilter filter
Definition: ttadata.h:40
int channels
number of audio channels
Definition: avcodec.h:2262
int32_t dl[MAX_ORDER]
Definition: ttadata.h:31
static enum AVSampleFormat sample_fmts[]
Definition: adpcmenc.c:701
#define av_freep(p)
static av_always_inline int64_t ff_samples_to_time_base(AVCodecContext *avctx, int64_t samples)
Rescale from sample rate to AVCodecContext.time_base.
Definition: internal.h:233
#define av_malloc_array(a, b)
TTARice rice
Definition: ttadata.h:41
void ff_tta_filter_init(TTAFilter *c, int32_t shift)
Definition: ttadata.c:48
This structure stores compressed data.
Definition: avcodec.h:1399
int nb_samples
number of audio samples (per channel) described by this frame
Definition: frame.h:225
uint32_t AVCRC
Definition: crc.h:34
int64_t pts
Presentation timestamp in AVStream->time_base units; the time at which the decompressed packet will b...
Definition: avcodec.h:1415
uint32_t sum0
Definition: ttadata.h:35
static int tta_encode_frame(AVCodecContext *avctx, AVPacket *avpkt, const AVFrame *frame, int *got_packet_ptr)
Definition: ttaenc.c:112
bitstream writer API