3 * Copyright (c) 2003 Michael Niedermayer
5 * This file is part of Libav.
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.
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.
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
28 #include "libavutil/common.h"
31 #include "mpeg12data.h"
37 #define ASV2_LEVEL_VLC_BITS 10
39 typedef struct ASV1Context{
40 AVCodecContext *avctx;
51 DECLARE_ALIGNED(16, DCTELEM, block)[6][64];
52 uint16_t intra_matrix[64];
53 int q_intra_matrix[64];
54 uint8_t *bitstream_buffer;
55 unsigned int bitstream_buffer_size;
58 static const uint8_t scantab[64]={
59 0x00,0x08,0x01,0x09,0x10,0x18,0x11,0x19,
60 0x02,0x0A,0x03,0x0B,0x12,0x1A,0x13,0x1B,
61 0x04,0x0C,0x05,0x0D,0x20,0x28,0x21,0x29,
62 0x06,0x0E,0x07,0x0F,0x14,0x1C,0x15,0x1D,
63 0x22,0x2A,0x23,0x2B,0x30,0x38,0x31,0x39,
64 0x16,0x1E,0x17,0x1F,0x24,0x2C,0x25,0x2D,
65 0x32,0x3A,0x33,0x3B,0x26,0x2E,0x27,0x2F,
66 0x34,0x3C,0x35,0x3D,0x36,0x3E,0x37,0x3F,
70 static const uint8_t ccp_tab[17][2]={
71 {0x2,2}, {0x7,5}, {0xB,5}, {0x3,5},
72 {0xD,5}, {0x5,5}, {0x9,5}, {0x1,5},
73 {0xE,5}, {0x6,5}, {0xA,5}, {0x2,5},
74 {0xC,5}, {0x4,5}, {0x8,5}, {0x3,2},
78 static const uint8_t level_tab[7][2]={
79 {3,4}, {3,3}, {3,2}, {0,3}, {2,2}, {2,3}, {2,4}
82 static const uint8_t dc_ccp_tab[8][2]={
83 {0x1,2}, {0xD,4}, {0xF,4}, {0xC,4},
84 {0x5,3}, {0xE,4}, {0x4,3}, {0x0,2},
87 static const uint8_t ac_ccp_tab[16][2]={
88 {0x00,2}, {0x3B,6}, {0x0A,4}, {0x3A,6},
89 {0x02,3}, {0x39,6}, {0x3C,6}, {0x38,6},
90 {0x03,3}, {0x3D,6}, {0x08,4}, {0x1F,5},
91 {0x09,4}, {0x0B,4}, {0x0D,4}, {0x0C,4},
94 static const uint8_t asv2_level_tab[63][2]={
95 {0x3F,10},{0x2F,10},{0x37,10},{0x27,10},{0x3B,10},{0x2B,10},{0x33,10},{0x23,10},
96 {0x3D,10},{0x2D,10},{0x35,10},{0x25,10},{0x39,10},{0x29,10},{0x31,10},{0x21,10},
97 {0x1F, 8},{0x17, 8},{0x1B, 8},{0x13, 8},{0x1D, 8},{0x15, 8},{0x19, 8},{0x11, 8},
98 {0x0F, 6},{0x0B, 6},{0x0D, 6},{0x09, 6},
104 {0x08, 6},{0x0C, 6},{0x0A, 6},{0x0E, 6},
105 {0x10, 8},{0x18, 8},{0x14, 8},{0x1C, 8},{0x12, 8},{0x1A, 8},{0x16, 8},{0x1E, 8},
106 {0x20,10},{0x30,10},{0x28,10},{0x38,10},{0x24,10},{0x34,10},{0x2C,10},{0x3C,10},
107 {0x22,10},{0x32,10},{0x2A,10},{0x3A,10},{0x26,10},{0x36,10},{0x2E,10},{0x3E,10},
112 static VLC level_vlc;
113 static VLC dc_ccp_vlc;
114 static VLC ac_ccp_vlc;
115 static VLC asv2_level_vlc;
117 static av_cold void init_vlcs(ASV1Context *a){
123 INIT_VLC_STATIC(&ccp_vlc, VLC_BITS, 17,
124 &ccp_tab[0][1], 2, 1,
125 &ccp_tab[0][0], 2, 1, 64);
126 INIT_VLC_STATIC(&dc_ccp_vlc, VLC_BITS, 8,
127 &dc_ccp_tab[0][1], 2, 1,
128 &dc_ccp_tab[0][0], 2, 1, 64);
129 INIT_VLC_STATIC(&ac_ccp_vlc, VLC_BITS, 16,
130 &ac_ccp_tab[0][1], 2, 1,
131 &ac_ccp_tab[0][0], 2, 1, 64);
132 INIT_VLC_STATIC(&level_vlc, VLC_BITS, 7,
133 &level_tab[0][1], 2, 1,
134 &level_tab[0][0], 2, 1, 64);
135 INIT_VLC_STATIC(&asv2_level_vlc, ASV2_LEVEL_VLC_BITS, 63,
136 &asv2_level_tab[0][1], 2, 1,
137 &asv2_level_tab[0][0], 2, 1, 1024);
141 //FIXME write a reversed bitstream reader to avoid the double reverse
142 static inline int asv2_get_bits(GetBitContext *gb, int n){
143 return av_reverse[ get_bits(gb, n) << (8-n) ];
146 static inline void asv2_put_bits(PutBitContext *pb, int n, int v){
147 put_bits(pb, n, av_reverse[ v << (8-n) ]);
150 static inline int asv1_get_level(GetBitContext *gb){
151 int code= get_vlc2(gb, level_vlc.table, VLC_BITS, 1);
153 if(code==3) return get_sbits(gb, 8);
154 else return code - 3;
157 static inline int asv2_get_level(GetBitContext *gb){
158 int code= get_vlc2(gb, asv2_level_vlc.table, ASV2_LEVEL_VLC_BITS, 1);
160 if(code==31) return (int8_t)asv2_get_bits(gb, 8);
161 else return code - 31;
164 static inline void asv1_put_level(PutBitContext *pb, int level){
165 unsigned int index= level + 3;
167 if(index <= 6) put_bits(pb, level_tab[index][1], level_tab[index][0]);
169 put_bits(pb, level_tab[3][1], level_tab[3][0]);
170 put_sbits(pb, 8, level);
174 static inline void asv2_put_level(PutBitContext *pb, int level){
175 unsigned int index= level + 31;
177 if(index <= 62) put_bits(pb, asv2_level_tab[index][1], asv2_level_tab[index][0]);
179 put_bits(pb, asv2_level_tab[31][1], asv2_level_tab[31][0]);
180 asv2_put_bits(pb, 8, level&0xFF);
184 static inline int asv1_decode_block(ASV1Context *a, DCTELEM block[64]){
187 block[0]= 8*get_bits(&a->gb, 8);
190 const int ccp= get_vlc2(&a->gb, ccp_vlc.table, VLC_BITS, 1);
194 if(ccp < 0 || i>=10){
195 av_log(a->avctx, AV_LOG_ERROR, "coded coeff pattern damaged\n");
199 if(ccp&8) block[a->scantable.permutated[4*i+0]]= (asv1_get_level(&a->gb) * a->intra_matrix[4*i+0])>>4;
200 if(ccp&4) block[a->scantable.permutated[4*i+1]]= (asv1_get_level(&a->gb) * a->intra_matrix[4*i+1])>>4;
201 if(ccp&2) block[a->scantable.permutated[4*i+2]]= (asv1_get_level(&a->gb) * a->intra_matrix[4*i+2])>>4;
202 if(ccp&1) block[a->scantable.permutated[4*i+3]]= (asv1_get_level(&a->gb) * a->intra_matrix[4*i+3])>>4;
209 static inline int asv2_decode_block(ASV1Context *a, DCTELEM block[64]){
212 count= asv2_get_bits(&a->gb, 4);
214 block[0]= 8*asv2_get_bits(&a->gb, 8);
216 ccp= get_vlc2(&a->gb, dc_ccp_vlc.table, VLC_BITS, 1);
218 if(ccp&4) block[a->scantable.permutated[1]]= (asv2_get_level(&a->gb) * a->intra_matrix[1])>>4;
219 if(ccp&2) block[a->scantable.permutated[2]]= (asv2_get_level(&a->gb) * a->intra_matrix[2])>>4;
220 if(ccp&1) block[a->scantable.permutated[3]]= (asv2_get_level(&a->gb) * a->intra_matrix[3])>>4;
223 for(i=1; i<count+1; i++){
224 const int ccp= get_vlc2(&a->gb, ac_ccp_vlc.table, VLC_BITS, 1);
227 if(ccp&8) block[a->scantable.permutated[4*i+0]]= (asv2_get_level(&a->gb) * a->intra_matrix[4*i+0])>>4;
228 if(ccp&4) block[a->scantable.permutated[4*i+1]]= (asv2_get_level(&a->gb) * a->intra_matrix[4*i+1])>>4;
229 if(ccp&2) block[a->scantable.permutated[4*i+2]]= (asv2_get_level(&a->gb) * a->intra_matrix[4*i+2])>>4;
230 if(ccp&1) block[a->scantable.permutated[4*i+3]]= (asv2_get_level(&a->gb) * a->intra_matrix[4*i+3])>>4;
237 static inline void asv1_encode_block(ASV1Context *a, DCTELEM block[64]){
241 put_bits(&a->pb, 8, (block[0] + 32)>>6);
245 const int index= scantab[4*i];
248 if( (block[index + 0] = (block[index + 0]*a->q_intra_matrix[index + 0] + (1<<15))>>16) ) ccp |= 8;
249 if( (block[index + 8] = (block[index + 8]*a->q_intra_matrix[index + 8] + (1<<15))>>16) ) ccp |= 4;
250 if( (block[index + 1] = (block[index + 1]*a->q_intra_matrix[index + 1] + (1<<15))>>16) ) ccp |= 2;
251 if( (block[index + 9] = (block[index + 9]*a->q_intra_matrix[index + 9] + (1<<15))>>16) ) ccp |= 1;
254 for(;nc_count; nc_count--)
255 put_bits(&a->pb, ccp_tab[0][1], ccp_tab[0][0]);
257 put_bits(&a->pb, ccp_tab[ccp][1], ccp_tab[ccp][0]);
259 if(ccp&8) asv1_put_level(&a->pb, block[index + 0]);
260 if(ccp&4) asv1_put_level(&a->pb, block[index + 8]);
261 if(ccp&2) asv1_put_level(&a->pb, block[index + 1]);
262 if(ccp&1) asv1_put_level(&a->pb, block[index + 9]);
267 put_bits(&a->pb, ccp_tab[16][1], ccp_tab[16][0]);
270 static inline void asv2_encode_block(ASV1Context *a, DCTELEM block[64]){
274 for(count=63; count>3; count--){
275 const int index= scantab[count];
277 if( (block[index]*a->q_intra_matrix[index] + (1<<15))>>16 )
283 asv2_put_bits(&a->pb, 4, count);
284 asv2_put_bits(&a->pb, 8, (block[0] + 32)>>6);
287 for(i=0; i<=count; i++){
288 const int index= scantab[4*i];
291 if( (block[index + 0] = (block[index + 0]*a->q_intra_matrix[index + 0] + (1<<15))>>16) ) ccp |= 8;
292 if( (block[index + 8] = (block[index + 8]*a->q_intra_matrix[index + 8] + (1<<15))>>16) ) ccp |= 4;
293 if( (block[index + 1] = (block[index + 1]*a->q_intra_matrix[index + 1] + (1<<15))>>16) ) ccp |= 2;
294 if( (block[index + 9] = (block[index + 9]*a->q_intra_matrix[index + 9] + (1<<15))>>16) ) ccp |= 1;
297 if(i) put_bits(&a->pb, ac_ccp_tab[ccp][1], ac_ccp_tab[ccp][0]);
298 else put_bits(&a->pb, dc_ccp_tab[ccp][1], dc_ccp_tab[ccp][0]);
301 if(ccp&8) asv2_put_level(&a->pb, block[index + 0]);
302 if(ccp&4) asv2_put_level(&a->pb, block[index + 8]);
303 if(ccp&2) asv2_put_level(&a->pb, block[index + 1]);
304 if(ccp&1) asv2_put_level(&a->pb, block[index + 9]);
309 static inline int decode_mb(ASV1Context *a, DCTELEM block[6][64]){
312 a->dsp.clear_blocks(block[0]);
314 if(a->avctx->codec_id == CODEC_ID_ASV1){
316 if( asv1_decode_block(a, block[i]) < 0)
321 if( asv2_decode_block(a, block[i]) < 0)
328 #define MAX_MB_SIZE (30*16*16*3/2/8)
330 static inline int encode_mb(ASV1Context *a, DCTELEM block[6][64]){
333 if (a->pb.buf_end - a->pb.buf - (put_bits_count(&a->pb)>>3) < MAX_MB_SIZE) {
334 av_log(a->avctx, AV_LOG_ERROR, "encoded frame too large\n");
338 if(a->avctx->codec_id == CODEC_ID_ASV1){
340 asv1_encode_block(a, block[i]);
343 asv2_encode_block(a, block[i]);
348 static inline void idct_put(ASV1Context *a, int mb_x, int mb_y){
349 DCTELEM (*block)[64]= a->block;
350 int linesize= a->picture.linesize[0];
352 uint8_t *dest_y = a->picture.data[0] + (mb_y * 16* linesize ) + mb_x * 16;
353 uint8_t *dest_cb = a->picture.data[1] + (mb_y * 8 * a->picture.linesize[1]) + mb_x * 8;
354 uint8_t *dest_cr = a->picture.data[2] + (mb_y * 8 * a->picture.linesize[2]) + mb_x * 8;
356 a->dsp.idct_put(dest_y , linesize, block[0]);
357 a->dsp.idct_put(dest_y + 8, linesize, block[1]);
358 a->dsp.idct_put(dest_y + 8*linesize , linesize, block[2]);
359 a->dsp.idct_put(dest_y + 8*linesize + 8, linesize, block[3]);
361 if(!(a->avctx->flags&CODEC_FLAG_GRAY)){
362 a->dsp.idct_put(dest_cb, a->picture.linesize[1], block[4]);
363 a->dsp.idct_put(dest_cr, a->picture.linesize[2], block[5]);
367 static inline void dct_get(ASV1Context *a, int mb_x, int mb_y){
368 DCTELEM (*block)[64]= a->block;
369 int linesize= a->picture.linesize[0];
372 uint8_t *ptr_y = a->picture.data[0] + (mb_y * 16* linesize ) + mb_x * 16;
373 uint8_t *ptr_cb = a->picture.data[1] + (mb_y * 8 * a->picture.linesize[1]) + mb_x * 8;
374 uint8_t *ptr_cr = a->picture.data[2] + (mb_y * 8 * a->picture.linesize[2]) + mb_x * 8;
376 a->dsp.get_pixels(block[0], ptr_y , linesize);
377 a->dsp.get_pixels(block[1], ptr_y + 8, linesize);
378 a->dsp.get_pixels(block[2], ptr_y + 8*linesize , linesize);
379 a->dsp.get_pixels(block[3], ptr_y + 8*linesize + 8, linesize);
381 a->dsp.fdct(block[i]);
383 if(!(a->avctx->flags&CODEC_FLAG_GRAY)){
384 a->dsp.get_pixels(block[4], ptr_cb, a->picture.linesize[1]);
385 a->dsp.get_pixels(block[5], ptr_cr, a->picture.linesize[2]);
387 a->dsp.fdct(block[i]);
391 static int decode_frame(AVCodecContext *avctx,
392 void *data, int *data_size,
395 const uint8_t *buf = avpkt->data;
396 int buf_size = avpkt->size;
397 ASV1Context * const a = avctx->priv_data;
398 AVFrame *picture = data;
399 AVFrame * const p= &a->picture;
403 avctx->release_buffer(avctx, p);
406 if(avctx->get_buffer(avctx, p) < 0){
407 av_log(avctx, AV_LOG_ERROR, "get_buffer() failed\n");
410 p->pict_type= AV_PICTURE_TYPE_I;
413 av_fast_padded_malloc(&a->bitstream_buffer, &a->bitstream_buffer_size,
415 if (!a->bitstream_buffer)
416 return AVERROR(ENOMEM);
418 if(avctx->codec_id == CODEC_ID_ASV1)
419 a->dsp.bswap_buf((uint32_t*)a->bitstream_buffer, (const uint32_t*)buf, buf_size/4);
422 for(i=0; i<buf_size; i++)
423 a->bitstream_buffer[i]= av_reverse[ buf[i] ];
426 init_get_bits(&a->gb, a->bitstream_buffer, buf_size*8);
428 for(mb_y=0; mb_y<a->mb_height2; mb_y++){
429 for(mb_x=0; mb_x<a->mb_width2; mb_x++){
430 if( decode_mb(a, a->block) <0)
433 idct_put(a, mb_x, mb_y);
437 if(a->mb_width2 != a->mb_width){
439 for(mb_y=0; mb_y<a->mb_height2; mb_y++){
440 if( decode_mb(a, a->block) <0)
443 idct_put(a, mb_x, mb_y);
447 if(a->mb_height2 != a->mb_height){
449 for(mb_x=0; mb_x<a->mb_width; mb_x++){
450 if( decode_mb(a, a->block) <0)
453 idct_put(a, mb_x, mb_y);
457 *picture = a->picture;
458 *data_size = sizeof(AVPicture);
462 return (get_bits_count(&a->gb)+31)/32*4;
465 #if CONFIG_ASV1_ENCODER || CONFIG_ASV2_ENCODER
466 static int encode_frame(AVCodecContext *avctx, AVPacket *pkt,
467 const AVFrame *pict, int *got_packet)
469 ASV1Context * const a = avctx->priv_data;
470 AVFrame * const p= &a->picture;
475 (ret = av_new_packet(pkt, a->mb_height*a->mb_width*MAX_MB_SIZE +
476 FF_MIN_BUFFER_SIZE)) < 0) {
477 av_log(avctx, AV_LOG_ERROR, "Error getting output packet.\n");
481 init_put_bits(&a->pb, pkt->data, pkt->size);
484 p->pict_type= AV_PICTURE_TYPE_I;
487 for(mb_y=0; mb_y<a->mb_height2; mb_y++){
488 for(mb_x=0; mb_x<a->mb_width2; mb_x++){
489 dct_get(a, mb_x, mb_y);
490 encode_mb(a, a->block);
494 if(a->mb_width2 != a->mb_width){
496 for(mb_y=0; mb_y<a->mb_height2; mb_y++){
497 dct_get(a, mb_x, mb_y);
498 encode_mb(a, a->block);
502 if(a->mb_height2 != a->mb_height){
504 for(mb_x=0; mb_x<a->mb_width; mb_x++){
505 dct_get(a, mb_x, mb_y);
506 encode_mb(a, a->block);
511 avpriv_align_put_bits(&a->pb);
512 while(put_bits_count(&a->pb)&31)
513 put_bits(&a->pb, 8, 0);
515 size= put_bits_count(&a->pb)/32;
517 if(avctx->codec_id == CODEC_ID_ASV1)
518 a->dsp.bswap_buf((uint32_t*)pkt->data, (uint32_t*)pkt->data, size);
521 for(i=0; i<4*size; i++)
522 pkt->data[i] = av_reverse[pkt->data[i]];
526 pkt->flags |= AV_PKT_FLAG_KEY;
531 #endif /* CONFIG_ASV1_ENCODER || CONFIG_ASV2_ENCODER */
533 static av_cold void common_init(AVCodecContext *avctx){
534 ASV1Context * const a = avctx->priv_data;
536 ff_dsputil_init(&a->dsp, avctx);
538 a->mb_width = (avctx->width + 15) / 16;
539 a->mb_height = (avctx->height + 15) / 16;
540 a->mb_width2 = (avctx->width + 0) / 16;
541 a->mb_height2 = (avctx->height + 0) / 16;
543 avctx->coded_frame= &a->picture;
547 static av_cold int decode_init(AVCodecContext *avctx){
548 ASV1Context * const a = avctx->priv_data;
549 AVFrame *p= &a->picture;
551 const int scale= avctx->codec_id == CODEC_ID_ASV1 ? 1 : 2;
555 ff_init_scantable(a->dsp.idct_permutation, &a->scantable, scantab);
556 avctx->pix_fmt= PIX_FMT_YUV420P;
558 a->inv_qscale= avctx->extradata[0];
559 if(a->inv_qscale == 0){
560 av_log(avctx, AV_LOG_ERROR, "illegal qscale 0\n");
561 if(avctx->codec_id == CODEC_ID_ASV1)
568 int index= scantab[i];
570 a->intra_matrix[i]= 64*scale*ff_mpeg1_default_intra_matrix[index] / a->inv_qscale;
573 p->qstride= a->mb_width;
574 p->qscale_table= av_malloc( p->qstride * a->mb_height);
575 p->quality= (32*scale + a->inv_qscale/2)/a->inv_qscale;
576 memset(p->qscale_table, p->quality, p->qstride*a->mb_height);
581 #if CONFIG_ASV1_ENCODER || CONFIG_ASV2_ENCODER
582 static av_cold int encode_init(AVCodecContext *avctx){
583 ASV1Context * const a = avctx->priv_data;
585 const int scale= avctx->codec_id == CODEC_ID_ASV1 ? 1 : 2;
589 if(avctx->global_quality == 0) avctx->global_quality= 4*FF_QUALITY_SCALE;
591 a->inv_qscale= (32*scale*FF_QUALITY_SCALE + avctx->global_quality/2) / avctx->global_quality;
593 avctx->extradata= av_mallocz(8);
594 avctx->extradata_size=8;
595 ((uint32_t*)avctx->extradata)[0]= av_le2ne32(a->inv_qscale);
596 ((uint32_t*)avctx->extradata)[1]= av_le2ne32(AV_RL32("ASUS"));
599 int q= 32*scale*ff_mpeg1_default_intra_matrix[i];
600 a->q_intra_matrix[i]= ((a->inv_qscale<<16) + q/2) / q;
605 #endif /* CONFIG_ASV1_ENCODER || CONFIG_ASV2_ENCODER */
607 static av_cold int decode_end(AVCodecContext *avctx){
608 ASV1Context * const a = avctx->priv_data;
610 av_freep(&a->bitstream_buffer);
611 av_freep(&a->picture.qscale_table);
612 a->bitstream_buffer_size=0;
614 if(a->picture.data[0])
615 avctx->release_buffer(avctx, &a->picture);
620 AVCodec ff_asv1_decoder = {
622 .type = AVMEDIA_TYPE_VIDEO,
624 .priv_data_size = sizeof(ASV1Context),
627 .decode = decode_frame,
628 .capabilities = CODEC_CAP_DR1,
629 .long_name = NULL_IF_CONFIG_SMALL("ASUS V1"),
632 AVCodec ff_asv2_decoder = {
634 .type = AVMEDIA_TYPE_VIDEO,
636 .priv_data_size = sizeof(ASV1Context),
639 .decode = decode_frame,
640 .capabilities = CODEC_CAP_DR1,
641 .long_name = NULL_IF_CONFIG_SMALL("ASUS V2"),
644 #if CONFIG_ASV1_ENCODER
645 AVCodec ff_asv1_encoder = {
647 .type = AVMEDIA_TYPE_VIDEO,
649 .priv_data_size = sizeof(ASV1Context),
651 .encode2 = encode_frame,
652 .pix_fmts = (const enum PixelFormat[]){ PIX_FMT_YUV420P, PIX_FMT_NONE },
653 .long_name = NULL_IF_CONFIG_SMALL("ASUS V1"),
657 #if CONFIG_ASV2_ENCODER
658 AVCodec ff_asv2_encoder = {
660 .type = AVMEDIA_TYPE_VIDEO,
662 .priv_data_size = sizeof(ASV1Context),
664 .encode2 = encode_frame,
665 .pix_fmts = (const enum PixelFormat[]){ PIX_FMT_YUV420P, PIX_FMT_NONE },
666 .long_name = NULL_IF_CONFIG_SMALL("ASUS V2"),