22 #ifndef AVCODEC_SNOW_H
23 #define AVCODEC_SNOW_H
36 #define QROOT (1<<QSHIFT)
37 #define LOSSLESS_QLOG -128
39 #define MAX_REF_FRAMES 8
41 #define LOG2_OBMC_MAX 8
42 #define OBMC_MAX (1<<(LOG2_OBMC_MAX))
57 .
color= {128,128,128},
65 #define LOG2_MB_SIZE 4
66 #define MB_SIZE (1<<LOG2_MB_SIZE)
67 #define ENCODER_EXTRA_BITS 4
149 #define QBIAS_SHIFT 3
156 #define ME_CACHE_SIZE 1024
168 extern const uint8_t *
const obmc_tab[4];
182 low[(*i)+1] = low[((*i)+1)>>1];
188 for (; (*i)>=0; (*i)-=2){
189 low[(*i)+1] = high[(*i)>>1];
190 low[*i] = low[(*i)>>1];
196 dst[i] = src[i] - ((mul * (ref[i] + ref[i + 1]) + add) >> shift);
199 if((width^lift_high)&1){
200 dst[w] = src[w] - ((mul * 2 * ref[w] + add) >> shift);
206 dst[i] = src[i] + ((ref[i] + ref[(i+1)]+
W_BO + 4 * src[i]) >>
W_BS);
210 dst[w] = src[w] + ((2 * ref[w] +
W_BO + 4 * src[w]) >>
W_BS);
225 int plane_index,
int w,
int h);
229 static inline void snow_set_blocks(
SnowContext *s,
int level,
int x,
int y,
int l,
int cb,
int cr,
int mx,
int my,
int ref,
int type){
232 const int index= (x + y*w) << rem_depth;
233 const int block_w= 1<<rem_depth;
246 for(j=0; j<block_w; j++){
247 for(i=0; i<block_w; i++){
261 (top ->
mx * scale[top ->
ref] + 128) >>8,
262 (tr ->
mx * scale[tr ->
ref] + 128) >>8);
264 (top ->
my * scale[top ->
ref] + 128) >>8,
265 (tr ->
my * scale[tr ->
ref] + 128) >>8);
279 static av_always_inline void add_yblock(
SnowContext *s,
int sliced,
slice_buffer *sb,
IDWTELEM *dst, uint8_t *dst8,
const uint8_t *obmc,
int src_x,
int src_y,
int b_w,
int b_h,
int w,
int h,
int dst_stride,
int src_stride,
int obmc_stride,
int b_x,
int b_y,
int add,
int offset_dst,
int plane_index){
282 const int b_stride= b_width;
296 }
else if(b_x + 1 >= b_width){
303 }
else if(b_y + 1 >= b_height){
311 if(!sliced && !offset_dst)
314 }
else if(src_x + b_w > w){
318 obmc -= src_y*obmc_stride;
320 if(!sliced && !offset_dst)
321 dst -= src_y*dst_stride;
323 }
else if(src_y + b_h> h){
327 if(b_w<=0 || b_h<=0)
return;
329 assert(src_stride > 2*
MB_SIZE + 5);
331 if(!sliced && offset_dst)
332 dst += src_x + src_y*dst_stride;
333 dst8+= src_x + src_y*src_stride;
336 ptmp= tmp + 3*tmp_step;
339 ff_snow_pred_block(s, block[0], tmp, src_stride, src_x, src_y, b_w, b_h, lt, plane_index, w, h);
346 ff_snow_pred_block(s, block[1], tmp, src_stride, src_x, src_y, b_w, b_h, rt, plane_index, w, h);
356 ff_snow_pred_block(s, block[2], tmp, src_stride, src_x, src_y, b_w, b_h, lb, plane_index, w, h);
367 ff_snow_pred_block(s, block[3], tmp, src_stride, src_x, src_y, b_w, b_h, rb, plane_index, w, h);
370 s->
dwt.
inner_add_yblock(obmc, obmc_stride, block, b_w, b_h, src_x,src_y, src_stride, sb, add, dst8);
372 for(y=0; y<b_h; y++){
374 const uint8_t *obmc1= obmc + y*obmc_stride;
375 const uint8_t *obmc2= obmc1+ (obmc_stride>>1);
376 const uint8_t *obmc3= obmc1+ obmc_stride*(obmc_stride>>1);
377 const uint8_t *
obmc4= obmc3+ (obmc_stride>>1);
378 for(x=0; x<b_w; x++){
379 int v= obmc1[x] * block[3][x + y*src_stride]
380 +obmc2[x] * block[2][x + y*src_stride]
381 +obmc3[x] * block[1][x + y*src_stride]
382 +obmc4[x] * block[0][x + y*src_stride];
389 v += dst[x + y*dst_stride];
391 if(v&(~255)) v= ~(v>>31);
392 dst8[x + y*src_stride] =
v;
394 dst[x + y*dst_stride] -=
v;
407 int block_w = plane_index ? block_size/2 : block_size;
409 const int obmc_stride= plane_index ? block_size : 2*block_size;
420 for(y=block_w*mb_y; y<
FFMIN(h,block_w*(mb_y+1)); y++){
424 if(v&(~255)) v= ~(v>>31);
425 dst8[x + y*ref_stride]=
v;
429 for(y=block_w*mb_y; y<
FFMIN(h,block_w*(mb_y+1)); y++){
439 for(mb_x=0; mb_x<=mb_w; mb_x++){
441 block_w*mb_x - block_w/2,
442 block_w*mb_y - block_w/2,
445 w, ref_stride, obmc_stride,
447 add, 1, plane_index);
454 for(mb_y=0; mb_y<=mb_h; mb_y++)
458 static inline void set_blocks(
SnowContext *s,
int level,
int x,
int y,
int l,
int cb,
int cr,
int mx,
int my,
int ref,
int type){
461 const int index= (x + y*w) << rem_depth;
462 const int block_w= 1<<rem_depth;
475 for(j=0; j<block_w; j++){
476 for(i=0; i<block_w; i++){
483 const int offset[3]= {
490 c->
src[0][i]= src [i];
491 c->
ref[0][i]= ref [i] + offset[i];
501 #define QEXPSHIFT (7-FRAC_BITS+8) //FIXME try to change this to 0
507 const int a=
FFABS(v);
508 const int e= av_log2(a);
509 const int el=
FFMIN(e, 10);
520 for(i=e-1; i>=el; i--){
521 put_rac(c, state+22+9, (a>>i)&1);
524 put_rac(c, state+22+i, (a>>i)&1);
528 put_rac(c, state+11 + el, v < 0);
545 for(i=e-1; i>=0; i--){
556 int r= log2>=0 ? 1<<log2 : 1;
569 for(i=log2-1; i>=0; i--){
570 put_rac(c, state+31-i, (v>>i)&1);
576 int r= log2>=0 ? 1<<log2 : 1;
581 while(
get_rac(c, state+4+log2)){
587 for(i=log2-1; i>=0; i--){
595 const int w= b->
width;
614 if(y && prev_xc->
x == 0){
626 if(prev_xc->
x == x + 1)
632 if(x>>1 > parent_xc->
x){
635 if(x>>1 == parent_xc->
x){
640 int context= av_log2(3*(l>>1) + (lt>>1) + (
t&~1) + (rt>>1) + (p>>1));
664 if(y) max_run=
FFMIN(run, prev_xc->
x - x - 2);
665 else max_run=
FFMIN(run, w-x-1);
667 max_run=
FFMIN(max_run, 2*parent_xc->
x - x - 1);
679 while(parent_xc->
x != parent->
width+1)
682 prev_parent_xc= parent_xc;
684 parent_xc= prev_parent_xc;