libavcodec/cabac.c
d592f67f
 /*
  * H.26L/H.264/AVC/JVT/14496-10/... encoder/decoder
  * Copyright (c) 2003 Michael Niedermayer <michaelni@gmx.at>
  *
b78e7197
  * This file is part of FFmpeg.
  *
  * FFmpeg is free software; you can redistribute it and/or
d592f67f
  * modify it under the terms of the GNU Lesser General Public
  * License as published by the Free Software Foundation; either
b78e7197
  * version 2.1 of the License, or (at your option) any later version.
d592f67f
  *
b78e7197
  * FFmpeg is distributed in the hope that it will be useful,
d592f67f
  * but WITHOUT ANY WARRANTY; without even the implied warranty of
  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
  * Lesser General Public License for more details.
  *
  * You should have received a copy of the GNU Lesser General Public
b78e7197
  * License along with FFmpeg; if not, write to the Free Software
5509bffa
  * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
d592f67f
  */
 
 /**
ba87f080
  * @file
d592f67f
  * Context Adaptive Binary Arithmetic Coder.
  */
 
 #include <string.h>
 
245976da
 #include "libavutil/common.h"
965fa6b0
 #include "libavutil/timer.h"
9106a698
 #include "get_bits.h"
d592f67f
 #include "cabac.h"
55b9ef18
 #include "cabac_functions.h"
d592f67f
 
092d1977
 #include "cabac_tablegen.h"
ec7eb896
 
d592f67f
 /**
  *
  * @param buf_size size of buf in bits
  */
 void ff_init_cabac_encoder(CABACContext *c, uint8_t *buf, int buf_size){
ed7debda
     init_put_bits(&c->pb, buf, buf_size);
d592f67f
 
     c->low= 0;
     c->range= 0x1FE;
     c->outstanding_count= 0;
     c->pb.bit_left++; //avoids firstBitFlag
 }
 
 /**
  *
  * @param buf_size size of buf in bits
  */
e96682e6
 void ff_init_cabac_decoder(CABACContext *c, const uint8_t *buf, int buf_size){
115329f1
     c->bytestream_start=
d592f67f
     c->bytestream= buf;
bba83349
     c->bytestream_end= buf + buf_size;
d592f67f
 
ec7eb896
 #if CABAC_BITS == 16
     c->low =  (*c->bytestream++)<<18;
     c->low+=  (*c->bytestream++)<<10;
 #else
     c->low =  (*c->bytestream++)<<10;
 #endif
     c->low+= ((*c->bytestream++)<<2) + 2;
f24a5159
     c->range= 0x1FE;
d592f67f
 }
 
0fe4b485
 void ff_init_cabac_states(void)
 {
cab8c5f8
     static int initialized = 0;
 
     if (initialized)
         return;
115329f1
 
092d1977
     cabac_tableinit();
cab8c5f8
 
     initialized = 1;
d592f67f
 }
 
07e4e3ea
 #ifdef TEST
d592f67f
 #define SIZE 10240
6041c217
 
294eaa26
 #include "libavutil/lfg.h"
6041c217
 #include "avcodec.h"
 
badb195d
 static inline void put_cabac_bit(CABACContext *c, int b){
     put_bits(&c->pb, 1, b);
     for(;c->outstanding_count; c->outstanding_count--){
         put_bits(&c->pb, 1, 1-b);
     }
 }
 
 static inline void renorm_cabac_encoder(CABACContext *c){
     while(c->range < 0x100){
         //FIXME optimize
         if(c->low<0x100){
             put_cabac_bit(c, 0);
         }else if(c->low<0x200){
             c->outstanding_count++;
             c->low -= 0x100;
         }else{
             put_cabac_bit(c, 1);
             c->low -= 0x200;
         }
 
         c->range+= c->range;
         c->low += c->low;
     }
 }
 
fd065d3e
 static void put_cabac(CABACContext *c, uint8_t * const state, int bit){
     int RangeLPS= ff_h264_lps_range[2*(c->range&0xC0) + *state];
 
     if(bit == ((*state)&1)){
         c->range -= RangeLPS;
b0b66269
         *state    = ff_h264_mlps_state[128 + *state];
fd065d3e
     }else{
         c->low += c->range - RangeLPS;
         c->range = RangeLPS;
ab6fbe82
         *state= ff_h264_mlps_state[127 - *state];
fd065d3e
     }
 
     renorm_cabac_encoder(c);
 }
 
 /**
  * @param bit 0 -> write zero bit, !=0 write one bit
  */
 static void put_cabac_bypass(CABACContext *c, int bit){
     c->low += c->low;
 
     if(bit){
         c->low += c->range;
     }
 //FIXME optimize
     if(c->low<0x200){
         put_cabac_bit(c, 0);
     }else if(c->low<0x400){
         c->outstanding_count++;
         c->low -= 0x200;
     }else{
         put_cabac_bit(c, 1);
         c->low -= 0x400;
     }
 }
 
 /**
  *
  * @return the number of bytes written
  */
 static int put_cabac_terminate(CABACContext *c, int bit){
     c->range -= 2;
 
     if(!bit){
         renorm_cabac_encoder(c);
     }else{
         c->low += c->range;
         c->range= 2;
 
         renorm_cabac_encoder(c);
 
dc7aecd8
         av_assert0(c->low <= 0x1FF);
fd065d3e
         put_cabac_bit(c, c->low>>9);
         put_bits(&c->pb, 2, ((c->low>>7)&3)|1);
 
         flush_put_bits(&c->pb); //FIXME FIXME FIXME XXX wrong
     }
 
     return (put_bits_count(&c->pb)+7)>>3;
 }
 
f8a80fd6
 int main(void){
d592f67f
     CABACContext c;
     uint8_t b[9*SIZE];
61ccfcc0
     uint8_t r[9*SIZE];
d592f67f
     int i;
61ccfcc0
     uint8_t state[10]= {0};
64bde197
     AVLFG prng;
115329f1
 
64bde197
     av_lfg_init(&prng, 1);
d592f67f
     ff_init_cabac_encoder(&c, b, SIZE);
bfb39023
     ff_init_cabac_states();
115329f1
 
d592f67f
     for(i=0; i<SIZE; i++){
1eb805ed
         if(2*i<SIZE) r[i] = av_lfg_get(&prng) % 7;
         else         r[i] = (i>>8)&1;
d592f67f
     }
115329f1
 
d592f67f
     for(i=0; i<SIZE; i++){
 START_TIMER
61ccfcc0
         put_cabac_bypass(&c, r[i]&1);
d592f67f
 STOP_TIMER("put_cabac_bypass")
     }
115329f1
 
d592f67f
     for(i=0; i<SIZE; i++){
 START_TIMER
61ccfcc0
         put_cabac(&c, state, r[i]&1);
d592f67f
 STOP_TIMER("put_cabac")
     }
 
     put_cabac_terminate(&c, 1);
115329f1
 
d592f67f
     ff_init_cabac_decoder(&c, b, SIZE);
115329f1
 
61ccfcc0
     memset(state, 0, sizeof(state));
115329f1
 
d592f67f
     for(i=0; i<SIZE; i++){
 START_TIMER
61ccfcc0
         if( (r[i]&1) != get_cabac_bypass(&c) )
6041c217
             av_log(NULL, AV_LOG_ERROR, "CABAC bypass failure at %d\n", i);
d592f67f
 STOP_TIMER("get_cabac_bypass")
     }
115329f1
 
d592f67f
     for(i=0; i<SIZE; i++){
 START_TIMER
0538b29a
         if( (r[i]&1) != get_cabac_noinline(&c, state) )
6041c217
             av_log(NULL, AV_LOG_ERROR, "CABAC failure at %d\n", i);
d592f67f
 STOP_TIMER("get_cabac")
     }
     if(!get_cabac_terminate(&c))
6041c217
         av_log(NULL, AV_LOG_ERROR, "where's the Terminator?\n");
115329f1
 
d592f67f
     return 0;
 }
 
07e4e3ea
 #endif /* TEST */