encodemv.c 6.41 KB
Newer Older
Jingning Han's avatar
Jingning Han committed
1
/*
2
 * Copyright (c) 2016, Alliance for Open Media. All rights reserved
Jingning Han's avatar
Jingning Han committed
3
 *
4 5 6 7 8 9
 * This source code is subject to the terms of the BSD 2 Clause License and
 * the Alliance for Open Media Patent License 1.0. If the BSD 2 Clause License
 * was not distributed with this source code in the LICENSE file, you can
 * obtain it at www.aomedia.org/license/software. If the Alliance for Open
 * Media Patent License 1.0 was not distributed with this source code in the
 * PATENTS file, you can obtain it at www.aomedia.org/license/patent.
Jingning Han's avatar
Jingning Han committed
10 11 12 13
 */

#include <math.h>

14 15
#include "av1/common/common.h"
#include "av1/common/entropymode.h"
Jingning Han's avatar
Jingning Han committed
16

17 18 19
#include "av1/encoder/cost.h"
#include "av1/encoder/encodemv.h"
#include "av1/encoder/subexp.h"
Jingning Han's avatar
Jingning Han committed
20

21 22 23 24 25 26 27 28 29 30
#include "aom_dsp/aom_dsp_common.h"

static struct av1_token mv_joint_encodings[MV_JOINTS];
static struct av1_token mv_class_encodings[MV_CLASSES];
static struct av1_token mv_fp_encodings[MV_FP_SIZE];

void av1_entropy_mv_init(void) {
  av1_tokens_from_tree(mv_joint_encodings, av1_mv_joint_tree);
  av1_tokens_from_tree(mv_class_encodings, av1_mv_class_tree);
  av1_tokens_from_tree(mv_fp_encodings, av1_mv_fp_tree);
Jingning Han's avatar
Jingning Han committed
31 32
}

33
static void encode_mv_component(aom_writer *w, int comp, nmv_component *mvcomp,
Alex Converse's avatar
Alex Converse committed
34
                                MvSubpelPrecision precision) {
Jingning Han's avatar
Jingning Han committed
35 36 37
  int offset;
  const int sign = comp < 0;
  const int mag = sign ? -comp : comp;
38
  const int mv_class = av1_get_mv_class(mag - 1, &offset);
39 40 41
  const int d = offset >> 3;         // int mv data
  const int fr = (offset >> 1) & 3;  // fractional mv data
  const int hp = offset & 1;         // high precision mv data
Jingning Han's avatar
Jingning Han committed
42 43 44

  assert(comp != 0);

45
  // Sign
46
  aom_write_symbol(w, sign, mvcomp->sign_cdf, 2);
Jingning Han's avatar
Jingning Han committed
47

48
  // Class
49
  aom_write_symbol(w, mv_class, mvcomp->classes_cdf, MV_CLASSES);
Jingning Han's avatar
Jingning Han committed
50 51 52

  // Integer bits
  if (mv_class == MV_CLASS_0) {
53
    aom_write_symbol(w, d, mvcomp->class0_cdf, CLASS0_SIZE);
Jingning Han's avatar
Jingning Han committed
54 55 56
  } else {
    int i;
    const int n = mv_class + CLASS0_BITS - 1;  // number of bits
57
    for (i = 0; i < n; ++i)
58
      aom_write_symbol(w, (d >> i) & 1, mvcomp->bits_cdf[i], 2);
Jingning Han's avatar
Jingning Han committed
59
  }
Alex Converse's avatar
Alex Converse committed
60
// Fractional bits
61
#if CONFIG_INTRABC || CONFIG_AMVR
Alex Converse's avatar
Alex Converse committed
62
  if (precision > MV_SUBPEL_NONE)
63
#endif  // CONFIG_INTRABC || CONFIG_AMVR
Alex Converse's avatar
Alex Converse committed
64
  {
clang-format's avatar
clang-format committed
65 66 67 68
    aom_write_symbol(
        w, fr,
        mv_class == MV_CLASS_0 ? mvcomp->class0_fp_cdf[d] : mvcomp->fp_cdf,
        MV_FP_SIZE);
Alex Converse's avatar
Alex Converse committed
69
  }
Jingning Han's avatar
Jingning Han committed
70 71

  // High precision bit
Alex Converse's avatar
Alex Converse committed
72
  if (precision > MV_SUBPEL_LOW_PRECISION)
73 74 75
    aom_write_symbol(
        w, hp, mv_class == MV_CLASS_0 ? mvcomp->class0_hp_cdf : mvcomp->hp_cdf,
        2);
Jingning Han's avatar
Jingning Han committed
76 77 78
}

static void build_nmv_component_cost_table(int *mvcost,
79
                                           const nmv_component *const mvcomp,
80
                                           MvSubpelPrecision precision) {
Jingning Han's avatar
Jingning Han committed
81 82 83 84 85 86
  int i, v;
  int sign_cost[2], class_cost[MV_CLASSES], class0_cost[CLASS0_SIZE];
  int bits_cost[MV_OFFSET_BITS][2];
  int class0_fp_cost[CLASS0_SIZE][MV_FP_SIZE], fp_cost[MV_FP_SIZE];
  int class0_hp_cost[2], hp_cost[2];

87 88 89
  av1_cost_tokens_from_cdf(sign_cost, mvcomp->sign_cdf, NULL);
  av1_cost_tokens_from_cdf(class_cost, mvcomp->classes_cdf, NULL);
  av1_cost_tokens_from_cdf(class0_cost, mvcomp->class0_cdf, NULL);
Jingning Han's avatar
Jingning Han committed
90
  for (i = 0; i < MV_OFFSET_BITS; ++i) {
91
    av1_cost_tokens_from_cdf(bits_cost[i], mvcomp->bits_cdf[i], NULL);
Jingning Han's avatar
Jingning Han committed
92 93 94
  }

  for (i = 0; i < CLASS0_SIZE; ++i)
95 96
    av1_cost_tokens_from_cdf(class0_fp_cost[i], mvcomp->class0_fp_cdf[i], NULL);
  av1_cost_tokens_from_cdf(fp_cost, mvcomp->fp_cdf, NULL);
Jingning Han's avatar
Jingning Han committed
97

98
  if (precision > MV_SUBPEL_LOW_PRECISION) {
99 100
    av1_cost_tokens_from_cdf(class0_hp_cost, mvcomp->class0_hp_cdf, NULL);
    av1_cost_tokens_from_cdf(hp_cost, mvcomp->hp_cdf, NULL);
Jingning Han's avatar
Jingning Han committed
101 102 103 104 105
  }
  mvcost[0] = 0;
  for (v = 1; v <= MV_MAX; ++v) {
    int z, c, o, d, e, f, cost = 0;
    z = v - 1;
106
    c = av1_get_mv_class(z, &o);
Jingning Han's avatar
Jingning Han committed
107
    cost += class_cost[c];
108 109 110
    d = (o >> 3);     /* int mv data */
    f = (o >> 1) & 3; /* fractional pel mv data */
    e = (o & 1);      /* high precision mv data */
Jingning Han's avatar
Jingning Han committed
111 112 113
    if (c == MV_CLASS_0) {
      cost += class0_cost[d];
    } else {
114
      const int b = c + CLASS0_BITS - 1; /* number of bits */
115
      for (i = 0; i < b; ++i) cost += bits_cost[i][((d >> i) & 1)];
Jingning Han's avatar
Jingning Han committed
116
    }
117
#if CONFIG_INTRABC || CONFIG_AMVR
118
    if (precision > MV_SUBPEL_NONE)
119
#endif  // CONFIG_INTRABC || CONFIG_AMVR
120
    {
Jingning Han's avatar
Jingning Han committed
121
      if (c == MV_CLASS_0) {
122
        cost += class0_fp_cost[d][f];
Jingning Han's avatar
Jingning Han committed
123
      } else {
124 125 126 127 128 129 130 131
        cost += fp_cost[f];
      }
      if (precision > MV_SUBPEL_LOW_PRECISION) {
        if (c == MV_CLASS_0) {
          cost += class0_hp_cost[e];
        } else {
          cost += hp_cost[e];
        }
Jingning Han's avatar
Jingning Han committed
132 133 134 135 136 137 138
      }
    }
    mvcost[v] = cost + sign_cost[0];
    mvcost[-v] = cost + sign_cost[1];
  }
}

139
void av1_encode_mv(AV1_COMP *cpi, aom_writer *w, const MV *mv, const MV *ref,
140
                   nmv_context *mvctx, int usehp) {
141
  const MV diff = { mv->row - ref->row, mv->col - ref->col };
142
  const MV_JOINT_TYPE j = av1_get_mv_joint(&diff);
143
#if CONFIG_AMVR
144
  if (cpi->common.cur_frame_force_integer_mv) {
145 146 147
    usehp = MV_SUBPEL_NONE;
  }
#endif
148
  aom_write_symbol(w, j, mvctx->joints_cdf, MV_JOINTS);
Jingning Han's avatar
Jingning Han committed
149 150 151 152 153 154 155 156 157
  if (mv_joint_vertical(j))
    encode_mv_component(w, diff.row, &mvctx->comps[0], usehp);

  if (mv_joint_horizontal(j))
    encode_mv_component(w, diff.col, &mvctx->comps[1], usehp);

  // If auto_mv_step_size is enabled then keep track of the largest
  // motion vector component used.
  if (cpi->sf.mv.auto_mv_step_size) {
158 159
    unsigned int maxv = AOMMAX(abs(mv->row), abs(mv->col)) >> 3;
    cpi->max_mv_magnitude = AOMMAX(maxv, cpi->max_mv_magnitude);
Jingning Han's avatar
Jingning Han committed
160 161 162
  }
}

163 164 165
#if CONFIG_INTRABC
void av1_encode_dv(aom_writer *w, const MV *mv, const MV *ref,
                   nmv_context *mvctx) {
166 167 168 169 170
  // DV and ref DV should not have sub-pel.
  assert((mv->col & 7) == 0);
  assert((mv->row & 7) == 0);
  assert((ref->col & 7) == 0);
  assert((ref->row & 7) == 0);
171 172 173
  const MV diff = { mv->row - ref->row, mv->col - ref->col };
  const MV_JOINT_TYPE j = av1_get_mv_joint(&diff);

174
  aom_write_symbol(w, j, mvctx->joints_cdf, MV_JOINTS);
175
  if (mv_joint_vertical(j))
Alex Converse's avatar
Alex Converse committed
176
    encode_mv_component(w, diff.row, &mvctx->comps[0], MV_SUBPEL_NONE);
177 178

  if (mv_joint_horizontal(j))
Alex Converse's avatar
Alex Converse committed
179
    encode_mv_component(w, diff.col, &mvctx->comps[1], MV_SUBPEL_NONE);
180 181 182
}
#endif  // CONFIG_INTRABC

183
void av1_build_nmv_cost_table(int *mvjoint, int *mvcost[2],
184 185
                              const nmv_context *ctx,
                              MvSubpelPrecision precision) {
186
  av1_cost_tokens_from_cdf(mvjoint, ctx->joints_cdf, NULL);
187 188
  build_nmv_component_cost_table(mvcost[0], &ctx->comps[0], precision);
  build_nmv_component_cost_table(mvcost[1], &ctx->comps[1], precision);
Jingning Han's avatar
Jingning Han committed
189
}