0e9e685989390ccd6966c42e554b4e7366a72b90
[ffmpeg.git] / 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     switch(avctx->codec->id) {
99     case CODEC_ID_PCM_U32LE:
100         ENCODE(uint32_t, le32, samples, dst, n, 0, 0x80000000)
101         break;
102     case CODEC_ID_PCM_U32BE:
103         ENCODE(uint32_t, be32, samples, dst, n, 0, 0x80000000)
104         break;
105     case CODEC_ID_PCM_S24LE:
106         ENCODE(int32_t, le24, samples, dst, n, 8, 0)
107         break;
108     case CODEC_ID_PCM_S24BE:
109         ENCODE(int32_t, be24, samples, dst, n, 8, 0)
110         break;
111     case CODEC_ID_PCM_U24LE:
112         ENCODE(uint32_t, le24, samples, dst, n, 8, 0x800000)
113         break;
114     case CODEC_ID_PCM_U24BE:
115         ENCODE(uint32_t, be24, samples, dst, n, 8, 0x800000)
116         break;
117     case CODEC_ID_PCM_S24DAUD:
118         for(;n>0;n--) {
119             uint32_t tmp = av_reverse[(*samples >> 8) & 0xff] +
120                            (av_reverse[*samples & 0xff] << 8);
121             tmp <<= 4; // sync flags would go here
122             bytestream_put_be24(&dst, tmp);
123             samples++;
124         }
125         break;
126     case CODEC_ID_PCM_U16LE:
127         ENCODE(uint16_t, le16, samples, dst, n, 0, 0x8000)
128         break;
129     case CODEC_ID_PCM_U16BE:
130         ENCODE(uint16_t, be16, samples, dst, n, 0, 0x8000)
131         break;
132     case CODEC_ID_PCM_S8:
133         srcu8= data;
134         for(;n>0;n--) {
135             v = *srcu8++;
136             *dst++ = v - 128;
137         }
138         break;
139 #if HAVE_BIGENDIAN
140     case CODEC_ID_PCM_F64LE:
141         ENCODE(int64_t, le64, samples, dst, n, 0, 0)
142         break;
143     case CODEC_ID_PCM_S32LE:
144     case CODEC_ID_PCM_F32LE:
145         ENCODE(int32_t, le32, samples, dst, n, 0, 0)
146         break;
147     case CODEC_ID_PCM_S16LE:
148         ENCODE(int16_t, le16, samples, dst, n, 0, 0)
149         break;
150     case CODEC_ID_PCM_F64BE:
151     case CODEC_ID_PCM_F32BE:
152     case CODEC_ID_PCM_S32BE:
153     case CODEC_ID_PCM_S16BE:
154 #else
155     case CODEC_ID_PCM_F64BE:
156         ENCODE(int64_t, be64, samples, dst, n, 0, 0)
157         break;
158     case CODEC_ID_PCM_F32BE:
159     case CODEC_ID_PCM_S32BE:
160         ENCODE(int32_t, be32, samples, dst, n, 0, 0)
161         break;
162     case CODEC_ID_PCM_S16BE:
163         ENCODE(int16_t, be16, samples, dst, n, 0, 0)
164         break;
165     case CODEC_ID_PCM_F64LE:
166     case CODEC_ID_PCM_F32LE:
167     case CODEC_ID_PCM_S32LE:
168     case CODEC_ID_PCM_S16LE:
169 #endif /* HAVE_BIGENDIAN */
170     case CODEC_ID_PCM_U8:
171         memcpy(dst, samples, n*sample_size);
172         dst += n*sample_size;
173         break;
174     case CODEC_ID_PCM_ALAW:
175         for(;n>0;n--) {
176             v = *samples++;
177             *dst++ = linear_to_alaw[(v + 32768) >> 2];
178         }
179         break;
180     case CODEC_ID_PCM_MULAW:
181         for(;n>0;n--) {
182             v = *samples++;
183             *dst++ = linear_to_ulaw[(v + 32768) >> 2];
184         }
185         break;
186     default:
187         return -1;
188     }
189     //avctx->frame_size = (dst - frame) / (sample_size * avctx->channels);
190
191     return dst - frame;
192 }
193
194 typedef struct PCMDecode {
195     short table[256];
196 } PCMDecode;
197
198 static av_cold int pcm_decode_init(AVCodecContext * avctx)
199 {
200     PCMDecode *s = avctx->priv_data;
201     int i;
202
203     if (avctx->channels <= 0 || avctx->channels > MAX_CHANNELS) {
204         av_log(avctx, AV_LOG_ERROR, "PCM channels out of bounds\n");
205         return AVERROR(EINVAL);
206     }
207
208     switch(avctx->codec->id) {
209     case CODEC_ID_PCM_ALAW:
210         for(i=0;i<256;i++)
211             s->table[i] = alaw2linear(i);
212         break;
213     case CODEC_ID_PCM_MULAW:
214         for(i=0;i<256;i++)
215             s->table[i] = ulaw2linear(i);
216         break;
217     default:
218         break;
219     }
220
221     avctx->sample_fmt = avctx->codec->sample_fmts[0];
222
223     if (avctx->sample_fmt == AV_SAMPLE_FMT_S32)
224         avctx->bits_per_raw_sample = av_get_bits_per_sample(avctx->codec->id);
225
226     return 0;
227 }
228
229 /**
230  * Read PCM samples macro
231  * @param size Data size of native machine format
232  * @param endian bytestream_get_xxx() endian suffix
233  * @param src Source pointer (variable name)
234  * @param dst Destination pointer (variable name)
235  * @param n Total number of samples (variable name)
236  * @param shift Bitshift (bits)
237  * @param offset Sample value offset
238  */
239 #define DECODE(size, endian, src, dst, n, shift, offset) \
240     for(;n>0;n--) { \
241         uint##size##_t v = bytestream_get_##endian(&src); \
242         AV_WN##size##A(dst, (v - offset) << shift); \
243         dst += size / 8; \
244     }
245
246 static int pcm_decode_frame(AVCodecContext *avctx,
247                             void *data, int *data_size,
248                             AVPacket *avpkt)
249 {
250     const uint8_t *src = avpkt->data;
251     int buf_size = avpkt->size;
252     PCMDecode *s = avctx->priv_data;
253     int sample_size, c, n, out_size;
254     uint8_t *samples;
255     int32_t *dst_int32_t;
256
257     samples = data;
258
259     sample_size = av_get_bits_per_sample(avctx->codec_id)/8;
260
261     /* av_get_bits_per_sample returns 0 for CODEC_ID_PCM_DVD */
262     if (CODEC_ID_PCM_DVD == avctx->codec_id) {
263         if (avctx->bits_per_coded_sample != 20 &&
264             avctx->bits_per_coded_sample != 24) {
265             av_log(avctx, AV_LOG_ERROR, "PCM DVD unsupported sample depth\n");
266             return AVERROR(EINVAL);
267         }
268         /* 2 samples are interleaved per block in PCM_DVD */
269         sample_size = avctx->bits_per_coded_sample * 2 / 8;
270     } else if (avctx->codec_id == CODEC_ID_PCM_LXF)
271         /* we process 40-bit blocks per channel for LXF */
272         sample_size = 5;
273
274     if (sample_size == 0) {
275         av_log(avctx, AV_LOG_ERROR, "Invalid sample_size\n");
276         return AVERROR(EINVAL);
277     }
278
279     n = avctx->channels * sample_size;
280
281     if(n && buf_size % n){
282         if (buf_size < n) {
283             av_log(avctx, AV_LOG_ERROR, "invalid PCM packet\n");
284             return -1;
285         }else
286             buf_size -= buf_size % n;
287     }
288
289     n = buf_size/sample_size;
290
291     out_size = n * av_get_bytes_per_sample(avctx->sample_fmt);
292     if (avctx->codec_id == CODEC_ID_PCM_DVD ||
293         avctx->codec_id == CODEC_ID_PCM_LXF)
294         out_size *= 2;
295     if (*data_size < out_size) {
296         av_log(avctx, AV_LOG_ERROR, "output buffer too small\n");
297         return AVERROR(EINVAL);
298     }
299
300     switch(avctx->codec->id) {
301     case CODEC_ID_PCM_U32LE:
302         DECODE(32, le32, src, samples, n, 0, 0x80000000)
303         break;
304     case CODEC_ID_PCM_U32BE:
305         DECODE(32, be32, src, samples, n, 0, 0x80000000)
306         break;
307     case CODEC_ID_PCM_S24LE:
308         DECODE(32, le24, src, samples, n, 8, 0)
309         break;
310     case CODEC_ID_PCM_S24BE:
311         DECODE(32, be24, src, samples, n, 8, 0)
312         break;
313     case CODEC_ID_PCM_U24LE:
314         DECODE(32, le24, src, samples, n, 8, 0x800000)
315         break;
316     case CODEC_ID_PCM_U24BE:
317         DECODE(32, be24, src, samples, n, 8, 0x800000)
318         break;
319     case CODEC_ID_PCM_S24DAUD:
320         for(;n>0;n--) {
321           uint32_t v = bytestream_get_be24(&src);
322           v >>= 4; // sync flags are here
323           AV_WN16A(samples, av_reverse[(v >> 8) & 0xff] +
324                            (av_reverse[v & 0xff] << 8));
325           samples += 2;
326         }
327         break;
328     case CODEC_ID_PCM_S16LE_PLANAR:
329     {
330         const uint8_t *src2[MAX_CHANNELS];
331         n /= avctx->channels;
332         for(c=0;c<avctx->channels;c++)
333             src2[c] = &src[c*n*2];
334         for(;n>0;n--)
335             for(c=0;c<avctx->channels;c++) {
336                 AV_WN16A(samples, bytestream_get_le16(&src2[c]));
337                 samples += 2;
338             }
339         break;
340     }
341     case CODEC_ID_PCM_U16LE:
342         DECODE(16, le16, src, samples, n, 0, 0x8000)
343         break;
344     case CODEC_ID_PCM_U16BE:
345         DECODE(16, be16, src, samples, n, 0, 0x8000)
346         break;
347     case CODEC_ID_PCM_S8:
348         for(;n>0;n--) {
349             *samples++ = *src++ + 128;
350         }
351         break;
352 #if HAVE_BIGENDIAN
353     case CODEC_ID_PCM_F64LE:
354         DECODE(64, le64, src, samples, n, 0, 0)
355         break;
356     case CODEC_ID_PCM_S32LE:
357     case CODEC_ID_PCM_F32LE:
358         DECODE(32, le32, src, samples, n, 0, 0)
359         break;
360     case CODEC_ID_PCM_S16LE:
361         DECODE(16, le16, src, samples, n, 0, 0)
362         break;
363     case CODEC_ID_PCM_F64BE:
364     case CODEC_ID_PCM_F32BE:
365     case CODEC_ID_PCM_S32BE:
366     case CODEC_ID_PCM_S16BE:
367 #else
368     case CODEC_ID_PCM_F64BE:
369         DECODE(64, be64, src, samples, n, 0, 0)
370         break;
371     case CODEC_ID_PCM_F32BE:
372     case CODEC_ID_PCM_S32BE:
373         DECODE(32, be32, src, samples, n, 0, 0)
374         break;
375     case CODEC_ID_PCM_S16BE:
376         DECODE(16, be16, src, samples, n, 0, 0)
377         break;
378     case CODEC_ID_PCM_F64LE:
379     case CODEC_ID_PCM_F32LE:
380     case CODEC_ID_PCM_S32LE:
381     case CODEC_ID_PCM_S16LE:
382 #endif /* HAVE_BIGENDIAN */
383     case CODEC_ID_PCM_U8:
384         memcpy(samples, src, n*sample_size);
385         break;
386     case CODEC_ID_PCM_ZORK:
387         for (; n > 0; n--) {
388             int v = *src++;
389             if (v < 128)
390                 v = 128 - v;
391             *samples++ = v;
392         }
393         break;
394     case CODEC_ID_PCM_ALAW:
395     case CODEC_ID_PCM_MULAW:
396         for(;n>0;n--) {
397             AV_WN16A(samples, s->table[*src++]);
398             samples += 2;
399         }
400         break;
401     case CODEC_ID_PCM_DVD:
402     {
403         const uint8_t *src8;
404         dst_int32_t = data;
405         n /= avctx->channels;
406         switch (avctx->bits_per_coded_sample) {
407         case 20:
408             while (n--) {
409                 c = avctx->channels;
410                 src8 = src + 4*c;
411                 while (c--) {
412                     *dst_int32_t++ = (bytestream_get_be16(&src) << 16) + ((*src8   &0xf0) << 8);
413                     *dst_int32_t++ = (bytestream_get_be16(&src) << 16) + ((*src8++ &0x0f) << 12);
414                 }
415                 src = src8;
416             }
417             break;
418         case 24:
419             while (n--) {
420                 c = avctx->channels;
421                 src8 = src + 4*c;
422                 while (c--) {
423                     *dst_int32_t++ = (bytestream_get_be16(&src) << 16) + ((*src8++) << 8);
424                     *dst_int32_t++ = (bytestream_get_be16(&src) << 16) + ((*src8++) << 8);
425                 }
426                 src = src8;
427             }
428             break;
429         }
430         break;
431     }
432     case CODEC_ID_PCM_LXF:
433     {
434         int i;
435         const uint8_t *src8;
436         dst_int32_t = data;
437         n /= avctx->channels;
438         //unpack and de-planerize
439         for (i = 0; i < n; i++) {
440             for (c = 0, src8 = src + i*5; c < avctx->channels; c++, src8 += n*5) {
441                 //extract low 20 bits and expand to 32 bits
442                 *dst_int32_t++ = (src8[2] << 28) | (src8[1] << 20) | (src8[0] << 12) |
443                                  ((src8[2] & 0xF) << 8) | src8[1];
444             }
445
446             for (c = 0, src8 = src + i*5; c < avctx->channels; c++, src8 += n*5) {
447                 //extract high 20 bits and expand to 32 bits
448                 *dst_int32_t++ = (src8[4] << 24) | (src8[3] << 16) |
449                                  ((src8[2] & 0xF0) << 8) | (src8[4] << 4) | (src8[3] >> 4);
450             }
451         }
452         break;
453     }
454     default:
455         return -1;
456     }
457     *data_size = out_size;
458     return buf_size;
459 }
460
461 #if CONFIG_ENCODERS
462 #define PCM_ENCODER(id_,sample_fmt_,name_,long_name_) \
463 AVCodec ff_ ## name_ ## _encoder = {            \
464     .name        = #name_,                      \
465     .type        = AVMEDIA_TYPE_AUDIO,          \
466     .id          = id_,                         \
467     .init        = pcm_encode_init,             \
468     .encode      = pcm_encode_frame,            \
469     .close       = pcm_encode_close,            \
470     .sample_fmts = (const enum AVSampleFormat[]){sample_fmt_,AV_SAMPLE_FMT_NONE}, \
471     .long_name = NULL_IF_CONFIG_SMALL(long_name_), \
472 }
473 #else
474 #define PCM_ENCODER(id,sample_fmt_,name,long_name_)
475 #endif
476
477 #if CONFIG_DECODERS
478 #define PCM_DECODER(id_,sample_fmt_,name_,long_name_)         \
479 AVCodec ff_ ## name_ ## _decoder = {            \
480     .name           = #name_,                   \
481     .type           = AVMEDIA_TYPE_AUDIO,       \
482     .id             = id_,                      \
483     .priv_data_size = sizeof(PCMDecode),        \
484     .init           = pcm_decode_init,          \
485     .decode         = pcm_decode_frame,         \
486     .sample_fmts = (const enum AVSampleFormat[]){sample_fmt_,AV_SAMPLE_FMT_NONE}, \
487     .long_name = NULL_IF_CONFIG_SMALL(long_name_), \
488 }
489 #else
490 #define PCM_DECODER(id,sample_fmt_,name,long_name_)
491 #endif
492
493 #define PCM_CODEC(id, sample_fmt_, name, long_name_)         \
494     PCM_ENCODER(id,sample_fmt_,name,long_name_); PCM_DECODER(id,sample_fmt_,name,long_name_)
495
496 /* Note: Do not forget to add new entries to the Makefile as well. */
497 PCM_CODEC  (CODEC_ID_PCM_ALAW,  AV_SAMPLE_FMT_S16, pcm_alaw, "PCM A-law");
498 PCM_DECODER(CODEC_ID_PCM_DVD,   AV_SAMPLE_FMT_S32, pcm_dvd, "PCM signed 20|24-bit big-endian");
499 PCM_CODEC  (CODEC_ID_PCM_F32BE, AV_SAMPLE_FMT_FLT, pcm_f32be, "PCM 32-bit floating point big-endian");
500 PCM_CODEC  (CODEC_ID_PCM_F32LE, AV_SAMPLE_FMT_FLT, pcm_f32le, "PCM 32-bit floating point little-endian");
501 PCM_CODEC  (CODEC_ID_PCM_F64BE, AV_SAMPLE_FMT_DBL, pcm_f64be, "PCM 64-bit floating point big-endian");
502 PCM_CODEC  (CODEC_ID_PCM_F64LE, AV_SAMPLE_FMT_DBL, pcm_f64le, "PCM 64-bit floating point little-endian");
503 PCM_DECODER(CODEC_ID_PCM_LXF,   AV_SAMPLE_FMT_S32, pcm_lxf, "PCM signed 20-bit little-endian planar");
504 PCM_CODEC  (CODEC_ID_PCM_MULAW, AV_SAMPLE_FMT_S16, pcm_mulaw, "PCM mu-law");
505 PCM_CODEC  (CODEC_ID_PCM_S8,    AV_SAMPLE_FMT_U8,  pcm_s8, "PCM signed 8-bit");
506 PCM_CODEC  (CODEC_ID_PCM_S16BE, AV_SAMPLE_FMT_S16, pcm_s16be, "PCM signed 16-bit big-endian");
507 PCM_CODEC  (CODEC_ID_PCM_S16LE, AV_SAMPLE_FMT_S16, pcm_s16le, "PCM signed 16-bit little-endian");
508 PCM_DECODER(CODEC_ID_PCM_S16LE_PLANAR, AV_SAMPLE_FMT_S16, pcm_s16le_planar, "PCM 16-bit little-endian planar");
509 PCM_CODEC  (CODEC_ID_PCM_S24BE, AV_SAMPLE_FMT_S32, pcm_s24be, "PCM signed 24-bit big-endian");
510 PCM_CODEC  (CODEC_ID_PCM_S24DAUD, AV_SAMPLE_FMT_S16,  pcm_s24daud, "PCM D-Cinema audio signed 24-bit");
511 PCM_CODEC  (CODEC_ID_PCM_S24LE, AV_SAMPLE_FMT_S32, pcm_s24le, "PCM signed 24-bit little-endian");
512 PCM_CODEC  (CODEC_ID_PCM_S32BE, AV_SAMPLE_FMT_S32, pcm_s32be, "PCM signed 32-bit big-endian");
513 PCM_CODEC  (CODEC_ID_PCM_S32LE, AV_SAMPLE_FMT_S32, pcm_s32le, "PCM signed 32-bit little-endian");
514 PCM_CODEC  (CODEC_ID_PCM_U8,    AV_SAMPLE_FMT_U8,  pcm_u8, "PCM unsigned 8-bit");
515 PCM_CODEC  (CODEC_ID_PCM_U16BE, AV_SAMPLE_FMT_S16, pcm_u16be, "PCM unsigned 16-bit big-endian");
516 PCM_CODEC  (CODEC_ID_PCM_U16LE, AV_SAMPLE_FMT_S16, pcm_u16le, "PCM unsigned 16-bit little-endian");
517 PCM_CODEC  (CODEC_ID_PCM_U24BE, AV_SAMPLE_FMT_S32, pcm_u24be, "PCM unsigned 24-bit big-endian");
518 PCM_CODEC  (CODEC_ID_PCM_U24LE, AV_SAMPLE_FMT_S32, pcm_u24le, "PCM unsigned 24-bit little-endian");
519 PCM_CODEC  (CODEC_ID_PCM_U32BE, AV_SAMPLE_FMT_S32, pcm_u32be, "PCM unsigned 32-bit big-endian");
520 PCM_CODEC  (CODEC_ID_PCM_U32LE, AV_SAMPLE_FMT_S32, pcm_u32le, "PCM unsigned 32-bit little-endian");
521 PCM_DECODER(CODEC_ID_PCM_ZORK,  AV_SAMPLE_FMT_U8,  pcm_zork,  "PCM Zork");