adpcm: Store the trellis nodes in a heap instead of a sorted array
[ffmpeg.git] / libavcodec / adpcm.c
1 /*
2  * ADPCM codecs
3  * Copyright (c) 2001-2003 The ffmpeg Project
4  *
5  * This file is part of FFmpeg.
6  *
7  * FFmpeg 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  * FFmpeg 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 FFmpeg; if not, write to the Free Software
19  * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
20  */
21 #include "avcodec.h"
22 #include "get_bits.h"
23 #include "put_bits.h"
24 #include "bytestream.h"
25
26 /**
27  * @file
28  * ADPCM codecs.
29  * First version by Francois Revol (revol@free.fr)
30  * Fringe ADPCM codecs (e.g., DK3, DK4, Westwood)
31  *   by Mike Melanson (melanson@pcisys.net)
32  * CD-ROM XA ADPCM codec by BERO
33  * EA ADPCM decoder by Robin Kay (komadori@myrealbox.com)
34  * EA ADPCM R1/R2/R3 decoder by Peter Ross (pross@xvid.org)
35  * EA IMA EACS decoder by Peter Ross (pross@xvid.org)
36  * EA IMA SEAD decoder by Peter Ross (pross@xvid.org)
37  * EA ADPCM XAS decoder by Peter Ross (pross@xvid.org)
38  * MAXIS EA ADPCM decoder by Robert Marston (rmarston@gmail.com)
39  * THP ADPCM decoder by Marco Gerards (mgerards@xs4all.nl)
40  *
41  * Features and limitations:
42  *
43  * Reference documents:
44  * http://www.pcisys.net/~melanson/codecs/simpleaudio.html
45  * http://www.geocities.com/SiliconValley/8682/aud3.txt
46  * http://openquicktime.sourceforge.net/plugins.htm
47  * XAnim sources (xa_codec.c) http://www.rasnaimaging.com/people/lapus/download.html
48  * http://www.cs.ucla.edu/~leec/mediabench/applications.html
49  * SoX source code http://home.sprynet.com/~cbagwell/sox.html
50  *
51  * CD-ROM XA:
52  * http://ku-www.ss.titech.ac.jp/~yatsushi/xaadpcm.html
53  * vagpack & depack http://homepages.compuserve.de/bITmASTER32/psx-index.html
54  * readstr http://www.geocities.co.jp/Playtown/2004/
55  */
56
57 #define BLKSIZE 1024
58
59 /* step_table[] and index_table[] are from the ADPCM reference source */
60 /* This is the index table: */
61 static const int index_table[16] = {
62     -1, -1, -1, -1, 2, 4, 6, 8,
63     -1, -1, -1, -1, 2, 4, 6, 8,
64 };
65
66 /**
67  * This is the step table. Note that many programs use slight deviations from
68  * this table, but such deviations are negligible:
69  */
70 static const int step_table[89] = {
71     7, 8, 9, 10, 11, 12, 13, 14, 16, 17,
72     19, 21, 23, 25, 28, 31, 34, 37, 41, 45,
73     50, 55, 60, 66, 73, 80, 88, 97, 107, 118,
74     130, 143, 157, 173, 190, 209, 230, 253, 279, 307,
75     337, 371, 408, 449, 494, 544, 598, 658, 724, 796,
76     876, 963, 1060, 1166, 1282, 1411, 1552, 1707, 1878, 2066,
77     2272, 2499, 2749, 3024, 3327, 3660, 4026, 4428, 4871, 5358,
78     5894, 6484, 7132, 7845, 8630, 9493, 10442, 11487, 12635, 13899,
79     15289, 16818, 18500, 20350, 22385, 24623, 27086, 29794, 32767
80 };
81
82 /* These are for MS-ADPCM */
83 /* AdaptationTable[], AdaptCoeff1[], and AdaptCoeff2[] are from libsndfile */
84 static const int AdaptationTable[] = {
85         230, 230, 230, 230, 307, 409, 512, 614,
86         768, 614, 512, 409, 307, 230, 230, 230
87 };
88
89 /** Divided by 4 to fit in 8-bit integers */
90 static const uint8_t AdaptCoeff1[] = {
91         64, 128, 0, 48, 60, 115, 98
92 };
93
94 /** Divided by 4 to fit in 8-bit integers */
95 static const int8_t AdaptCoeff2[] = {
96         0, -64, 0, 16, 0, -52, -58
97 };
98
99 /* These are for CD-ROM XA ADPCM */
100 static const int xa_adpcm_table[5][2] = {
101    {   0,   0 },
102    {  60,   0 },
103    { 115, -52 },
104    {  98, -55 },
105    { 122, -60 }
106 };
107
108 static const int ea_adpcm_table[] = {
109     0, 240, 460, 392, 0, 0, -208, -220, 0, 1,
110     3, 4, 7, 8, 10, 11, 0, -1, -3, -4
111 };
112
113 // padded to zero where table size is less then 16
114 static const int swf_index_tables[4][16] = {
115     /*2*/ { -1, 2 },
116     /*3*/ { -1, -1, 2, 4 },
117     /*4*/ { -1, -1, -1, -1, 2, 4, 6, 8 },
118     /*5*/ { -1, -1, -1, -1, -1, -1, -1, -1, 1, 2, 4, 6, 8, 10, 13, 16 }
119 };
120
121 static const int yamaha_indexscale[] = {
122     230, 230, 230, 230, 307, 409, 512, 614,
123     230, 230, 230, 230, 307, 409, 512, 614
124 };
125
126 static const int yamaha_difflookup[] = {
127     1, 3, 5, 7, 9, 11, 13, 15,
128     -1, -3, -5, -7, -9, -11, -13, -15
129 };
130
131 /* end of tables */
132
133 typedef struct ADPCMChannelStatus {
134     int predictor;
135     short int step_index;
136     int step;
137     /* for encoding */
138     int prev_sample;
139
140     /* MS version */
141     short sample1;
142     short sample2;
143     int coeff1;
144     int coeff2;
145     int idelta;
146 } ADPCMChannelStatus;
147
148 typedef struct TrellisPath {
149     int nibble;
150     int prev;
151 } TrellisPath;
152
153 typedef struct TrellisNode {
154     uint32_t ssd;
155     int path;
156     int sample1;
157     int sample2;
158     int step;
159 } TrellisNode;
160
161 typedef struct ADPCMContext {
162     ADPCMChannelStatus status[6];
163     TrellisPath *paths;
164     TrellisNode *node_buf;
165     TrellisNode **nodep_buf;
166 } ADPCMContext;
167
168 #define FREEZE_INTERVAL 128
169
170 /* XXX: implement encoding */
171
172 #if CONFIG_ENCODERS
173 static av_cold int adpcm_encode_init(AVCodecContext *avctx)
174 {
175     ADPCMContext *s = avctx->priv_data;
176     uint8_t *extradata;
177     int i;
178     if (avctx->channels > 2)
179         return -1; /* only stereo or mono =) */
180
181     if(avctx->trellis && (unsigned)avctx->trellis > 16U){
182         av_log(avctx, AV_LOG_ERROR, "invalid trellis size\n");
183         return -1;
184     }
185
186     if (avctx->trellis) {
187         int frontier = 1 << avctx->trellis;
188         int max_paths =  frontier * FREEZE_INTERVAL;
189         FF_ALLOC_OR_GOTO(avctx, s->paths,     max_paths * sizeof(*s->paths), error);
190         FF_ALLOC_OR_GOTO(avctx, s->node_buf,  2 * frontier * sizeof(*s->node_buf), error);
191         FF_ALLOC_OR_GOTO(avctx, s->nodep_buf, 2 * frontier * sizeof(*s->nodep_buf), error);
192     }
193
194     switch(avctx->codec->id) {
195     case CODEC_ID_ADPCM_IMA_WAV:
196         avctx->frame_size = (BLKSIZE - 4 * avctx->channels) * 8 / (4 * avctx->channels) + 1; /* each 16 bits sample gives one nibble */
197                                                              /* and we have 4 bytes per channel overhead */
198         avctx->block_align = BLKSIZE;
199         /* seems frame_size isn't taken into account... have to buffer the samples :-( */
200         break;
201     case CODEC_ID_ADPCM_IMA_QT:
202         avctx->frame_size = 64;
203         avctx->block_align = 34 * avctx->channels;
204         break;
205     case CODEC_ID_ADPCM_MS:
206         avctx->frame_size = (BLKSIZE - 7 * avctx->channels) * 2 / avctx->channels + 2; /* each 16 bits sample gives one nibble */
207                                                              /* and we have 7 bytes per channel overhead */
208         avctx->block_align = BLKSIZE;
209         avctx->extradata_size = 32;
210         extradata = avctx->extradata = av_malloc(avctx->extradata_size);
211         if (!extradata)
212             return AVERROR(ENOMEM);
213         bytestream_put_le16(&extradata, avctx->frame_size);
214         bytestream_put_le16(&extradata, 7); /* wNumCoef */
215         for (i = 0; i < 7; i++) {
216             bytestream_put_le16(&extradata, AdaptCoeff1[i] * 4);
217             bytestream_put_le16(&extradata, AdaptCoeff2[i] * 4);
218         }
219         break;
220     case CODEC_ID_ADPCM_YAMAHA:
221         avctx->frame_size = BLKSIZE * avctx->channels;
222         avctx->block_align = BLKSIZE;
223         break;
224     case CODEC_ID_ADPCM_SWF:
225         if (avctx->sample_rate != 11025 &&
226             avctx->sample_rate != 22050 &&
227             avctx->sample_rate != 44100) {
228             av_log(avctx, AV_LOG_ERROR, "Sample rate must be 11025, 22050 or 44100\n");
229             goto error;
230         }
231         avctx->frame_size = 512 * (avctx->sample_rate / 11025);
232         break;
233     default:
234         goto error;
235     }
236
237     avctx->coded_frame= avcodec_alloc_frame();
238     avctx->coded_frame->key_frame= 1;
239
240     return 0;
241 error:
242     av_freep(&s->paths);
243     av_freep(&s->node_buf);
244     av_freep(&s->nodep_buf);
245     return -1;
246 }
247
248 static av_cold int adpcm_encode_close(AVCodecContext *avctx)
249 {
250     ADPCMContext *s = avctx->priv_data;
251     av_freep(&avctx->coded_frame);
252     av_freep(&s->paths);
253     av_freep(&s->node_buf);
254     av_freep(&s->nodep_buf);
255
256     return 0;
257 }
258
259
260 static inline unsigned char adpcm_ima_compress_sample(ADPCMChannelStatus *c, short sample)
261 {
262     int delta = sample - c->prev_sample;
263     int nibble = FFMIN(7, abs(delta)*4/step_table[c->step_index]) + (delta<0)*8;
264     c->prev_sample += ((step_table[c->step_index] * yamaha_difflookup[nibble]) / 8);
265     c->prev_sample = av_clip_int16(c->prev_sample);
266     c->step_index = av_clip(c->step_index + index_table[nibble], 0, 88);
267     return nibble;
268 }
269
270 static inline unsigned char adpcm_ms_compress_sample(ADPCMChannelStatus *c, short sample)
271 {
272     int predictor, nibble, bias;
273
274     predictor = (((c->sample1) * (c->coeff1)) + ((c->sample2) * (c->coeff2))) / 64;
275
276     nibble= sample - predictor;
277     if(nibble>=0) bias= c->idelta/2;
278     else          bias=-c->idelta/2;
279
280     nibble= (nibble + bias) / c->idelta;
281     nibble= av_clip(nibble, -8, 7)&0x0F;
282
283     predictor += (signed)((nibble & 0x08)?(nibble - 0x10):(nibble)) * c->idelta;
284
285     c->sample2 = c->sample1;
286     c->sample1 = av_clip_int16(predictor);
287
288     c->idelta = (AdaptationTable[(int)nibble] * c->idelta) >> 8;
289     if (c->idelta < 16) c->idelta = 16;
290
291     return nibble;
292 }
293
294 static inline unsigned char adpcm_yamaha_compress_sample(ADPCMChannelStatus *c, short sample)
295 {
296     int nibble, delta;
297
298     if(!c->step) {
299         c->predictor = 0;
300         c->step = 127;
301     }
302
303     delta = sample - c->predictor;
304
305     nibble = FFMIN(7, abs(delta)*4/c->step) + (delta<0)*8;
306
307     c->predictor += ((c->step * yamaha_difflookup[nibble]) / 8);
308     c->predictor = av_clip_int16(c->predictor);
309     c->step = (c->step * yamaha_indexscale[nibble]) >> 8;
310     c->step = av_clip(c->step, 127, 24567);
311
312     return nibble;
313 }
314
315 static void adpcm_compress_trellis(AVCodecContext *avctx, const short *samples,
316                                    uint8_t *dst, ADPCMChannelStatus *c, int n)
317 {
318     //FIXME 6% faster if frontier is a compile-time constant
319     ADPCMContext *s = avctx->priv_data;
320     const int frontier = 1 << avctx->trellis;
321     const int stride = avctx->channels;
322     const int version = avctx->codec->id;
323     TrellisPath *paths = s->paths, *p;
324     TrellisNode *node_buf = s->node_buf;
325     TrellisNode **nodep_buf = s->nodep_buf;
326     TrellisNode **nodes = nodep_buf; // nodes[] is always sorted by .ssd
327     TrellisNode **nodes_next = nodep_buf + frontier;
328     int pathn = 0, froze = -1, i, j, k;
329
330     memset(nodep_buf, 0, 2 * frontier * sizeof(*nodep_buf));
331     nodes[0] = node_buf + frontier;
332     nodes[0]->ssd = 0;
333     nodes[0]->path = 0;
334     nodes[0]->step = c->step_index;
335     nodes[0]->sample1 = c->sample1;
336     nodes[0]->sample2 = c->sample2;
337     if((version == CODEC_ID_ADPCM_IMA_WAV) || (version == CODEC_ID_ADPCM_IMA_QT) || (version == CODEC_ID_ADPCM_SWF))
338         nodes[0]->sample1 = c->prev_sample;
339     if(version == CODEC_ID_ADPCM_MS)
340         nodes[0]->step = c->idelta;
341     if(version == CODEC_ID_ADPCM_YAMAHA) {
342         if(c->step == 0) {
343             nodes[0]->step = 127;
344             nodes[0]->sample1 = 0;
345         } else {
346             nodes[0]->step = c->step;
347             nodes[0]->sample1 = c->predictor;
348         }
349     }
350
351     for(i=0; i<n; i++) {
352         TrellisNode *t = node_buf + frontier*(i&1);
353         TrellisNode **u;
354         int sample = samples[i*stride];
355         int heap_pos = 0;
356         memset(nodes_next, 0, frontier*sizeof(TrellisNode*));
357         for(j=0; j<frontier && nodes[j]; j++) {
358             // higher j have higher ssd already, so they're unlikely to use a suboptimal next sample too
359             const int range = (j < frontier/2) ? 1 : 0;
360             const int step = nodes[j]->step;
361             int nidx;
362             if(version == CODEC_ID_ADPCM_MS) {
363                 const int predictor = ((nodes[j]->sample1 * c->coeff1) + (nodes[j]->sample2 * c->coeff2)) / 64;
364                 const int div = (sample - predictor) / step;
365                 const int nmin = av_clip(div-range, -8, 6);
366                 const int nmax = av_clip(div+range, -7, 7);
367                 for(nidx=nmin; nidx<=nmax; nidx++) {
368                     const int nibble = nidx & 0xf;
369                     int dec_sample = predictor + nidx * step;
370 #define STORE_NODE(NAME, STEP_INDEX)\
371                     int d;\
372                     uint32_t ssd;\
373                     int pos;\
374                     TrellisNode *u;\
375                     dec_sample = av_clip_int16(dec_sample);\
376                     d = sample - dec_sample;\
377                     ssd = nodes[j]->ssd + d*d;\
378                     /* Collapse any two states with the same previous sample value. \
379                      * One could also distinguish states by step and by 2nd to last
380                      * sample, but the effects of that are negligible. */\
381                     for(k=0; k<frontier && nodes_next[k]; k++) {\
382                         if(dec_sample == nodes_next[k]->sample1) {\
383                             assert(ssd >= nodes_next[k]->ssd);\
384                             goto next_##NAME;\
385                         }\
386                     }\
387                     if (heap_pos < frontier) {\
388                         pos = heap_pos++;\
389                     } else {\
390                         /* Find the largest node in the heap, which is one \
391                          * of the leaf nodes. */\
392                         int maxpos = 0;\
393                         uint32_t max_ssd = 0;\
394                         for (k = frontier >> 1; k < frontier; k++) {\
395                             if (nodes_next[k]->ssd > max_ssd) {\
396                                 maxpos = k;\
397                                 max_ssd = nodes_next[k]->ssd;\
398                             }\
399                         }\
400                         pos = maxpos;\
401                         if (ssd > max_ssd)\
402                             goto next_##NAME;\
403                     }\
404                     u = nodes_next[pos];\
405                             if(!u) {\
406                                 assert(pathn < FREEZE_INTERVAL<<avctx->trellis);\
407                                 u = t++;\
408                                 nodes_next[pos] = u;\
409                                 u->path = pathn++;\
410                             }\
411                             u->ssd = ssd;\
412                             u->step = STEP_INDEX;\
413                             u->sample2 = nodes[j]->sample1;\
414                             u->sample1 = dec_sample;\
415                             paths[u->path].nibble = nibble;\
416                             paths[u->path].prev = nodes[j]->path;\
417                     /* Sift the newly inserted node down in the heap to \
418                      * restore the heap property. */\
419                     while (pos > 0) {\
420                         int parent = (pos - 1) >> 1;\
421                         if (nodes_next[parent]->ssd <= ssd)\
422                             break;\
423                         FFSWAP(TrellisNode*, nodes_next[parent], nodes_next[pos]);\
424                         pos = parent;\
425                     }\
426                     next_##NAME:;
427                     STORE_NODE(ms, FFMAX(16, (AdaptationTable[nibble] * step) >> 8));
428                 }
429             } else if((version == CODEC_ID_ADPCM_IMA_WAV)|| (version == CODEC_ID_ADPCM_IMA_QT)|| (version == CODEC_ID_ADPCM_SWF)) {
430 #define LOOP_NODES(NAME, STEP_TABLE, STEP_INDEX)\
431                 const int predictor = nodes[j]->sample1;\
432                 const int div = (sample - predictor) * 4 / STEP_TABLE;\
433                 int nmin = av_clip(div-range, -7, 6);\
434                 int nmax = av_clip(div+range, -6, 7);\
435                 if(nmin<=0) nmin--; /* distinguish -0 from +0 */\
436                 if(nmax<0) nmax--;\
437                 for(nidx=nmin; nidx<=nmax; nidx++) {\
438                     const int nibble = nidx<0 ? 7-nidx : nidx;\
439                     int dec_sample = predictor + (STEP_TABLE * yamaha_difflookup[nibble]) / 8;\
440                     STORE_NODE(NAME, STEP_INDEX);\
441                 }
442                 LOOP_NODES(ima, step_table[step], av_clip(step + index_table[nibble], 0, 88));
443             } else { //CODEC_ID_ADPCM_YAMAHA
444                 LOOP_NODES(yamaha, step, av_clip((step * yamaha_indexscale[nibble]) >> 8, 127, 24567));
445 #undef LOOP_NODES
446 #undef STORE_NODE
447             }
448         }
449
450         u = nodes;
451         nodes = nodes_next;
452         nodes_next = u;
453
454         // prevent overflow
455         if(nodes[0]->ssd > (1<<28)) {
456             for(j=1; j<frontier && nodes[j]; j++)
457                 nodes[j]->ssd -= nodes[0]->ssd;
458             nodes[0]->ssd = 0;
459         }
460
461         // merge old paths to save memory
462         if(i == froze + FREEZE_INTERVAL) {
463             p = &paths[nodes[0]->path];
464             for(k=i; k>froze; k--) {
465                 dst[k] = p->nibble;
466                 p = &paths[p->prev];
467             }
468             froze = i;
469             pathn = 0;
470             // other nodes might use paths that don't coincide with the frozen one.
471             // checking which nodes do so is too slow, so just kill them all.
472             // this also slightly improves quality, but I don't know why.
473             memset(nodes+1, 0, (frontier-1)*sizeof(TrellisNode*));
474         }
475     }
476
477     p = &paths[nodes[0]->path];
478     for(i=n-1; i>froze; i--) {
479         dst[i] = p->nibble;
480         p = &paths[p->prev];
481     }
482
483     c->predictor = nodes[0]->sample1;
484     c->sample1 = nodes[0]->sample1;
485     c->sample2 = nodes[0]->sample2;
486     c->step_index = nodes[0]->step;
487     c->step = nodes[0]->step;
488     c->idelta = nodes[0]->step;
489 }
490
491 static int adpcm_encode_frame(AVCodecContext *avctx,
492                             unsigned char *frame, int buf_size, void *data)
493 {
494     int n, i, st;
495     short *samples;
496     unsigned char *dst;
497     ADPCMContext *c = avctx->priv_data;
498     uint8_t *buf;
499
500     dst = frame;
501     samples = (short *)data;
502     st= avctx->channels == 2;
503 /*    n = (BLKSIZE - 4 * avctx->channels) / (2 * 8 * avctx->channels); */
504
505     switch(avctx->codec->id) {
506     case CODEC_ID_ADPCM_IMA_WAV:
507         n = avctx->frame_size / 8;
508             c->status[0].prev_sample = (signed short)samples[0]; /* XXX */
509 /*            c->status[0].step_index = 0; *//* XXX: not sure how to init the state machine */
510             bytestream_put_le16(&dst, c->status[0].prev_sample);
511             *dst++ = (unsigned char)c->status[0].step_index;
512             *dst++ = 0; /* unknown */
513             samples++;
514             if (avctx->channels == 2) {
515                 c->status[1].prev_sample = (signed short)samples[0];
516 /*                c->status[1].step_index = 0; */
517                 bytestream_put_le16(&dst, c->status[1].prev_sample);
518                 *dst++ = (unsigned char)c->status[1].step_index;
519                 *dst++ = 0;
520                 samples++;
521             }
522
523             /* stereo: 4 bytes (8 samples) for left, 4 bytes for right, 4 bytes left, ... */
524             if(avctx->trellis > 0) {
525                 FF_ALLOC_OR_GOTO(avctx, buf, 2*n*8, error);
526                 adpcm_compress_trellis(avctx, samples, buf, &c->status[0], n*8);
527                 if(avctx->channels == 2)
528                     adpcm_compress_trellis(avctx, samples+1, buf + n*8, &c->status[1], n*8);
529                 for(i=0; i<n; i++) {
530                     *dst++ = buf[8*i+0] | (buf[8*i+1] << 4);
531                     *dst++ = buf[8*i+2] | (buf[8*i+3] << 4);
532                     *dst++ = buf[8*i+4] | (buf[8*i+5] << 4);
533                     *dst++ = buf[8*i+6] | (buf[8*i+7] << 4);
534                     if (avctx->channels == 2) {
535                         uint8_t *buf1 = buf + n*8;
536                         *dst++ = buf1[8*i+0] | (buf1[8*i+1] << 4);
537                         *dst++ = buf1[8*i+2] | (buf1[8*i+3] << 4);
538                         *dst++ = buf1[8*i+4] | (buf1[8*i+5] << 4);
539                         *dst++ = buf1[8*i+6] | (buf1[8*i+7] << 4);
540                     }
541                 }
542                 av_free(buf);
543             } else
544             for (; n>0; n--) {
545                 *dst = adpcm_ima_compress_sample(&c->status[0], samples[0]);
546                 *dst |= adpcm_ima_compress_sample(&c->status[0], samples[avctx->channels]) << 4;
547                 dst++;
548                 *dst = adpcm_ima_compress_sample(&c->status[0], samples[avctx->channels * 2]);
549                 *dst |= adpcm_ima_compress_sample(&c->status[0], samples[avctx->channels * 3]) << 4;
550                 dst++;
551                 *dst = adpcm_ima_compress_sample(&c->status[0], samples[avctx->channels * 4]);
552                 *dst |= adpcm_ima_compress_sample(&c->status[0], samples[avctx->channels * 5]) << 4;
553                 dst++;
554                 *dst = adpcm_ima_compress_sample(&c->status[0], samples[avctx->channels * 6]);
555                 *dst |= adpcm_ima_compress_sample(&c->status[0], samples[avctx->channels * 7]) << 4;
556                 dst++;
557                 /* right channel */
558                 if (avctx->channels == 2) {
559                     *dst = adpcm_ima_compress_sample(&c->status[1], samples[1]);
560                     *dst |= adpcm_ima_compress_sample(&c->status[1], samples[3]) << 4;
561                     dst++;
562                     *dst = adpcm_ima_compress_sample(&c->status[1], samples[5]);
563                     *dst |= adpcm_ima_compress_sample(&c->status[1], samples[7]) << 4;
564                     dst++;
565                     *dst = adpcm_ima_compress_sample(&c->status[1], samples[9]);
566                     *dst |= adpcm_ima_compress_sample(&c->status[1], samples[11]) << 4;
567                     dst++;
568                     *dst = adpcm_ima_compress_sample(&c->status[1], samples[13]);
569                     *dst |= adpcm_ima_compress_sample(&c->status[1], samples[15]) << 4;
570                     dst++;
571                 }
572                 samples += 8 * avctx->channels;
573             }
574         break;
575     case CODEC_ID_ADPCM_IMA_QT:
576     {
577         int ch, i;
578         PutBitContext pb;
579         init_put_bits(&pb, dst, buf_size*8);
580
581         for(ch=0; ch<avctx->channels; ch++){
582             put_bits(&pb, 9, (c->status[ch].prev_sample + 0x10000) >> 7);
583             put_bits(&pb, 7, c->status[ch].step_index);
584             if(avctx->trellis > 0) {
585                 uint8_t buf[64];
586                 adpcm_compress_trellis(avctx, samples+ch, buf, &c->status[ch], 64);
587                 for(i=0; i<64; i++)
588                     put_bits(&pb, 4, buf[i^1]);
589                 c->status[ch].prev_sample = c->status[ch].predictor & ~0x7F;
590             } else {
591                 for (i=0; i<64; i+=2){
592                     int t1, t2;
593                     t1 = adpcm_ima_compress_sample(&c->status[ch], samples[avctx->channels*(i+0)+ch]);
594                     t2 = adpcm_ima_compress_sample(&c->status[ch], samples[avctx->channels*(i+1)+ch]);
595                     put_bits(&pb, 4, t2);
596                     put_bits(&pb, 4, t1);
597                 }
598                 c->status[ch].prev_sample &= ~0x7F;
599             }
600         }
601
602         flush_put_bits(&pb);
603         dst += put_bits_count(&pb)>>3;
604         break;
605     }
606     case CODEC_ID_ADPCM_SWF:
607     {
608         int i;
609         PutBitContext pb;
610         init_put_bits(&pb, dst, buf_size*8);
611
612         n = avctx->frame_size-1;
613
614         //Store AdpcmCodeSize
615         put_bits(&pb, 2, 2);                //Set 4bits flash adpcm format
616
617         //Init the encoder state
618         for(i=0; i<avctx->channels; i++){
619             c->status[i].step_index = av_clip(c->status[i].step_index, 0, 63); // clip step so it fits 6 bits
620             put_sbits(&pb, 16, samples[i]);
621             put_bits(&pb, 6, c->status[i].step_index);
622             c->status[i].prev_sample = (signed short)samples[i];
623         }
624
625         if(avctx->trellis > 0) {
626             FF_ALLOC_OR_GOTO(avctx, buf, 2*n, error);
627             adpcm_compress_trellis(avctx, samples+2, buf, &c->status[0], n);
628             if (avctx->channels == 2)
629                 adpcm_compress_trellis(avctx, samples+3, buf+n, &c->status[1], n);
630             for(i=0; i<n; i++) {
631                 put_bits(&pb, 4, buf[i]);
632                 if (avctx->channels == 2)
633                     put_bits(&pb, 4, buf[n+i]);
634             }
635             av_free(buf);
636         } else {
637             for (i=1; i<avctx->frame_size; i++) {
638                 put_bits(&pb, 4, adpcm_ima_compress_sample(&c->status[0], samples[avctx->channels*i]));
639                 if (avctx->channels == 2)
640                     put_bits(&pb, 4, adpcm_ima_compress_sample(&c->status[1], samples[2*i+1]));
641             }
642         }
643         flush_put_bits(&pb);
644         dst += put_bits_count(&pb)>>3;
645         break;
646     }
647     case CODEC_ID_ADPCM_MS:
648         for(i=0; i<avctx->channels; i++){
649             int predictor=0;
650
651             *dst++ = predictor;
652             c->status[i].coeff1 = AdaptCoeff1[predictor];
653             c->status[i].coeff2 = AdaptCoeff2[predictor];
654         }
655         for(i=0; i<avctx->channels; i++){
656             if (c->status[i].idelta < 16)
657                 c->status[i].idelta = 16;
658
659             bytestream_put_le16(&dst, c->status[i].idelta);
660         }
661         for(i=0; i<avctx->channels; i++){
662             c->status[i].sample2= *samples++;
663         }
664         for(i=0; i<avctx->channels; i++){
665             c->status[i].sample1= *samples++;
666
667             bytestream_put_le16(&dst, c->status[i].sample1);
668         }
669         for(i=0; i<avctx->channels; i++)
670             bytestream_put_le16(&dst, c->status[i].sample2);
671
672         if(avctx->trellis > 0) {
673             int n = avctx->block_align - 7*avctx->channels;
674             FF_ALLOC_OR_GOTO(avctx, buf, 2*n, error);
675             if(avctx->channels == 1) {
676                 adpcm_compress_trellis(avctx, samples, buf, &c->status[0], n);
677                 for(i=0; i<n; i+=2)
678                     *dst++ = (buf[i] << 4) | buf[i+1];
679             } else {
680                 adpcm_compress_trellis(avctx, samples, buf, &c->status[0], n);
681                 adpcm_compress_trellis(avctx, samples+1, buf+n, &c->status[1], n);
682                 for(i=0; i<n; i++)
683                     *dst++ = (buf[i] << 4) | buf[n+i];
684             }
685             av_free(buf);
686         } else
687         for(i=7*avctx->channels; i<avctx->block_align; i++) {
688             int nibble;
689             nibble = adpcm_ms_compress_sample(&c->status[ 0], *samples++)<<4;
690             nibble|= adpcm_ms_compress_sample(&c->status[st], *samples++);
691             *dst++ = nibble;
692         }
693         break;
694     case CODEC_ID_ADPCM_YAMAHA:
695         n = avctx->frame_size / 2;
696         if(avctx->trellis > 0) {
697             FF_ALLOC_OR_GOTO(avctx, buf, 2*n*2, error);
698             n *= 2;
699             if(avctx->channels == 1) {
700                 adpcm_compress_trellis(avctx, samples, buf, &c->status[0], n);
701                 for(i=0; i<n; i+=2)
702                     *dst++ = buf[i] | (buf[i+1] << 4);
703             } else {
704                 adpcm_compress_trellis(avctx, samples, buf, &c->status[0], n);
705                 adpcm_compress_trellis(avctx, samples+1, buf+n, &c->status[1], n);
706                 for(i=0; i<n; i++)
707                     *dst++ = buf[i] | (buf[n+i] << 4);
708             }
709             av_free(buf);
710         } else
711             for (n *= avctx->channels; n>0; n--) {
712                 int nibble;
713                 nibble  = adpcm_yamaha_compress_sample(&c->status[ 0], *samples++);
714                 nibble |= adpcm_yamaha_compress_sample(&c->status[st], *samples++) << 4;
715                 *dst++ = nibble;
716             }
717         break;
718     default:
719     error:
720         return -1;
721     }
722     return dst - frame;
723 }
724 #endif //CONFIG_ENCODERS
725
726 static av_cold int adpcm_decode_init(AVCodecContext * avctx)
727 {
728     ADPCMContext *c = avctx->priv_data;
729     unsigned int max_channels = 2;
730
731     switch(avctx->codec->id) {
732     case CODEC_ID_ADPCM_EA_R1:
733     case CODEC_ID_ADPCM_EA_R2:
734     case CODEC_ID_ADPCM_EA_R3:
735         max_channels = 6;
736         break;
737     }
738     if(avctx->channels > max_channels){
739         return -1;
740     }
741
742     switch(avctx->codec->id) {
743     case CODEC_ID_ADPCM_CT:
744         c->status[0].step = c->status[1].step = 511;
745         break;
746     case CODEC_ID_ADPCM_IMA_WAV:
747         if (avctx->bits_per_coded_sample != 4) {
748             av_log(avctx, AV_LOG_ERROR, "Only 4-bit ADPCM IMA WAV files are supported\n");
749             return -1;
750         }
751         break;
752     case CODEC_ID_ADPCM_IMA_WS:
753         if (avctx->extradata && avctx->extradata_size == 2 * 4) {
754             c->status[0].predictor = AV_RL32(avctx->extradata);
755             c->status[1].predictor = AV_RL32(avctx->extradata + 4);
756         }
757         break;
758     default:
759         break;
760     }
761     avctx->sample_fmt = AV_SAMPLE_FMT_S16;
762     return 0;
763 }
764
765 static inline short adpcm_ima_expand_nibble(ADPCMChannelStatus *c, char nibble, int shift)
766 {
767     int step_index;
768     int predictor;
769     int sign, delta, diff, step;
770
771     step = step_table[c->step_index];
772     step_index = c->step_index + index_table[(unsigned)nibble];
773     if (step_index < 0) step_index = 0;
774     else if (step_index > 88) step_index = 88;
775
776     sign = nibble & 8;
777     delta = nibble & 7;
778     /* perform direct multiplication instead of series of jumps proposed by
779      * the reference ADPCM implementation since modern CPUs can do the mults
780      * quickly enough */
781     diff = ((2 * delta + 1) * step) >> shift;
782     predictor = c->predictor;
783     if (sign) predictor -= diff;
784     else predictor += diff;
785
786     c->predictor = av_clip_int16(predictor);
787     c->step_index = step_index;
788
789     return (short)c->predictor;
790 }
791
792 static inline short adpcm_ms_expand_nibble(ADPCMChannelStatus *c, char nibble)
793 {
794     int predictor;
795
796     predictor = (((c->sample1) * (c->coeff1)) + ((c->sample2) * (c->coeff2))) / 64;
797     predictor += (signed)((nibble & 0x08)?(nibble - 0x10):(nibble)) * c->idelta;
798
799     c->sample2 = c->sample1;
800     c->sample1 = av_clip_int16(predictor);
801     c->idelta = (AdaptationTable[(int)nibble] * c->idelta) >> 8;
802     if (c->idelta < 16) c->idelta = 16;
803
804     return c->sample1;
805 }
806
807 static inline short adpcm_ct_expand_nibble(ADPCMChannelStatus *c, char nibble)
808 {
809     int sign, delta, diff;
810     int new_step;
811
812     sign = nibble & 8;
813     delta = nibble & 7;
814     /* perform direct multiplication instead of series of jumps proposed by
815      * the reference ADPCM implementation since modern CPUs can do the mults
816      * quickly enough */
817     diff = ((2 * delta + 1) * c->step) >> 3;
818     /* predictor update is not so trivial: predictor is multiplied on 254/256 before updating */
819     c->predictor = ((c->predictor * 254) >> 8) + (sign ? -diff : diff);
820     c->predictor = av_clip_int16(c->predictor);
821     /* calculate new step and clamp it to range 511..32767 */
822     new_step = (AdaptationTable[nibble & 7] * c->step) >> 8;
823     c->step = av_clip(new_step, 511, 32767);
824
825     return (short)c->predictor;
826 }
827
828 static inline short adpcm_sbpro_expand_nibble(ADPCMChannelStatus *c, char nibble, int size, int shift)
829 {
830     int sign, delta, diff;
831
832     sign = nibble & (1<<(size-1));
833     delta = nibble & ((1<<(size-1))-1);
834     diff = delta << (7 + c->step + shift);
835
836     /* clamp result */
837     c->predictor = av_clip(c->predictor + (sign ? -diff : diff), -16384,16256);
838
839     /* calculate new step */
840     if (delta >= (2*size - 3) && c->step < 3)
841         c->step++;
842     else if (delta == 0 && c->step > 0)
843         c->step--;
844
845     return (short) c->predictor;
846 }
847
848 static inline short adpcm_yamaha_expand_nibble(ADPCMChannelStatus *c, unsigned char nibble)
849 {
850     if(!c->step) {
851         c->predictor = 0;
852         c->step = 127;
853     }
854
855     c->predictor += (c->step * yamaha_difflookup[nibble]) / 8;
856     c->predictor = av_clip_int16(c->predictor);
857     c->step = (c->step * yamaha_indexscale[nibble]) >> 8;
858     c->step = av_clip(c->step, 127, 24567);
859     return c->predictor;
860 }
861
862 static void xa_decode(short *out, const unsigned char *in,
863     ADPCMChannelStatus *left, ADPCMChannelStatus *right, int inc)
864 {
865     int i, j;
866     int shift,filter,f0,f1;
867     int s_1,s_2;
868     int d,s,t;
869
870     for(i=0;i<4;i++) {
871
872         shift  = 12 - (in[4+i*2] & 15);
873         filter = in[4+i*2] >> 4;
874         f0 = xa_adpcm_table[filter][0];
875         f1 = xa_adpcm_table[filter][1];
876
877         s_1 = left->sample1;
878         s_2 = left->sample2;
879
880         for(j=0;j<28;j++) {
881             d = in[16+i+j*4];
882
883             t = (signed char)(d<<4)>>4;
884             s = ( t<<shift ) + ((s_1*f0 + s_2*f1+32)>>6);
885             s_2 = s_1;
886             s_1 = av_clip_int16(s);
887             *out = s_1;
888             out += inc;
889         }
890
891         if (inc==2) { /* stereo */
892             left->sample1 = s_1;
893             left->sample2 = s_2;
894             s_1 = right->sample1;
895             s_2 = right->sample2;
896             out = out + 1 - 28*2;
897         }
898
899         shift  = 12 - (in[5+i*2] & 15);
900         filter = in[5+i*2] >> 4;
901
902         f0 = xa_adpcm_table[filter][0];
903         f1 = xa_adpcm_table[filter][1];
904
905         for(j=0;j<28;j++) {
906             d = in[16+i+j*4];
907
908             t = (signed char)d >> 4;
909             s = ( t<<shift ) + ((s_1*f0 + s_2*f1+32)>>6);
910             s_2 = s_1;
911             s_1 = av_clip_int16(s);
912             *out = s_1;
913             out += inc;
914         }
915
916         if (inc==2) { /* stereo */
917             right->sample1 = s_1;
918             right->sample2 = s_2;
919             out -= 1;
920         } else {
921             left->sample1 = s_1;
922             left->sample2 = s_2;
923         }
924     }
925 }
926
927
928 /* DK3 ADPCM support macro */
929 #define DK3_GET_NEXT_NIBBLE() \
930     if (decode_top_nibble_next) \
931     { \
932         nibble = last_byte >> 4; \
933         decode_top_nibble_next = 0; \
934     } \
935     else \
936     { \
937         last_byte = *src++; \
938         if (src >= buf + buf_size) break; \
939         nibble = last_byte & 0x0F; \
940         decode_top_nibble_next = 1; \
941     }
942
943 static int adpcm_decode_frame(AVCodecContext *avctx,
944                             void *data, int *data_size,
945                             AVPacket *avpkt)
946 {
947     const uint8_t *buf = avpkt->data;
948     int buf_size = avpkt->size;
949     ADPCMContext *c = avctx->priv_data;
950     ADPCMChannelStatus *cs;
951     int n, m, channel, i;
952     int block_predictor[2];
953     short *samples;
954     short *samples_end;
955     const uint8_t *src;
956     int st; /* stereo */
957
958     /* DK3 ADPCM accounting variables */
959     unsigned char last_byte = 0;
960     unsigned char nibble;
961     int decode_top_nibble_next = 0;
962     int diff_channel;
963
964     /* EA ADPCM state variables */
965     uint32_t samples_in_chunk;
966     int32_t previous_left_sample, previous_right_sample;
967     int32_t current_left_sample, current_right_sample;
968     int32_t next_left_sample, next_right_sample;
969     int32_t coeff1l, coeff2l, coeff1r, coeff2r;
970     uint8_t shift_left, shift_right;
971     int count1, count2;
972     int coeff[2][2], shift[2];//used in EA MAXIS ADPCM
973
974     if (!buf_size)
975         return 0;
976
977     //should protect all 4bit ADPCM variants
978     //8 is needed for CODEC_ID_ADPCM_IMA_WAV with 2 channels
979     //
980     if(*data_size/4 < buf_size + 8)
981         return -1;
982
983     samples = data;
984     samples_end= samples + *data_size/2;
985     *data_size= 0;
986     src = buf;
987
988     st = avctx->channels == 2 ? 1 : 0;
989
990     switch(avctx->codec->id) {
991     case CODEC_ID_ADPCM_IMA_QT:
992         n = buf_size - 2*avctx->channels;
993         for (channel = 0; channel < avctx->channels; channel++) {
994             cs = &(c->status[channel]);
995             /* (pppppp) (piiiiiii) */
996
997             /* Bits 15-7 are the _top_ 9 bits of the 16-bit initial predictor value */
998             cs->predictor = (*src++) << 8;
999             cs->predictor |= (*src & 0x80);
1000             cs->predictor &= 0xFF80;
1001
1002             /* sign extension */
1003             if(cs->predictor & 0x8000)
1004                 cs->predictor -= 0x10000;
1005
1006             cs->predictor = av_clip_int16(cs->predictor);
1007
1008             cs->step_index = (*src++) & 0x7F;
1009
1010             if (cs->step_index > 88){
1011                 av_log(avctx, AV_LOG_ERROR, "ERROR: step_index = %i\n", cs->step_index);
1012                 cs->step_index = 88;
1013             }
1014
1015             cs->step = step_table[cs->step_index];
1016
1017             samples = (short*)data + channel;
1018
1019             for(m=32; n>0 && m>0; n--, m--) { /* in QuickTime, IMA is encoded by chuncks of 34 bytes (=64 samples) */
1020                 *samples = adpcm_ima_expand_nibble(cs, src[0] & 0x0F, 3);
1021                 samples += avctx->channels;
1022                 *samples = adpcm_ima_expand_nibble(cs, src[0] >> 4  , 3);
1023                 samples += avctx->channels;
1024                 src ++;
1025             }
1026         }
1027         if (st)
1028             samples--;
1029         break;
1030     case CODEC_ID_ADPCM_IMA_WAV:
1031         if (avctx->block_align != 0 && buf_size > avctx->block_align)
1032             buf_size = avctx->block_align;
1033
1034 //        samples_per_block= (block_align-4*chanels)*8 / (bits_per_sample * chanels) + 1;
1035
1036         for(i=0; i<avctx->channels; i++){
1037             cs = &(c->status[i]);
1038             cs->predictor = *samples++ = (int16_t)bytestream_get_le16(&src);
1039
1040             cs->step_index = *src++;
1041             if (cs->step_index > 88){
1042                 av_log(avctx, AV_LOG_ERROR, "ERROR: step_index = %i\n", cs->step_index);
1043                 cs->step_index = 88;
1044             }
1045             if (*src++) av_log(avctx, AV_LOG_ERROR, "unused byte should be null but is %d!!\n", src[-1]); /* unused */
1046         }
1047
1048         while(src < buf + buf_size){
1049             for(m=0; m<4; m++){
1050                 for(i=0; i<=st; i++)
1051                     *samples++ = adpcm_ima_expand_nibble(&c->status[i], src[4*i] & 0x0F, 3);
1052                 for(i=0; i<=st; i++)
1053                     *samples++ = adpcm_ima_expand_nibble(&c->status[i], src[4*i] >> 4  , 3);
1054                 src++;
1055             }
1056             src += 4*st;
1057         }
1058         break;
1059     case CODEC_ID_ADPCM_4XM:
1060         cs = &(c->status[0]);
1061         c->status[0].predictor= (int16_t)bytestream_get_le16(&src);
1062         if(st){
1063             c->status[1].predictor= (int16_t)bytestream_get_le16(&src);
1064         }
1065         c->status[0].step_index= (int16_t)bytestream_get_le16(&src);
1066         if(st){
1067             c->status[1].step_index= (int16_t)bytestream_get_le16(&src);
1068         }
1069         if (cs->step_index < 0) cs->step_index = 0;
1070         if (cs->step_index > 88) cs->step_index = 88;
1071
1072         m= (buf_size - (src - buf))>>st;
1073         for(i=0; i<m; i++) {
1074             *samples++ = adpcm_ima_expand_nibble(&c->status[0], src[i] & 0x0F, 4);
1075             if (st)
1076                 *samples++ = adpcm_ima_expand_nibble(&c->status[1], src[i+m] & 0x0F, 4);
1077             *samples++ = adpcm_ima_expand_nibble(&c->status[0], src[i] >> 4, 4);
1078             if (st)
1079                 *samples++ = adpcm_ima_expand_nibble(&c->status[1], src[i+m] >> 4, 4);
1080         }
1081
1082         src += m<<st;
1083
1084         break;
1085     case CODEC_ID_ADPCM_MS:
1086         if (avctx->block_align != 0 && buf_size > avctx->block_align)
1087             buf_size = avctx->block_align;
1088         n = buf_size - 7 * avctx->channels;
1089         if (n < 0)
1090             return -1;
1091         block_predictor[0] = av_clip(*src++, 0, 6);
1092         block_predictor[1] = 0;
1093         if (st)
1094             block_predictor[1] = av_clip(*src++, 0, 6);
1095         c->status[0].idelta = (int16_t)bytestream_get_le16(&src);
1096         if (st){
1097             c->status[1].idelta = (int16_t)bytestream_get_le16(&src);
1098         }
1099         c->status[0].coeff1 = AdaptCoeff1[block_predictor[0]];
1100         c->status[0].coeff2 = AdaptCoeff2[block_predictor[0]];
1101         c->status[1].coeff1 = AdaptCoeff1[block_predictor[1]];
1102         c->status[1].coeff2 = AdaptCoeff2[block_predictor[1]];
1103
1104         c->status[0].sample1 = bytestream_get_le16(&src);
1105         if (st) c->status[1].sample1 = bytestream_get_le16(&src);
1106         c->status[0].sample2 = bytestream_get_le16(&src);
1107         if (st) c->status[1].sample2 = bytestream_get_le16(&src);
1108
1109         *samples++ = c->status[0].sample2;
1110         if (st) *samples++ = c->status[1].sample2;
1111         *samples++ = c->status[0].sample1;
1112         if (st) *samples++ = c->status[1].sample1;
1113         for(;n>0;n--) {
1114             *samples++ = adpcm_ms_expand_nibble(&c->status[0 ], src[0] >> 4  );
1115             *samples++ = adpcm_ms_expand_nibble(&c->status[st], src[0] & 0x0F);
1116             src ++;
1117         }
1118         break;
1119     case CODEC_ID_ADPCM_IMA_DK4:
1120         if (avctx->block_align != 0 && buf_size > avctx->block_align)
1121             buf_size = avctx->block_align;
1122
1123         c->status[0].predictor  = (int16_t)bytestream_get_le16(&src);
1124         c->status[0].step_index = *src++;
1125         src++;
1126         *samples++ = c->status[0].predictor;
1127         if (st) {
1128             c->status[1].predictor  = (int16_t)bytestream_get_le16(&src);
1129             c->status[1].step_index = *src++;
1130             src++;
1131             *samples++ = c->status[1].predictor;
1132         }
1133         while (src < buf + buf_size) {
1134
1135             /* take care of the top nibble (always left or mono channel) */
1136             *samples++ = adpcm_ima_expand_nibble(&c->status[0],
1137                 src[0] >> 4, 3);
1138
1139             /* take care of the bottom nibble, which is right sample for
1140              * stereo, or another mono sample */
1141             if (st)
1142                 *samples++ = adpcm_ima_expand_nibble(&c->status[1],
1143                     src[0] & 0x0F, 3);
1144             else
1145                 *samples++ = adpcm_ima_expand_nibble(&c->status[0],
1146                     src[0] & 0x0F, 3);
1147
1148             src++;
1149         }
1150         break;
1151     case CODEC_ID_ADPCM_IMA_DK3:
1152         if (avctx->block_align != 0 && buf_size > avctx->block_align)
1153             buf_size = avctx->block_align;
1154
1155         if(buf_size + 16 > (samples_end - samples)*3/8)
1156             return -1;
1157
1158         c->status[0].predictor  = (int16_t)AV_RL16(src + 10);
1159         c->status[1].predictor  = (int16_t)AV_RL16(src + 12);
1160         c->status[0].step_index = src[14];
1161         c->status[1].step_index = src[15];
1162         /* sign extend the predictors */
1163         src += 16;
1164         diff_channel = c->status[1].predictor;
1165
1166         /* the DK3_GET_NEXT_NIBBLE macro issues the break statement when
1167          * the buffer is consumed */
1168         while (1) {
1169
1170             /* for this algorithm, c->status[0] is the sum channel and
1171              * c->status[1] is the diff channel */
1172
1173             /* process the first predictor of the sum channel */
1174             DK3_GET_NEXT_NIBBLE();
1175             adpcm_ima_expand_nibble(&c->status[0], nibble, 3);
1176
1177             /* process the diff channel predictor */
1178             DK3_GET_NEXT_NIBBLE();
1179             adpcm_ima_expand_nibble(&c->status[1], nibble, 3);
1180
1181             /* process the first pair of stereo PCM samples */
1182             diff_channel = (diff_channel + c->status[1].predictor) / 2;
1183             *samples++ = c->status[0].predictor + c->status[1].predictor;
1184             *samples++ = c->status[0].predictor - c->status[1].predictor;
1185
1186             /* process the second predictor of the sum channel */
1187             DK3_GET_NEXT_NIBBLE();
1188             adpcm_ima_expand_nibble(&c->status[0], nibble, 3);
1189
1190             /* process the second pair of stereo PCM samples */
1191             diff_channel = (diff_channel + c->status[1].predictor) / 2;
1192             *samples++ = c->status[0].predictor + c->status[1].predictor;
1193             *samples++ = c->status[0].predictor - c->status[1].predictor;
1194         }
1195         break;
1196     case CODEC_ID_ADPCM_IMA_ISS:
1197         c->status[0].predictor  = (int16_t)AV_RL16(src + 0);
1198         c->status[0].step_index = src[2];
1199         src += 4;
1200         if(st) {
1201             c->status[1].predictor  = (int16_t)AV_RL16(src + 0);
1202             c->status[1].step_index = src[2];
1203             src += 4;
1204         }
1205
1206         while (src < buf + buf_size) {
1207
1208             if (st) {
1209                 *samples++ = adpcm_ima_expand_nibble(&c->status[0],
1210                     src[0] >> 4  , 3);
1211                 *samples++ = adpcm_ima_expand_nibble(&c->status[1],
1212                     src[0] & 0x0F, 3);
1213             } else {
1214                 *samples++ = adpcm_ima_expand_nibble(&c->status[0],
1215                     src[0] & 0x0F, 3);
1216                 *samples++ = adpcm_ima_expand_nibble(&c->status[0],
1217                     src[0] >> 4  , 3);
1218             }
1219
1220             src++;
1221         }
1222         break;
1223     case CODEC_ID_ADPCM_IMA_WS:
1224         /* no per-block initialization; just start decoding the data */
1225         while (src < buf + buf_size) {
1226
1227             if (st) {
1228                 *samples++ = adpcm_ima_expand_nibble(&c->status[0],
1229                     src[0] >> 4  , 3);
1230                 *samples++ = adpcm_ima_expand_nibble(&c->status[1],
1231                     src[0] & 0x0F, 3);
1232             } else {
1233                 *samples++ = adpcm_ima_expand_nibble(&c->status[0],
1234                     src[0] >> 4  , 3);
1235                 *samples++ = adpcm_ima_expand_nibble(&c->status[0],
1236                     src[0] & 0x0F, 3);
1237             }
1238
1239             src++;
1240         }
1241         break;
1242     case CODEC_ID_ADPCM_XA:
1243         while (buf_size >= 128) {
1244             xa_decode(samples, src, &c->status[0], &c->status[1],
1245                 avctx->channels);
1246             src += 128;
1247             samples += 28 * 8;
1248             buf_size -= 128;
1249         }
1250         break;
1251     case CODEC_ID_ADPCM_IMA_EA_EACS:
1252         samples_in_chunk = bytestream_get_le32(&src) >> (1-st);
1253
1254         if (samples_in_chunk > buf_size-4-(8<<st)) {
1255             src += buf_size - 4;
1256             break;
1257         }
1258
1259         for (i=0; i<=st; i++)
1260             c->status[i].step_index = bytestream_get_le32(&src);
1261         for (i=0; i<=st; i++)
1262             c->status[i].predictor  = bytestream_get_le32(&src);
1263
1264         for (; samples_in_chunk; samples_in_chunk--, src++) {
1265             *samples++ = adpcm_ima_expand_nibble(&c->status[0],  *src>>4,   3);
1266             *samples++ = adpcm_ima_expand_nibble(&c->status[st], *src&0x0F, 3);
1267         }
1268         break;
1269     case CODEC_ID_ADPCM_IMA_EA_SEAD:
1270         for (; src < buf+buf_size; src++) {
1271             *samples++ = adpcm_ima_expand_nibble(&c->status[0], src[0] >> 4, 6);
1272             *samples++ = adpcm_ima_expand_nibble(&c->status[st],src[0]&0x0F, 6);
1273         }
1274         break;
1275     case CODEC_ID_ADPCM_EA:
1276         if (buf_size < 4 || AV_RL32(src) >= ((buf_size - 12) * 2)) {
1277             src += buf_size;
1278             break;
1279         }
1280         samples_in_chunk = AV_RL32(src);
1281         src += 4;
1282         current_left_sample   = (int16_t)bytestream_get_le16(&src);
1283         previous_left_sample  = (int16_t)bytestream_get_le16(&src);
1284         current_right_sample  = (int16_t)bytestream_get_le16(&src);
1285         previous_right_sample = (int16_t)bytestream_get_le16(&src);
1286
1287         for (count1 = 0; count1 < samples_in_chunk/28;count1++) {
1288             coeff1l = ea_adpcm_table[ *src >> 4       ];
1289             coeff2l = ea_adpcm_table[(*src >> 4  ) + 4];
1290             coeff1r = ea_adpcm_table[*src & 0x0F];
1291             coeff2r = ea_adpcm_table[(*src & 0x0F) + 4];
1292             src++;
1293
1294             shift_left  = (*src >> 4  ) + 8;
1295             shift_right = (*src & 0x0F) + 8;
1296             src++;
1297
1298             for (count2 = 0; count2 < 28; count2++) {
1299                 next_left_sample  = (int32_t)((*src & 0xF0) << 24) >> shift_left;
1300                 next_right_sample = (int32_t)((*src & 0x0F) << 28) >> shift_right;
1301                 src++;
1302
1303                 next_left_sample = (next_left_sample +
1304                     (current_left_sample * coeff1l) +
1305                     (previous_left_sample * coeff2l) + 0x80) >> 8;
1306                 next_right_sample = (next_right_sample +
1307                     (current_right_sample * coeff1r) +
1308                     (previous_right_sample * coeff2r) + 0x80) >> 8;
1309
1310                 previous_left_sample = current_left_sample;
1311                 current_left_sample = av_clip_int16(next_left_sample);
1312                 previous_right_sample = current_right_sample;
1313                 current_right_sample = av_clip_int16(next_right_sample);
1314                 *samples++ = (unsigned short)current_left_sample;
1315                 *samples++ = (unsigned short)current_right_sample;
1316             }
1317         }
1318
1319         if (src - buf == buf_size - 2)
1320             src += 2; // Skip terminating 0x0000
1321
1322         break;
1323     case CODEC_ID_ADPCM_EA_MAXIS_XA:
1324         for(channel = 0; channel < avctx->channels; channel++) {
1325             for (i=0; i<2; i++)
1326                 coeff[channel][i] = ea_adpcm_table[(*src >> 4) + 4*i];
1327             shift[channel] = (*src & 0x0F) + 8;
1328             src++;
1329         }
1330         for (count1 = 0; count1 < (buf_size - avctx->channels) / avctx->channels; count1++) {
1331             for(i = 4; i >= 0; i-=4) { /* Pairwise samples LL RR (st) or LL LL (mono) */
1332                 for(channel = 0; channel < avctx->channels; channel++) {
1333                     int32_t sample = (int32_t)(((*(src+channel) >> i) & 0x0F) << 0x1C) >> shift[channel];
1334                     sample = (sample +
1335                              c->status[channel].sample1 * coeff[channel][0] +
1336                              c->status[channel].sample2 * coeff[channel][1] + 0x80) >> 8;
1337                     c->status[channel].sample2 = c->status[channel].sample1;
1338                     c->status[channel].sample1 = av_clip_int16(sample);
1339                     *samples++ = c->status[channel].sample1;
1340                 }
1341             }
1342             src+=avctx->channels;
1343         }
1344         break;
1345     case CODEC_ID_ADPCM_EA_R1:
1346     case CODEC_ID_ADPCM_EA_R2:
1347     case CODEC_ID_ADPCM_EA_R3: {
1348         /* channel numbering
1349            2chan: 0=fl, 1=fr
1350            4chan: 0=fl, 1=rl, 2=fr, 3=rr
1351            6chan: 0=fl, 1=c,  2=fr, 3=rl,  4=rr, 5=sub */
1352         const int big_endian = avctx->codec->id == CODEC_ID_ADPCM_EA_R3;
1353         int32_t previous_sample, current_sample, next_sample;
1354         int32_t coeff1, coeff2;
1355         uint8_t shift;
1356         unsigned int channel;
1357         uint16_t *samplesC;
1358         const uint8_t *srcC;
1359         const uint8_t *src_end = buf + buf_size;
1360
1361         samples_in_chunk = (big_endian ? bytestream_get_be32(&src)
1362                                        : bytestream_get_le32(&src)) / 28;
1363         if (samples_in_chunk > UINT32_MAX/(28*avctx->channels) ||
1364             28*samples_in_chunk*avctx->channels > samples_end-samples) {
1365             src += buf_size - 4;
1366             break;
1367         }
1368
1369         for (channel=0; channel<avctx->channels; channel++) {
1370             int32_t offset = (big_endian ? bytestream_get_be32(&src)
1371                                          : bytestream_get_le32(&src))
1372                            + (avctx->channels-channel-1) * 4;
1373
1374             if ((offset < 0) || (offset >= src_end - src - 4)) break;
1375             srcC  = src + offset;
1376             samplesC = samples + channel;
1377
1378             if (avctx->codec->id == CODEC_ID_ADPCM_EA_R1) {
1379                 current_sample  = (int16_t)bytestream_get_le16(&srcC);
1380                 previous_sample = (int16_t)bytestream_get_le16(&srcC);
1381             } else {
1382                 current_sample  = c->status[channel].predictor;
1383                 previous_sample = c->status[channel].prev_sample;
1384             }
1385
1386             for (count1=0; count1<samples_in_chunk; count1++) {
1387                 if (*srcC == 0xEE) {  /* only seen in R2 and R3 */
1388                     srcC++;
1389                     if (srcC > src_end - 30*2) break;
1390                     current_sample  = (int16_t)bytestream_get_be16(&srcC);
1391                     previous_sample = (int16_t)bytestream_get_be16(&srcC);
1392
1393                     for (count2=0; count2<28; count2++) {
1394                         *samplesC = (int16_t)bytestream_get_be16(&srcC);
1395                         samplesC += avctx->channels;
1396                     }
1397                 } else {
1398                     coeff1 = ea_adpcm_table[ *srcC>>4     ];
1399                     coeff2 = ea_adpcm_table[(*srcC>>4) + 4];
1400                     shift = (*srcC++ & 0x0F) + 8;
1401
1402                     if (srcC > src_end - 14) break;
1403                     for (count2=0; count2<28; count2++) {
1404                         if (count2 & 1)
1405                             next_sample = (int32_t)((*srcC++ & 0x0F) << 28) >> shift;
1406                         else
1407                             next_sample = (int32_t)((*srcC   & 0xF0) << 24) >> shift;
1408
1409                         next_sample += (current_sample  * coeff1) +
1410                                        (previous_sample * coeff2);
1411                         next_sample = av_clip_int16(next_sample >> 8);
1412
1413                         previous_sample = current_sample;
1414                         current_sample  = next_sample;
1415                         *samplesC = current_sample;
1416                         samplesC += avctx->channels;
1417                     }
1418                 }
1419             }
1420
1421             if (avctx->codec->id != CODEC_ID_ADPCM_EA_R1) {
1422                 c->status[channel].predictor   = current_sample;
1423                 c->status[channel].prev_sample = previous_sample;
1424             }
1425         }
1426
1427         src = src + buf_size - (4 + 4*avctx->channels);
1428         samples += 28 * samples_in_chunk * avctx->channels;
1429         break;
1430     }
1431     case CODEC_ID_ADPCM_EA_XAS:
1432         if (samples_end-samples < 32*4*avctx->channels
1433             || buf_size < (4+15)*4*avctx->channels) {
1434             src += buf_size;
1435             break;
1436         }
1437         for (channel=0; channel<avctx->channels; channel++) {
1438             int coeff[2][4], shift[4];
1439             short *s2, *s = &samples[channel];
1440             for (n=0; n<4; n++, s+=32*avctx->channels) {
1441                 for (i=0; i<2; i++)
1442                     coeff[i][n] = ea_adpcm_table[(src[0]&0x0F)+4*i];
1443                 shift[n] = (src[2]&0x0F) + 8;
1444                 for (s2=s, i=0; i<2; i++, src+=2, s2+=avctx->channels)
1445                     s2[0] = (src[0]&0xF0) + (src[1]<<8);
1446             }
1447
1448             for (m=2; m<32; m+=2) {
1449                 s = &samples[m*avctx->channels + channel];
1450                 for (n=0; n<4; n++, src++, s+=32*avctx->channels) {
1451                     for (s2=s, i=0; i<8; i+=4, s2+=avctx->channels) {
1452                         int level = (int32_t)((*src & (0xF0>>i)) << (24+i)) >> shift[n];
1453                         int pred  = s2[-1*avctx->channels] * coeff[0][n]
1454                                   + s2[-2*avctx->channels] * coeff[1][n];
1455                         s2[0] = av_clip_int16((level + pred + 0x80) >> 8);
1456                     }
1457                 }
1458             }
1459         }
1460         samples += 32*4*avctx->channels;
1461         break;
1462     case CODEC_ID_ADPCM_IMA_AMV:
1463     case CODEC_ID_ADPCM_IMA_SMJPEG:
1464         c->status[0].predictor = (int16_t)bytestream_get_le16(&src);
1465         c->status[0].step_index = bytestream_get_le16(&src);
1466
1467         if (avctx->codec->id == CODEC_ID_ADPCM_IMA_AMV)
1468             src+=4;
1469
1470         while (src < buf + buf_size) {
1471             char hi, lo;
1472             lo = *src & 0x0F;
1473             hi = *src >> 4;
1474
1475             if (avctx->codec->id == CODEC_ID_ADPCM_IMA_AMV)
1476                 FFSWAP(char, hi, lo);
1477
1478             *samples++ = adpcm_ima_expand_nibble(&c->status[0],
1479                 lo, 3);
1480             *samples++ = adpcm_ima_expand_nibble(&c->status[0],
1481                 hi, 3);
1482             src++;
1483         }
1484         break;
1485     case CODEC_ID_ADPCM_CT:
1486         while (src < buf + buf_size) {
1487             if (st) {
1488                 *samples++ = adpcm_ct_expand_nibble(&c->status[0],
1489                     src[0] >> 4);
1490                 *samples++ = adpcm_ct_expand_nibble(&c->status[1],
1491                     src[0] & 0x0F);
1492             } else {
1493                 *samples++ = adpcm_ct_expand_nibble(&c->status[0],
1494                     src[0] >> 4);
1495                 *samples++ = adpcm_ct_expand_nibble(&c->status[0],
1496                     src[0] & 0x0F);
1497             }
1498             src++;
1499         }
1500         break;
1501     case CODEC_ID_ADPCM_SBPRO_4:
1502     case CODEC_ID_ADPCM_SBPRO_3:
1503     case CODEC_ID_ADPCM_SBPRO_2:
1504         if (!c->status[0].step_index) {
1505             /* the first byte is a raw sample */
1506             *samples++ = 128 * (*src++ - 0x80);
1507             if (st)
1508               *samples++ = 128 * (*src++ - 0x80);
1509             c->status[0].step_index = 1;
1510         }
1511         if (avctx->codec->id == CODEC_ID_ADPCM_SBPRO_4) {
1512             while (src < buf + buf_size) {
1513                 *samples++ = adpcm_sbpro_expand_nibble(&c->status[0],
1514                     src[0] >> 4, 4, 0);
1515                 *samples++ = adpcm_sbpro_expand_nibble(&c->status[st],
1516                     src[0] & 0x0F, 4, 0);
1517                 src++;
1518             }
1519         } else if (avctx->codec->id == CODEC_ID_ADPCM_SBPRO_3) {
1520             while (src < buf + buf_size && samples + 2 < samples_end) {
1521                 *samples++ = adpcm_sbpro_expand_nibble(&c->status[0],
1522                      src[0] >> 5        , 3, 0);
1523                 *samples++ = adpcm_sbpro_expand_nibble(&c->status[0],
1524                     (src[0] >> 2) & 0x07, 3, 0);
1525                 *samples++ = adpcm_sbpro_expand_nibble(&c->status[0],
1526                     src[0] & 0x03, 2, 0);
1527                 src++;
1528             }
1529         } else {
1530             while (src < buf + buf_size && samples + 3 < samples_end) {
1531                 *samples++ = adpcm_sbpro_expand_nibble(&c->status[0],
1532                      src[0] >> 6        , 2, 2);
1533                 *samples++ = adpcm_sbpro_expand_nibble(&c->status[st],
1534                     (src[0] >> 4) & 0x03, 2, 2);
1535                 *samples++ = adpcm_sbpro_expand_nibble(&c->status[0],
1536                     (src[0] >> 2) & 0x03, 2, 2);
1537                 *samples++ = adpcm_sbpro_expand_nibble(&c->status[st],
1538                     src[0] & 0x03, 2, 2);
1539                 src++;
1540             }
1541         }
1542         break;
1543     case CODEC_ID_ADPCM_SWF:
1544     {
1545         GetBitContext gb;
1546         const int *table;
1547         int k0, signmask, nb_bits, count;
1548         int size = buf_size*8;
1549
1550         init_get_bits(&gb, buf, size);
1551
1552         //read bits & initial values
1553         nb_bits = get_bits(&gb, 2)+2;
1554         //av_log(NULL,AV_LOG_INFO,"nb_bits: %d\n", nb_bits);
1555         table = swf_index_tables[nb_bits-2];
1556         k0 = 1 << (nb_bits-2);
1557         signmask = 1 << (nb_bits-1);
1558
1559         while (get_bits_count(&gb) <= size - 22*avctx->channels) {
1560             for (i = 0; i < avctx->channels; i++) {
1561                 *samples++ = c->status[i].predictor = get_sbits(&gb, 16);
1562                 c->status[i].step_index = get_bits(&gb, 6);
1563             }
1564
1565             for (count = 0; get_bits_count(&gb) <= size - nb_bits*avctx->channels && count < 4095; count++) {
1566                 int i;
1567
1568                 for (i = 0; i < avctx->channels; i++) {
1569                     // similar to IMA adpcm
1570                     int delta = get_bits(&gb, nb_bits);
1571                     int step = step_table[c->status[i].step_index];
1572                     long vpdiff = 0; // vpdiff = (delta+0.5)*step/4
1573                     int k = k0;
1574
1575                     do {
1576                         if (delta & k)
1577                             vpdiff += step;
1578                         step >>= 1;
1579                         k >>= 1;
1580                     } while(k);
1581                     vpdiff += step;
1582
1583                     if (delta & signmask)
1584                         c->status[i].predictor -= vpdiff;
1585                     else
1586                         c->status[i].predictor += vpdiff;
1587
1588                     c->status[i].step_index += table[delta & (~signmask)];
1589
1590                     c->status[i].step_index = av_clip(c->status[i].step_index, 0, 88);
1591                     c->status[i].predictor = av_clip_int16(c->status[i].predictor);
1592
1593                     *samples++ = c->status[i].predictor;
1594                     if (samples >= samples_end) {
1595                         av_log(avctx, AV_LOG_ERROR, "allocated output buffer is too small\n");
1596                         return -1;
1597                     }
1598                 }
1599             }
1600         }
1601         src += buf_size;
1602         break;
1603     }
1604     case CODEC_ID_ADPCM_YAMAHA:
1605         while (src < buf + buf_size) {
1606             if (st) {
1607                 *samples++ = adpcm_yamaha_expand_nibble(&c->status[0],
1608                         src[0] & 0x0F);
1609                 *samples++ = adpcm_yamaha_expand_nibble(&c->status[1],
1610                         src[0] >> 4  );
1611             } else {
1612                 *samples++ = adpcm_yamaha_expand_nibble(&c->status[0],
1613                         src[0] & 0x0F);
1614                 *samples++ = adpcm_yamaha_expand_nibble(&c->status[0],
1615                         src[0] >> 4  );
1616             }
1617             src++;
1618         }
1619         break;
1620     case CODEC_ID_ADPCM_THP:
1621     {
1622         int table[2][16];
1623         unsigned int samplecnt;
1624         int prev[2][2];
1625         int ch;
1626
1627         if (buf_size < 80) {
1628             av_log(avctx, AV_LOG_ERROR, "frame too small\n");
1629             return -1;
1630         }
1631
1632         src+=4;
1633         samplecnt = bytestream_get_be32(&src);
1634
1635         for (i = 0; i < 32; i++)
1636             table[0][i] = (int16_t)bytestream_get_be16(&src);
1637
1638         /* Initialize the previous sample.  */
1639         for (i = 0; i < 4; i++)
1640             prev[0][i] = (int16_t)bytestream_get_be16(&src);
1641
1642         if (samplecnt >= (samples_end - samples) /  (st + 1)) {
1643             av_log(avctx, AV_LOG_ERROR, "allocated output buffer is too small\n");
1644             return -1;
1645         }
1646
1647         for (ch = 0; ch <= st; ch++) {
1648             samples = (unsigned short *) data + ch;
1649
1650             /* Read in every sample for this channel.  */
1651             for (i = 0; i < samplecnt / 14; i++) {
1652                 int index = (*src >> 4) & 7;
1653                 unsigned int exp = 28 - (*src++ & 15);
1654                 int factor1 = table[ch][index * 2];
1655                 int factor2 = table[ch][index * 2 + 1];
1656
1657                 /* Decode 14 samples.  */
1658                 for (n = 0; n < 14; n++) {
1659                     int32_t sampledat;
1660                     if(n&1) sampledat=  *src++    <<28;
1661                     else    sampledat= (*src&0xF0)<<24;
1662
1663                     sampledat = ((prev[ch][0]*factor1
1664                                 + prev[ch][1]*factor2) >> 11) + (sampledat>>exp);
1665                     *samples = av_clip_int16(sampledat);
1666                     prev[ch][1] = prev[ch][0];
1667                     prev[ch][0] = *samples++;
1668
1669                     /* In case of stereo, skip one sample, this sample
1670                        is for the other channel.  */
1671                     samples += st;
1672                 }
1673             }
1674         }
1675
1676         /* In the previous loop, in case stereo is used, samples is
1677            increased exactly one time too often.  */
1678         samples -= st;
1679         break;
1680     }
1681
1682     default:
1683         return -1;
1684     }
1685     *data_size = (uint8_t *)samples - (uint8_t *)data;
1686     return src - buf;
1687 }
1688
1689
1690
1691 #if CONFIG_ENCODERS
1692 #define ADPCM_ENCODER(id,name,long_name_)       \
1693 AVCodec name ## _encoder = {                    \
1694     #name,                                      \
1695     AVMEDIA_TYPE_AUDIO,                         \
1696     id,                                         \
1697     sizeof(ADPCMContext),                       \
1698     adpcm_encode_init,                          \
1699     adpcm_encode_frame,                         \
1700     adpcm_encode_close,                         \
1701     NULL,                                       \
1702     .sample_fmts = (const enum AVSampleFormat[]){AV_SAMPLE_FMT_S16,AV_SAMPLE_FMT_NONE}, \
1703     .long_name = NULL_IF_CONFIG_SMALL(long_name_), \
1704 };
1705 #else
1706 #define ADPCM_ENCODER(id,name,long_name_)
1707 #endif
1708
1709 #if CONFIG_DECODERS
1710 #define ADPCM_DECODER(id,name,long_name_)       \
1711 AVCodec name ## _decoder = {                    \
1712     #name,                                      \
1713     AVMEDIA_TYPE_AUDIO,                         \
1714     id,                                         \
1715     sizeof(ADPCMContext),                       \
1716     adpcm_decode_init,                          \
1717     NULL,                                       \
1718     NULL,                                       \
1719     adpcm_decode_frame,                         \
1720     .long_name = NULL_IF_CONFIG_SMALL(long_name_), \
1721 };
1722 #else
1723 #define ADPCM_DECODER(id,name,long_name_)
1724 #endif
1725
1726 #define ADPCM_CODEC(id,name,long_name_)         \
1727     ADPCM_ENCODER(id,name,long_name_) ADPCM_DECODER(id,name,long_name_)
1728
1729 /* Note: Do not forget to add new entries to the Makefile as well. */
1730 ADPCM_DECODER(CODEC_ID_ADPCM_4XM, adpcm_4xm, "ADPCM 4X Movie");
1731 ADPCM_DECODER(CODEC_ID_ADPCM_CT, adpcm_ct, "ADPCM Creative Technology");
1732 ADPCM_DECODER(CODEC_ID_ADPCM_EA, adpcm_ea, "ADPCM Electronic Arts");
1733 ADPCM_DECODER(CODEC_ID_ADPCM_EA_MAXIS_XA, adpcm_ea_maxis_xa, "ADPCM Electronic Arts Maxis CDROM XA");
1734 ADPCM_DECODER(CODEC_ID_ADPCM_EA_R1, adpcm_ea_r1, "ADPCM Electronic Arts R1");
1735 ADPCM_DECODER(CODEC_ID_ADPCM_EA_R2, adpcm_ea_r2, "ADPCM Electronic Arts R2");
1736 ADPCM_DECODER(CODEC_ID_ADPCM_EA_R3, adpcm_ea_r3, "ADPCM Electronic Arts R3");
1737 ADPCM_DECODER(CODEC_ID_ADPCM_EA_XAS, adpcm_ea_xas, "ADPCM Electronic Arts XAS");
1738 ADPCM_DECODER(CODEC_ID_ADPCM_IMA_AMV, adpcm_ima_amv, "ADPCM IMA AMV");
1739 ADPCM_DECODER(CODEC_ID_ADPCM_IMA_DK3, adpcm_ima_dk3, "ADPCM IMA Duck DK3");
1740 ADPCM_DECODER(CODEC_ID_ADPCM_IMA_DK4, adpcm_ima_dk4, "ADPCM IMA Duck DK4");
1741 ADPCM_DECODER(CODEC_ID_ADPCM_IMA_EA_EACS, adpcm_ima_ea_eacs, "ADPCM IMA Electronic Arts EACS");
1742 ADPCM_DECODER(CODEC_ID_ADPCM_IMA_EA_SEAD, adpcm_ima_ea_sead, "ADPCM IMA Electronic Arts SEAD");
1743 ADPCM_DECODER(CODEC_ID_ADPCM_IMA_ISS, adpcm_ima_iss, "ADPCM IMA Funcom ISS");
1744 ADPCM_CODEC  (CODEC_ID_ADPCM_IMA_QT, adpcm_ima_qt, "ADPCM IMA QuickTime");
1745 ADPCM_DECODER(CODEC_ID_ADPCM_IMA_SMJPEG, adpcm_ima_smjpeg, "ADPCM IMA Loki SDL MJPEG");
1746 ADPCM_CODEC  (CODEC_ID_ADPCM_IMA_WAV, adpcm_ima_wav, "ADPCM IMA WAV");
1747 ADPCM_DECODER(CODEC_ID_ADPCM_IMA_WS, adpcm_ima_ws, "ADPCM IMA Westwood");
1748 ADPCM_CODEC  (CODEC_ID_ADPCM_MS, adpcm_ms, "ADPCM Microsoft");
1749 ADPCM_DECODER(CODEC_ID_ADPCM_SBPRO_2, adpcm_sbpro_2, "ADPCM Sound Blaster Pro 2-bit");
1750 ADPCM_DECODER(CODEC_ID_ADPCM_SBPRO_3, adpcm_sbpro_3, "ADPCM Sound Blaster Pro 2.6-bit");
1751 ADPCM_DECODER(CODEC_ID_ADPCM_SBPRO_4, adpcm_sbpro_4, "ADPCM Sound Blaster Pro 4-bit");
1752 ADPCM_CODEC  (CODEC_ID_ADPCM_SWF, adpcm_swf, "ADPCM Shockwave Flash");
1753 ADPCM_DECODER(CODEC_ID_ADPCM_THP, adpcm_thp, "ADPCM Nintendo Gamecube THP");
1754 ADPCM_DECODER(CODEC_ID_ADPCM_XA, adpcm_xa, "ADPCM CDROM XA");
1755 ADPCM_CODEC  (CODEC_ID_ADPCM_YAMAHA, adpcm_yamaha, "ADPCM Yamaha");