]> git.sesse.net Git - ffmpeg/blob - libavcodec/pcm.c
158366c19abaab32659bfc6d1965a2fa7977e268
[ffmpeg] / libavcodec / pcm.c
1 /*
2  * PCM codecs
3  * Copyright (c) 2001 Fabrice Bellard
4  *
5  * This file is part of Libav.
6  *
7  * Libav is free software; you can redistribute it and/or
8  * modify it under the terms of the GNU Lesser General Public
9  * License as published by the Free Software Foundation; either
10  * version 2.1 of the License, or (at your option) any later version.
11  *
12  * Libav is distributed in the hope that it will be useful,
13  * but WITHOUT ANY WARRANTY; without even the implied warranty of
14  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
15  * Lesser General Public License for more details.
16  *
17  * You should have received a copy of the GNU Lesser General Public
18  * License along with Libav; if not, write to the Free Software
19  * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
20  */
21
22 /**
23  * @file
24  * PCM codecs
25  */
26
27 #include "avcodec.h"
28 #include "libavutil/common.h" /* for av_reverse */
29 #include "bytestream.h"
30 #include "pcm_tablegen.h"
31
32 #define MAX_CHANNELS 64
33
34 static av_cold int pcm_encode_init(AVCodecContext *avctx)
35 {
36     avctx->frame_size = 1;
37     switch(avctx->codec->id) {
38     case CODEC_ID_PCM_ALAW:
39         pcm_alaw_tableinit();
40         break;
41     case CODEC_ID_PCM_MULAW:
42         pcm_ulaw_tableinit();
43         break;
44     default:
45         break;
46     }
47
48     avctx->bits_per_coded_sample = av_get_bits_per_sample(avctx->codec->id);
49     avctx->block_align = avctx->channels * avctx->bits_per_coded_sample/8;
50     avctx->coded_frame= avcodec_alloc_frame();
51     avctx->coded_frame->key_frame= 1;
52
53     return 0;
54 }
55
56 static av_cold int pcm_encode_close(AVCodecContext *avctx)
57 {
58     av_freep(&avctx->coded_frame);
59
60     return 0;
61 }
62
63 /**
64  * Write PCM samples macro
65  * @param type Datatype of native machine format
66  * @param endian bytestream_put_xxx() suffix
67  * @param src Source pointer (variable name)
68  * @param dst Destination pointer (variable name)
69  * @param n Total number of samples (variable name)
70  * @param shift Bitshift (bits)
71  * @param offset Sample value offset
72  */
73 #define ENCODE(type, endian, src, dst, n, shift, offset) \
74     samples_##type = (const type*) src; \
75     for(;n>0;n--) { \
76         register type v = (*samples_##type++ >> shift) + offset; \
77         bytestream_put_##endian(&dst, v); \
78     }
79
80 static int pcm_encode_frame(AVCodecContext *avctx,
81                             unsigned char *frame, int buf_size, void *data)
82 {
83     int n, sample_size, v;
84     const short *samples;
85     unsigned char *dst;
86     const uint8_t *srcu8;
87     const int16_t *samples_int16_t;
88     const int32_t *samples_int32_t;
89     const int64_t *samples_int64_t;
90     const uint16_t *samples_uint16_t;
91     const uint32_t *samples_uint32_t;
92
93     sample_size = av_get_bits_per_sample(avctx->codec->id)/8;
94     n = buf_size / sample_size;
95     samples = data;
96     dst = frame;
97
98     if (avctx->sample_fmt!=avctx->codec->sample_fmts[0]) {
99         av_log(avctx, AV_LOG_ERROR, "invalid sample_fmt\n");
100         return -1;
101     }
102
103     switch(avctx->codec->id) {
104     case CODEC_ID_PCM_U32LE:
105         ENCODE(uint32_t, le32, samples, dst, n, 0, 0x80000000)
106         break;
107     case CODEC_ID_PCM_U32BE:
108         ENCODE(uint32_t, be32, samples, dst, n, 0, 0x80000000)
109         break;
110     case CODEC_ID_PCM_S24LE:
111         ENCODE(int32_t, le24, samples, dst, n, 8, 0)
112         break;
113     case CODEC_ID_PCM_S24BE:
114         ENCODE(int32_t, be24, samples, dst, n, 8, 0)
115         break;
116     case CODEC_ID_PCM_U24LE:
117         ENCODE(uint32_t, le24, samples, dst, n, 8, 0x800000)
118         break;
119     case CODEC_ID_PCM_U24BE:
120         ENCODE(uint32_t, be24, samples, dst, n, 8, 0x800000)
121         break;
122     case CODEC_ID_PCM_S24DAUD:
123         for(;n>0;n--) {
124             uint32_t tmp = av_reverse[(*samples >> 8) & 0xff] +
125                            (av_reverse[*samples & 0xff] << 8);
126             tmp <<= 4; // sync flags would go here
127             bytestream_put_be24(&dst, tmp);
128             samples++;
129         }
130         break;
131     case CODEC_ID_PCM_U16LE:
132         ENCODE(uint16_t, le16, samples, dst, n, 0, 0x8000)
133         break;
134     case CODEC_ID_PCM_U16BE:
135         ENCODE(uint16_t, be16, samples, dst, n, 0, 0x8000)
136         break;
137     case CODEC_ID_PCM_S8:
138         srcu8= data;
139         for(;n>0;n--) {
140             v = *srcu8++;
141             *dst++ = v - 128;
142         }
143         break;
144 #if HAVE_BIGENDIAN
145     case CODEC_ID_PCM_F64LE:
146         ENCODE(int64_t, le64, samples, dst, n, 0, 0)
147         break;
148     case CODEC_ID_PCM_S32LE:
149     case CODEC_ID_PCM_F32LE:
150         ENCODE(int32_t, le32, samples, dst, n, 0, 0)
151         break;
152     case CODEC_ID_PCM_S16LE:
153         ENCODE(int16_t, le16, samples, dst, n, 0, 0)
154         break;
155     case CODEC_ID_PCM_F64BE:
156     case CODEC_ID_PCM_F32BE:
157     case CODEC_ID_PCM_S32BE:
158     case CODEC_ID_PCM_S16BE:
159 #else
160     case CODEC_ID_PCM_F64BE:
161         ENCODE(int64_t, be64, samples, dst, n, 0, 0)
162         break;
163     case CODEC_ID_PCM_F32BE:
164     case CODEC_ID_PCM_S32BE:
165         ENCODE(int32_t, be32, samples, dst, n, 0, 0)
166         break;
167     case CODEC_ID_PCM_S16BE:
168         ENCODE(int16_t, be16, samples, dst, n, 0, 0)
169         break;
170     case CODEC_ID_PCM_F64LE:
171     case CODEC_ID_PCM_F32LE:
172     case CODEC_ID_PCM_S32LE:
173     case CODEC_ID_PCM_S16LE:
174 #endif /* HAVE_BIGENDIAN */
175     case CODEC_ID_PCM_U8:
176         memcpy(dst, samples, n*sample_size);
177         dst += n*sample_size;
178         break;
179     case CODEC_ID_PCM_ZORK:
180         for(;n>0;n--) {
181             v= *samples++ >> 8;
182             if(v<0)   v = -v;
183             else      v+= 128;
184             *dst++ = v;
185         }
186         break;
187     case CODEC_ID_PCM_ALAW:
188         for(;n>0;n--) {
189             v = *samples++;
190             *dst++ = linear_to_alaw[(v + 32768) >> 2];
191         }
192         break;
193     case CODEC_ID_PCM_MULAW:
194         for(;n>0;n--) {
195             v = *samples++;
196             *dst++ = linear_to_ulaw[(v + 32768) >> 2];
197         }
198         break;
199     default:
200         return -1;
201     }
202     //avctx->frame_size = (dst - frame) / (sample_size * avctx->channels);
203
204     return dst - frame;
205 }
206
207 typedef struct PCMDecode {
208     short table[256];
209 } PCMDecode;
210
211 static av_cold int pcm_decode_init(AVCodecContext * avctx)
212 {
213     PCMDecode *s = avctx->priv_data;
214     int i;
215
216     switch(avctx->codec->id) {
217     case CODEC_ID_PCM_ALAW:
218         for(i=0;i<256;i++)
219             s->table[i] = alaw2linear(i);
220         break;
221     case CODEC_ID_PCM_MULAW:
222         for(i=0;i<256;i++)
223             s->table[i] = ulaw2linear(i);
224         break;
225     default:
226         break;
227     }
228
229     avctx->sample_fmt = avctx->codec->sample_fmts[0];
230
231     if (avctx->sample_fmt == AV_SAMPLE_FMT_S32)
232         avctx->bits_per_raw_sample = av_get_bits_per_sample(avctx->codec->id);
233
234     return 0;
235 }
236
237 /**
238  * Read PCM samples macro
239  * @param size Data size of native machine format
240  * @param endian bytestream_get_xxx() endian suffix
241  * @param src Source pointer (variable name)
242  * @param dst Destination pointer (variable name)
243  * @param n Total number of samples (variable name)
244  * @param shift Bitshift (bits)
245  * @param offset Sample value offset
246  */
247 #define DECODE(size, endian, src, dst, n, shift, offset) \
248     for(;n>0;n--) { \
249         uint##size##_t v = bytestream_get_##endian(&src); \
250         AV_WN##size##A(dst, (v - offset) << shift); \
251         dst += size / 8; \
252     }
253
254 static int pcm_decode_frame(AVCodecContext *avctx,
255                             void *data, int *data_size,
256                             AVPacket *avpkt)
257 {
258     const uint8_t *buf = avpkt->data;
259     int buf_size = avpkt->size;
260     PCMDecode *s = avctx->priv_data;
261     int sample_size, c, n, i;
262     uint8_t *samples;
263     const uint8_t *src, *src8, *src2[MAX_CHANNELS];
264     int32_t *dst_int32_t;
265
266     samples = data;
267     src = buf;
268
269     if (avctx->sample_fmt!=avctx->codec->sample_fmts[0]) {
270         av_log(avctx, AV_LOG_ERROR, "invalid sample_fmt\n");
271         return -1;
272     }
273
274     if(avctx->channels <= 0 || avctx->channels > MAX_CHANNELS){
275         av_log(avctx, AV_LOG_ERROR, "PCM channels out of bounds\n");
276         return -1;
277     }
278
279     sample_size = av_get_bits_per_sample(avctx->codec_id)/8;
280
281     /* av_get_bits_per_sample returns 0 for CODEC_ID_PCM_DVD */
282     if (CODEC_ID_PCM_DVD == avctx->codec_id) {
283         if (avctx->bits_per_coded_sample != 20 &&
284             avctx->bits_per_coded_sample != 24) {
285             av_log(avctx, AV_LOG_ERROR, "PCM DVD unsupported sample depth\n");
286             return AVERROR(EINVAL);
287         }
288         /* 2 samples are interleaved per block in PCM_DVD */
289         sample_size = avctx->bits_per_coded_sample * 2 / 8;
290     } else if (avctx->codec_id == CODEC_ID_PCM_LXF)
291         /* we process 40-bit blocks per channel for LXF */
292         sample_size = 5;
293
294     if (sample_size == 0) {
295         av_log(avctx, AV_LOG_ERROR, "Invalid sample_size\n");
296         return AVERROR(EINVAL);
297     }
298
299     n = avctx->channels * sample_size;
300
301     if(n && buf_size % n){
302         if (buf_size < n) {
303             av_log(avctx, AV_LOG_ERROR, "invalid PCM packet\n");
304             return -1;
305         }else
306             buf_size -= buf_size % n;
307     }
308
309     buf_size= FFMIN(buf_size, *data_size/2);
310
311     n = buf_size/sample_size;
312
313     switch(avctx->codec->id) {
314     case CODEC_ID_PCM_U32LE:
315         DECODE(32, le32, src, samples, n, 0, 0x80000000)
316         break;
317     case CODEC_ID_PCM_U32BE:
318         DECODE(32, be32, src, samples, n, 0, 0x80000000)
319         break;
320     case CODEC_ID_PCM_S24LE:
321         DECODE(32, le24, src, samples, n, 8, 0)
322         break;
323     case CODEC_ID_PCM_S24BE:
324         DECODE(32, be24, src, samples, n, 8, 0)
325         break;
326     case CODEC_ID_PCM_U24LE:
327         DECODE(32, le24, src, samples, n, 8, 0x800000)
328         break;
329     case CODEC_ID_PCM_U24BE:
330         DECODE(32, be24, src, samples, n, 8, 0x800000)
331         break;
332     case CODEC_ID_PCM_S24DAUD:
333         for(;n>0;n--) {
334           uint32_t v = bytestream_get_be24(&src);
335           v >>= 4; // sync flags are here
336           AV_WN16A(samples, av_reverse[(v >> 8) & 0xff] +
337                            (av_reverse[v & 0xff] << 8));
338           samples += 2;
339         }
340         break;
341     case CODEC_ID_PCM_S16LE_PLANAR:
342         n /= avctx->channels;
343         for(c=0;c<avctx->channels;c++)
344             src2[c] = &src[c*n*2];
345         for(;n>0;n--)
346             for(c=0;c<avctx->channels;c++) {
347                 AV_WN16A(samples, bytestream_get_le16(&src2[c]));
348                 samples += 2;
349             }
350         src = src2[avctx->channels-1];
351         break;
352     case CODEC_ID_PCM_U16LE:
353         DECODE(16, le16, src, samples, n, 0, 0x8000)
354         break;
355     case CODEC_ID_PCM_U16BE:
356         DECODE(16, be16, src, samples, n, 0, 0x8000)
357         break;
358     case CODEC_ID_PCM_S8:
359         for(;n>0;n--) {
360             *samples++ = *src++ + 128;
361         }
362         break;
363 #if HAVE_BIGENDIAN
364     case CODEC_ID_PCM_F64LE:
365         DECODE(64, le64, src, samples, n, 0, 0)
366         break;
367     case CODEC_ID_PCM_S32LE:
368     case CODEC_ID_PCM_F32LE:
369         DECODE(32, le32, src, samples, n, 0, 0)
370         break;
371     case CODEC_ID_PCM_S16LE:
372         DECODE(16, le16, src, samples, n, 0, 0)
373         break;
374     case CODEC_ID_PCM_F64BE:
375     case CODEC_ID_PCM_F32BE:
376     case CODEC_ID_PCM_S32BE:
377     case CODEC_ID_PCM_S16BE:
378 #else
379     case CODEC_ID_PCM_F64BE:
380         DECODE(64, be64, src, samples, n, 0, 0)
381         break;
382     case CODEC_ID_PCM_F32BE:
383     case CODEC_ID_PCM_S32BE:
384         DECODE(32, be32, src, samples, n, 0, 0)
385         break;
386     case CODEC_ID_PCM_S16BE:
387         DECODE(16, be16, src, samples, n, 0, 0)
388         break;
389     case CODEC_ID_PCM_F64LE:
390     case CODEC_ID_PCM_F32LE:
391     case CODEC_ID_PCM_S32LE:
392     case CODEC_ID_PCM_S16LE:
393 #endif /* HAVE_BIGENDIAN */
394     case CODEC_ID_PCM_U8:
395         memcpy(samples, src, n*sample_size);
396         src += n*sample_size;
397         samples += n * sample_size;
398         break;
399     case CODEC_ID_PCM_ZORK:
400         for(;n>0;n--) {
401             int x= *src++;
402             if(x&128) x-= 128;
403             else      x = -x;
404             AV_WN16A(samples, x << 8);
405             samples += 2;
406         }
407         break;
408     case CODEC_ID_PCM_ALAW:
409     case CODEC_ID_PCM_MULAW:
410         for(;n>0;n--) {
411             AV_WN16A(samples, s->table[*src++]);
412             samples += 2;
413         }
414         break;
415     case CODEC_ID_PCM_DVD:
416         dst_int32_t = data;
417         n /= avctx->channels;
418         switch (avctx->bits_per_coded_sample) {
419         case 20:
420             while (n--) {
421                 c = avctx->channels;
422                 src8 = src + 4*c;
423                 while (c--) {
424                     *dst_int32_t++ = (bytestream_get_be16(&src) << 16) + ((*src8   &0xf0) << 8);
425                     *dst_int32_t++ = (bytestream_get_be16(&src) << 16) + ((*src8++ &0x0f) << 12);
426                 }
427                 src = src8;
428             }
429             break;
430         case 24:
431             while (n--) {
432                 c = avctx->channels;
433                 src8 = src + 4*c;
434                 while (c--) {
435                     *dst_int32_t++ = (bytestream_get_be16(&src) << 16) + ((*src8++) << 8);
436                     *dst_int32_t++ = (bytestream_get_be16(&src) << 16) + ((*src8++) << 8);
437                 }
438                 src = src8;
439             }
440             break;
441         }
442         samples = (uint8_t *) dst_int32_t;
443         break;
444     case CODEC_ID_PCM_LXF:
445         dst_int32_t = data;
446         n /= avctx->channels;
447         //unpack and de-planerize
448         for (i = 0; i < n; i++) {
449             for (c = 0, src8 = src + i*5; c < avctx->channels; c++, src8 += n*5) {
450                 //extract low 20 bits and expand to 32 bits
451                 *dst_int32_t++ = (src8[2] << 28) | (src8[1] << 20) | (src8[0] << 12) |
452                                  ((src8[2] & 0xF) << 8) | src8[1];
453             }
454
455             for (c = 0, src8 = src + i*5; c < avctx->channels; c++, src8 += n*5) {
456                 //extract high 20 bits and expand to 32 bits
457                 *dst_int32_t++ = (src8[4] << 24) | (src8[3] << 16) |
458                                  ((src8[2] & 0xF0) << 8) | (src8[4] << 4) | (src8[3] >> 4);
459             }
460         }
461         src += n * avctx->channels * 5;
462         samples = (uint8_t *) dst_int32_t;
463         break;
464     default:
465         return -1;
466     }
467     *data_size = samples - (uint8_t *)data;
468     return src - buf;
469 }
470
471 #if CONFIG_ENCODERS
472 #define PCM_ENCODER(id_,sample_fmt_,name_,long_name_) \
473 AVCodec ff_ ## name_ ## _encoder = {            \
474     .name        = #name_,                      \
475     .type        = AVMEDIA_TYPE_AUDIO,          \
476     .id          = id_,                         \
477     .init        = pcm_encode_init,             \
478     .encode      = pcm_encode_frame,            \
479     .close       = pcm_encode_close,            \
480     .sample_fmts = (const enum AVSampleFormat[]){sample_fmt_,AV_SAMPLE_FMT_NONE}, \
481     .long_name = NULL_IF_CONFIG_SMALL(long_name_), \
482 }
483 #else
484 #define PCM_ENCODER(id,sample_fmt_,name,long_name_)
485 #endif
486
487 #if CONFIG_DECODERS
488 #define PCM_DECODER(id_,sample_fmt_,name_,long_name_)         \
489 AVCodec ff_ ## name_ ## _decoder = {            \
490     .name           = #name_,                   \
491     .type           = AVMEDIA_TYPE_AUDIO,       \
492     .id             = id_,                      \
493     .priv_data_size = sizeof(PCMDecode),        \
494     .init           = pcm_decode_init,          \
495     .decode         = pcm_decode_frame,         \
496     .sample_fmts = (const enum AVSampleFormat[]){sample_fmt_,AV_SAMPLE_FMT_NONE}, \
497     .long_name = NULL_IF_CONFIG_SMALL(long_name_), \
498 }
499 #else
500 #define PCM_DECODER(id,sample_fmt_,name,long_name_)
501 #endif
502
503 #define PCM_CODEC(id, sample_fmt_, name, long_name_)         \
504     PCM_ENCODER(id,sample_fmt_,name,long_name_); PCM_DECODER(id,sample_fmt_,name,long_name_)
505
506 /* Note: Do not forget to add new entries to the Makefile as well. */
507 PCM_CODEC  (CODEC_ID_PCM_ALAW,  AV_SAMPLE_FMT_S16, pcm_alaw, "PCM A-law");
508 PCM_DECODER(CODEC_ID_PCM_DVD,   AV_SAMPLE_FMT_S32, pcm_dvd, "PCM signed 20|24-bit big-endian");
509 PCM_CODEC  (CODEC_ID_PCM_F32BE, AV_SAMPLE_FMT_FLT, pcm_f32be, "PCM 32-bit floating point big-endian");
510 PCM_CODEC  (CODEC_ID_PCM_F32LE, AV_SAMPLE_FMT_FLT, pcm_f32le, "PCM 32-bit floating point little-endian");
511 PCM_CODEC  (CODEC_ID_PCM_F64BE, AV_SAMPLE_FMT_DBL, pcm_f64be, "PCM 64-bit floating point big-endian");
512 PCM_CODEC  (CODEC_ID_PCM_F64LE, AV_SAMPLE_FMT_DBL, pcm_f64le, "PCM 64-bit floating point little-endian");
513 PCM_DECODER(CODEC_ID_PCM_LXF,   AV_SAMPLE_FMT_S32, pcm_lxf, "PCM signed 20-bit little-endian planar");
514 PCM_CODEC  (CODEC_ID_PCM_MULAW, AV_SAMPLE_FMT_S16, pcm_mulaw, "PCM mu-law");
515 PCM_CODEC  (CODEC_ID_PCM_S8,    AV_SAMPLE_FMT_U8,  pcm_s8, "PCM signed 8-bit");
516 PCM_CODEC  (CODEC_ID_PCM_S16BE, AV_SAMPLE_FMT_S16, pcm_s16be, "PCM signed 16-bit big-endian");
517 PCM_CODEC  (CODEC_ID_PCM_S16LE, AV_SAMPLE_FMT_S16, pcm_s16le, "PCM signed 16-bit little-endian");
518 PCM_DECODER(CODEC_ID_PCM_S16LE_PLANAR, AV_SAMPLE_FMT_S16, pcm_s16le_planar, "PCM 16-bit little-endian planar");
519 PCM_CODEC  (CODEC_ID_PCM_S24BE, AV_SAMPLE_FMT_S32, pcm_s24be, "PCM signed 24-bit big-endian");
520 PCM_CODEC  (CODEC_ID_PCM_S24DAUD, AV_SAMPLE_FMT_S16,  pcm_s24daud, "PCM D-Cinema audio signed 24-bit");
521 PCM_CODEC  (CODEC_ID_PCM_S24LE, AV_SAMPLE_FMT_S32, pcm_s24le, "PCM signed 24-bit little-endian");
522 PCM_CODEC  (CODEC_ID_PCM_S32BE, AV_SAMPLE_FMT_S32, pcm_s32be, "PCM signed 32-bit big-endian");
523 PCM_CODEC  (CODEC_ID_PCM_S32LE, AV_SAMPLE_FMT_S32, pcm_s32le, "PCM signed 32-bit little-endian");
524 PCM_CODEC  (CODEC_ID_PCM_U8,    AV_SAMPLE_FMT_U8,  pcm_u8, "PCM unsigned 8-bit");
525 PCM_CODEC  (CODEC_ID_PCM_U16BE, AV_SAMPLE_FMT_S16, pcm_u16be, "PCM unsigned 16-bit big-endian");
526 PCM_CODEC  (CODEC_ID_PCM_U16LE, AV_SAMPLE_FMT_S16, pcm_u16le, "PCM unsigned 16-bit little-endian");
527 PCM_CODEC  (CODEC_ID_PCM_U24BE, AV_SAMPLE_FMT_S32, pcm_u24be, "PCM unsigned 24-bit big-endian");
528 PCM_CODEC  (CODEC_ID_PCM_U24LE, AV_SAMPLE_FMT_S32, pcm_u24le, "PCM unsigned 24-bit little-endian");
529 PCM_CODEC  (CODEC_ID_PCM_U32BE, AV_SAMPLE_FMT_S32, pcm_u32be, "PCM unsigned 32-bit big-endian");
530 PCM_CODEC  (CODEC_ID_PCM_U32LE, AV_SAMPLE_FMT_S32, pcm_u32le, "PCM unsigned 32-bit little-endian");
531 PCM_CODEC  (CODEC_ID_PCM_ZORK,  AV_SAMPLE_FMT_S16, pcm_zork, "PCM Zork");