hevc_mvs.c 25.9 KB
Newer Older
1
/*
2
 * HEVC video decoder
3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23
 *
 * Copyright (C) 2012 - 2013 Guillaume Martres
 * Copyright (C) 2013 Anand Meher Kotra
 *
 * 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
 */

24
#include "hevc.h"
25
#include "hevcdec.h"
26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41

static const uint8_t l0_l1_cand_idx[12][2] = {
    { 0, 1, },
    { 1, 0, },
    { 0, 2, },
    { 2, 0, },
    { 1, 2, },
    { 2, 1, },
    { 0, 3, },
    { 3, 0, },
    { 1, 3, },
    { 3, 1, },
    { 2, 3, },
    { 3, 2, },
};

42 43
void ff_hevc_set_neighbour_available(HEVCContext *s, int x0, int y0,
                                     int nPbW, int nPbH)
44
{
45
    HEVCLocalContext *lc = s->HEVClc;
46 47
    int x0b = av_mod_uintp2(x0, s->ps.sps->log2_ctb_size);
    int y0b = av_mod_uintp2(y0, s->ps.sps->log2_ctb_size);
48 49 50 51 52

    lc->na.cand_up       = (lc->ctb_up_flag   || y0b);
    lc->na.cand_left     = (lc->ctb_left_flag || x0b);
    lc->na.cand_up_left  = (!x0b && !y0b) ? lc->ctb_up_left_flag : lc->na.cand_left && lc->na.cand_up;
    lc->na.cand_up_right_sap =
53
            ((x0b + nPbW) == (1 << s->ps.sps->log2_ctb_size)) ?
54 55
                    lc->ctb_up_right_flag && !y0b : lc->na.cand_up;
    lc->na.cand_up_right =
56
            lc->na.cand_up_right_sap
57 58 59 60 61 62 63
                     && (x0 + nPbW) < lc->end_of_tiles_x;
    lc->na.cand_bottom_left = ((y0 + nPbH) >= lc->end_of_tiles_y) ? 0 : lc->na.cand_left;
}

/*
 * 6.4.1 Derivation process for z-scan order block availability
 */
Mickaël Raulet's avatar
Mickaël Raulet committed
64
static av_always_inline int z_scan_block_avail(HEVCContext *s, int xCurr, int yCurr,
65 66 67
                              int xN, int yN)
{
#define MIN_TB_ADDR_ZS(x, y)                                            \
68
    s->ps.pps->min_tb_addr_zs[(y) * (s->ps.sps->tb_mask+2) + (x)]
69

70 71 72 73
    int xCurr_ctb = xCurr >> s->ps.sps->log2_ctb_size;
    int yCurr_ctb = yCurr >> s->ps.sps->log2_ctb_size;
    int xN_ctb    = xN    >> s->ps.sps->log2_ctb_size;
    int yN_ctb    = yN    >> s->ps.sps->log2_ctb_size;
74 75 76
    if( yN_ctb < yCurr_ctb || xN_ctb < xCurr_ctb )
        return 1;
    else {
77 78 79 80
        int Curr = MIN_TB_ADDR_ZS((xCurr >> s->ps.sps->log2_min_tb_size) & s->ps.sps->tb_mask,
                (yCurr >> s->ps.sps->log2_min_tb_size) & s->ps.sps->tb_mask);
        int N    = MIN_TB_ADDR_ZS((xN >> s->ps.sps->log2_min_tb_size) & s->ps.sps->tb_mask,
                (yN >> s->ps.sps->log2_min_tb_size) & s->ps.sps->tb_mask);
81 82
        return N <= Curr;
    }
83 84
}

85
//check if the two luma locations belong to the same motion estimation region
Mickaël Raulet's avatar
Mickaël Raulet committed
86
static av_always_inline int is_diff_mer(HEVCContext *s, int xN, int yN, int xP, int yP)
87
{
88
    uint8_t plevel = s->ps.pps->log2_parallel_merge_level;
89 90 91 92 93

    return xN >> plevel == xP >> plevel &&
           yN >> plevel == yP >> plevel;
}

94
#define MATCH_MV(x) (AV_RN32A(&A.x) == AV_RN32A(&B.x))
95 96 97
#define MATCH(x) (A.x == B.x)

// check if the mv's and refidx are the same between A and B
Mickaël Raulet's avatar
Mickaël Raulet committed
98
static av_always_inline int compare_mv_ref_idx(struct MvField A, struct MvField B)
99
{
100 101 102 103
    int a_pf = A.pred_flag;
    int b_pf = B.pred_flag;
    if (a_pf == b_pf) {
        if (a_pf == PF_BI) {
104 105
            return MATCH(ref_idx[0]) && MATCH_MV(mv[0]) &&
                   MATCH(ref_idx[1]) && MATCH_MV(mv[1]);
106
        } else if (a_pf == PF_L0) {
107
            return MATCH(ref_idx[0]) && MATCH_MV(mv[0]);
108
        } else if (a_pf == PF_L1) {
109
            return MATCH(ref_idx[1]) && MATCH_MV(mv[1]);
110 111
        }
    }
112 113 114 115 116 117 118
    return 0;
}

static av_always_inline void mv_scale(Mv *dst, Mv *src, int td, int tb)
{
    int tx, scale_factor;

119 120
    td = av_clip_int8(td);
    tb = av_clip_int8(tb);
121
    tx = (0x4000 + abs(td / 2)) / td;
122
    scale_factor = av_clip_intp2((tb * tx + 32) >> 6, 12);
123 124 125 126
    dst->x = av_clip_int16((scale_factor * src->x + 127 +
                           (scale_factor * src->x < 0)) >> 8);
    dst->y = av_clip_int16((scale_factor * src->y + 127 +
                           (scale_factor * src->y < 0)) >> 8);
127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146
}

static int check_mvset(Mv *mvLXCol, Mv *mvCol,
                       int colPic, int poc,
                       RefPicList *refPicList, int X, int refIdxLx,
                       RefPicList *refPicList_col, int listCol, int refidxCol)
{
    int cur_lt = refPicList[X].isLongTerm[refIdxLx];
    int col_lt = refPicList_col[listCol].isLongTerm[refidxCol];
    int col_poc_diff, cur_poc_diff;

    if (cur_lt != col_lt) {
        mvLXCol->x = 0;
        mvLXCol->y = 0;
        return 0;
    }

    col_poc_diff = colPic - refPicList_col[listCol].list[refidxCol];
    cur_poc_diff = poc    - refPicList[X].list[refIdxLx];

147
    if (cur_lt || col_poc_diff == cur_poc_diff || !col_poc_diff) {
148 149 150 151 152 153 154 155
        mvLXCol->x = mvCol->x;
        mvLXCol->y = mvCol->y;
    } else {
        mv_scale(mvLXCol, mvCol, col_poc_diff, cur_poc_diff);
    }
    return 1;
}

156 157 158 159 160
#define CHECK_MVSET(l)                                          \
    check_mvset(mvLXCol, temp_col.mv + l,                       \
                colPic, s->poc,                                 \
                refPicList, X, refIdxLx,                        \
                refPicList_col, L ## l, temp_col.ref_idx[l])
161 162 163

// derive the motion vectors section 8.5.3.1.8
static int derive_temporal_colocated_mvs(HEVCContext *s, MvField temp_col,
164 165
                                         int refIdxLx, Mv *mvLXCol, int X,
                                         int colPic, RefPicList *refPicList_col)
166 167 168
{
    RefPicList *refPicList = s->ref->refPicList;

169
    if (temp_col.pred_flag == PF_INTRA)
170 171
        return 0;

172
    if (!(temp_col.pred_flag & PF_L0))
173
        return CHECK_MVSET(1);
174
    else if (temp_col.pred_flag == PF_L0)
175
        return CHECK_MVSET(0);
176
    else if (temp_col.pred_flag == PF_BI) {
177
        int check_diffpicount = 0;
178 179 180 181 182 183 184 185
        int i, j;
        for (j = 0; j < 2; j++) {
            for (i = 0; i < refPicList[j].nb_refs; i++) {
                if (refPicList[j].list[i] > s->poc) {
                    check_diffpicount++;
                    break;
                }
            }
186
        }
187 188 189 190 191 192
        if (!check_diffpicount) {
            if (X==0)
                return CHECK_MVSET(0);
            else
                return CHECK_MVSET(1);
        } else {
193 194 195 196 197 198 199 200 201 202
            if (s->sh.collocated_list == L1)
                return CHECK_MVSET(0);
            else
                return CHECK_MVSET(1);
        }
    }

    return 0;
}

203
#define TAB_MVF(x, y)                                                   \
204
    tab_mvf[(y) * min_pu_width + x]
205

206
#define TAB_MVF_PU(v)                                                   \
207 208
    TAB_MVF(((x ## v) >> s->ps.sps->log2_min_pu_size),                     \
            ((y ## v) >> s->ps.sps->log2_min_pu_size))
209

210 211 212
#define DERIVE_TEMPORAL_COLOCATED_MVS                                   \
    derive_temporal_colocated_mvs(s, temp_col,                          \
                                  refIdxLx, mvLXCol, X, colPic,         \
213
                                  ff_hevc_get_ref_list(s, ref, x, y))
214 215 216 217 218 219

/*
 * 8.5.3.1.7  temporal luma motion vector prediction
 */
static int temporal_luma_motion_vector(HEVCContext *s, int x0, int y0,
                                       int nPbW, int nPbH, int refIdxLx,
220
                                       Mv *mvLXCol, int X)
221 222 223
{
    MvField *tab_mvf;
    MvField temp_col;
224
    int x, y, x_pu, y_pu;
225
    int min_pu_width = s->ps.sps->min_pu_width;
226 227 228 229 230
    int availableFlagLXCol = 0;
    int colPic;

    HEVCFrame *ref = s->ref->collocated_ref;

231 232
    if (!ref) {
        memset(mvLXCol, 0, sizeof(*mvLXCol));
233
        return 0;
234
    }
235 236 237 238 239

    tab_mvf = ref->tab_mvf;
    colPic  = ref->poc;

    //bottom right collocated motion vector
240 241
    x = x0 + nPbW;
    y = y0 + nPbH;
242 243

    if (tab_mvf &&
244 245 246
        (y0 >> s->ps.sps->log2_ctb_size) == (y >> s->ps.sps->log2_ctb_size) &&
        y < s->ps.sps->height &&
        x < s->ps.sps->width) {
247 248
        x                 &= ~15;
        y                 &= ~15;
249 250
        if (s->threads_type == FF_THREAD_FRAME)
            ff_thread_await_progress(&ref->tf, y, 0);
251 252
        x_pu               = x >> s->ps.sps->log2_min_pu_size;
        y_pu               = y >> s->ps.sps->log2_min_pu_size;
253
        temp_col           = TAB_MVF(x_pu, y_pu);
254
        availableFlagLXCol = DERIVE_TEMPORAL_COLOCATED_MVS;
255 256 257
    }

    // derive center collocated motion vector
258
    if (tab_mvf && !availableFlagLXCol) {
259 260
        x                  = x0 + (nPbW >> 1);
        y                  = y0 + (nPbH >> 1);
261 262
        x                 &= ~15;
        y                 &= ~15;
263 264
        if (s->threads_type == FF_THREAD_FRAME)
            ff_thread_await_progress(&ref->tf, y, 0);
265 266
        x_pu               = x >> s->ps.sps->log2_min_pu_size;
        y_pu               = y >> s->ps.sps->log2_min_pu_size;
267
        temp_col           = TAB_MVF(x_pu, y_pu);
268
        availableFlagLXCol = DERIVE_TEMPORAL_COLOCATED_MVS;
269 270 271 272
    }
    return availableFlagLXCol;
}

273
#define AVAILABLE(cand, v)                                      \
274
    (cand && !(TAB_MVF_PU(v).pred_flag == PF_INTRA))
275

276
#define PRED_BLOCK_AVAILABLE(v)                                 \
277
    z_scan_block_avail(s, x0, y0, x ## v, y ## v)
278

279
#define COMPARE_MV_REFIDX(a, b)                                 \
Mickaël Raulet's avatar
Mickaël Raulet committed
280
    compare_mv_ref_idx(TAB_MVF_PU(a), TAB_MVF_PU(b))
281 282 283 284 285

/*
 * 8.5.3.1.2  Derivation process for spatial merging candidates
 */
static void derive_spatial_merge_candidates(HEVCContext *s, int x0, int y0,
286 287
                                            int nPbW, int nPbH,
                                            int log2_cb_size,
288
                                            int singleMCLFlag, int part_idx,
289
                                            int merge_idx,
290 291
                                            struct MvField mergecandlist[])
{
292
    HEVCLocalContext *lc   = s->HEVClc;
293
    RefPicList *refPicList = s->ref->refPicList;
294 295
    MvField *tab_mvf       = s->ref->tab_mvf;

296
    const int min_pu_width = s->ps.sps->min_pu_width;
297 298 299 300 301 302 303 304 305 306 307 308 309 310 311 312 313 314 315 316 317 318

    const int cand_bottom_left = lc->na.cand_bottom_left;
    const int cand_left        = lc->na.cand_left;
    const int cand_up_left     = lc->na.cand_up_left;
    const int cand_up          = lc->na.cand_up;
    const int cand_up_right    = lc->na.cand_up_right_sap;

    const int xA1    = x0 - 1;
    const int yA1    = y0 + nPbH - 1;

    const int xB1    = x0 + nPbW - 1;
    const int yB1    = y0 - 1;

    const int xB0    = x0 + nPbW;
    const int yB0    = y0 - 1;

    const int xA0    = x0 - 1;
    const int yA0    = y0 + nPbH;

    const int xB2    = x0 - 1;
    const int yB2    = y0 - 1;

319
    const int nb_refs = (s->sh.slice_type == HEVC_SLICE_P) ?
320
                        s->sh.nb_refs[0] : FFMIN(s->sh.nb_refs[0], s->sh.nb_refs[1]);
321

322
    int zero_idx = 0;
323

324 325
    int nb_merge_cand = 0;
    int nb_orig_merge_cand = 0;
326

327 328 329
    int is_available_a0;
    int is_available_a1;
    int is_available_b0;
330
    int is_available_b1;
331
    int is_available_b2;
332 333 334 335 336 337


    if (!singleMCLFlag && part_idx == 1 &&
        (lc->cu.part_mode == PART_Nx2N ||
         lc->cu.part_mode == PART_nLx2N ||
         lc->cu.part_mode == PART_nRx2N) ||
Mickaël Raulet's avatar
Mickaël Raulet committed
338
        is_diff_mer(s, xA1, yA1, x0, y0)) {
339
        is_available_a1 = 0;
340 341 342 343
    } else {
        is_available_a1 = AVAILABLE(cand_left, A1);
        if (is_available_a1) {
            mergecandlist[nb_merge_cand] = TAB_MVF_PU(A1);
344 345
            if (merge_idx == 0)
                return;
346 347
            nb_merge_cand++;
        }
348 349 350 351 352 353
    }

    if (!singleMCLFlag && part_idx == 1 &&
        (lc->cu.part_mode == PART_2NxN ||
         lc->cu.part_mode == PART_2NxnU ||
         lc->cu.part_mode == PART_2NxnD) ||
Mickaël Raulet's avatar
Mickaël Raulet committed
354
        is_diff_mer(s, xB1, yB1, x0, y0)) {
355
        is_available_b1 = 0;
356 357 358 359 360
    } else {
        is_available_b1 = AVAILABLE(cand_up, B1);
        if (is_available_b1 &&
            !(is_available_a1 && COMPARE_MV_REFIDX(B1, A1))) {
            mergecandlist[nb_merge_cand] = TAB_MVF_PU(B1);
361 362
            if (merge_idx == nb_merge_cand)
                return;
363 364
            nb_merge_cand++;
        }
365 366 367
    }

    // above right spatial merge candidate
368
    is_available_b0 = AVAILABLE(cand_up_right, B0) &&
369
                      xB0 < s->ps.sps->width &&
370
                      PRED_BLOCK_AVAILABLE(B0) &&
Mickaël Raulet's avatar
Mickaël Raulet committed
371
                      !is_diff_mer(s, xB0, yB0, x0, y0);
372

373 374
    if (is_available_b0 &&
        !(is_available_b1 && COMPARE_MV_REFIDX(B0, B1))) {
375
        mergecandlist[nb_merge_cand] = TAB_MVF_PU(B0);
376 377
        if (merge_idx == nb_merge_cand)
            return;
378 379
        nb_merge_cand++;
    }
380 381

    // left bottom spatial merge candidate
382
    is_available_a0 = AVAILABLE(cand_bottom_left, A0) &&
383
                      yA0 < s->ps.sps->height &&
384
                      PRED_BLOCK_AVAILABLE(A0) &&
Mickaël Raulet's avatar
Mickaël Raulet committed
385
                      !is_diff_mer(s, xA0, yA0, x0, y0);
386

387 388
    if (is_available_a0 &&
        !(is_available_a1 && COMPARE_MV_REFIDX(A0, A1))) {
389
        mergecandlist[nb_merge_cand] = TAB_MVF_PU(A0);
390 391
        if (merge_idx == nb_merge_cand)
            return;
392 393
        nb_merge_cand++;
    }
394 395

    // above left spatial merge candidate
396
    is_available_b2 = AVAILABLE(cand_up_left, B2) &&
Mickaël Raulet's avatar
Mickaël Raulet committed
397
                      !is_diff_mer(s, xB2, yB2, x0, y0);
398

399 400 401 402
    if (is_available_b2 &&
        !(is_available_a1 && COMPARE_MV_REFIDX(B2, A1)) &&
        !(is_available_b1 && COMPARE_MV_REFIDX(B2, B1)) &&
        nb_merge_cand != 4) {
403
        mergecandlist[nb_merge_cand] = TAB_MVF_PU(B2);
404 405
        if (merge_idx == nb_merge_cand)
            return;
406 407
        nb_merge_cand++;
    }
408 409

    // temporal motion vector candidate
410 411
    if (s->sh.slice_temporal_mvp_enabled_flag &&
        nb_merge_cand < s->sh.max_num_merge_cand) {
412
        Mv mv_l0_col = { 0 }, mv_l1_col = { 0 };
413 414
        int available_l0 = temporal_luma_motion_vector(s, x0, y0, nPbW, nPbH,
                                                       0, &mv_l0_col, 0);
415
        int available_l1 = (s->sh.slice_type == HEVC_SLICE_B) ?
416 417 418 419
                           temporal_luma_motion_vector(s, x0, y0, nPbW, nPbH,
                                                       0, &mv_l1_col, 1) : 0;

        if (available_l0 || available_l1) {
420
            mergecandlist[nb_merge_cand].pred_flag = available_l0 + (available_l1 << 1);
421 422 423 424
            AV_ZERO16(mergecandlist[nb_merge_cand].ref_idx);
            mergecandlist[nb_merge_cand].mv[0]      = mv_l0_col;
            mergecandlist[nb_merge_cand].mv[1]      = mv_l1_col;

425 426
            if (merge_idx == nb_merge_cand)
                return;
427
            nb_merge_cand++;
428 429 430
        }
    }

431
    nb_orig_merge_cand = nb_merge_cand;
432 433

    // combined bi-predictive merge candidates  (applies for B slices)
434
    if (s->sh.slice_type == HEVC_SLICE_B && nb_orig_merge_cand > 1 &&
435 436 437 438 439 440 441 442 443 444
        nb_orig_merge_cand < s->sh.max_num_merge_cand) {
        int comb_idx = 0;

        for (comb_idx = 0; nb_merge_cand < s->sh.max_num_merge_cand &&
                           comb_idx < nb_orig_merge_cand * (nb_orig_merge_cand - 1); comb_idx++) {
            int l0_cand_idx = l0_l1_cand_idx[comb_idx][0];
            int l1_cand_idx = l0_l1_cand_idx[comb_idx][1];
            MvField l0_cand = mergecandlist[l0_cand_idx];
            MvField l1_cand = mergecandlist[l1_cand_idx];

445
            if ((l0_cand.pred_flag & PF_L0) && (l1_cand.pred_flag & PF_L1) &&
446 447
                (refPicList[0].list[l0_cand.ref_idx[0]] !=
                 refPicList[1].list[l1_cand.ref_idx[1]] ||
448
                 AV_RN32A(&l0_cand.mv[0]) != AV_RN32A(&l1_cand.mv[1]))) {
449 450
                mergecandlist[nb_merge_cand].ref_idx[0]   = l0_cand.ref_idx[0];
                mergecandlist[nb_merge_cand].ref_idx[1]   = l1_cand.ref_idx[1];
451
                mergecandlist[nb_merge_cand].pred_flag    = PF_BI;
452 453
                AV_COPY32(&mergecandlist[nb_merge_cand].mv[0], &l0_cand.mv[0]);
                AV_COPY32(&mergecandlist[nb_merge_cand].mv[1], &l1_cand.mv[1]);
454 455
                if (merge_idx == nb_merge_cand)
                    return;
456
                nb_merge_cand++;
457 458 459 460
            }
        }
    }

461 462
    // append Zero motion vector candidates
    while (nb_merge_cand < s->sh.max_num_merge_cand) {
James Almer's avatar
James Almer committed
463
        mergecandlist[nb_merge_cand].pred_flag    = PF_L0 + ((s->sh.slice_type == HEVC_SLICE_B) << 1);
464 465
        AV_ZERO32(mergecandlist[nb_merge_cand].mv + 0);
        AV_ZERO32(mergecandlist[nb_merge_cand].mv + 1);
466 467
        mergecandlist[nb_merge_cand].ref_idx[0]   = zero_idx < nb_refs ? zero_idx : 0;
        mergecandlist[nb_merge_cand].ref_idx[1]   = zero_idx < nb_refs ? zero_idx : 0;
468

469 470
        if (merge_idx == nb_merge_cand)
            return;
471 472
        nb_merge_cand++;
        zero_idx++;
473 474 475 476 477 478 479 480 481 482 483 484
    }
}

/*
 * 8.5.3.1.1 Derivation process of luma Mvs for merge mode
 */
void ff_hevc_luma_mv_merge_mode(HEVCContext *s, int x0, int y0, int nPbW,
                                int nPbH, int log2_cb_size, int part_idx,
                                int merge_idx, MvField *mv)
{
    int singleMCLFlag = 0;
    int nCS = 1 << log2_cb_size;
485
    MvField mergecand_list[MRG_MAX_NUM_CANDS];
486 487
    int nPbW2 = nPbW;
    int nPbH2 = nPbH;
488
    HEVCLocalContext *lc = s->HEVClc;
489

490
    if (s->ps.pps->log2_parallel_merge_level > 2 && nCS == 8) {
491
        singleMCLFlag = 1;
492 493 494 495 496
        x0            = lc->cu.x;
        y0            = lc->cu.y;
        nPbW          = nCS;
        nPbH          = nCS;
        part_idx      = 0;
497 498 499 500
    }

    ff_hevc_set_neighbour_available(s, x0, y0, nPbW, nPbH);
    derive_spatial_merge_candidates(s, x0, y0, nPbW, nPbH, log2_cb_size,
501 502
                                    singleMCLFlag, part_idx,
                                    merge_idx, mergecand_list);
503

504
    if (mergecand_list[merge_idx].pred_flag == PF_BI &&
505
        (nPbW2 + nPbH2) == 12) {
506
        mergecand_list[merge_idx].pred_flag = PF_L0;
507 508 509 510 511
    }

    *mv = mergecand_list[merge_idx];
}

512
static av_always_inline void dist_scale(HEVCContext *s, Mv *mv,
513
                                        int min_pu_width, int x, int y,
514 515 516
                                        int elist, int ref_idx_curr, int ref_idx)
{
    RefPicList *refPicList = s->ref->refPicList;
517 518 519
    MvField *tab_mvf       = s->ref->tab_mvf;
    int ref_pic_elist      = refPicList[elist].list[TAB_MVF(x, y).ref_idx[elist]];
    int ref_pic_curr       = refPicList[ref_idx_curr].list[ref_idx];
520

Luca Barbato's avatar
Luca Barbato committed
521 522 523 524 525 526
    if (ref_pic_elist != ref_pic_curr) {
        int poc_diff = s->poc - ref_pic_elist;
        if (!poc_diff)
            poc_diff = 1;
        mv_scale(mv, mv, poc_diff, s->poc - ref_pic_curr);
    }
527 528 529 530 531 532
}

static int mv_mp_mode_mx(HEVCContext *s, int x, int y, int pred_flag_index,
                         Mv *mv, int ref_idx_curr, int ref_idx)
{
    MvField *tab_mvf = s->ref->tab_mvf;
533
    int min_pu_width = s->ps.sps->min_pu_width;
534 535 536

    RefPicList *refPicList = s->ref->refPicList;

537
    if (((TAB_MVF(x, y).pred_flag) & (1 << pred_flag_index)) &&
538 539 540 541 542 543 544 545 546 547 548
        refPicList[pred_flag_index].list[TAB_MVF(x, y).ref_idx[pred_flag_index]] == refPicList[ref_idx_curr].list[ref_idx]) {
        *mv = TAB_MVF(x, y).mv[pred_flag_index];
        return 1;
    }
    return 0;
}

static int mv_mp_mode_mx_lt(HEVCContext *s, int x, int y, int pred_flag_index,
                            Mv *mv, int ref_idx_curr, int ref_idx)
{
    MvField *tab_mvf = s->ref->tab_mvf;
549
    int min_pu_width = s->ps.sps->min_pu_width;
550 551 552

    RefPicList *refPicList = s->ref->refPicList;

553 554
    if ((TAB_MVF(x, y).pred_flag) & (1 << pred_flag_index)) {
        int currIsLongTerm     = refPicList[ref_idx_curr].isLongTerm[ref_idx];
555

556 557 558 559 560 561 562 563 564 565
        int colIsLongTerm =
            refPicList[pred_flag_index].isLongTerm[(TAB_MVF(x, y).ref_idx[pred_flag_index])];

        if (colIsLongTerm == currIsLongTerm) {
            *mv = TAB_MVF(x, y).mv[pred_flag_index];
            if (!currIsLongTerm)
                dist_scale(s, mv, min_pu_width, x, y,
                           pred_flag_index, ref_idx_curr, ref_idx);
            return 1;
        }
566 567 568 569
    }
    return 0;
}

570
#define MP_MX(v, pred, mx)                                      \
571
    mv_mp_mode_mx(s,                                            \
572 573
                  (x ## v) >> s->ps.sps->log2_min_pu_size,         \
                  (y ## v) >> s->ps.sps->log2_min_pu_size,         \
574
                  pred, &mx, ref_idx_curr, ref_idx)
575

576
#define MP_MX_LT(v, pred, mx)                                   \
577
    mv_mp_mode_mx_lt(s,                                         \
578 579
                     (x ## v) >> s->ps.sps->log2_min_pu_size,      \
                     (y ## v) >> s->ps.sps->log2_min_pu_size,      \
580
                     pred, &mx, ref_idx_curr, ref_idx)
581 582 583 584 585 586

void ff_hevc_luma_mv_mvp_mode(HEVCContext *s, int x0, int y0, int nPbW,
                              int nPbH, int log2_cb_size, int part_idx,
                              int merge_idx, MvField *mv,
                              int mvp_lx_flag, int LX)
{
587
    HEVCLocalContext *lc = s->HEVClc;
588 589
    MvField *tab_mvf = s->ref->tab_mvf;
    int isScaledFlag_L0 = 0;
590 591
    int availableFlagLXA0 = 1;
    int availableFlagLXB0 = 1;
592
    int numMVPCandLX = 0;
593
    int min_pu_width = s->ps.sps->min_pu_width;
594 595 596 597 598 599 600 601

    int xA0, yA0;
    int is_available_a0;
    int xA1, yA1;
    int is_available_a1;
    int xB0, yB0;
    int is_available_b0;
    int xB1, yB1;
602
    int is_available_b1;
603
    int xB2, yB2;
604 605
    int is_available_b2;

606
    Mv mvpcand_list[2] = { { 0 } };
607 608
    Mv mxA;
    Mv mxB;
609
    int ref_idx_curr;
610 611 612 613
    int ref_idx = 0;
    int pred_flag_index_l0;
    int pred_flag_index_l1;

614 615 616 617 618
    const int cand_bottom_left = lc->na.cand_bottom_left;
    const int cand_left        = lc->na.cand_left;
    const int cand_up_left     = lc->na.cand_up_left;
    const int cand_up          = lc->na.cand_up;
    const int cand_up_right    = lc->na.cand_up_right_sap;
619 620 621 622 623 624 625 626 627
    ref_idx_curr       = LX;
    ref_idx            = mv->ref_idx[LX];
    pred_flag_index_l0 = LX;
    pred_flag_index_l1 = !LX;

    // left bottom spatial candidate
    xA0 = x0 - 1;
    yA0 = y0 + nPbH;

628
    is_available_a0 = AVAILABLE(cand_bottom_left, A0) &&
629
                      yA0 < s->ps.sps->height &&
630
                      PRED_BLOCK_AVAILABLE(A0);
631 632

    //left spatial merge candidate
633 634
    xA1    = x0 - 1;
    yA1    = y0 + nPbH - 1;
635 636

    is_available_a1 = AVAILABLE(cand_left, A1);
637
    if (is_available_a0 || is_available_a1)
638 639 640
        isScaledFlag_L0 = 1;

    if (is_available_a0) {
641 642 643 644 645 646
        if (MP_MX(A0, pred_flag_index_l0, mxA)) {
            goto b_candidates;
        }
        if (MP_MX(A0, pred_flag_index_l1, mxA)) {
            goto b_candidates;
        }
647 648
    }

649 650 651 652 653 654 655
    if (is_available_a1) {
        if (MP_MX(A1, pred_flag_index_l0, mxA)) {
            goto b_candidates;
        }
        if (MP_MX(A1, pred_flag_index_l1, mxA)) {
            goto b_candidates;
        }
656 657
    }

658 659 660 661 662 663 664
    if (is_available_a0) {
        if (MP_MX_LT(A0, pred_flag_index_l0, mxA)) {
            goto b_candidates;
        }
        if (MP_MX_LT(A0, pred_flag_index_l1, mxA)) {
            goto b_candidates;
        }
665 666
    }

667 668 669 670 671 672 673
    if (is_available_a1) {
        if (MP_MX_LT(A1, pred_flag_index_l0, mxA)) {
            goto b_candidates;
        }
        if (MP_MX_LT(A1, pred_flag_index_l1, mxA)) {
            goto b_candidates;
        }
674
    }
675
    availableFlagLXA0 = 0;
676

677
b_candidates:
678 679
    // B candidates
    // above right spatial merge candidate
680 681
    xB0    = x0 + nPbW;
    yB0    = y0 - 1;
682

683
    is_available_b0 =  AVAILABLE(cand_up_right, B0) &&
684
                       xB0 < s->ps.sps->width &&
685
                       PRED_BLOCK_AVAILABLE(B0);
686

687 688 689 690 691 692 693 694 695 696 697
    // above spatial merge candidate
    xB1    = x0 + nPbW - 1;
    yB1    = y0 - 1;
    is_available_b1 = AVAILABLE(cand_up, B1);

    // above left spatial merge candidate
    xB2 = x0 - 1;
    yB2 = y0 - 1;
    is_available_b2 = AVAILABLE(cand_up_left, B2);

    // above right spatial merge candidate
698
    if (is_available_b0) {
699 700 701 702 703 704
        if (MP_MX(B0, pred_flag_index_l0, mxB)) {
            goto scalef;
        }
        if (MP_MX(B0, pred_flag_index_l1, mxB)) {
            goto scalef;
        }
705 706
    }

707 708 709 710 711 712 713
    // above spatial merge candidate
    if (is_available_b1) {
        if (MP_MX(B1, pred_flag_index_l0, mxB)) {
            goto scalef;
        }
        if (MP_MX(B1, pred_flag_index_l1, mxB)) {
            goto scalef;
714 715 716
        }
    }

717 718 719 720 721 722 723
    // above left spatial merge candidate
    if (is_available_b2) {
        if (MP_MX(B2, pred_flag_index_l0, mxB)) {
            goto scalef;
        }
        if (MP_MX(B2, pred_flag_index_l1, mxB)) {
            goto scalef;
724 725
        }
    }
726
    availableFlagLXB0 = 0;
727

728 729
scalef:
    if (!isScaledFlag_L0) {
730 731 732 733 734 735 736 737 738 739 740 741 742 743 744 745 746 747 748 749 750 751 752 753 754 755
        if (availableFlagLXB0) {
            availableFlagLXA0 = 1;
            mxA = mxB;
        }
        availableFlagLXB0 = 0;

        // XB0 and L1
        if (is_available_b0) {
            availableFlagLXB0 = MP_MX_LT(B0, pred_flag_index_l0, mxB);
            if (!availableFlagLXB0)
                availableFlagLXB0 = MP_MX_LT(B0, pred_flag_index_l1, mxB);
        }

        if (is_available_b1 && !availableFlagLXB0) {
            availableFlagLXB0 = MP_MX_LT(B1, pred_flag_index_l0, mxB);
            if (!availableFlagLXB0)
                availableFlagLXB0 = MP_MX_LT(B1, pred_flag_index_l1, mxB);
        }

        if (is_available_b2 && !availableFlagLXB0) {
            availableFlagLXB0 = MP_MX_LT(B2, pred_flag_index_l0, mxB);
            if (!availableFlagLXB0)
                availableFlagLXB0 = MP_MX_LT(B2, pred_flag_index_l1, mxB);
        }
    }

756 757
    if (availableFlagLXA0)
        mvpcand_list[numMVPCandLX++] = mxA;
758

759 760
    if (availableFlagLXB0 && (!availableFlagLXA0 || mxA.x != mxB.x || mxA.y != mxB.y))
        mvpcand_list[numMVPCandLX++] = mxB;
761

762
    //temporal motion vector prediction candidate
763 764
    if (numMVPCandLX < 2 && s->sh.slice_temporal_mvp_enabled_flag &&
        mvp_lx_flag == numMVPCandLX) {
765 766
        Mv mv_col;
        int available_col = temporal_luma_motion_vector(s, x0, y0, nPbW,
767 768
                                                        nPbH, ref_idx,
                                                        &mv_col, LX);
769 770
        if (available_col)
            mvpcand_list[numMVPCandLX++] = mv_col;
771 772
    }

773
    mv->mv[LX] = mvpcand_list[mvp_lx_flag];
774
}