libavcodec/vp9mvs.c
1c9f4b50
 /*
  * VP9 compatible video decoder
  *
  * Copyright (C) 2013 Ronald S. Bultje <rsbultje gmail com>
  * Copyright (C) 2013 Clément Bœsch <u pkh me>
  *
  * This file is part of FFmpeg.
  *
  * FFmpeg is free software; you can redistribute it and/or
  * modify it under the terms of the GNU Lesser General Public
  * License as published by the Free Software Foundation; either
  * version 2.1 of the License, or (at your option) any later version.
  *
  * FFmpeg is distributed in the hope that it will be useful,
  * 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
  * License along with FFmpeg; if not, write to the Free Software
  * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
  */
 
 #include "internal.h"
 #include "vp56.h"
 #include "vp9.h"
 #include "vp9data.h"
f8c01994
 #include "vp9dec.h"
1c9f4b50
 
 static av_always_inline void clamp_mv(VP56mv *dst, const VP56mv *src,
e59da0f7
                                       VP9TileData *td)
1c9f4b50
 {
e59da0f7
     dst->x = av_clip(src->x, td->min_mv.x, td->max_mv.x);
     dst->y = av_clip(src->y, td->min_mv.y, td->max_mv.y);
1c9f4b50
 }
 
e59da0f7
 static void find_ref_mvs(VP9TileData *td,
1c9f4b50
                          VP56mv *pmv, int ref, int z, int idx, int sb)
 {
     static const int8_t mv_ref_blk_off[N_BS_SIZES][8][2] = {
875f6955
         [BS_64x64] = { {  3, -1 }, { -1,  3 }, {  4, -1 }, { -1,  4 },
                        { -1, -1 }, {  0, -1 }, { -1,  0 }, {  6, -1 } },
         [BS_64x32] = { {  0, -1 }, { -1,  0 }, {  4, -1 }, { -1,  2 },
                        { -1, -1 }, {  0, -3 }, { -3,  0 }, {  2, -1 } },
         [BS_32x64] = { { -1,  0 }, {  0, -1 }, { -1,  4 }, {  2, -1 },
                        { -1, -1 }, { -3,  0 }, {  0, -3 }, { -1,  2 } },
         [BS_32x32] = { {  1, -1 }, { -1,  1 }, {  2, -1 }, { -1,  2 },
                        { -1, -1 }, {  0, -3 }, { -3,  0 }, { -3, -3 } },
         [BS_32x16] = { {  0, -1 }, { -1,  0 }, {  2, -1 }, { -1, -1 },
                        { -1,  1 }, {  0, -3 }, { -3,  0 }, { -3, -3 } },
         [BS_16x32] = { { -1,  0 }, {  0, -1 }, { -1,  2 }, { -1, -1 },
                        {  1, -1 }, { -3,  0 }, {  0, -3 }, { -3, -3 } },
         [BS_16x16] = { {  0, -1 }, { -1,  0 }, {  1, -1 }, { -1,  1 },
                        { -1, -1 }, {  0, -3 }, { -3,  0 }, { -3, -3 } },
         [BS_16x8]  = { {  0, -1 }, { -1,  0 }, {  1, -1 }, { -1, -1 },
                        {  0, -2 }, { -2,  0 }, { -2, -1 }, { -1, -2 } },
         [BS_8x16]  = { { -1,  0 }, {  0, -1 }, { -1,  1 }, { -1, -1 },
                        { -2,  0 }, {  0, -2 }, { -1, -2 }, { -2, -1 } },
         [BS_8x8]   = { {  0, -1 }, { -1,  0 }, { -1, -1 }, {  0, -2 },
                        { -2,  0 }, { -1, -2 }, { -2, -1 }, { -2, -2 } },
         [BS_8x4]   = { {  0, -1 }, { -1,  0 }, { -1, -1 }, {  0, -2 },
                        { -2,  0 }, { -1, -2 }, { -2, -1 }, { -2, -2 } },
         [BS_4x8]   = { {  0, -1 }, { -1,  0 }, { -1, -1 }, {  0, -2 },
                        { -2,  0 }, { -1, -2 }, { -2, -1 }, { -2, -2 } },
         [BS_4x4]   = { {  0, -1 }, { -1,  0 }, { -1, -1 }, {  0, -2 },
                        { -2,  0 }, { -1, -2 }, { -2, -1 }, { -2, -2 } },
1c9f4b50
     };
e59da0f7
     VP9Context *s = td->s;
     VP9Block *b = td->b;
     int row = td->row, col = td->col, row7 = td->row7;
1c9f4b50
     const int8_t (*p)[2] = mv_ref_blk_off[b->bs];
 #define INVALID_MV 0x80008000U
     uint32_t mem = INVALID_MV, mem_sub8x8 = INVALID_MV;
     int i;
 
875f6955
 #define RETURN_DIRECT_MV(mv)                    \
     do {                                        \
         uint32_t m = AV_RN32A(&mv);             \
         if (!idx) {                             \
             AV_WN32A(pmv, m);                   \
             return;                             \
         } else if (mem == INVALID_MV) {         \
             mem = m;                            \
         } else if (m != mem) {                  \
             AV_WN32A(pmv, m);                   \
             return;                             \
         }                                       \
1c9f4b50
     } while (0)
 
     if (sb >= 0) {
         if (sb == 2 || sb == 1) {
             RETURN_DIRECT_MV(b->mv[0][z]);
         } else if (sb == 3) {
             RETURN_DIRECT_MV(b->mv[2][z]);
             RETURN_DIRECT_MV(b->mv[1][z]);
             RETURN_DIRECT_MV(b->mv[0][z]);
         }
 
875f6955
 #define RETURN_MV(mv)                                                  \
     do {                                                               \
         if (sb > 0) {                                                  \
             VP56mv tmp;                                                \
             uint32_t m;                                                \
             av_assert2(idx == 1);                                      \
             av_assert2(mem != INVALID_MV);                             \
             if (mem_sub8x8 == INVALID_MV) {                            \
e59da0f7
                 clamp_mv(&tmp, &mv, td);                               \
875f6955
                 m = AV_RN32A(&tmp);                                    \
                 if (m != mem) {                                        \
                     AV_WN32A(pmv, m);                                  \
                     return;                                            \
                 }                                                      \
                 mem_sub8x8 = AV_RN32A(&mv);                            \
             } else if (mem_sub8x8 != AV_RN32A(&mv)) {                  \
e59da0f7
                 clamp_mv(&tmp, &mv, td);                               \
875f6955
                 m = AV_RN32A(&tmp);                                    \
                 if (m != mem) {                                        \
                     AV_WN32A(pmv, m);                                  \
                 } else {                                               \
1c9f4b50
                     /* BUG I'm pretty sure this isn't the intention */ \
875f6955
                     AV_WN32A(pmv, 0);                                  \
                 }                                                      \
                 return;                                                \
             }                                                          \
         } else {                                                       \
             uint32_t m = AV_RN32A(&mv);                                \
             if (!idx) {                                                \
e59da0f7
                 clamp_mv(pmv, &mv, td);                                \
875f6955
                 return;                                                \
             } else if (mem == INVALID_MV) {                            \
                 mem = m;                                               \
             } else if (m != mem) {                                     \
e59da0f7
                 clamp_mv(pmv, &mv, td);                                \
875f6955
                 return;                                                \
             }                                                          \
         }                                                              \
1c9f4b50
     } while (0)
 
         if (row > 0) {
37814a21
             VP9mvrefPair *mv = &s->s.frames[CUR_FRAME].mv[(row - 1) * s->sb_cols * 8 + col];
875f6955
             if (mv->ref[0] == ref)
1c9f4b50
                 RETURN_MV(s->above_mv_ctx[2 * col + (sb & 1)][0]);
875f6955
             else if (mv->ref[1] == ref)
1c9f4b50
                 RETURN_MV(s->above_mv_ctx[2 * col + (sb & 1)][1]);
         }
e59da0f7
         if (col > td->tile_col_start) {
37814a21
             VP9mvrefPair *mv = &s->s.frames[CUR_FRAME].mv[row * s->sb_cols * 8 + col - 1];
875f6955
             if (mv->ref[0] == ref)
e59da0f7
                 RETURN_MV(td->left_mv_ctx[2 * row7 + (sb >> 1)][0]);
875f6955
             else if (mv->ref[1] == ref)
e59da0f7
                 RETURN_MV(td->left_mv_ctx[2 * row7 + (sb >> 1)][1]);
1c9f4b50
         }
         i = 2;
     } else {
         i = 0;
     }
 
875f6955
     // previously coded MVs in this neighborhood, using same reference frame
1c9f4b50
     for (; i < 8; i++) {
         int c = p[i][0] + col, r = p[i][1] + row;
 
e59da0f7
         if (c >= td->tile_col_start && c < s->cols &&
875f6955
             r >= 0 && r < s->rows) {
37814a21
             VP9mvrefPair *mv = &s->s.frames[CUR_FRAME].mv[r * s->sb_cols * 8 + c];
1c9f4b50
 
875f6955
             if (mv->ref[0] == ref)
1c9f4b50
                 RETURN_MV(mv->mv[0]);
875f6955
             else if (mv->ref[1] == ref)
1c9f4b50
                 RETURN_MV(mv->mv[1]);
         }
     }
 
     // MV at this position in previous frame, using same reference frame
     if (s->s.h.use_last_frame_mvs) {
37814a21
         VP9mvrefPair *mv = &s->s.frames[REF_FRAME_MVPAIR].mv[row * s->sb_cols * 8 + col];
1c9f4b50
 
         if (!s->s.frames[REF_FRAME_MVPAIR].uses_2pass)
             ff_thread_await_progress(&s->s.frames[REF_FRAME_MVPAIR].tf, row >> 3, 0);
875f6955
         if (mv->ref[0] == ref)
1c9f4b50
             RETURN_MV(mv->mv[0]);
875f6955
         else if (mv->ref[1] == ref)
1c9f4b50
             RETURN_MV(mv->mv[1]);
     }
 
875f6955
 #define RETURN_SCALE_MV(mv, scale)              \
     do {                                        \
         if (scale) {                            \
             VP56mv mv_temp = { -mv.x, -mv.y };  \
             RETURN_MV(mv_temp);                 \
         } else {                                \
             RETURN_MV(mv);                      \
         }                                       \
1c9f4b50
     } while (0)
 
875f6955
     // previously coded MVs in this neighborhood, using different reference frame
1c9f4b50
     for (i = 0; i < 8; i++) {
         int c = p[i][0] + col, r = p[i][1] + row;
 
e59da0f7
         if (c >= td->tile_col_start && c < s->cols && r >= 0 && r < s->rows) {
37814a21
             VP9mvrefPair *mv = &s->s.frames[CUR_FRAME].mv[r * s->sb_cols * 8 + c];
1c9f4b50
 
875f6955
             if (mv->ref[0] != ref && mv->ref[0] >= 0)
                 RETURN_SCALE_MV(mv->mv[0],
                                 s->s.h.signbias[mv->ref[0]] != s->s.h.signbias[ref]);
1c9f4b50
             if (mv->ref[1] != ref && mv->ref[1] >= 0 &&
                 // BUG - libvpx has this condition regardless of whether
                 // we used the first ref MV and pre-scaling
                 AV_RN32A(&mv->mv[0]) != AV_RN32A(&mv->mv[1])) {
                 RETURN_SCALE_MV(mv->mv[1], s->s.h.signbias[mv->ref[1]] != s->s.h.signbias[ref]);
             }
         }
     }
 
     // MV at this position in previous frame, using different reference frame
     if (s->s.h.use_last_frame_mvs) {
37814a21
         VP9mvrefPair *mv = &s->s.frames[REF_FRAME_MVPAIR].mv[row * s->sb_cols * 8 + col];
1c9f4b50
 
         // no need to await_progress, because we already did that above
875f6955
         if (mv->ref[0] != ref && mv->ref[0] >= 0)
1c9f4b50
             RETURN_SCALE_MV(mv->mv[0], s->s.h.signbias[mv->ref[0]] != s->s.h.signbias[ref]);
         if (mv->ref[1] != ref && mv->ref[1] >= 0 &&
             // BUG - libvpx has this condition regardless of whether
             // we used the first ref MV and pre-scaling
             AV_RN32A(&mv->mv[0]) != AV_RN32A(&mv->mv[1])) {
             RETURN_SCALE_MV(mv->mv[1], s->s.h.signbias[mv->ref[1]] != s->s.h.signbias[ref]);
         }
     }
 
     AV_ZERO32(pmv);
e59da0f7
     clamp_mv(pmv, pmv, td);
1c9f4b50
 #undef INVALID_MV
 #undef RETURN_MV
 #undef RETURN_SCALE_MV
 }
 
e59da0f7
 static av_always_inline int read_mv_component(VP9TileData *td, int idx, int hp)
1c9f4b50
 {
e59da0f7
     VP9Context *s = td->s;
     int bit, sign = vp56_rac_get_prob(td->c, s->prob.p.mv_comp[idx].sign);
     int n, c = vp8_rac_get_tree(td->c, ff_vp9_mv_class_tree,
1c9f4b50
                                 s->prob.p.mv_comp[idx].classes);
 
e59da0f7
     td->counts.mv_comp[idx].sign[sign]++;
     td->counts.mv_comp[idx].classes[c]++;
1c9f4b50
     if (c) {
         int m;
 
         for (n = 0, m = 0; m < c; m++) {
e59da0f7
             bit = vp56_rac_get_prob(td->c, s->prob.p.mv_comp[idx].bits[m]);
1c9f4b50
             n |= bit << m;
e59da0f7
             td->counts.mv_comp[idx].bits[m][bit]++;
1c9f4b50
         }
         n <<= 3;
e59da0f7
         bit = vp8_rac_get_tree(td->c, ff_vp9_mv_fp_tree,
875f6955
                                s->prob.p.mv_comp[idx].fp);
         n  |= bit << 1;
e59da0f7
         td->counts.mv_comp[idx].fp[bit]++;
1c9f4b50
         if (hp) {
e59da0f7
             bit = vp56_rac_get_prob(td->c, s->prob.p.mv_comp[idx].hp);
             td->counts.mv_comp[idx].hp[bit]++;
1c9f4b50
             n |= bit;
         } else {
             n |= 1;
             // bug in libvpx - we count for bw entropy purposes even if the
             // bit wasn't coded
e59da0f7
             td->counts.mv_comp[idx].hp[1]++;
1c9f4b50
         }
         n += 8 << c;
     } else {
e59da0f7
         n = vp56_rac_get_prob(td->c, s->prob.p.mv_comp[idx].class0);
         td->counts.mv_comp[idx].class0[n]++;
         bit = vp8_rac_get_tree(td->c, ff_vp9_mv_fp_tree,
1c9f4b50
                                s->prob.p.mv_comp[idx].class0_fp[n]);
e59da0f7
         td->counts.mv_comp[idx].class0_fp[n][bit]++;
1c9f4b50
         n = (n << 3) | (bit << 1);
         if (hp) {
e59da0f7
             bit = vp56_rac_get_prob(td->c, s->prob.p.mv_comp[idx].class0_hp);
             td->counts.mv_comp[idx].class0_hp[bit]++;
1c9f4b50
             n |= bit;
         } else {
             n |= 1;
             // bug in libvpx - we count for bw entropy purposes even if the
             // bit wasn't coded
e59da0f7
             td->counts.mv_comp[idx].class0_hp[1]++;
1c9f4b50
         }
     }
 
     return sign ? -(n + 1) : (n + 1);
 }
 
e59da0f7
 void ff_vp9_fill_mv(VP9TileData *td, VP56mv *mv, int mode, int sb)
1c9f4b50
 {
e59da0f7
     VP9Context *s = td->s;
     VP9Block *b = td->b;
1c9f4b50
 
     if (mode == ZEROMV) {
         AV_ZERO64(mv);
     } else {
         int hp;
 
         // FIXME cache this value and reuse for other subblocks
e59da0f7
         find_ref_mvs(td, &mv[0], b->ref[0], 0, mode == NEARMV,
1c9f4b50
                      mode == NEWMV ? -1 : sb);
         // FIXME maybe move this code into find_ref_mvs()
         if ((mode == NEWMV || sb == -1) &&
875f6955
             !(hp = s->s.h.highprecisionmvs &&
               abs(mv[0].x) < 64 && abs(mv[0].y) < 64)) {
1c9f4b50
             if (mv[0].y & 1) {
                 if (mv[0].y < 0)
                     mv[0].y++;
                 else
                     mv[0].y--;
             }
             if (mv[0].x & 1) {
                 if (mv[0].x < 0)
                     mv[0].x++;
                 else
                     mv[0].x--;
             }
         }
         if (mode == NEWMV) {
e59da0f7
             enum MVJoint j = vp8_rac_get_tree(td->c, ff_vp9_mv_joint_tree,
1c9f4b50
                                               s->prob.p.mv_joint);
 
e59da0f7
             td->counts.mv_joint[j]++;
1c9f4b50
             if (j >= MV_JOINT_V)
e59da0f7
                 mv[0].y += read_mv_component(td, 0, hp);
1c9f4b50
             if (j & 1)
e59da0f7
                 mv[0].x += read_mv_component(td, 1, hp);
1c9f4b50
         }
 
         if (b->comp) {
             // FIXME cache this value and reuse for other subblocks
e59da0f7
             find_ref_mvs(td, &mv[1], b->ref[1], 1, mode == NEARMV,
1c9f4b50
                          mode == NEWMV ? -1 : sb);
             if ((mode == NEWMV || sb == -1) &&
875f6955
                 !(hp = s->s.h.highprecisionmvs &&
                   abs(mv[1].x) < 64 && abs(mv[1].y) < 64)) {
1c9f4b50
                 if (mv[1].y & 1) {
                     if (mv[1].y < 0)
                         mv[1].y++;
                     else
                         mv[1].y--;
                 }
                 if (mv[1].x & 1) {
                     if (mv[1].x < 0)
                         mv[1].x++;
                     else
                         mv[1].x--;
                 }
             }
             if (mode == NEWMV) {
e59da0f7
                 enum MVJoint j = vp8_rac_get_tree(td->c, ff_vp9_mv_joint_tree,
1c9f4b50
                                                   s->prob.p.mv_joint);
 
e59da0f7
                 td->counts.mv_joint[j]++;
1c9f4b50
                 if (j >= MV_JOINT_V)
e59da0f7
                     mv[1].y += read_mv_component(td, 0, hp);
1c9f4b50
                 if (j & 1)
e59da0f7
                     mv[1].x += read_mv_component(td, 1, hp);
1c9f4b50
             }
         }
     }
 }