hevc/mvs: simplifying derive_spatial_merge_candidates function
[ffmpeg.git] / libavcodec / hevc_mvs.c
1 /*
2  * HEVC video decoder
3  *
4  * Copyright (C) 2012 - 2013 Guillaume Martres
5  * Copyright (C) 2013 Anand Meher Kotra
6  *
7  * This file is part of FFmpeg.
8  *
9  * FFmpeg is free software; you can redistribute it and/or
10  * modify it under the terms of the GNU Lesser General Public
11  * License as published by the Free Software Foundation; either
12  * version 2.1 of the License, or (at your option) any later version.
13  *
14  * FFmpeg is distributed in the hope that it will be useful,
15  * but WITHOUT ANY WARRANTY; without even the implied warranty of
16  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
17  * Lesser General Public License for more details.
18  *
19  * You should have received a copy of the GNU Lesser General Public
20  * License along with FFmpeg; if not, write to the Free Software
21  * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
22  */
23
24 #include "hevc.h"
25
26 static const uint8_t l0_l1_cand_idx[12][2] = {
27     { 0, 1, },
28     { 1, 0, },
29     { 0, 2, },
30     { 2, 0, },
31     { 1, 2, },
32     { 2, 1, },
33     { 0, 3, },
34     { 3, 0, },
35     { 1, 3, },
36     { 3, 1, },
37     { 2, 3, },
38     { 3, 2, },
39 };
40
41 void ff_hevc_set_neighbour_available(HEVCContext *s, int x0, int y0,
42                                      int nPbW, int nPbH)
43 {
44     HEVCLocalContext *lc = s->HEVClc;
45     int x0b = x0 & ((1 << s->sps->log2_ctb_size) - 1);
46     int y0b = y0 & ((1 << s->sps->log2_ctb_size) - 1);
47
48     lc->na.cand_up       = (lc->ctb_up_flag   || y0b);
49     lc->na.cand_left     = (lc->ctb_left_flag || x0b);
50     lc->na.cand_up_left  = (!x0b && !y0b) ? lc->ctb_up_left_flag : lc->na.cand_left && lc->na.cand_up;
51     lc->na.cand_up_right_sap =
52             ((x0b + nPbW) == (1 << s->sps->log2_ctb_size)) ?
53                     lc->ctb_up_right_flag && !y0b : lc->na.cand_up;
54     lc->na.cand_up_right =
55             ((x0b + nPbW) == (1 << s->sps->log2_ctb_size) ?
56                     lc->ctb_up_right_flag && !y0b : lc->na.cand_up )
57                      && (x0 + nPbW) < lc->end_of_tiles_x;
58     lc->na.cand_bottom_left = ((y0 + nPbH) >= lc->end_of_tiles_y) ? 0 : lc->na.cand_left;
59 }
60
61 /*
62  * 6.4.1 Derivation process for z-scan order block availability
63  */
64 static int z_scan_block_avail(HEVCContext *s, int xCurr, int yCurr,
65                               int xN, int yN)
66 {
67 #define MIN_TB_ADDR_ZS(x, y)                                            \
68     s->pps->min_tb_addr_zs[(y) * (s->sps->tb_mask+2) + (x)]
69
70     int xCurr_ctb = xCurr >> s->sps->log2_ctb_size;
71     int yCurr_ctb = yCurr >> s->sps->log2_ctb_size;
72     int xN_ctb    = xN    >> s->sps->log2_ctb_size;
73     int yN_ctb    = yN    >> s->sps->log2_ctb_size;
74
75     if (xN < 0 || yN < 0 ||
76         xN >= s->sps->width ||
77         yN >= s->sps->height)
78         return 0;
79
80     if( yN_ctb < yCurr_ctb || xN_ctb < xCurr_ctb )
81         return 1;
82     else {
83         int Curr = MIN_TB_ADDR_ZS((xCurr >> s->sps->log2_min_tb_size) & s->sps->tb_mask,
84                 (yCurr >> s->sps->log2_min_tb_size) & s->sps->tb_mask);
85         int N    = MIN_TB_ADDR_ZS((xN >> s->sps->log2_min_tb_size) & s->sps->tb_mask,
86                 (yN >> s->sps->log2_min_tb_size) & s->sps->tb_mask);
87         return N <= Curr;
88     }
89 }
90
91 //check if the two luma locations belong to the same mostion estimation region
92 static int isDiffMER(HEVCContext *s, int xN, int yN, int xP, int yP)
93 {
94     uint8_t plevel = s->pps->log2_parallel_merge_level;
95
96     return xN >> plevel == xP >> plevel &&
97            yN >> plevel == yP >> plevel;
98 }
99
100 #define MATCH_MV(x) (AV_RN32A(&A.x) == AV_RN32A(&B.x))
101 #define MATCH(x) (A.x == B.x)
102
103 // check if the mv's and refidx are the same between A and B
104 static int compareMVrefidx(struct MvField A, struct MvField B)
105 {
106     int a_pf = A.pred_flag;
107     int b_pf = B.pred_flag;
108     if (a_pf == b_pf) {
109         if (a_pf == PF_BI) {
110             return MATCH(ref_idx[0]) && MATCH_MV(mv[0]) &&
111                    MATCH(ref_idx[1]) && MATCH_MV(mv[1]);
112         } else if (a_pf == PF_L0) {
113             return MATCH(ref_idx[0]) && MATCH_MV(mv[0]);
114         } else if (a_pf == PF_L1) {
115             return MATCH(ref_idx[1]) && MATCH_MV(mv[1]);
116         }
117     }
118     return 0;
119 }
120
121 static av_always_inline void mv_scale(Mv *dst, Mv *src, int td, int tb)
122 {
123     int tx, scale_factor;
124
125     td = av_clip_int8(td);
126     tb = av_clip_int8(tb);
127     tx = (0x4000 + abs(td / 2)) / td;
128     scale_factor = av_clip((tb * tx + 32) >> 6, -4096, 4095);
129     dst->x = av_clip_int16((scale_factor * src->x + 127 +
130                            (scale_factor * src->x < 0)) >> 8);
131     dst->y = av_clip_int16((scale_factor * src->y + 127 +
132                            (scale_factor * src->y < 0)) >> 8);
133 }
134
135 static int check_mvset(Mv *mvLXCol, Mv *mvCol,
136                        int colPic, int poc,
137                        RefPicList *refPicList, int X, int refIdxLx,
138                        RefPicList *refPicList_col, int listCol, int refidxCol)
139 {
140     int cur_lt = refPicList[X].isLongTerm[refIdxLx];
141     int col_lt = refPicList_col[listCol].isLongTerm[refidxCol];
142     int col_poc_diff, cur_poc_diff;
143
144     if (cur_lt != col_lt) {
145         mvLXCol->x = 0;
146         mvLXCol->y = 0;
147         return 0;
148     }
149
150     col_poc_diff = colPic - refPicList_col[listCol].list[refidxCol];
151     cur_poc_diff = poc    - refPicList[X].list[refIdxLx];
152
153     if (cur_lt || col_poc_diff == cur_poc_diff || !col_poc_diff) {
154         mvLXCol->x = mvCol->x;
155         mvLXCol->y = mvCol->y;
156     } else {
157         mv_scale(mvLXCol, mvCol, col_poc_diff, cur_poc_diff);
158     }
159     return 1;
160 }
161
162 #define CHECK_MVSET(l)                                          \
163     check_mvset(mvLXCol, temp_col.mv + l,                       \
164                 colPic, s->poc,                                 \
165                 refPicList, X, refIdxLx,                        \
166                 refPicList_col, L ## l, temp_col.ref_idx[l])
167
168 // derive the motion vectors section 8.5.3.1.8
169 static int derive_temporal_colocated_mvs(HEVCContext *s, MvField temp_col,
170                                          int refIdxLx, Mv *mvLXCol, int X,
171                                          int colPic, RefPicList *refPicList_col)
172 {
173     RefPicList *refPicList = s->ref->refPicList;
174
175     if (temp_col.pred_flag == PF_INTRA)
176         return 0;
177
178     if (!(temp_col.pred_flag & PF_L0))
179         return CHECK_MVSET(1);
180     else if (temp_col.pred_flag == PF_L0)
181         return CHECK_MVSET(0);
182     else if (temp_col.pred_flag == PF_BI) {
183         int check_diffpicount = 0;
184         int i, j;
185         for (j = 0; j < 2; j++) {
186             for (i = 0; i < refPicList[j].nb_refs; i++) {
187                 if (refPicList[j].list[i] > s->poc) {
188                     check_diffpicount++;
189                     break;
190                 }
191             }
192         }
193         if (!check_diffpicount) {
194             if (X==0)
195                 return CHECK_MVSET(0);
196             else
197                 return CHECK_MVSET(1);
198         } else {
199             if (s->sh.collocated_list == L1)
200                 return CHECK_MVSET(0);
201             else
202                 return CHECK_MVSET(1);
203         }
204     }
205
206     return 0;
207 }
208
209 #define TAB_MVF(x, y)                                                   \
210     tab_mvf[(y) * min_pu_width + x]
211
212 #define TAB_MVF_PU(v)                                                   \
213     TAB_MVF(x ## v ## _pu, y ## v ## _pu)
214
215 #define DERIVE_TEMPORAL_COLOCATED_MVS                                   \
216     derive_temporal_colocated_mvs(s, temp_col,                          \
217                                   refIdxLx, mvLXCol, X, colPic,         \
218                                   ff_hevc_get_ref_list(s, ref, x, y))
219
220 /*
221  * 8.5.3.1.7  temporal luma motion vector prediction
222  */
223 static int temporal_luma_motion_vector(HEVCContext *s, int x0, int y0,
224                                        int nPbW, int nPbH, int refIdxLx,
225                                        Mv *mvLXCol, int X)
226 {
227     MvField *tab_mvf;
228     MvField temp_col;
229     int x, y, x_pu, y_pu;
230     int min_pu_width = s->sps->min_pu_width;
231     int availableFlagLXCol = 0;
232     int colPic;
233
234     HEVCFrame *ref = s->ref->collocated_ref;
235
236     if (!ref)
237         return 0;
238
239     tab_mvf = ref->tab_mvf;
240     colPic  = ref->poc;
241
242     //bottom right collocated motion vector
243     x = x0 + nPbW;
244     y = y0 + nPbH;
245
246     if (s->threads_type == FF_THREAD_FRAME )
247         ff_thread_await_progress(&ref->tf, y, 0);
248     if (tab_mvf &&
249         (y0 >> s->sps->log2_ctb_size) == (y >> s->sps->log2_ctb_size) &&
250         y < s->sps->height &&
251         x < s->sps->width) {
252         x                 &= -16;
253         y                 &= -16;
254         x_pu               = x >> s->sps->log2_min_pu_size;
255         y_pu               = y >> s->sps->log2_min_pu_size;
256         temp_col           = TAB_MVF(x_pu, y_pu);
257         availableFlagLXCol = DERIVE_TEMPORAL_COLOCATED_MVS;
258     }
259
260     // derive center collocated motion vector
261     if (tab_mvf && !availableFlagLXCol) {
262         x                  = x0 + (nPbW >> 1);
263         y                  = y0 + (nPbH >> 1);
264         x                 &= -16;
265         y                 &= -16;
266         x_pu               = x >> s->sps->log2_min_pu_size;
267         y_pu               = y >> s->sps->log2_min_pu_size;
268         temp_col           = TAB_MVF(x_pu, y_pu);
269         availableFlagLXCol = DERIVE_TEMPORAL_COLOCATED_MVS;
270     }
271     return availableFlagLXCol;
272 }
273
274 #define AVAILABLE(cand, v)                                      \
275     (cand && !(TAB_MVF_PU(v).pred_flag == PF_INTRA))
276
277 #define PRED_BLOCK_AVAILABLE(v)                                 \
278     z_scan_block_avail(s, x0, y0, x ## v, y ## v)
279
280 #define COMPARE_MV_REFIDX(a, b)                                 \
281     compareMVrefidx(TAB_MVF_PU(a), TAB_MVF_PU(b))
282
283 /*
284  * 8.5.3.1.2  Derivation process for spatial merging candidates
285  */
286 static void derive_spatial_merge_candidates(HEVCContext *s, int x0, int y0,
287                                             int nPbW, int nPbH,
288                                             int log2_cb_size,
289                                             int singleMCLFlag, int part_idx,
290                                             int merge_idx,
291                                             struct MvField mergecandlist[])
292 {
293     HEVCLocalContext *lc   = s->HEVClc;
294     RefPicList *refPicList = s->ref->refPicList;
295     MvField *tab_mvf       = s->ref->tab_mvf;
296
297     const int min_pu_width = s->sps->min_pu_width;
298
299     const int cand_bottom_left = lc->na.cand_bottom_left;
300     const int cand_left        = lc->na.cand_left;
301     const int cand_up_left     = lc->na.cand_up_left;
302     const int cand_up          = lc->na.cand_up;
303     const int cand_up_right    = lc->na.cand_up_right_sap;
304
305     const int xA1    = x0 - 1;
306     const int yA1    = y0 + nPbH - 1;
307     const int xA1_pu = xA1 >> s->sps->log2_min_pu_size;
308     const int yA1_pu = yA1 >> s->sps->log2_min_pu_size;
309
310     const int xB1    = x0 + nPbW - 1;
311     const int yB1    = y0 - 1;
312     const int xB1_pu = xB1 >> s->sps->log2_min_pu_size;
313     const int yB1_pu = yB1 >> s->sps->log2_min_pu_size;
314
315     const int xB0    = x0 + nPbW;
316     const int yB0    = y0 - 1;
317     const int xB0_pu = xB0 >> s->sps->log2_min_pu_size;
318     const int yB0_pu = yB0 >> s->sps->log2_min_pu_size;
319
320     const int xA0    = x0 - 1;
321     const int yA0    = y0 + nPbH;
322     const int xA0_pu = xA0 >> s->sps->log2_min_pu_size;
323     const int yA0_pu = yA0 >> s->sps->log2_min_pu_size;
324
325     const int xB2    = x0 - 1;
326     const int yB2    = y0 - 1;
327     const int xB2_pu = xB2 >> s->sps->log2_min_pu_size;
328     const int yB2_pu = yB2 >> s->sps->log2_min_pu_size;
329
330     const int nb_refs = (s->sh.slice_type == P_SLICE) ?
331                         s->sh.nb_refs[0] : FFMIN(s->sh.nb_refs[0], s->sh.nb_refs[1]);
332
333     int zero_idx = 0;
334
335     int nb_merge_cand = 0;
336     int nb_orig_merge_cand = 0;
337
338     int is_available_a0;
339     int is_available_a1;
340     int is_available_b0;
341     int is_available_b1;
342     int is_available_b2;
343
344
345     if (!singleMCLFlag && part_idx == 1 &&
346         (lc->cu.part_mode == PART_Nx2N ||
347          lc->cu.part_mode == PART_nLx2N ||
348          lc->cu.part_mode == PART_nRx2N) ||
349         isDiffMER(s, xA1, yA1, x0, y0)) {
350         is_available_a1 = 0;
351     } else {
352         is_available_a1 = AVAILABLE(cand_left, A1);
353         if (is_available_a1) {
354             mergecandlist[nb_merge_cand] = TAB_MVF_PU(A1);
355             if (merge_idx == 0) return;
356             nb_merge_cand++;
357         }
358     }
359
360     if (!singleMCLFlag && part_idx == 1 &&
361         (lc->cu.part_mode == PART_2NxN ||
362          lc->cu.part_mode == PART_2NxnU ||
363          lc->cu.part_mode == PART_2NxnD) ||
364         isDiffMER(s, xB1, yB1, x0, y0)) {
365         is_available_b1 = 0;
366     } else {
367         is_available_b1 = AVAILABLE(cand_up, B1);
368         if (is_available_b1 &&
369             !(is_available_a1 && COMPARE_MV_REFIDX(B1, A1))) {
370             mergecandlist[nb_merge_cand] = TAB_MVF_PU(B1);
371             if (merge_idx == nb_merge_cand) return;
372             nb_merge_cand++;
373         }
374     }
375
376     // above right spatial merge candidate
377     is_available_b0 = AVAILABLE(cand_up_right, B0) &&
378                       xB0 < s->sps->width &&
379                       PRED_BLOCK_AVAILABLE(B0) &&
380                       !isDiffMER(s, xB0, yB0, x0, y0);
381
382     if (is_available_b0 &&
383         !(is_available_b1 && COMPARE_MV_REFIDX(B0, B1))) {
384         mergecandlist[nb_merge_cand] = TAB_MVF_PU(B0);
385         if (merge_idx == nb_merge_cand) return;
386         nb_merge_cand++;
387     }
388
389     // left bottom spatial merge candidate
390     is_available_a0 = AVAILABLE(cand_bottom_left, A0) &&
391                       yA0 < s->sps->height &&
392                       PRED_BLOCK_AVAILABLE(A0) &&
393                       !isDiffMER(s, xA0, yA0, x0, y0);
394
395     if (is_available_a0 &&
396         !(is_available_a1 && COMPARE_MV_REFIDX(A0, A1))) {
397         mergecandlist[nb_merge_cand] = TAB_MVF_PU(A0);
398         if (merge_idx == nb_merge_cand) return;
399         nb_merge_cand++;
400     }
401
402     // above left spatial merge candidate
403     is_available_b2 = AVAILABLE(cand_up_left, B2) &&
404                       !isDiffMER(s, xB2, yB2, x0, y0);
405
406     if (is_available_b2 &&
407         !(is_available_a1 && COMPARE_MV_REFIDX(B2, A1)) &&
408         !(is_available_b1 && COMPARE_MV_REFIDX(B2, B1)) &&
409         nb_merge_cand != 4) {
410         mergecandlist[nb_merge_cand] = TAB_MVF_PU(B2);
411         if (merge_idx == nb_merge_cand) return;
412         nb_merge_cand++;
413     }
414
415     // temporal motion vector candidate
416     if (s->sh.slice_temporal_mvp_enabled_flag &&
417         nb_merge_cand < s->sh.max_num_merge_cand) {
418         Mv mv_l0_col, mv_l1_col;
419         int available_l0 = temporal_luma_motion_vector(s, x0, y0, nPbW, nPbH,
420                                                        0, &mv_l0_col, 0);
421         int available_l1 = (s->sh.slice_type == B_SLICE) ?
422                            temporal_luma_motion_vector(s, x0, y0, nPbW, nPbH,
423                                                        0, &mv_l1_col, 1) : 0;
424
425         if (available_l0 || available_l1) {
426             mergecandlist[nb_merge_cand].pred_flag = available_l0 + (available_l1 << 1);
427             if (available_l0) {
428                 mergecandlist[nb_merge_cand].mv[0]      = mv_l0_col;
429                 mergecandlist[nb_merge_cand].ref_idx[0] = 0;
430             }
431             if (available_l1) {
432                 mergecandlist[nb_merge_cand].mv[1]      = mv_l1_col;
433                 mergecandlist[nb_merge_cand].ref_idx[1] = 0;
434             }
435             if (merge_idx == nb_merge_cand) return;
436             nb_merge_cand++;
437         }
438     }
439
440     nb_orig_merge_cand = nb_merge_cand;
441
442     // combined bi-predictive merge candidates  (applies for B slices)
443     if (s->sh.slice_type == B_SLICE && nb_orig_merge_cand > 1 &&
444         nb_orig_merge_cand < s->sh.max_num_merge_cand) {
445         int comb_idx = 0;
446
447         for (comb_idx = 0; nb_merge_cand < s->sh.max_num_merge_cand &&
448                            comb_idx < nb_orig_merge_cand * (nb_orig_merge_cand - 1); comb_idx++) {
449             int l0_cand_idx = l0_l1_cand_idx[comb_idx][0];
450             int l1_cand_idx = l0_l1_cand_idx[comb_idx][1];
451             MvField l0_cand = mergecandlist[l0_cand_idx];
452             MvField l1_cand = mergecandlist[l1_cand_idx];
453
454             if ((l0_cand.pred_flag & PF_L0) && (l1_cand.pred_flag & PF_L1) &&
455                 (refPicList[0].list[l0_cand.ref_idx[0]] !=
456                  refPicList[1].list[l1_cand.ref_idx[1]] ||
457                  AV_RN32A(&l0_cand.mv[0]) != AV_RN32A(&l1_cand.mv[1]))) {
458                 mergecandlist[nb_merge_cand].ref_idx[0]   = l0_cand.ref_idx[0];
459                 mergecandlist[nb_merge_cand].ref_idx[1]   = l1_cand.ref_idx[1];
460                 mergecandlist[nb_merge_cand].pred_flag    = PF_BI;
461                 AV_COPY32(&mergecandlist[nb_merge_cand].mv[0], &l0_cand.mv[0]);
462                 AV_COPY32(&mergecandlist[nb_merge_cand].mv[1], &l1_cand.mv[1]);
463                 if (merge_idx == nb_merge_cand) return;
464                 nb_merge_cand++;
465             }
466         }
467     }
468
469     // append Zero motion vector candidates
470     while (nb_merge_cand < s->sh.max_num_merge_cand) {
471         mergecandlist[nb_merge_cand].pred_flag    = PF_L0 + ((s->sh.slice_type == B_SLICE) << 1);
472         AV_ZERO32(mergecandlist[nb_merge_cand].mv+0);
473         AV_ZERO32(mergecandlist[nb_merge_cand].mv+1);
474         mergecandlist[nb_merge_cand].ref_idx[0]   = zero_idx < nb_refs ? zero_idx : 0;
475         mergecandlist[nb_merge_cand].ref_idx[1]   = zero_idx < nb_refs ? zero_idx : 0;
476
477         if (merge_idx == nb_merge_cand) return;
478         nb_merge_cand++;
479         zero_idx++;
480     }
481 }
482
483 /*
484  * 8.5.3.1.1 Derivation process of luma Mvs for merge mode
485  */
486 void ff_hevc_luma_mv_merge_mode(HEVCContext *s, int x0, int y0, int nPbW,
487                                 int nPbH, int log2_cb_size, int part_idx,
488                                 int merge_idx, MvField *mv)
489 {
490     int singleMCLFlag = 0;
491     int nCS = 1 << log2_cb_size;
492     LOCAL_ALIGNED(4, MvField, mergecand_list, [MRG_MAX_NUM_CANDS]);
493     int nPbW2 = nPbW;
494     int nPbH2 = nPbH;
495     HEVCLocalContext *lc = s->HEVClc;
496
497     if (s->pps->log2_parallel_merge_level > 2 && nCS == 8) {
498         singleMCLFlag = 1;
499         x0            = lc->cu.x;
500         y0            = lc->cu.y;
501         nPbW          = nCS;
502         nPbH          = nCS;
503         part_idx      = 0;
504     }
505
506     ff_hevc_set_neighbour_available(s, x0, y0, nPbW, nPbH);
507     derive_spatial_merge_candidates(s, x0, y0, nPbW, nPbH, log2_cb_size,
508                                     singleMCLFlag, part_idx,
509                                     merge_idx, mergecand_list);
510
511     if (mergecand_list[merge_idx].pred_flag == PF_BI &&
512         (nPbW2 + nPbH2) == 12) {
513         mergecand_list[merge_idx].pred_flag = PF_L0;
514     }
515
516     *mv = mergecand_list[merge_idx];
517 }
518
519 static av_always_inline void dist_scale(HEVCContext *s, Mv *mv,
520                                         int min_pu_width, int x, int y,
521                                         int elist, int ref_idx_curr, int ref_idx)
522 {
523     RefPicList *refPicList = s->ref->refPicList;
524     MvField *tab_mvf       = s->ref->tab_mvf;
525     int ref_pic_elist      = refPicList[elist].list[TAB_MVF(x, y).ref_idx[elist]];
526     int ref_pic_curr       = refPicList[ref_idx_curr].list[ref_idx];
527
528     if (ref_pic_elist != ref_pic_curr) {
529         int poc_diff = s->poc - ref_pic_elist;
530         if (!poc_diff)
531             poc_diff = 1;
532         mv_scale(mv, mv, poc_diff, s->poc - ref_pic_curr);
533     }
534 }
535
536 static int mv_mp_mode_mx(HEVCContext *s, int x, int y, int pred_flag_index,
537                          Mv *mv, int ref_idx_curr, int ref_idx)
538 {
539     MvField *tab_mvf = s->ref->tab_mvf;
540     int min_pu_width = s->sps->min_pu_width;
541
542     RefPicList *refPicList = s->ref->refPicList;
543
544     if (((TAB_MVF(x, y).pred_flag) & (1 << pred_flag_index)) &&
545         refPicList[pred_flag_index].list[TAB_MVF(x, y).ref_idx[pred_flag_index]] == refPicList[ref_idx_curr].list[ref_idx]) {
546         *mv = TAB_MVF(x, y).mv[pred_flag_index];
547         return 1;
548     }
549     return 0;
550 }
551
552 static int mv_mp_mode_mx_lt(HEVCContext *s, int x, int y, int pred_flag_index,
553                             Mv *mv, int ref_idx_curr, int ref_idx)
554 {
555     MvField *tab_mvf = s->ref->tab_mvf;
556     int min_pu_width = s->sps->min_pu_width;
557
558     RefPicList *refPicList = s->ref->refPicList;
559
560     if ((TAB_MVF(x, y).pred_flag) & (1 << pred_flag_index)) {
561         int currIsLongTerm     = refPicList[ref_idx_curr].isLongTerm[ref_idx];
562
563         int colIsLongTerm =
564             refPicList[pred_flag_index].isLongTerm[(TAB_MVF(x, y).ref_idx[pred_flag_index])];
565
566         if (colIsLongTerm == currIsLongTerm) {
567             *mv = TAB_MVF(x, y).mv[pred_flag_index];
568             if (!currIsLongTerm)
569                 dist_scale(s, mv, min_pu_width, x, y,
570                            pred_flag_index, ref_idx_curr, ref_idx);
571             return 1;
572         }
573     }
574     return 0;
575 }
576
577 #define MP_MX(v, pred, mx)                                      \
578     mv_mp_mode_mx(s, x ## v ## _pu, y ## v ## _pu, pred,        \
579                   &mx, ref_idx_curr, ref_idx)
580
581 #define MP_MX_LT(v, pred, mx)                                   \
582     mv_mp_mode_mx_lt(s, x ## v ## _pu, y ## v ## _pu, pred,     \
583                      &mx, ref_idx_curr, ref_idx)
584
585 void ff_hevc_luma_mv_mvp_mode(HEVCContext *s, int x0, int y0, int nPbW,
586                               int nPbH, int log2_cb_size, int part_idx,
587                               int merge_idx, MvField *mv,
588                               int mvp_lx_flag, int LX)
589 {
590     HEVCLocalContext *lc = s->HEVClc;
591     MvField *tab_mvf = s->ref->tab_mvf;
592     int isScaledFlag_L0 = 0;
593     int availableFlagLXA0 = 0;
594     int availableFlagLXB0 = 0;
595     int numMVPCandLX = 0;
596     int min_pu_width = s->sps->min_pu_width;
597
598     int xA0, yA0;
599     int xA0_pu, yA0_pu;
600     int is_available_a0;
601
602     int xA1, yA1;
603     int xA1_pu, yA1_pu;
604     int is_available_a1;
605
606     int xB0, yB0;
607     int xB0_pu, yB0_pu;
608     int is_available_b0;
609
610     int xB1, yB1;
611     int xB1_pu = 0, yB1_pu = 0;
612     int is_available_b1 = 0;
613
614     int xB2, yB2;
615     int xB2_pu = 0, yB2_pu = 0;
616     int is_available_b2 = 0;
617     Mv mvpcand_list[2] = { { 0 } };
618     Mv mxA;
619     Mv mxB;
620     int ref_idx_curr = 0;
621     int ref_idx = 0;
622     int pred_flag_index_l0;
623     int pred_flag_index_l1;
624     int x0b = x0 & ((1 << s->sps->log2_ctb_size) - 1);
625     int y0b = y0 & ((1 << s->sps->log2_ctb_size) - 1);
626
627     int cand_up = (lc->ctb_up_flag || y0b);
628     int cand_left = (lc->ctb_left_flag || x0b);
629     int cand_up_left =
630             (!x0b && !y0b) ? lc->ctb_up_left_flag : cand_left && cand_up;
631     int cand_up_right =
632             (x0b + nPbW == (1 << s->sps->log2_ctb_size) ||
633              x0  + nPbW >= lc->end_of_tiles_x) ? lc->ctb_up_right_flag && !y0b
634                                                : cand_up;
635     int cand_bottom_left = (y0 + nPbH >= lc->end_of_tiles_y) ? 0 : cand_left;
636
637     ref_idx_curr       = LX;
638     ref_idx            = mv->ref_idx[LX];
639     pred_flag_index_l0 = LX;
640     pred_flag_index_l1 = !LX;
641
642     // left bottom spatial candidate
643     xA0 = x0 - 1;
644     yA0 = y0 + nPbH;
645     xA0_pu = xA0 >> s->sps->log2_min_pu_size;
646     yA0_pu = yA0 >> s->sps->log2_min_pu_size;
647
648     is_available_a0 = AVAILABLE(cand_bottom_left, A0) &&
649                       yA0 < s->sps->height &&
650                       PRED_BLOCK_AVAILABLE(A0);
651
652     //left spatial merge candidate
653     xA1    = x0 - 1;
654     yA1    = y0 + nPbH - 1;
655     xA1_pu = xA1 >> s->sps->log2_min_pu_size;
656     yA1_pu = yA1 >> s->sps->log2_min_pu_size;
657
658     is_available_a1 = AVAILABLE(cand_left, A1);
659     if (is_available_a0 || is_available_a1)
660         isScaledFlag_L0 = 1;
661
662     if (is_available_a0) {
663         availableFlagLXA0 = MP_MX(A0, pred_flag_index_l0, mxA);
664         if (!availableFlagLXA0)
665             availableFlagLXA0 = MP_MX(A0, pred_flag_index_l1, mxA);
666     }
667
668     if (is_available_a1 && !availableFlagLXA0) {
669         availableFlagLXA0 = MP_MX(A1, pred_flag_index_l0, mxA);
670         if (!availableFlagLXA0)
671             availableFlagLXA0 = MP_MX(A1, pred_flag_index_l1, mxA);
672     }
673
674     if (is_available_a0 && !availableFlagLXA0) {
675         availableFlagLXA0 = MP_MX_LT(A0, pred_flag_index_l0, mxA);
676         if (!availableFlagLXA0)
677             availableFlagLXA0 = MP_MX_LT(A0, pred_flag_index_l1, mxA);
678     }
679
680     if (is_available_a1 && !availableFlagLXA0) {
681         availableFlagLXA0 = MP_MX_LT(A1, pred_flag_index_l0, mxA);
682         if (!availableFlagLXA0)
683             availableFlagLXA0 = MP_MX_LT(A1, pred_flag_index_l1, mxA);
684     }
685
686     if(availableFlagLXA0 && !mvp_lx_flag) {
687         mv->mv[LX] = mxA;
688         return;
689     }
690
691     // B candidates
692     // above right spatial merge candidate
693     xB0    = x0 + nPbW;
694     yB0    = y0 - 1;
695     xB0_pu = xB0 >> s->sps->log2_min_pu_size;
696     yB0_pu = yB0 >> s->sps->log2_min_pu_size;
697
698     is_available_b0 = PRED_BLOCK_AVAILABLE(B0) && AVAILABLE(cand_up_right, B0);
699
700     if (is_available_b0) {
701         availableFlagLXB0 = MP_MX(B0, pred_flag_index_l0, mxB);
702         if (!availableFlagLXB0)
703             availableFlagLXB0 = MP_MX(B0, pred_flag_index_l1, mxB);
704     }
705
706     if (!availableFlagLXB0) {
707         // above spatial merge candidate
708         xB1    = x0 + nPbW - 1;
709         yB1    = y0 - 1;
710         xB1_pu = xB1 >> s->sps->log2_min_pu_size;
711         yB1_pu = yB1 >> s->sps->log2_min_pu_size;
712
713         is_available_b1 = AVAILABLE(cand_up, B1);
714
715         if (is_available_b1) {
716             availableFlagLXB0 = MP_MX(B1, pred_flag_index_l0, mxB);
717             if (!availableFlagLXB0)
718                 availableFlagLXB0 = MP_MX(B1, pred_flag_index_l1, mxB);
719         }
720     }
721
722     if (!availableFlagLXB0) {
723         // above left spatial merge candidate
724         xB2 = x0 - 1;
725         yB2 = y0 - 1;
726         xB2_pu = xB2 >> s->sps->log2_min_pu_size;
727         yB2_pu = yB2 >> s->sps->log2_min_pu_size;
728         is_available_b2 = AVAILABLE(cand_up_left, B2);
729
730         if (is_available_b2) {
731             availableFlagLXB0 = MP_MX(B2, pred_flag_index_l0, mxB);
732             if (!availableFlagLXB0)
733                 availableFlagLXB0 = MP_MX(B2, pred_flag_index_l1, mxB);
734         }
735     }
736
737     if (isScaledFlag_L0 == 0) {
738         if (availableFlagLXB0) {
739             availableFlagLXA0 = 1;
740             mxA = mxB;
741         }
742         availableFlagLXB0 = 0;
743
744         // XB0 and L1
745         if (is_available_b0) {
746             availableFlagLXB0 = MP_MX_LT(B0, pred_flag_index_l0, mxB);
747             if (!availableFlagLXB0)
748                 availableFlagLXB0 = MP_MX_LT(B0, pred_flag_index_l1, mxB);
749         }
750
751         if (is_available_b1 && !availableFlagLXB0) {
752             availableFlagLXB0 = MP_MX_LT(B1, pred_flag_index_l0, mxB);
753             if (!availableFlagLXB0)
754                 availableFlagLXB0 = MP_MX_LT(B1, pred_flag_index_l1, mxB);
755         }
756
757         if (is_available_b2 && !availableFlagLXB0) {
758             availableFlagLXB0 = MP_MX_LT(B2, pred_flag_index_l0, mxB);
759             if (!availableFlagLXB0)
760                 availableFlagLXB0 = MP_MX_LT(B2, pred_flag_index_l1, mxB);
761         }
762     }
763
764     if (availableFlagLXA0)
765         mvpcand_list[numMVPCandLX++] = mxA;
766
767     if (availableFlagLXB0 && (!availableFlagLXA0 || mxA.x != mxB.x || mxA.y != mxB.y))
768         mvpcand_list[numMVPCandLX++] = mxB;
769
770     //temporal motion vector prediction candidate
771     if (numMVPCandLX < 2 && s->sh.slice_temporal_mvp_enabled_flag) {
772         Mv mv_col;
773         int available_col = temporal_luma_motion_vector(s, x0, y0, nPbW,
774                                                         nPbH, ref_idx,
775                                                         &mv_col, LX);
776         if (available_col)
777             mvpcand_list[numMVPCandLX++] = mv_col;
778     }
779
780     mv->mv[LX] = mvpcand_list[mvp_lx_flag];
781 }