hevc: Use parsed VUI colorimetry in avcodec
[ffmpeg.git] / libavcodec / hevc_ps.c
1 /*
2  * HEVC Parameter Set decoding
3  *
4  * Copyright (C) 2012 - 2103 Guillaume Martres
5  * Copyright (C) 2012 - 2103 Mickael Raulet
6  * Copyright (C) 2012 - 2013 Gildas Cocherel
7  * Copyright (C) 2013 Vittorio Giovara
8  *
9  * This file is part of Libav.
10  *
11  * Libav is free software; you can redistribute it and/or
12  * modify it under the terms of the GNU Lesser General Public
13  * License as published by the Free Software Foundation; either
14  * version 2.1 of the License, or (at your option) any later version.
15  *
16  * Libav is distributed in the hope that it will be useful,
17  * but WITHOUT ANY WARRANTY; without even the implied warranty of
18  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
19  * Lesser General Public License for more details.
20  *
21  * You should have received a copy of the GNU Lesser General Public
22  * License along with Libav; if not, write to the Free Software
23  * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
24  */
25
26 #include "libavutil/imgutils.h"
27
28 #include "golomb.h"
29 #include "hevc.h"
30
31 static const uint8_t default_scaling_list_intra[] = {
32     16, 16, 16, 16, 17, 18, 21, 24,
33     16, 16, 16, 16, 17, 19, 22, 25,
34     16, 16, 17, 18, 20, 22, 25, 29,
35     16, 16, 18, 21, 24, 27, 31, 36,
36     17, 17, 20, 24, 30, 35, 41, 47,
37     18, 19, 22, 27, 35, 44, 54, 65,
38     21, 22, 25, 31, 41, 54, 70, 88,
39     24, 25, 29, 36, 47, 65, 88, 115
40 };
41
42 static const uint8_t default_scaling_list_inter[] = {
43     16, 16, 16, 16, 17, 18, 20, 24,
44     16, 16, 16, 17, 18, 20, 24, 25,
45     16, 16, 17, 18, 20, 24, 25, 28,
46     16, 17, 18, 20, 24, 25, 28, 33,
47     17, 18, 20, 24, 25, 28, 33, 41,
48     18, 20, 24, 25, 28, 33, 41, 54,
49     20, 24, 25, 28, 33, 41, 54, 71,
50     24, 25, 28, 33, 41, 54, 71, 91
51 };
52
53 static const AVRational vui_sar[] = {
54     {  0,   1 },
55     {  1,   1 },
56     { 12,  11 },
57     { 10,  11 },
58     { 16,  11 },
59     { 40,  33 },
60     { 24,  11 },
61     { 20,  11 },
62     { 32,  11 },
63     { 80,  33 },
64     { 18,  11 },
65     { 15,  11 },
66     { 64,  33 },
67     { 160, 99 },
68     {  4,   3 },
69     {  3,   2 },
70     {  2,   1 },
71 };
72
73 int ff_hevc_decode_short_term_rps(HEVCContext *s, ShortTermRPS *rps,
74                                   const HEVCSPS *sps, int is_slice_header)
75 {
76     HEVCLocalContext *lc = &s->HEVClc;
77     uint8_t rps_predict = 0;
78     int delta_poc;
79     int k0 = 0;
80     int k1 = 0;
81     int k  = 0;
82     int i;
83
84     GetBitContext *gb = &lc->gb;
85
86     if (rps != sps->st_rps && sps->nb_st_rps)
87         rps_predict = get_bits1(gb);
88
89     if (rps_predict) {
90         const ShortTermRPS *rps_ridx;
91         int delta_rps, abs_delta_rps;
92         uint8_t use_delta_flag = 0;
93         uint8_t delta_rps_sign;
94
95         if (is_slice_header) {
96             int delta_idx = get_ue_golomb_long(gb) + 1;
97             if (delta_idx > sps->nb_st_rps) {
98                 av_log(s->avctx, AV_LOG_ERROR,
99                        "Invalid value of delta_idx in slice header RPS: %d > %d.\n",
100                        delta_idx, sps->nb_st_rps);
101                 return AVERROR_INVALIDDATA;
102             }
103             rps_ridx = &sps->st_rps[sps->nb_st_rps - delta_idx];
104         } else
105             rps_ridx = &sps->st_rps[rps - sps->st_rps - 1];
106
107         delta_rps_sign = get_bits1(gb);
108         abs_delta_rps  = get_ue_golomb_long(gb) + 1;
109         delta_rps      = (1 - (delta_rps_sign << 1)) * abs_delta_rps;
110         for (i = 0; i <= rps_ridx->num_delta_pocs; i++) {
111             int used = rps->used[k] = get_bits1(gb);
112
113             if (!used)
114                 use_delta_flag = get_bits1(gb);
115
116             if (used || use_delta_flag) {
117                 if (i < rps_ridx->num_delta_pocs)
118                     delta_poc = delta_rps + rps_ridx->delta_poc[i];
119                 else
120                     delta_poc = delta_rps;
121                 rps->delta_poc[k] = delta_poc;
122                 if (delta_poc < 0)
123                     k0++;
124                 else
125                     k1++;
126                 k++;
127             }
128         }
129
130         rps->num_delta_pocs    = k;
131         rps->num_negative_pics = k0;
132         // sort in increasing order (smallest first)
133         if (rps->num_delta_pocs != 0) {
134             int used, tmp;
135             for (i = 1; i < rps->num_delta_pocs; i++) {
136                 delta_poc = rps->delta_poc[i];
137                 used      = rps->used[i];
138                 for (k = i - 1; k >= 0; k--) {
139                     tmp = rps->delta_poc[k];
140                     if (delta_poc < tmp) {
141                         rps->delta_poc[k + 1] = tmp;
142                         rps->used[k + 1]      = rps->used[k];
143                         rps->delta_poc[k]     = delta_poc;
144                         rps->used[k]          = used;
145                     }
146                 }
147             }
148         }
149         if ((rps->num_negative_pics >> 1) != 0) {
150             int used;
151             k = rps->num_negative_pics - 1;
152             // flip the negative values to largest first
153             for (i = 0; i < rps->num_negative_pics >> 1; i++) {
154                 delta_poc         = rps->delta_poc[i];
155                 used              = rps->used[i];
156                 rps->delta_poc[i] = rps->delta_poc[k];
157                 rps->used[i]      = rps->used[k];
158                 rps->delta_poc[k] = delta_poc;
159                 rps->used[k]      = used;
160                 k--;
161             }
162         }
163     } else {
164         unsigned int prev, nb_positive_pics;
165         rps->num_negative_pics = get_ue_golomb_long(gb);
166         nb_positive_pics       = get_ue_golomb_long(gb);
167
168         if (rps->num_negative_pics >= MAX_REFS ||
169             nb_positive_pics >= MAX_REFS) {
170             av_log(s->avctx, AV_LOG_ERROR, "Too many refs in a short term RPS.\n");
171             return AVERROR_INVALIDDATA;
172         }
173
174         rps->num_delta_pocs = rps->num_negative_pics + nb_positive_pics;
175         if (rps->num_delta_pocs) {
176             prev = 0;
177             for (i = 0; i < rps->num_negative_pics; i++) {
178                 delta_poc = get_ue_golomb_long(gb) + 1;
179                 prev -= delta_poc;
180                 rps->delta_poc[i] = prev;
181                 rps->used[i]      = get_bits1(gb);
182             }
183             prev = 0;
184             for (i = 0; i < nb_positive_pics; i++) {
185                 delta_poc = get_ue_golomb_long(gb) + 1;
186                 prev += delta_poc;
187                 rps->delta_poc[rps->num_negative_pics + i] = prev;
188                 rps->used[rps->num_negative_pics + i]      = get_bits1(gb);
189             }
190         }
191     }
192     return 0;
193 }
194
195 static int decode_profile_tier_level(HEVCLocalContext *lc, PTL *ptl,
196                                      int max_num_sub_layers)
197 {
198     int i, j;
199     GetBitContext *gb = &lc->gb;
200
201     ptl->general_profile_space = get_bits(gb, 2);
202     ptl->general_tier_flag     = get_bits1(gb);
203     ptl->general_profile_idc   = get_bits(gb, 5);
204     for (i = 0; i < 32; i++)
205         ptl->general_profile_compatibility_flag[i] = get_bits1(gb);
206     skip_bits1(gb); // general_progressive_source_flag
207     skip_bits1(gb); // general_interlaced_source_flag
208     skip_bits1(gb); // general_non_packed_constraint_flag
209     skip_bits1(gb); // general_frame_only_constraint_flag
210     if (get_bits(gb, 16) != 0) // XXX_reserved_zero_44bits[0..15]
211         return -1;
212     if (get_bits(gb, 16) != 0) // XXX_reserved_zero_44bits[16..31]
213         return -1;
214     if (get_bits(gb, 12) != 0) // XXX_reserved_zero_44bits[32..43]
215         return -1;
216
217     ptl->general_level_idc = get_bits(gb, 8);
218     for (i = 0; i < max_num_sub_layers - 1; i++) {
219         ptl->sub_layer_profile_present_flag[i] = get_bits1(gb);
220         ptl->sub_layer_level_present_flag[i]   = get_bits1(gb);
221     }
222     if (max_num_sub_layers - 1 > 0)
223         for (i = max_num_sub_layers - 1; i < 8; i++)
224             skip_bits(gb, 2);  // reserved_zero_2bits[i]
225     for (i = 0; i < max_num_sub_layers - 1; i++) {
226         if (ptl->sub_layer_profile_present_flag[i]) {
227             ptl->sub_layer_profile_space[i] = get_bits(gb, 2);
228             ptl->sub_layer_tier_flag[i]     = get_bits(gb, 1);
229             ptl->sub_layer_profile_idc[i]   = get_bits(gb, 5);
230             for (j = 0; j < 32; j++)
231                 ptl->sub_layer_profile_compatibility_flags[i][j] = get_bits1(gb);
232             skip_bits1(gb); // sub_layer_progressive_source_flag
233             skip_bits1(gb); // sub_layer_interlaced_source_flag
234             skip_bits1(gb); // sub_layer_non_packed_constraint_flag
235             skip_bits1(gb); // sub_layer_frame_only_constraint_flag
236
237             if (get_bits(gb, 16) != 0) // sub_layer_reserved_zero_44bits[0..15]
238                 return -1;
239             if (get_bits(gb, 16) != 0) // sub_layer_reserved_zero_44bits[16..31]
240                 return -1;
241             if (get_bits(gb, 12) != 0) // sub_layer_reserved_zero_44bits[32..43]
242                 return -1;
243         }
244         if (ptl->sub_layer_level_present_flag[i])
245             ptl->sub_layer_level_idc[i] = get_bits(gb, 8);
246     }
247     return 0;
248 }
249
250 static void decode_sublayer_hrd(HEVCContext *s, int nb_cpb,
251                                 int subpic_params_present)
252 {
253     GetBitContext *gb = &s->HEVClc.gb;
254     int i;
255
256     for (i = 0; i < nb_cpb; i++) {
257         get_ue_golomb_long(gb); // bit_rate_value_minus1
258         get_ue_golomb_long(gb); // cpb_size_value_minus1
259
260         if (subpic_params_present) {
261             get_ue_golomb_long(gb); // cpb_size_du_value_minus1
262             get_ue_golomb_long(gb); // bit_rate_du_value_minus1
263         }
264         skip_bits1(gb); // cbr_flag
265     }
266 }
267
268 static void decode_hrd(HEVCContext *s, int common_inf_present,
269                        int max_sublayers)
270 {
271     GetBitContext *gb = &s->HEVClc.gb;
272     int nal_params_present = 0, vcl_params_present = 0;
273     int subpic_params_present = 0;
274     int i;
275
276     if (common_inf_present) {
277         nal_params_present = get_bits1(gb);
278         vcl_params_present = get_bits1(gb);
279
280         if (nal_params_present || vcl_params_present) {
281             subpic_params_present = get_bits1(gb);
282
283             if (subpic_params_present) {
284                 skip_bits(gb, 8); // tick_divisor_minus2
285                 skip_bits(gb, 5); // du_cpb_removal_delay_increment_length_minus1
286                 skip_bits(gb, 1); // sub_pic_cpb_params_in_pic_timing_sei_flag
287                 skip_bits(gb, 5); // dpb_output_delay_du_length_minus1
288             }
289
290             skip_bits(gb, 4); // bit_rate_scale
291             skip_bits(gb, 4); // cpb_size_scale
292
293             if (subpic_params_present)
294                 skip_bits(gb, 4);  // cpb_size_du_scale
295
296             skip_bits(gb, 5); // initial_cpb_removal_delay_length_minus1
297             skip_bits(gb, 5); // au_cpb_removal_delay_length_minus1
298             skip_bits(gb, 5); // dpb_output_delay_length_minus1
299         }
300     }
301
302     for (i = 0; i < max_sublayers; i++) {
303         int low_delay = 0;
304         int nb_cpb = 1;
305         int fixed_rate = get_bits1(gb);
306
307         if (!fixed_rate)
308             fixed_rate = get_bits1(gb);
309
310         if (fixed_rate)
311             get_ue_golomb_long(gb);  // elemental_duration_in_tc_minus1
312         else
313             low_delay = get_bits1(gb);
314
315         if (!low_delay)
316             nb_cpb = get_ue_golomb_long(gb) + 1;
317
318         if (nal_params_present)
319             decode_sublayer_hrd(s, nb_cpb, subpic_params_present);
320         if (vcl_params_present)
321             decode_sublayer_hrd(s, nb_cpb, subpic_params_present);
322     }
323 }
324
325 int ff_hevc_decode_nal_vps(HEVCContext *s)
326 {
327     int i,j;
328     GetBitContext *gb = &s->HEVClc.gb;
329     int vps_id = 0;
330     HEVCVPS *vps;
331
332     av_log(s->avctx, AV_LOG_DEBUG, "Decoding VPS\n");
333
334     vps = av_mallocz(sizeof(*vps));
335     if (!vps)
336         return AVERROR(ENOMEM);
337
338     vps_id = get_bits(gb, 4);
339     if (vps_id >= MAX_VPS_COUNT) {
340         av_log(s->avctx, AV_LOG_ERROR, "VPS id out of range: %d\n", vps_id);
341         goto err;
342     }
343
344     if (get_bits(gb, 2) != 3) { // vps_reserved_three_2bits
345         av_log(s->avctx, AV_LOG_ERROR, "vps_reserved_three_2bits is not three\n");
346         goto err;
347     }
348
349     vps->vps_max_layers               = get_bits(gb, 6) + 1;
350     vps->vps_max_sub_layers           = get_bits(gb, 3) + 1;
351     vps->vps_temporal_id_nesting_flag = get_bits1(gb);
352
353     if (get_bits(gb, 16) != 0xffff) { // vps_reserved_ffff_16bits
354         av_log(s->avctx, AV_LOG_ERROR, "vps_reserved_ffff_16bits is not 0xffff\n");
355         goto err;
356     }
357
358     if (vps->vps_max_sub_layers > MAX_SUB_LAYERS) {
359         av_log(s->avctx, AV_LOG_ERROR, "vps_max_sub_layers out of range: %d\n",
360                vps->vps_max_sub_layers);
361         goto err;
362     }
363
364     if (decode_profile_tier_level(&s->HEVClc, &vps->ptl, vps->vps_max_sub_layers) < 0) {
365         av_log(s->avctx, AV_LOG_ERROR, "Error decoding profile tier level.\n");
366         goto err;
367     }
368     vps->vps_sub_layer_ordering_info_present_flag = get_bits1(gb);
369
370     i = vps->vps_sub_layer_ordering_info_present_flag ? 0 : vps->vps_max_sub_layers - 1;
371     for (; i < vps->vps_max_sub_layers; i++) {
372         vps->vps_max_dec_pic_buffering[i] = get_ue_golomb_long(gb) + 1;
373         vps->vps_num_reorder_pics[i]      = get_ue_golomb_long(gb);
374         vps->vps_max_latency_increase[i]  = get_ue_golomb_long(gb) - 1;
375
376         if (vps->vps_max_dec_pic_buffering[i] > MAX_DPB_SIZE) {
377             av_log(s->avctx, AV_LOG_ERROR, "vps_max_dec_pic_buffering_minus1 out of range: %d\n",
378                    vps->vps_max_dec_pic_buffering[i] - 1);
379             goto err;
380         }
381         if (vps->vps_num_reorder_pics[i] > vps->vps_max_dec_pic_buffering[i] - 1) {
382             av_log(s->avctx, AV_LOG_ERROR, "vps_max_num_reorder_pics out of range: %d\n",
383                    vps->vps_num_reorder_pics[i]);
384             goto err;
385         }
386     }
387
388     vps->vps_max_layer_id   = get_bits(gb, 6);
389     vps->vps_num_layer_sets = get_ue_golomb_long(gb) + 1;
390     for (i = 1; i < vps->vps_num_layer_sets; i++)
391         for (j = 0; j <= vps->vps_max_layer_id; j++)
392             skip_bits(gb, 1);  // layer_id_included_flag[i][j]
393
394     vps->vps_timing_info_present_flag = get_bits1(gb);
395     if (vps->vps_timing_info_present_flag) {
396         vps->vps_num_units_in_tick               = get_bits_long(gb, 32);
397         vps->vps_time_scale                      = get_bits_long(gb, 32);
398         vps->vps_poc_proportional_to_timing_flag = get_bits1(gb);
399         if (vps->vps_poc_proportional_to_timing_flag)
400             vps->vps_num_ticks_poc_diff_one = get_ue_golomb_long(gb) + 1;
401         vps->vps_num_hrd_parameters = get_ue_golomb_long(gb);
402         for (i = 0; i < vps->vps_num_hrd_parameters; i++) {
403             int common_inf_present = 1;
404
405             get_ue_golomb_long(gb); // hrd_layer_set_idx
406             if (i)
407                 common_inf_present = get_bits1(gb);
408             decode_hrd(s, common_inf_present, vps->vps_max_sub_layers);
409         }
410     }
411     get_bits1(gb); /* vps_extension_flag */
412
413     av_free(s->vps_list[vps_id]);
414     s->vps_list[vps_id] = vps;
415     return 0;
416
417 err:
418     av_free(vps);
419     return AVERROR_INVALIDDATA;
420 }
421
422 static void decode_vui(HEVCContext *s, HEVCSPS *sps)
423 {
424     VUI *vui          = &sps->vui;
425     GetBitContext *gb = &s->HEVClc.gb;
426     int sar_present;
427
428     av_log(s->avctx, AV_LOG_DEBUG, "Decoding VUI\n");
429
430     sar_present = get_bits1(gb);
431     if (sar_present) {
432         uint8_t sar_idx = get_bits(gb, 8);
433         if (sar_idx < FF_ARRAY_ELEMS(vui_sar))
434             vui->sar = vui_sar[sar_idx];
435         else if (sar_idx == 255) {
436             vui->sar.num = get_bits(gb, 16);
437             vui->sar.den = get_bits(gb, 16);
438         } else
439             av_log(s->avctx, AV_LOG_WARNING,
440                    "Unknown SAR index: %u.\n", sar_idx);
441     }
442
443     vui->overscan_info_present_flag = get_bits1(gb);
444     if (vui->overscan_info_present_flag)
445         vui->overscan_appropriate_flag = get_bits1(gb);
446
447     vui->video_signal_type_present_flag = get_bits1(gb);
448     if (vui->video_signal_type_present_flag) {
449         vui->video_format                    = get_bits(gb, 3);
450         vui->video_full_range_flag           = get_bits1(gb);
451         vui->colour_description_present_flag = get_bits1(gb);
452         if (vui->colour_description_present_flag) {
453             vui->colour_primaries        = get_bits(gb, 8);
454             vui->transfer_characteristic = get_bits(gb, 8);
455             vui->matrix_coeffs           = get_bits(gb, 8);
456
457             // Set invalid values to "unspecified"
458             if (vui->colour_primaries >= AVCOL_PRI_NB)
459                 vui->colour_primaries = AVCOL_PRI_UNSPECIFIED;
460             if (vui->transfer_characteristic >= AVCOL_TRC_NB)
461                 vui->transfer_characteristic = AVCOL_TRC_UNSPECIFIED;
462             if (vui->matrix_coeffs >= AVCOL_SPC_NB)
463                 vui->matrix_coeffs = AVCOL_SPC_UNSPECIFIED;
464         }
465     }
466
467     vui->chroma_loc_info_present_flag = get_bits1(gb);
468     if (vui->chroma_loc_info_present_flag) {
469         vui->chroma_sample_loc_type_top_field    = get_ue_golomb_long(gb);
470         vui->chroma_sample_loc_type_bottom_field = get_ue_golomb_long(gb);
471     }
472
473     vui->neutra_chroma_indication_flag = get_bits1(gb);
474     vui->field_seq_flag                = get_bits1(gb);
475     vui->frame_field_info_present_flag = get_bits1(gb);
476
477     vui->default_display_window_flag = get_bits1(gb);
478     if (vui->default_display_window_flag) {
479         //TODO: * 2 is only valid for 420
480         vui->def_disp_win.left_offset   = get_ue_golomb_long(gb) * 2;
481         vui->def_disp_win.right_offset  = get_ue_golomb_long(gb) * 2;
482         vui->def_disp_win.top_offset    = get_ue_golomb_long(gb) * 2;
483         vui->def_disp_win.bottom_offset = get_ue_golomb_long(gb) * 2;
484
485         if (s->apply_defdispwin &&
486             s->avctx->flags2 & CODEC_FLAG2_IGNORE_CROP) {
487             av_log(s->avctx, AV_LOG_DEBUG,
488                    "discarding vui default display window, "
489                    "original values are l:%u r:%u t:%u b:%u\n",
490                    vui->def_disp_win.left_offset,
491                    vui->def_disp_win.right_offset,
492                    vui->def_disp_win.top_offset,
493                    vui->def_disp_win.bottom_offset);
494
495             vui->def_disp_win.left_offset   =
496             vui->def_disp_win.right_offset  =
497             vui->def_disp_win.top_offset    =
498             vui->def_disp_win.bottom_offset = 0;
499         }
500     }
501
502     vui->vui_timing_info_present_flag = get_bits1(gb);
503     if (vui->vui_timing_info_present_flag) {
504         vui->vui_num_units_in_tick               = get_bits(gb, 32);
505         vui->vui_time_scale                      = get_bits(gb, 32);
506         vui->vui_poc_proportional_to_timing_flag = get_bits1(gb);
507         if (vui->vui_poc_proportional_to_timing_flag)
508             vui->vui_num_ticks_poc_diff_one_minus1 = get_ue_golomb_long(gb);
509         vui->vui_hrd_parameters_present_flag = get_bits1(gb);
510         if (vui->vui_hrd_parameters_present_flag)
511             decode_hrd(s, 1, sps->max_sub_layers);
512     }
513
514     vui->bitstream_restriction_flag = get_bits1(gb);
515     if (vui->bitstream_restriction_flag) {
516         vui->tiles_fixed_structure_flag              = get_bits1(gb);
517         vui->motion_vectors_over_pic_boundaries_flag = get_bits1(gb);
518         vui->restricted_ref_pic_lists_flag           = get_bits1(gb);
519         vui->min_spatial_segmentation_idc            = get_ue_golomb_long(gb);
520         vui->max_bytes_per_pic_denom                 = get_ue_golomb_long(gb);
521         vui->max_bits_per_min_cu_denom               = get_ue_golomb_long(gb);
522         vui->log2_max_mv_length_horizontal           = get_ue_golomb_long(gb);
523         vui->log2_max_mv_length_vertical             = get_ue_golomb_long(gb);
524     }
525 }
526
527 static void set_default_scaling_list_data(ScalingList *sl)
528 {
529     int matrixId;
530
531     for (matrixId = 0; matrixId < 6; matrixId++) {
532         // 4x4 default is 16
533         memset(sl->sl[0][matrixId], 16, 16);
534         sl->sl_dc[0][matrixId] = 16; // default for 16x16
535         sl->sl_dc[1][matrixId] = 16; // default for 32x32
536     }
537     memcpy(sl->sl[1][0], default_scaling_list_intra, 64);
538     memcpy(sl->sl[1][1], default_scaling_list_intra, 64);
539     memcpy(sl->sl[1][2], default_scaling_list_intra, 64);
540     memcpy(sl->sl[1][3], default_scaling_list_inter, 64);
541     memcpy(sl->sl[1][4], default_scaling_list_inter, 64);
542     memcpy(sl->sl[1][5], default_scaling_list_inter, 64);
543     memcpy(sl->sl[2][0], default_scaling_list_intra, 64);
544     memcpy(sl->sl[2][1], default_scaling_list_intra, 64);
545     memcpy(sl->sl[2][2], default_scaling_list_intra, 64);
546     memcpy(sl->sl[2][3], default_scaling_list_inter, 64);
547     memcpy(sl->sl[2][4], default_scaling_list_inter, 64);
548     memcpy(sl->sl[2][5], default_scaling_list_inter, 64);
549     memcpy(sl->sl[3][0], default_scaling_list_intra, 64);
550     memcpy(sl->sl[3][1], default_scaling_list_inter, 64);
551 }
552
553 static int scaling_list_data(HEVCContext *s, ScalingList *sl)
554 {
555     GetBitContext *gb = &s->HEVClc.gb;
556     uint8_t scaling_list_pred_mode_flag[4][6];
557     int32_t scaling_list_dc_coef[2][6];
558     int size_id, matrix_id, i, pos, delta;
559
560     for (size_id = 0; size_id < 4; size_id++)
561         for (matrix_id = 0; matrix_id < (size_id == 3 ? 2 : 6); matrix_id++) {
562             scaling_list_pred_mode_flag[size_id][matrix_id] = get_bits1(gb);
563             if (!scaling_list_pred_mode_flag[size_id][matrix_id]) {
564                 delta = get_ue_golomb_long(gb);
565                 /* Only need to handle non-zero delta. Zero means default,
566                  * which should already be in the arrays. */
567                 if (delta) {
568                     // Copy from previous array.
569                     if (matrix_id - delta < 0) {
570                         av_log(s->avctx, AV_LOG_ERROR,
571                                "Invalid delta in scaling list data: %d.\n", delta);
572                         return AVERROR_INVALIDDATA;
573                     }
574
575                     memcpy(sl->sl[size_id][matrix_id],
576                            sl->sl[size_id][matrix_id - delta],
577                            size_id > 0 ? 64 : 16);
578                     if (size_id > 1)
579                         sl->sl_dc[size_id - 2][matrix_id] = sl->sl_dc[size_id - 2][matrix_id - delta];
580                 }
581             } else {
582                 int next_coef, coef_num;
583                 int32_t scaling_list_delta_coef;
584
585                 next_coef = 8;
586                 coef_num  = FFMIN(64, 1 << (4 + (size_id << 1)));
587                 if (size_id > 1) {
588                     scaling_list_dc_coef[size_id - 2][matrix_id] = get_se_golomb(gb) + 8;
589                     next_coef = scaling_list_dc_coef[size_id - 2][matrix_id];
590                     sl->sl_dc[size_id - 2][matrix_id] = next_coef;
591                 }
592                 for (i = 0; i < coef_num; i++) {
593                     if (size_id == 0)
594                         pos = 4 * ff_hevc_diag_scan4x4_y[i] +
595                                   ff_hevc_diag_scan4x4_x[i];
596                     else
597                         pos = 8 * ff_hevc_diag_scan8x8_y[i] +
598                                   ff_hevc_diag_scan8x8_x[i];
599
600                     scaling_list_delta_coef = get_se_golomb(gb);
601                     next_coef = (next_coef + scaling_list_delta_coef + 256) % 256;
602                     sl->sl[size_id][matrix_id][pos] = next_coef;
603                 }
604             }
605         }
606
607     return 0;
608 }
609
610 int ff_hevc_decode_nal_sps(HEVCContext *s)
611 {
612     const AVPixFmtDescriptor *desc;
613     GetBitContext *gb = &s->HEVClc.gb;
614     int ret    = 0;
615     int sps_id = 0;
616     int log2_diff_max_min_transform_block_size;
617     int bit_depth_chroma, start, vui_present, sublayer_ordering_info;
618     int i;
619
620     HEVCSPS *sps;
621     AVBufferRef *sps_buf = av_buffer_allocz(sizeof(*sps));
622
623     if (!sps_buf)
624         return AVERROR(ENOMEM);
625     sps = (HEVCSPS*)sps_buf->data;
626
627     av_log(s->avctx, AV_LOG_DEBUG, "Decoding SPS\n");
628
629     // Coded parameters
630
631     sps->vps_id = get_bits(gb, 4);
632     if (sps->vps_id >= MAX_VPS_COUNT) {
633         av_log(s->avctx, AV_LOG_ERROR, "VPS id out of range: %d\n", sps->vps_id);
634         ret = AVERROR_INVALIDDATA;
635         goto err;
636     }
637
638     sps->max_sub_layers = get_bits(gb, 3) + 1;
639     if (sps->max_sub_layers > MAX_SUB_LAYERS) {
640         av_log(s->avctx, AV_LOG_ERROR, "sps_max_sub_layers out of range: %d\n",
641                sps->max_sub_layers);
642         ret = AVERROR_INVALIDDATA;
643         goto err;
644     }
645
646     skip_bits1(gb); // temporal_id_nesting_flag
647     if (decode_profile_tier_level(&s->HEVClc, &sps->ptl,
648                                   sps->max_sub_layers) < 0) {
649         av_log(s->avctx, AV_LOG_ERROR, "error decoding profile tier level\n");
650         ret = AVERROR_INVALIDDATA;
651         goto err;
652     }
653     sps_id = get_ue_golomb_long(gb);
654     if (sps_id >= MAX_SPS_COUNT) {
655         av_log(s->avctx, AV_LOG_ERROR, "SPS id out of range: %d\n", sps_id);
656         ret = AVERROR_INVALIDDATA;
657         goto err;
658     }
659
660     sps->chroma_format_idc = get_ue_golomb_long(gb);
661     if (sps->chroma_format_idc != 1) {
662         avpriv_report_missing_feature(s->avctx, "chroma_format_idc != 1\n");
663         ret = AVERROR_PATCHWELCOME;
664         goto err;
665     }
666
667     if (sps->chroma_format_idc == 3)
668         sps->separate_colour_plane_flag = get_bits1(gb);
669
670     sps->width  = get_ue_golomb_long(gb);
671     sps->height = get_ue_golomb_long(gb);
672     if ((ret = av_image_check_size(sps->width,
673                                    sps->height, 0, s->avctx)) < 0)
674         goto err;
675
676     if (get_bits1(gb)) { // pic_conformance_flag
677         //TODO: * 2 is only valid for 420
678         sps->pic_conf_win.left_offset   = get_ue_golomb_long(gb) * 2;
679         sps->pic_conf_win.right_offset  = get_ue_golomb_long(gb) * 2;
680         sps->pic_conf_win.top_offset    = get_ue_golomb_long(gb) * 2;
681         sps->pic_conf_win.bottom_offset = get_ue_golomb_long(gb) * 2;
682
683         if (s->avctx->flags2 & CODEC_FLAG2_IGNORE_CROP) {
684             av_log(s->avctx, AV_LOG_DEBUG,
685                    "discarding sps conformance window, "
686                    "original values are l:%u r:%u t:%u b:%u\n",
687                    sps->pic_conf_win.left_offset,
688                    sps->pic_conf_win.right_offset,
689                    sps->pic_conf_win.top_offset,
690                    sps->pic_conf_win.bottom_offset);
691
692             sps->pic_conf_win.left_offset   =
693             sps->pic_conf_win.right_offset  =
694             sps->pic_conf_win.top_offset    =
695             sps->pic_conf_win.bottom_offset = 0;
696         }
697         sps->output_window = sps->pic_conf_win;
698     }
699
700     sps->bit_depth   = get_ue_golomb_long(gb) + 8;
701     bit_depth_chroma = get_ue_golomb_long(gb) + 8;
702     if (bit_depth_chroma != sps->bit_depth) {
703         av_log(s->avctx, AV_LOG_ERROR,
704                "Luma bit depth (%d) is different from chroma bit depth (%d), "
705                "this is unsupported.\n",
706                sps->bit_depth, bit_depth_chroma);
707         ret = AVERROR_INVALIDDATA;
708         goto err;
709     }
710
711     if (sps->chroma_format_idc == 1) {
712         switch (sps->bit_depth) {
713         case 8:  sps->pix_fmt = AV_PIX_FMT_YUV420P;   break;
714         case 9:  sps->pix_fmt = AV_PIX_FMT_YUV420P9;  break;
715         case 10: sps->pix_fmt = AV_PIX_FMT_YUV420P10; break;
716         default:
717             av_log(s->avctx, AV_LOG_ERROR, "Unsupported bit depth: %d\n",
718                    sps->bit_depth);
719             ret = AVERROR_PATCHWELCOME;
720             goto err;
721         }
722     } else {
723         av_log(s->avctx, AV_LOG_ERROR,
724                "non-4:2:0 support is currently unspecified.\n");
725         return AVERROR_PATCHWELCOME;
726     }
727
728     desc = av_pix_fmt_desc_get(sps->pix_fmt);
729     if (!desc) {
730         ret = AVERROR(EINVAL);
731         goto err;
732     }
733
734     sps->hshift[0] = sps->vshift[0] = 0;
735     sps->hshift[2] = sps->hshift[1] = desc->log2_chroma_w;
736     sps->vshift[2] = sps->vshift[1] = desc->log2_chroma_h;
737
738     sps->pixel_shift = sps->bit_depth > 8;
739
740     sps->log2_max_poc_lsb = get_ue_golomb_long(gb) + 4;
741     if (sps->log2_max_poc_lsb > 16) {
742         av_log(s->avctx, AV_LOG_ERROR, "log2_max_pic_order_cnt_lsb_minus4 out range: %d\n",
743                sps->log2_max_poc_lsb - 4);
744         ret = AVERROR_INVALIDDATA;
745         goto err;
746     }
747
748     sublayer_ordering_info = get_bits1(gb);
749     start = sublayer_ordering_info ? 0 : sps->max_sub_layers - 1;
750     for (i = start; i < sps->max_sub_layers; i++) {
751         sps->temporal_layer[i].max_dec_pic_buffering = get_ue_golomb_long(gb) + 1;
752         sps->temporal_layer[i].num_reorder_pics      = get_ue_golomb_long(gb);
753         sps->temporal_layer[i].max_latency_increase  = get_ue_golomb_long(gb) - 1;
754         if (sps->temporal_layer[i].max_dec_pic_buffering > MAX_DPB_SIZE) {
755             av_log(s->avctx, AV_LOG_ERROR, "sps_max_dec_pic_buffering_minus1 out of range: %d\n",
756                    sps->temporal_layer[i].max_dec_pic_buffering - 1);
757             ret = AVERROR_INVALIDDATA;
758             goto err;
759         }
760         if (sps->temporal_layer[i].num_reorder_pics > sps->temporal_layer[i].max_dec_pic_buffering - 1) {
761             av_log(s->avctx, AV_LOG_ERROR, "sps_max_num_reorder_pics out of range: %d\n",
762                    sps->temporal_layer[i].num_reorder_pics);
763             ret = AVERROR_INVALIDDATA;
764             goto err;
765         }
766     }
767
768     if (!sublayer_ordering_info) {
769         for (i = 0; i < start; i++) {
770             sps->temporal_layer[i].max_dec_pic_buffering = sps->temporal_layer[start].max_dec_pic_buffering;
771             sps->temporal_layer[i].num_reorder_pics      = sps->temporal_layer[start].num_reorder_pics;
772             sps->temporal_layer[i].max_latency_increase  = sps->temporal_layer[start].max_latency_increase;
773         }
774     }
775
776     sps->log2_min_cb_size                    = get_ue_golomb_long(gb) + 3;
777     sps->log2_diff_max_min_coding_block_size = get_ue_golomb_long(gb);
778     sps->log2_min_tb_size                    = get_ue_golomb_long(gb) + 2;
779     log2_diff_max_min_transform_block_size   = get_ue_golomb_long(gb);
780     sps->log2_max_trafo_size                 = log2_diff_max_min_transform_block_size +
781                                                sps->log2_min_tb_size;
782
783     if (sps->log2_min_tb_size >= sps->log2_min_cb_size) {
784         av_log(s->avctx, AV_LOG_ERROR, "Invalid value for log2_min_tb_size");
785         ret = AVERROR_INVALIDDATA;
786         goto err;
787     }
788     sps->max_transform_hierarchy_depth_inter = get_ue_golomb_long(gb);
789     sps->max_transform_hierarchy_depth_intra = get_ue_golomb_long(gb);
790
791     sps->scaling_list_enable_flag = get_bits1(gb);
792     if (sps->scaling_list_enable_flag) {
793         set_default_scaling_list_data(&sps->scaling_list);
794
795         if (get_bits1(gb)) {
796             ret = scaling_list_data(s, &sps->scaling_list);
797             if (ret < 0)
798                 goto err;
799         }
800     }
801
802     sps->amp_enabled_flag = get_bits1(gb);
803     sps->sao_enabled      = get_bits1(gb);
804
805     sps->pcm_enabled_flag = get_bits1(gb);
806     if (sps->pcm_enabled_flag) {
807         int pcm_bit_depth_chroma;
808         sps->pcm.bit_depth   = get_bits(gb, 4) + 1;
809         pcm_bit_depth_chroma = get_bits(gb, 4) + 1;
810         if (pcm_bit_depth_chroma != sps->pcm.bit_depth) {
811             av_log(s->avctx, AV_LOG_ERROR,
812                    "PCM Luma bit depth (%d) is different from PCM chroma"
813                    "bit depth (%d), this is unsupported.\n",
814                    sps->pcm.bit_depth, pcm_bit_depth_chroma);
815             ret = AVERROR_INVALIDDATA;
816             goto err;
817         }
818
819         sps->pcm.log2_min_pcm_cb_size = get_ue_golomb_long(gb) + 3;
820         sps->pcm.log2_max_pcm_cb_size = sps->pcm.log2_min_pcm_cb_size +
821                                         get_ue_golomb_long(gb);
822         if (sps->pcm.bit_depth > sps->bit_depth) {
823             av_log(s->avctx, AV_LOG_ERROR,
824                    "PCM bit depth (%d) is greater than normal bit depth (%d)\n",
825                    sps->pcm.bit_depth, sps->bit_depth);
826             ret = AVERROR_INVALIDDATA;
827             goto err;
828         }
829
830         sps->pcm.loop_filter_disable_flag = get_bits1(gb);
831     }
832
833     sps->nb_st_rps = get_ue_golomb_long(gb);
834     if (sps->nb_st_rps > MAX_SHORT_TERM_RPS_COUNT) {
835         av_log(s->avctx, AV_LOG_ERROR, "Too many short term RPS: %d.\n",
836                sps->nb_st_rps);
837         ret = AVERROR_INVALIDDATA;
838         goto err;
839     }
840     for (i = 0; i < sps->nb_st_rps; i++) {
841         if ((ret = ff_hevc_decode_short_term_rps(s, &sps->st_rps[i],
842                                                  sps, 0)) < 0)
843             goto err;
844     }
845
846     sps->long_term_ref_pics_present_flag = get_bits1(gb);
847     if (sps->long_term_ref_pics_present_flag) {
848         sps->num_long_term_ref_pics_sps = get_ue_golomb_long(gb);
849         for (i = 0; i < sps->num_long_term_ref_pics_sps; i++) {
850             sps->lt_ref_pic_poc_lsb_sps[i]       = get_bits(gb, sps->log2_max_poc_lsb);
851             sps->used_by_curr_pic_lt_sps_flag[i] = get_bits1(gb);
852         }
853     }
854
855     sps->sps_temporal_mvp_enabled_flag          = get_bits1(gb);
856     sps->sps_strong_intra_smoothing_enable_flag = get_bits1(gb);
857     sps->vui.sar = (AVRational){0, 1};
858     vui_present = get_bits1(gb);
859     if (vui_present)
860         decode_vui(s, sps);
861     skip_bits1(gb); // sps_extension_flag
862
863     if (s->apply_defdispwin) {
864         sps->output_window.left_offset   += sps->vui.def_disp_win.left_offset;
865         sps->output_window.right_offset  += sps->vui.def_disp_win.right_offset;
866         sps->output_window.top_offset    += sps->vui.def_disp_win.top_offset;
867         sps->output_window.bottom_offset += sps->vui.def_disp_win.bottom_offset;
868     }
869     if (sps->output_window.left_offset & (0x1F >> (sps->pixel_shift)) &&
870         !(s->avctx->flags & CODEC_FLAG_UNALIGNED)) {
871         sps->output_window.left_offset &= ~(0x1F >> (sps->pixel_shift));
872         av_log(s->avctx, AV_LOG_WARNING, "Reducing left output window to %d "
873                "chroma samples to preserve alignment.\n",
874                sps->output_window.left_offset);
875     }
876     sps->output_width  = sps->width -
877                          (sps->output_window.left_offset + sps->output_window.right_offset);
878     sps->output_height = sps->height -
879                          (sps->output_window.top_offset + sps->output_window.bottom_offset);
880     if (sps->output_width <= 0 || sps->output_height <= 0) {
881         av_log(s->avctx, AV_LOG_WARNING, "Invalid visible frame dimensions: %dx%d.\n",
882                sps->output_width, sps->output_height);
883         if (s->avctx->err_recognition & AV_EF_EXPLODE) {
884             ret = AVERROR_INVALIDDATA;
885             goto err;
886         }
887         av_log(s->avctx, AV_LOG_WARNING,
888                "Displaying the whole video surface.\n");
889         sps->pic_conf_win.left_offset   =
890         sps->pic_conf_win.right_offset  =
891         sps->pic_conf_win.top_offset    =
892         sps->pic_conf_win.bottom_offset = 0;
893         sps->output_width               = sps->width;
894         sps->output_height              = sps->height;
895     }
896
897     // Inferred parameters
898     sps->log2_ctb_size = sps->log2_min_cb_size +
899                          sps->log2_diff_max_min_coding_block_size;
900     sps->log2_min_pu_size = sps->log2_min_cb_size - 1;
901
902     sps->ctb_width  = (sps->width  + (1 << sps->log2_ctb_size) - 1) >> sps->log2_ctb_size;
903     sps->ctb_height = (sps->height + (1 << sps->log2_ctb_size) - 1) >> sps->log2_ctb_size;
904     sps->ctb_size   = sps->ctb_width * sps->ctb_height;
905
906     sps->min_cb_width  = sps->width  >> sps->log2_min_cb_size;
907     sps->min_cb_height = sps->height >> sps->log2_min_cb_size;
908     sps->min_tb_width  = sps->width  >> sps->log2_min_tb_size;
909     sps->min_tb_height = sps->height >> sps->log2_min_tb_size;
910     sps->min_pu_width  = sps->width  >> sps->log2_min_pu_size;
911     sps->min_pu_height = sps->height >> sps->log2_min_pu_size;
912
913     sps->qp_bd_offset = 6 * (sps->bit_depth - 8);
914
915     if (sps->width  & ((1 << sps->log2_min_cb_size) - 1) ||
916         sps->height & ((1 << sps->log2_min_cb_size) - 1)) {
917         av_log(s->avctx, AV_LOG_ERROR, "Invalid coded frame dimensions.\n");
918         goto err;
919     }
920
921     if (sps->log2_ctb_size > MAX_LOG2_CTB_SIZE) {
922         av_log(s->avctx, AV_LOG_ERROR, "CTB size out of range: 2^%d\n", sps->log2_ctb_size);
923         goto err;
924     }
925     if (sps->max_transform_hierarchy_depth_inter > sps->log2_ctb_size - sps->log2_min_tb_size) {
926         av_log(s->avctx, AV_LOG_ERROR, "max_transform_hierarchy_depth_inter out of range: %d\n",
927                sps->max_transform_hierarchy_depth_inter);
928         goto err;
929     }
930     if (sps->max_transform_hierarchy_depth_intra > sps->log2_ctb_size - sps->log2_min_tb_size) {
931         av_log(s->avctx, AV_LOG_ERROR, "max_transform_hierarchy_depth_intra out of range: %d\n",
932                sps->max_transform_hierarchy_depth_intra);
933         goto err;
934     }
935     if (sps->log2_max_trafo_size > FFMIN(sps->log2_ctb_size, 5)) {
936         av_log(s->avctx, AV_LOG_ERROR,
937                "max transform block size out of range: %d\n",
938                sps->log2_max_trafo_size);
939         goto err;
940     }
941
942     if (s->avctx->debug & FF_DEBUG_BITSTREAM) {
943         av_log(s->avctx, AV_LOG_DEBUG,
944                "Parsed SPS: id %d; coded wxh: %dx%d; "
945                "cropped wxh: %dx%d; pix_fmt: %s.\n",
946                sps_id, sps->width, sps->height,
947                sps->output_width, sps->output_height,
948                av_get_pix_fmt_name(sps->pix_fmt));
949     }
950
951     /* check if this is a repeat of an already parsed SPS, then keep the
952      * original one.
953      * otherwise drop all PPSes that depend on it */
954     if (s->sps_list[sps_id] &&
955         !memcmp(s->sps_list[sps_id]->data, sps_buf->data, sps_buf->size)) {
956         av_buffer_unref(&sps_buf);
957     } else {
958         for (i = 0; i < FF_ARRAY_ELEMS(s->pps_list); i++) {
959             if (s->pps_list[i] && ((HEVCPPS*)s->pps_list[i]->data)->sps_id == sps_id)
960                 av_buffer_unref(&s->pps_list[i]);
961         }
962         av_buffer_unref(&s->sps_list[sps_id]);
963         s->sps_list[sps_id] = sps_buf;
964     }
965
966     return 0;
967
968 err:
969     av_buffer_unref(&sps_buf);
970     return ret;
971 }
972
973 static void hevc_pps_free(void *opaque, uint8_t *data)
974 {
975     HEVCPPS *pps = (HEVCPPS*)data;
976
977     av_freep(&pps->column_width);
978     av_freep(&pps->row_height);
979     av_freep(&pps->col_bd);
980     av_freep(&pps->row_bd);
981     av_freep(&pps->col_idxX);
982     av_freep(&pps->ctb_addr_rs_to_ts);
983     av_freep(&pps->ctb_addr_ts_to_rs);
984     av_freep(&pps->tile_pos_rs);
985     av_freep(&pps->tile_id);
986     av_freep(&pps->min_cb_addr_zs);
987     av_freep(&pps->min_tb_addr_zs);
988
989     av_freep(&pps);
990 }
991
992 int ff_hevc_decode_nal_pps(HEVCContext *s)
993 {
994     GetBitContext *gb = &s->HEVClc.gb;
995     HEVCSPS      *sps = NULL;
996     int pic_area_in_ctbs, pic_area_in_min_cbs, pic_area_in_min_tbs;
997     int log2_diff_ctb_min_tb_size;
998     int i, j, x, y, ctb_addr_rs, tile_id;
999     int ret    = 0;
1000     int pps_id = 0;
1001
1002     AVBufferRef *pps_buf;
1003     HEVCPPS *pps = av_mallocz(sizeof(*pps));
1004
1005     if (!pps)
1006         return AVERROR(ENOMEM);
1007
1008     pps_buf = av_buffer_create((uint8_t *)pps, sizeof(*pps),
1009                                hevc_pps_free, NULL, 0);
1010     if (!pps_buf) {
1011         av_freep(&pps);
1012         return AVERROR(ENOMEM);
1013     }
1014
1015     av_log(s->avctx, AV_LOG_DEBUG, "Decoding PPS\n");
1016
1017     // Default values
1018     pps->loop_filter_across_tiles_enabled_flag = 1;
1019     pps->num_tile_columns                      = 1;
1020     pps->num_tile_rows                         = 1;
1021     pps->uniform_spacing_flag                  = 1;
1022     pps->disable_dbf                           = 0;
1023     pps->beta_offset                           = 0;
1024     pps->tc_offset                             = 0;
1025
1026     // Coded parameters
1027     pps_id = get_ue_golomb_long(gb);
1028     if (pps_id >= MAX_PPS_COUNT) {
1029         av_log(s->avctx, AV_LOG_ERROR, "PPS id out of range: %d\n", pps_id);
1030         ret = AVERROR_INVALIDDATA;
1031         goto err;
1032     }
1033     pps->sps_id = get_ue_golomb_long(gb);
1034     if (pps->sps_id >= MAX_SPS_COUNT) {
1035         av_log(s->avctx, AV_LOG_ERROR, "SPS id out of range: %d\n", pps->sps_id);
1036         ret = AVERROR_INVALIDDATA;
1037         goto err;
1038     }
1039     if (!s->sps_list[pps->sps_id]) {
1040         av_log(s->avctx, AV_LOG_ERROR, "SPS does not exist \n");
1041         ret = AVERROR_INVALIDDATA;
1042         goto err;
1043     }
1044     sps = (HEVCSPS *)s->sps_list[pps->sps_id]->data;
1045
1046     pps->dependent_slice_segments_enabled_flag = get_bits1(gb);
1047     pps->output_flag_present_flag              = get_bits1(gb);
1048     pps->num_extra_slice_header_bits           = get_bits(gb, 3);
1049
1050     pps->sign_data_hiding_flag = get_bits1(gb);
1051
1052     pps->cabac_init_present_flag = get_bits1(gb);
1053
1054     pps->num_ref_idx_l0_default_active = get_ue_golomb_long(gb) + 1;
1055     pps->num_ref_idx_l1_default_active = get_ue_golomb_long(gb) + 1;
1056
1057     pps->pic_init_qp_minus26 = get_se_golomb(gb);
1058
1059     pps->constrained_intra_pred_flag = get_bits1(gb);
1060     pps->transform_skip_enabled_flag = get_bits1(gb);
1061
1062     pps->cu_qp_delta_enabled_flag = get_bits1(gb);
1063     pps->diff_cu_qp_delta_depth   = 0;
1064     if (pps->cu_qp_delta_enabled_flag)
1065         pps->diff_cu_qp_delta_depth = get_ue_golomb_long(gb);
1066
1067     pps->cb_qp_offset = get_se_golomb(gb);
1068     if (pps->cb_qp_offset < -12 || pps->cb_qp_offset > 12) {
1069         av_log(s->avctx, AV_LOG_ERROR, "pps_cb_qp_offset out of range: %d\n",
1070                pps->cb_qp_offset);
1071         ret = AVERROR_INVALIDDATA;
1072         goto err;
1073     }
1074     pps->cr_qp_offset = get_se_golomb(gb);
1075     if (pps->cr_qp_offset < -12 || pps->cr_qp_offset > 12) {
1076         av_log(s->avctx, AV_LOG_ERROR, "pps_cr_qp_offset out of range: %d\n",
1077                pps->cr_qp_offset);
1078         ret = AVERROR_INVALIDDATA;
1079         goto err;
1080     }
1081     pps->pic_slice_level_chroma_qp_offsets_present_flag = get_bits1(gb);
1082
1083     pps->weighted_pred_flag   = get_bits1(gb);
1084     pps->weighted_bipred_flag = get_bits1(gb);
1085
1086     pps->transquant_bypass_enable_flag    = get_bits1(gb);
1087     pps->tiles_enabled_flag               = get_bits1(gb);
1088     pps->entropy_coding_sync_enabled_flag = get_bits1(gb);
1089
1090     if (pps->tiles_enabled_flag) {
1091         pps->num_tile_columns = get_ue_golomb_long(gb) + 1;
1092         pps->num_tile_rows    = get_ue_golomb_long(gb) + 1;
1093         if (pps->num_tile_columns == 0 ||
1094             pps->num_tile_columns >= sps->width) {
1095             av_log(s->avctx, AV_LOG_ERROR, "num_tile_columns_minus1 out of range: %d\n",
1096                    pps->num_tile_columns - 1);
1097             ret = AVERROR_INVALIDDATA;
1098             goto err;
1099         }
1100         if (pps->num_tile_rows == 0 ||
1101             pps->num_tile_rows >= sps->height) {
1102             av_log(s->avctx, AV_LOG_ERROR, "num_tile_rows_minus1 out of range: %d\n",
1103                    pps->num_tile_rows - 1);
1104             ret = AVERROR_INVALIDDATA;
1105             goto err;
1106         }
1107
1108         pps->column_width = av_malloc_array(pps->num_tile_columns, sizeof(*pps->column_width));
1109         pps->row_height   = av_malloc_array(pps->num_tile_rows,    sizeof(*pps->row_height));
1110         if (!pps->column_width || !pps->row_height) {
1111             ret = AVERROR(ENOMEM);
1112             goto err;
1113         }
1114
1115         pps->uniform_spacing_flag = get_bits1(gb);
1116         if (!pps->uniform_spacing_flag) {
1117             int sum = 0;
1118             for (i = 0; i < pps->num_tile_columns - 1; i++) {
1119                 pps->column_width[i] = get_ue_golomb_long(gb) + 1;
1120                 sum                 += pps->column_width[i];
1121             }
1122             if (sum >= sps->ctb_width) {
1123                 av_log(s->avctx, AV_LOG_ERROR, "Invalid tile widths.\n");
1124                 ret = AVERROR_INVALIDDATA;
1125                 goto err;
1126             }
1127             pps->column_width[pps->num_tile_columns - 1] = sps->ctb_width - sum;
1128
1129             sum = 0;
1130             for (i = 0; i < pps->num_tile_rows - 1; i++) {
1131                 pps->row_height[i] = get_ue_golomb_long(gb) + 1;
1132                 sum               += pps->row_height[i];
1133             }
1134             if (sum >= sps->ctb_height) {
1135                 av_log(s->avctx, AV_LOG_ERROR, "Invalid tile heights.\n");
1136                 ret = AVERROR_INVALIDDATA;
1137                 goto err;
1138             }
1139             pps->row_height[pps->num_tile_rows - 1] = sps->ctb_height - sum;
1140         }
1141         pps->loop_filter_across_tiles_enabled_flag = get_bits1(gb);
1142     }
1143
1144     pps->seq_loop_filter_across_slices_enabled_flag = get_bits1(gb);
1145
1146     pps->deblocking_filter_control_present_flag = get_bits1(gb);
1147     if (pps->deblocking_filter_control_present_flag) {
1148         pps->deblocking_filter_override_enabled_flag = get_bits1(gb);
1149         pps->disable_dbf                             = get_bits1(gb);
1150         if (!pps->disable_dbf) {
1151             pps->beta_offset = get_se_golomb(gb) * 2;
1152             pps->tc_offset = get_se_golomb(gb) * 2;
1153             if (pps->beta_offset/2 < -6 || pps->beta_offset/2 > 6) {
1154                 av_log(s->avctx, AV_LOG_ERROR, "pps_beta_offset_div2 out of range: %d\n",
1155                        pps->beta_offset/2);
1156                 ret = AVERROR_INVALIDDATA;
1157                 goto err;
1158             }
1159             if (pps->tc_offset/2 < -6 || pps->tc_offset/2 > 6) {
1160                 av_log(s->avctx, AV_LOG_ERROR, "pps_tc_offset_div2 out of range: %d\n",
1161                        pps->tc_offset/2);
1162                 ret = AVERROR_INVALIDDATA;
1163                 goto err;
1164             }
1165         }
1166     }
1167
1168     pps->scaling_list_data_present_flag = get_bits1(gb);
1169     if (pps->scaling_list_data_present_flag) {
1170         set_default_scaling_list_data(&pps->scaling_list);
1171         ret = scaling_list_data(s, &pps->scaling_list);
1172         if (ret < 0)
1173             goto err;
1174     }
1175     pps->lists_modification_present_flag = get_bits1(gb);
1176     pps->log2_parallel_merge_level       = get_ue_golomb_long(gb) + 2;
1177     if (pps->log2_parallel_merge_level > sps->log2_ctb_size) {
1178         av_log(s->avctx, AV_LOG_ERROR, "log2_parallel_merge_level_minus2 out of range: %d\n",
1179                pps->log2_parallel_merge_level - 2);
1180         ret = AVERROR_INVALIDDATA;
1181         goto err;
1182     }
1183
1184     pps->slice_header_extension_present_flag = get_bits1(gb);
1185     skip_bits1(gb);     // pps_extension_flag
1186
1187     // Inferred parameters
1188     pps->col_bd   = av_malloc_array(pps->num_tile_columns + 1, sizeof(*pps->col_bd));
1189     pps->row_bd   = av_malloc_array(pps->num_tile_rows + 1,    sizeof(*pps->row_bd));
1190     pps->col_idxX = av_malloc_array(sps->ctb_width,    sizeof(*pps->col_idxX));
1191     if (!pps->col_bd || !pps->row_bd || !pps->col_idxX) {
1192         ret = AVERROR(ENOMEM);
1193         goto err;
1194     }
1195
1196     if (pps->uniform_spacing_flag) {
1197         if (!pps->column_width) {
1198             pps->column_width = av_malloc_array(pps->num_tile_columns, sizeof(*pps->column_width));
1199             pps->row_height   = av_malloc_array(pps->num_tile_rows,    sizeof(*pps->row_height));
1200         }
1201         if (!pps->column_width || !pps->row_height) {
1202             ret = AVERROR(ENOMEM);
1203             goto err;
1204         }
1205
1206         for (i = 0; i < pps->num_tile_columns; i++) {
1207             pps->column_width[i] = ((i + 1) * sps->ctb_width) / pps->num_tile_columns -
1208                                    (i * sps->ctb_width) / pps->num_tile_columns;
1209         }
1210
1211         for (i = 0; i < pps->num_tile_rows; i++) {
1212             pps->row_height[i] = ((i + 1) * sps->ctb_height) / pps->num_tile_rows -
1213                                  (i * sps->ctb_height) / pps->num_tile_rows;
1214         }
1215     }
1216
1217     pps->col_bd[0] = 0;
1218     for (i = 0; i < pps->num_tile_columns; i++)
1219         pps->col_bd[i + 1] = pps->col_bd[i] + pps->column_width[i];
1220
1221     pps->row_bd[0] = 0;
1222     for (i = 0; i < pps->num_tile_rows; i++)
1223         pps->row_bd[i + 1] = pps->row_bd[i] + pps->row_height[i];
1224
1225     for (i = 0, j = 0; i < sps->ctb_width; i++) {
1226         if (i > pps->col_bd[j])
1227             j++;
1228         pps->col_idxX[i] = j;
1229     }
1230
1231     /**
1232      * 6.5
1233      */
1234     pic_area_in_ctbs     = sps->ctb_width    * sps->ctb_height;
1235     pic_area_in_min_cbs  = sps->min_cb_width * sps->min_cb_height;
1236     pic_area_in_min_tbs  = sps->min_tb_width * sps->min_tb_height;
1237
1238     pps->ctb_addr_rs_to_ts = av_malloc_array(pic_area_in_ctbs,    sizeof(*pps->ctb_addr_rs_to_ts));
1239     pps->ctb_addr_ts_to_rs = av_malloc_array(pic_area_in_ctbs,    sizeof(*pps->ctb_addr_ts_to_rs));
1240     pps->tile_id           = av_malloc_array(pic_area_in_ctbs,    sizeof(*pps->tile_id));
1241     pps->min_cb_addr_zs    = av_malloc_array(pic_area_in_min_cbs, sizeof(*pps->min_cb_addr_zs));
1242     pps->min_tb_addr_zs    = av_malloc_array(pic_area_in_min_tbs, sizeof(*pps->min_tb_addr_zs));
1243     if (!pps->ctb_addr_rs_to_ts || !pps->ctb_addr_ts_to_rs ||
1244         !pps->tile_id || !pps->min_cb_addr_zs || !pps->min_tb_addr_zs) {
1245         ret = AVERROR(ENOMEM);
1246         goto err;
1247     }
1248
1249     for (ctb_addr_rs = 0; ctb_addr_rs < pic_area_in_ctbs; ctb_addr_rs++) {
1250         int tb_x   = ctb_addr_rs % sps->ctb_width;
1251         int tb_y   = ctb_addr_rs / sps->ctb_width;
1252         int tile_x = 0;
1253         int tile_y = 0;
1254         int val    = 0;
1255
1256         for (i = 0; i < pps->num_tile_columns; i++) {
1257             if (tb_x < pps->col_bd[i + 1]) {
1258                 tile_x = i;
1259                 break;
1260             }
1261         }
1262
1263         for (i = 0; i < pps->num_tile_rows; i++) {
1264             if (tb_y < pps->row_bd[i + 1]) {
1265                 tile_y = i;
1266                 break;
1267             }
1268         }
1269
1270         for (i = 0; i < tile_x; i++)
1271             val += pps->row_height[tile_y] * pps->column_width[i];
1272         for (i = 0; i < tile_y; i++)
1273             val += sps->ctb_width * pps->row_height[i];
1274
1275         val += (tb_y - pps->row_bd[tile_y]) * pps->column_width[tile_x] +
1276                tb_x - pps->col_bd[tile_x];
1277
1278         pps->ctb_addr_rs_to_ts[ctb_addr_rs] = val;
1279         pps->ctb_addr_ts_to_rs[val]         = ctb_addr_rs;
1280     }
1281
1282     for (j = 0, tile_id = 0; j < pps->num_tile_rows; j++)
1283         for (i = 0; i < pps->num_tile_columns; i++, tile_id++)
1284             for (y = pps->row_bd[j]; y < pps->row_bd[j + 1]; y++)
1285                 for (x = pps->col_bd[i]; x < pps->col_bd[i + 1]; x++)
1286                     pps->tile_id[pps->ctb_addr_rs_to_ts[y * sps->ctb_width + x]] = tile_id;
1287
1288     pps->tile_pos_rs = av_malloc_array(tile_id, sizeof(*pps->tile_pos_rs));
1289     if (!pps->tile_pos_rs) {
1290         ret = AVERROR(ENOMEM);
1291         goto err;
1292     }
1293
1294     for (j = 0; j < pps->num_tile_rows; j++)
1295         for (i = 0; i < pps->num_tile_columns; i++)
1296             pps->tile_pos_rs[j * pps->num_tile_columns + i] = pps->row_bd[j] * sps->ctb_width + pps->col_bd[i];
1297
1298     for (y = 0; y < sps->min_cb_height; y++) {
1299         for (x = 0; x < sps->min_cb_width; x++) {
1300             int tb_x        = x >> sps->log2_diff_max_min_coding_block_size;
1301             int tb_y        = y >> sps->log2_diff_max_min_coding_block_size;
1302             int ctb_addr_rs = sps->ctb_width * tb_y + tb_x;
1303             int val         = pps->ctb_addr_rs_to_ts[ctb_addr_rs] <<
1304                               (sps->log2_diff_max_min_coding_block_size * 2);
1305             for (i = 0; i < sps->log2_diff_max_min_coding_block_size; i++) {
1306                 int m = 1 << i;
1307                 val += (m & x ? m * m : 0) + (m & y ? 2 * m * m : 0);
1308             }
1309             pps->min_cb_addr_zs[y * sps->min_cb_width + x] = val;
1310         }
1311     }
1312
1313     log2_diff_ctb_min_tb_size = sps->log2_ctb_size - sps->log2_min_tb_size;
1314     for (y = 0; y < sps->min_tb_height; y++) {
1315         for (x = 0; x < sps->min_tb_width; x++) {
1316             int tb_x        = x >> log2_diff_ctb_min_tb_size;
1317             int tb_y        = y >> log2_diff_ctb_min_tb_size;
1318             int ctb_addr_rs = sps->ctb_width * tb_y + tb_x;
1319             int val         = pps->ctb_addr_rs_to_ts[ctb_addr_rs] <<
1320                               (log2_diff_ctb_min_tb_size * 2);
1321             for (i = 0; i < log2_diff_ctb_min_tb_size; i++) {
1322                 int m = 1 << i;
1323                 val += (m & x ? m * m : 0) + (m & y ? 2 * m * m : 0);
1324             }
1325             pps->min_tb_addr_zs[y * sps->min_tb_width + x] = val;
1326         }
1327     }
1328
1329     av_buffer_unref(&s->pps_list[pps_id]);
1330     s->pps_list[pps_id] = pps_buf;
1331
1332     return 0;
1333
1334 err:
1335     av_buffer_unref(&pps_buf);
1336     return ret;
1337 }