vp9_findnearmv.c 11.6 KB
Newer Older
John Koleszar's avatar
John Koleszar committed
1
/*
2
 *  Copyright (c) 2010 The WebM project authors. All Rights Reserved.
John Koleszar's avatar
John Koleszar committed
3
 *
4
 *  Use of this source code is governed by a BSD-style license
5
6
 *  that can be found in the LICENSE file in the root of the source
 *  tree. An additional intellectual property rights grant can be found
7
 *  in the file PATENTS.  All contributing project authors may
8
 *  be found in the AUTHORS file in the root of the source tree.
John Koleszar's avatar
John Koleszar committed
9
10
11
 */


12
#include "vp9/common/vp9_findnearmv.h"
13
14
#include "vp9/common/vp9_sadmxn.h"
#include "vp9/common/vp9_subpelvar.h"
15
#include <limits.h>
John Koleszar's avatar
John Koleszar committed
16

17
const uint8_t vp9_mbsplit_offset[4][16] = {
John Koleszar's avatar
John Koleszar committed
18
19
20
21
  { 0,  8,  0,  0,  0,  0,  0,  0,  0,  0,   0,  0,  0,  0,  0,  0},
  { 0,  2,  0,  0,  0,  0,  0,  0,  0,  0,   0,  0,  0,  0,  0,  0},
  { 0,  2,  8, 10,  0,  0,  0,  0,  0,  0,   0,  0,  0,  0,  0,  0},
  { 0,  1,  2,  3,  4,  5,  6,  7,  8,  9,  10, 11, 12, 13, 14, 15}
22
23
};

24
static void lower_mv_precision(int_mv *mv, int usehp)
25
{
26
  if (!usehp || !vp9_use_nmv_hp(&mv->as_mv)) {
27
28
29
30
31
    if (mv->as_mv.row & 1)
      mv->as_mv.row += (mv->as_mv.row > 0 ? -1 : 1);
    if (mv->as_mv.col & 1)
      mv->as_mv.col += (mv->as_mv.col > 0 ? -1 : 1);
  }
32
33
}

34
vp9_prob *vp9_mv_ref_probs(VP9_COMMON *pc,
Paul Wilkins's avatar
Paul Wilkins committed
35
36
37
38
39
40
                           vp9_prob p[4], const int context
                          ) {
  p[0] = pc->fc.vp9_mode_contexts[context][0];
  p[1] = pc->fc.vp9_mode_contexts[context][1];
  p[2] = pc->fc.vp9_mode_contexts[context][2];
  p[3] = pc->fc.vp9_mode_contexts[context][3];
John Koleszar's avatar
John Koleszar committed
41
  return p;
John Koleszar's avatar
John Koleszar committed
42
}
43

44
#define SP(x) (((x) & 7) << 1)
45
unsigned int vp9_sad3x16_c(const uint8_t *src_ptr,
46
                           int  src_stride,
47
                           const uint8_t *ref_ptr,
Yunqing Wang's avatar
Yunqing Wang committed
48
                           int  ref_stride) {
49
50
  return sad_mx_n_c(src_ptr, src_stride, ref_ptr, ref_stride, 3, 16);
}
51
unsigned int vp9_sad16x3_c(const uint8_t *src_ptr,
52
                           int  src_stride,
53
                           const uint8_t *ref_ptr,
Yunqing Wang's avatar
Yunqing Wang committed
54
                           int  ref_stride) {
55
56
57
  return sad_mx_n_c(src_ptr, src_stride, ref_ptr, ref_stride, 16, 3);
}

Yaowu Xu's avatar
Yaowu Xu committed
58

59
unsigned int vp9_variance2x16_c(const uint8_t *src_ptr,
60
                                const int  source_stride,
61
                                const uint8_t *ref_ptr,
62
63
                                const int  recon_stride,
                                unsigned int *sse) {
Yaowu Xu's avatar
Yaowu Xu committed
64
65
66
  int sum;
  variance(src_ptr, source_stride, ref_ptr, recon_stride, 2, 16, sse, &sum);
  return (*sse - (((unsigned int)sum * sum) >> 5));
67
68
}

69
unsigned int vp9_variance16x2_c(const uint8_t *src_ptr,
70
                                const int  source_stride,
71
                                const uint8_t *ref_ptr,
72
73
                                const int  recon_stride,
                                unsigned int *sse) {
Yaowu Xu's avatar
Yaowu Xu committed
74
75
76
  int sum;
  variance(src_ptr, source_stride, ref_ptr, recon_stride, 16, 2, sse, &sum);
  return (*sse - (((unsigned int)sum * sum) >> 5));
77
78
}

79
unsigned int vp9_sub_pixel_variance16x2_c(const uint8_t *src_ptr,
80
81
82
                                          const int  src_pixels_per_line,
                                          const int  xoffset,
                                          const int  yoffset,
83
                                          const uint8_t *dst_ptr,
84
85
                                          const int dst_pixels_per_line,
                                          unsigned int *sse) {
86
87
88
  uint16_t FData3[16 * 3];  // Temp data buffer used in filtering
  uint8_t temp2[2 * 16];
  const int16_t *HFilter, *VFilter;
89

90
91
  HFilter = VP9_BILINEAR_FILTERS_2TAP(xoffset);
  VFilter = VP9_BILINEAR_FILTERS_2TAP(yoffset);
92
93
94
95
96
97
98
99

  var_filter_block2d_bil_first_pass(src_ptr, FData3,
                                    src_pixels_per_line, 1, 3, 16, HFilter);
  var_filter_block2d_bil_second_pass(FData3, temp2, 16, 16, 2, 16, VFilter);

  return vp9_variance16x2_c(temp2, 16, dst_ptr, dst_pixels_per_line, sse);
}

100
unsigned int vp9_sub_pixel_variance2x16_c(const uint8_t *src_ptr,
101
102
103
                                          const int  src_pixels_per_line,
                                          const int  xoffset,
                                          const int  yoffset,
104
                                          const uint8_t *dst_ptr,
105
106
                                          const int dst_pixels_per_line,
                                          unsigned int *sse) {
107
108
109
  uint16_t FData3[2 * 17];  // Temp data buffer used in filtering
  uint8_t temp2[2 * 16];
  const int16_t *HFilter, *VFilter;
110

111
112
  HFilter = VP9_BILINEAR_FILTERS_2TAP(xoffset);
  VFilter = VP9_BILINEAR_FILTERS_2TAP(yoffset);
113
114
115
116
117
118
119
120

  var_filter_block2d_bil_first_pass(src_ptr, FData3,
                                    src_pixels_per_line, 1, 17, 2, HFilter);
  var_filter_block2d_bil_second_pass(FData3, temp2, 2, 2, 16, 2, VFilter);

  return vp9_variance2x16_c(temp2, 2, dst_ptr, dst_pixels_per_line, sse);
}

121
122
123
124
/* check a list of motion vectors by sad score using a number rows of pixels
 * above and a number cols of pixels in the left to select the one with best
 * score to use as ref motion vector
 */
125
void vp9_find_best_ref_mvs(MACROBLOCKD *xd,
126
                           uint8_t *ref_y_buffer,
Paul Wilkins's avatar
Paul Wilkins committed
127
                           int ref_y_stride,
Paul Wilkins's avatar
Paul Wilkins committed
128
                           int_mv *mvlist,
Paul Wilkins's avatar
Paul Wilkins committed
129
130
131
                           int_mv *nearest,
                           int_mv *near) {
  int i, j;
132
133
  uint8_t *above_src;
  uint8_t *above_ref;
134
135
#if !CONFIG_ABOVESPREFMV
  uint8_t *left_src;
136
  uint8_t *left_ref;
137
#endif
Yaowu Xu's avatar
Yaowu Xu committed
138
139
  unsigned int score;
  unsigned int sse;
Paul Wilkins's avatar
Paul Wilkins committed
140
141
  unsigned int ref_scores[MAX_MV_REF_CANDIDATES] = {0};
  int_mv sorted_mvs[MAX_MV_REF_CANDIDATES];
Paul Wilkins's avatar
Paul Wilkins committed
142
143
  int zero_seen = FALSE;

144
  if (ref_y_buffer) {
Paul Wilkins's avatar
Paul Wilkins committed
145

146
147
148
149
150
151
    // Default all to 0,0 if nothing else available
    nearest->as_int = near->as_int = 0;
    vpx_memset(sorted_mvs, 0, sizeof(sorted_mvs));

    above_src = xd->dst.y_buffer - xd->dst.y_stride * 2;
    above_ref = ref_y_buffer - ref_y_stride * 2;
152
#if CONFIG_ABOVESPREFMV
153
154
    above_src -= 4;
    above_ref -= 4;
155
#else
156
157
    left_src  = xd->dst.y_buffer - 2;
    left_ref  = ref_y_buffer - 2;
158
#endif
Paul Wilkins's avatar
Paul Wilkins committed
159

160
161
162
163
164
    // Limit search to the predicted best few candidates
    for (i = 0; i < MAX_MV_REF_CANDIDATES; ++i) {
      int_mv this_mv;
      int offset = 0;
      int row_offset, col_offset;
Paul Wilkins's avatar
Paul Wilkins committed
165

166
      this_mv.as_int = mvlist[i].as_int;
Paul Wilkins's avatar
Paul Wilkins committed
167

168
169
170
171
      // If we see a 0,0 vector for a second time we have reached the end of
      // the list of valid candidate vectors.
      if (!this_mv.as_int && zero_seen)
        break;
Yaowu Xu's avatar
Yaowu Xu committed
172

173
      zero_seen = zero_seen || !this_mv.as_int;
Paul Wilkins's avatar
Paul Wilkins committed
174

175
#if !CONFIG_ABOVESPREFMV
176
177
178
179
180
      clamp_mv(&this_mv,
               xd->mb_to_left_edge - LEFT_TOP_MARGIN + 24,
               xd->mb_to_right_edge + RIGHT_BOTTOM_MARGIN,
               xd->mb_to_top_edge - LEFT_TOP_MARGIN + 24,
               xd->mb_to_bottom_edge + RIGHT_BOTTOM_MARGIN);
181
#else
182
183
184
185
186
      clamp_mv(&this_mv,
               xd->mb_to_left_edge - LEFT_TOP_MARGIN + 32,
               xd->mb_to_right_edge + RIGHT_BOTTOM_MARGIN,
               xd->mb_to_top_edge - LEFT_TOP_MARGIN + 24,
               xd->mb_to_bottom_edge + RIGHT_BOTTOM_MARGIN);
187
#endif
Yaowu Xu's avatar
Yaowu Xu committed
188

189
190
191
192
      row_offset = this_mv.as_mv.row >> 3;
      col_offset = this_mv.as_mv.col >> 3;
      offset = ref_y_stride * row_offset + col_offset;
      score = 0;
193
#if !CONFIG_ABOVESPREFMV
194
      if (xd->up_available) {
195
196
197
#else
      if (xd->up_available && xd->left_available) {
#endif
198
        vp9_sub_pixel_variance16x2(above_ref + offset, ref_y_stride,
199
200
                                   SP(this_mv.as_mv.col),
                                   SP(this_mv.as_mv.row),
201
                                   above_src, xd->dst.y_stride, &sse);
Yaowu Xu's avatar
Yaowu Xu committed
202
        score += sse;
203
204
        if (xd->mode_info_context->mbmi.sb_type >= BLOCK_SIZE_SB32X32) {
          vp9_sub_pixel_variance16x2(above_ref + offset + 16,
Yaowu Xu's avatar
Yaowu Xu committed
205
                                     ref_y_stride,
206
207
                                     SP(this_mv.as_mv.col),
                                     SP(this_mv.as_mv.row),
208
209
210
211
212
                                     above_src + 16, xd->dst.y_stride, &sse);
          score += sse;
        }
        if (xd->mode_info_context->mbmi.sb_type >= BLOCK_SIZE_SB64X64) {
          vp9_sub_pixel_variance16x2(above_ref + offset + 32,
Ronald S. Bultje's avatar
Ronald S. Bultje committed
213
214
215
                                     ref_y_stride,
                                     SP(this_mv.as_mv.col),
                                     SP(this_mv.as_mv.row),
216
217
218
                                     above_src + 32, xd->dst.y_stride, &sse);
          score += sse;
          vp9_sub_pixel_variance16x2(above_ref + offset + 48,
Ronald S. Bultje's avatar
Ronald S. Bultje committed
219
220
221
                                     ref_y_stride,
                                     SP(this_mv.as_mv.col),
                                     SP(this_mv.as_mv.row),
222
223
224
225
226
227
228
229
230
231
                                     above_src + 48, xd->dst.y_stride, &sse);
          score += sse;
        }
      }
#if !CONFIG_ABOVESPREFMV
      if (xd->left_available) {
        vp9_sub_pixel_variance2x16_c(left_ref + offset, ref_y_stride,
                                     SP(this_mv.as_mv.col),
                                     SP(this_mv.as_mv.row),
                                     left_src, xd->dst.y_stride, &sse);
Ronald S. Bultje's avatar
Ronald S. Bultje committed
232
        score += sse;
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
        if (xd->mode_info_context->mbmi.sb_type >= BLOCK_SIZE_SB32X32) {
          vp9_sub_pixel_variance2x16_c(left_ref + offset + ref_y_stride * 16,
                                       ref_y_stride,
                                       SP(this_mv.as_mv.col),
                                       SP(this_mv.as_mv.row),
                                       left_src + xd->dst.y_stride * 16,
                                       xd->dst.y_stride, &sse);
          score += sse;
        }
        if (xd->mode_info_context->mbmi.sb_type >= BLOCK_SIZE_SB64X64) {
          vp9_sub_pixel_variance2x16_c(left_ref + offset + ref_y_stride * 32,
                                     ref_y_stride,
                                       SP(this_mv.as_mv.col),
                                       SP(this_mv.as_mv.row),
                                       left_src + xd->dst.y_stride * 32,
                                       xd->dst.y_stride, &sse);
          score += sse;
          vp9_sub_pixel_variance2x16_c(left_ref + offset + ref_y_stride * 48,
                                       ref_y_stride,
                                       SP(this_mv.as_mv.col),
                                       SP(this_mv.as_mv.row),
                                       left_src + xd->dst.y_stride * 48,
                                       xd->dst.y_stride, &sse);
          score += sse;
        }
Ronald S. Bultje's avatar
Ronald S. Bultje committed
258
      }
259
#endif
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
      // Add the entry to our list and then resort the list on score.
      ref_scores[i] = score;
      sorted_mvs[i].as_int = this_mv.as_int;
      j = i;
      while (j > 0) {
        if (ref_scores[j] < ref_scores[j-1]) {
          ref_scores[j] = ref_scores[j-1];
          sorted_mvs[j].as_int = sorted_mvs[j-1].as_int;
          ref_scores[j-1] = score;
          sorted_mvs[j-1].as_int = this_mv.as_int;
          j--;
        } else {
          break;
        }
      }
Paul Wilkins's avatar
Paul Wilkins committed
275
    }
276
277
  } else {
    vpx_memcpy(sorted_mvs, mvlist, sizeof(sorted_mvs));
Paul Wilkins's avatar
Paul Wilkins committed
278
279
  }

280
  // Make sure all the candidates are properly clamped etc
Paul Wilkins's avatar
Paul Wilkins committed
281
  for (i = 0; i < MAX_MV_REF_CANDIDATES; ++i) {
282
    lower_mv_precision(&sorted_mvs[i], xd->allow_high_precision_mv);
283
    clamp_mv2(&sorted_mvs[i], xd);
284
285
  }

Paul Wilkins's avatar
Paul Wilkins committed
286
287
288
289
290
291
292
293
  // Nearest may be a 0,0 or non zero vector and now matches the chosen
  // "best reference". This has advantages when it is used as part of a
  // compound predictor as it means a non zero vector can be paired using
  // this mode with a 0 vector. The Near vector is still forced to be a
  // non zero candidate if one is avaialble.
  nearest->as_int = sorted_mvs[0].as_int;
  if ( sorted_mvs[1].as_int ) {
    near->as_int = sorted_mvs[1].as_int;
Paul Wilkins's avatar
Paul Wilkins committed
294
  } else {
Paul Wilkins's avatar
Paul Wilkins committed
295
    near->as_int = sorted_mvs[2].as_int;
Paul Wilkins's avatar
Paul Wilkins committed
296
297
  }

Paul Wilkins's avatar
Paul Wilkins committed
298
299
  // Copy back the re-ordered mv list
  vpx_memcpy(mvlist, sorted_mvs, sizeof(sorted_mvs));
Paul Wilkins's avatar
Paul Wilkins committed
300
}