22 #ifndef AVCODEC_SNOW_H
23 #define AVCODEC_SNOW_H
35 #define FF_MPV_OFFSET(x) (offsetof(MpegEncContext, x) + offsetof(SnowContext, m))
43 #define QROOT (1<<QSHIFT)
44 #define LOSSLESS_QLOG -128
46 #define MAX_REF_FRAMES 8
48 #define LOG2_OBMC_MAX 8
49 #define OBMC_MAX (1<<(LOG2_OBMC_MAX))
64 .
color= {128,128,128},
72 #define LOG2_MB_SIZE 4
73 #define MB_SIZE (1<<LOG2_MB_SIZE)
74 #define ENCODER_EXTRA_BITS 4
164 #define QBIAS_SHIFT 3
172 #define ME_CACHE_SIZE 1024
199 low[(*i)+1] = low[((*i)+1)>>1];
205 for (; (*i)>=0; (*i)-=2){
206 low[(*i)+1] = high[(*i)>>1];
207 low[*i] = low[(*i)>>1];
213 dst[i] = src[i] - ((mul * (ref[i] + ref[i + 1]) + add) >>
shift);
216 if((width^lift_high)&1){
217 dst[w] = src[w] - ((mul * 2 * ref[w] + add) >> shift);
223 dst[i] = src[i] + ((ref[i] + ref[(i+1)]+
W_BO + 4 * src[i]) >>
W_BS);
227 dst[w] = src[w] + ((2 * ref[w] +
W_BO + 4 * src[w]) >>
W_BS);
242 int plane_index,
int w,
int h);
247 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){
250 const int index= (x + y*w) << rem_depth;
251 const int block_w= 1<<rem_depth;
264 for(j=0; j<block_w; j++){
265 for(i=0; i<block_w; i++){
279 (top ->
mx * scale[top ->
ref] + 128) >>8,
280 (tr ->
mx * scale[tr ->
ref] + 128) >>8);
282 (top ->
my * scale[top ->
ref] + 128) >>8,
283 (tr ->
my * scale[tr ->
ref] + 128) >>8);
297 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){
300 const int b_stride= b_width;
314 }
else if(b_x + 1 >= b_width){
321 }
else if(b_y + 1 >= b_height){
329 if(!sliced && !offset_dst)
337 obmc -= src_y*obmc_stride;
339 if(!sliced && !offset_dst)
340 dst -= src_y*dst_stride;
347 if(b_w<=0 || b_h<=0)
return;
351 if(!sliced && offset_dst)
352 dst += src_x + src_y*dst_stride;
353 dst8+= src_x + src_y*src_stride;
356 ptmp= tmp + 3*tmp_step;
359 ff_snow_pred_block(s, block[0], tmp, src_stride, src_x, src_y, b_w, b_h, lt, plane_index, w, h);
366 ff_snow_pred_block(s, block[1], tmp, src_stride, src_x, src_y, b_w, b_h, rt, plane_index, w, h);
376 ff_snow_pred_block(s, block[2], tmp, src_stride, src_x, src_y, b_w, b_h, lb, plane_index, w, h);
387 ff_snow_pred_block(s, block[3], tmp, src_stride, src_x, src_y, b_w, b_h, rb, plane_index, w, h);
390 s->
dwt.
inner_add_yblock(obmc, obmc_stride, block, b_w, b_h, src_x,src_y, src_stride, sb, add, dst8);
392 for(y=0; y<b_h; y++){
394 const uint8_t *obmc1= obmc + y*obmc_stride;
395 const uint8_t *obmc2= obmc1+ (obmc_stride>>1);
396 const uint8_t *obmc3= obmc1+ obmc_stride*(obmc_stride>>1);
398 for(x=0; x<b_w; x++){
399 int v= obmc1[x] * block[3][x + y*src_stride]
400 +obmc2[x] * block[2][x + y*src_stride]
401 +obmc3[x] * block[1][x + y*src_stride]
402 +obmc4[x] * block[0][x + y*src_stride];
409 v += dst[x + y*dst_stride];
411 if(v&(~255)) v= ~(v>>31);
412 dst8[x + y*src_stride] =
v;
414 dst[x + y*dst_stride] -=
v;
427 int block_w = plane_index ? block_size>>s->
chroma_h_shift : block_size;
428 int block_h = plane_index ? block_size>>s->
chroma_v_shift : block_size;
430 const int obmc_stride= plane_index ? (2*block_size)>>s->
chroma_h_shift : 2*block_size;
441 for(y=block_h*mb_y; y<
FFMIN(h,block_h*(mb_y+1)); y++){
445 if(v&(~255)) v= ~(v>>31);
446 dst8[x + y*ref_stride]=
v;
450 for(y=block_h*mb_y; y<
FFMIN(h,block_h*(mb_y+1)); y++){
460 for(mb_x=0; mb_x<=mb_w; mb_x++){
462 block_w*mb_x - block_w/2,
463 block_h*mb_y - block_h/2,
466 w, ref_stride, obmc_stride,
468 add, 1, plane_index);
475 for(mb_y=0; mb_y<=mb_h; mb_y++)
479 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){
482 const int index= (x + y*w) << rem_depth;
483 const int block_w= 1<<rem_depth;
484 const int block_h= 1<<rem_depth;
497 for(j=0; j<block_h; j++){
498 for(i=0; i<block_w; i++){
513 c->
src[0][i]= src [i];
514 c->
ref[0][i]= ref [i] + offset[i];
524 #define QEXPSHIFT (7-FRAC_BITS+8) //FIXME try to change this to 0
532 const int el=
FFMIN(e, 10);
543 for(i=e-1; i>=el; i--){
544 put_rac(c, state+22+9, (a>>i)&1);
547 put_rac(c, state+22+i, (a>>i)&1);
551 put_rac(c, state+11 + el, v < 0);
568 for(i=e-1; i>=0; i--){
579 int r= log2>=0 ? 1<<log2 : 1;
592 for(i=log2-1; i>=0; i--){
593 put_rac(c, state+31-i, (v>>i)&1);
599 int r= log2>=0 ? 1<<log2 : 1;
604 while(log2<28 &&
get_rac(c, state+4+log2)){
610 for(i=log2-1; i>=0; i--){
618 const int w= b->
width;
637 if(y && prev_xc->
x == 0){
649 if(prev_xc->
x == x + 1)
655 if(x>>1 > parent_xc->
x){
658 if(x>>1 == parent_xc->
x){
663 int context=
av_log2(3*(l>>1) + (lt>>1) + (t&~1) + (rt>>1) + (p>>1));
669 if ((uint16_t)v != v) {
682 if ((uint16_t)v != v) {
694 if(y) max_run=
FFMIN(run, prev_xc->
x - x - 2);
695 else max_run=
FFMIN(run, w-x-1);
697 max_run=
FFMIN(max_run, 2*parent_xc->
x - x - 1);
711 while(parent_xc->
x != parent->
width+1)
714 prev_parent_xc= parent_xc;
716 parent_xc= prev_parent_xc;