Statistics
| Branch: | Revision:

ffmpeg / libavcodec / vp56.c @ 68a4d349

History | View | Annotate | Download (23 KB)

1 5ce117c3 Aurelien Jacobs
/**
2 bad5537e Diego Biurrun
 * @file libavcodec/vp56.c
3 5ce117c3 Aurelien Jacobs
 * VP5 and VP6 compatible video decoder (common features)
4
 *
5
 * Copyright (C) 2006  Aurelien Jacobs <aurel@gnuage.org>
6
 *
7 b78e7197 Diego Biurrun
 * This file is part of FFmpeg.
8
 *
9
 * FFmpeg is free software; you can redistribute it and/or
10 5ce117c3 Aurelien Jacobs
 * 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 b78e7197 Diego Biurrun
 * FFmpeg is distributed in the hope that it will be useful,
15 5ce117c3 Aurelien Jacobs
 * 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 b78e7197 Diego Biurrun
 * License along with FFmpeg; if not, write to the Free Software
21 e5a389a1 Diego Biurrun
 * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
22 5ce117c3 Aurelien Jacobs
 */
23
24
#include "avcodec.h"
25 91fc2cf1 Aurelien Jacobs
#include "bytestream.h"
26 5ce117c3 Aurelien Jacobs
27
#include "vp56.h"
28
#include "vp56data.h"
29
30
31 3d52bca6 Aurelien Jacobs
void vp56_init_dequant(VP56Context *s, int quantizer)
32 5ce117c3 Aurelien Jacobs
{
33
    s->quantizer = quantizer;
34
    s->dequant_dc = vp56_dc_dequant[quantizer] << 2;
35
    s->dequant_ac = vp56_ac_dequant[quantizer] << 2;
36
}
37
38 3d52bca6 Aurelien Jacobs
static int vp56_get_vectors_predictors(VP56Context *s, int row, int col,
39
                                       VP56Frame ref_frame)
40 5ce117c3 Aurelien Jacobs
{
41
    int nb_pred = 0;
42 3d52bca6 Aurelien Jacobs
    VP56mv vect[2] = {{0,0}, {0,0}};
43 5ce117c3 Aurelien Jacobs
    int pos, offset;
44 3d52bca6 Aurelien Jacobs
    VP56mv mvp;
45 5ce117c3 Aurelien Jacobs
46
    for (pos=0; pos<12; pos++) {
47
        mvp.x = col + vp56_candidate_predictor_pos[pos][0];
48
        mvp.y = row + vp56_candidate_predictor_pos[pos][1];
49
        if (mvp.x < 0 || mvp.x >= s->mb_width ||
50
            mvp.y < 0 || mvp.y >= s->mb_height)
51
            continue;
52
        offset = mvp.x + s->mb_width*mvp.y;
53
54
        if (vp56_reference_frame[s->macroblocks[offset].type] != ref_frame)
55
            continue;
56 d120e402 Aurelien Jacobs
        if ((s->macroblocks[offset].mv.x == vect[0].x &&
57
             s->macroblocks[offset].mv.y == vect[0].y) ||
58 5ce117c3 Aurelien Jacobs
            (s->macroblocks[offset].mv.x == 0 &&
59
             s->macroblocks[offset].mv.y == 0))
60
            continue;
61
62 d120e402 Aurelien Jacobs
        vect[nb_pred++] = s->macroblocks[offset].mv;
63 5ce117c3 Aurelien Jacobs
        if (nb_pred > 1) {
64
            nb_pred = -1;
65
            break;
66
        }
67
        s->vector_candidate_pos = pos;
68
    }
69
70 d120e402 Aurelien Jacobs
    s->vector_candidate[0] = vect[0];
71
    s->vector_candidate[1] = vect[1];
72 5ce117c3 Aurelien Jacobs
73
    return nb_pred+1;
74
}
75
76 3d52bca6 Aurelien Jacobs
static void vp56_parse_mb_type_models(VP56Context *s)
77 5ce117c3 Aurelien Jacobs
{
78 3d52bca6 Aurelien Jacobs
    VP56RangeCoder *c = &s->c;
79 d887151d Aurelien Jacobs
    VP56Model *model = s->modelp;
80 5ce117c3 Aurelien Jacobs
    int i, ctx, type;
81
82
    for (ctx=0; ctx<3; ctx++) {
83
        if (vp56_rac_get_prob(c, 174)) {
84
            int idx = vp56_rac_gets(c, 4);
85 247df384 Aurelien Jacobs
            memcpy(model->mb_types_stats[ctx],
86
                   vp56_pre_def_mb_type_stats[idx][ctx],
87
                   sizeof(model->mb_types_stats[ctx]));
88 5ce117c3 Aurelien Jacobs
        }
89
        if (vp56_rac_get_prob(c, 254)) {
90
            for (type=0; type<10; type++) {
91
                for(i=0; i<2; i++) {
92
                    if (vp56_rac_get_prob(c, 205)) {
93
                        int delta, sign = vp56_rac_get(c);
94
95
                        delta = vp56_rac_get_tree(c, vp56_pmbtm_tree,
96
                                                  vp56_mb_type_model_model);
97
                        if (!delta)
98
                            delta = 4 * vp56_rac_gets(c, 7);
99 247df384 Aurelien Jacobs
                        model->mb_types_stats[ctx][type][i] += (delta ^ -sign) + sign;
100 5ce117c3 Aurelien Jacobs
                    }
101
                }
102
            }
103
        }
104
    }
105
106
    /* compute MB type probability tables based on previous MB type */
107
    for (ctx=0; ctx<3; ctx++) {
108
        int p[10];
109
110
        for (type=0; type<10; type++)
111 247df384 Aurelien Jacobs
            p[type] = 100 * model->mb_types_stats[ctx][type][1];
112 5ce117c3 Aurelien Jacobs
113
        for (type=0; type<10; type++) {
114
            int p02, p34, p0234, p17, p56, p89, p5689, p156789;
115
116
            /* conservative MB type probability */
117 247df384 Aurelien Jacobs
            model->mb_type[ctx][type][0] = 255 - (255 * model->mb_types_stats[ctx][type][0]) / (1 + model->mb_types_stats[ctx][type][0] + model->mb_types_stats[ctx][type][1]);
118 5ce117c3 Aurelien Jacobs
119
            p[type] = 0;    /* same MB type => weight is null */
120
121
            /* binary tree parsing probabilities */
122
            p02 = p[0] + p[2];
123
            p34 = p[3] + p[4];
124
            p0234 = p02 + p34;
125
            p17 = p[1] + p[7];
126
            p56 = p[5] + p[6];
127
            p89 = p[8] + p[9];
128
            p5689 = p56 + p89;
129
            p156789 = p17 + p5689;
130
131 247df384 Aurelien Jacobs
            model->mb_type[ctx][type][1] = 1 + 255 * p0234/(1+p0234+p156789);
132
            model->mb_type[ctx][type][2] = 1 + 255 * p02  / (1+p0234);
133
            model->mb_type[ctx][type][3] = 1 + 255 * p17  / (1+p156789);
134
            model->mb_type[ctx][type][4] = 1 + 255 * p[0] / (1+p02);
135
            model->mb_type[ctx][type][5] = 1 + 255 * p[3] / (1+p34);
136
            model->mb_type[ctx][type][6] = 1 + 255 * p[1] / (1+p17);
137
            model->mb_type[ctx][type][7] = 1 + 255 * p56  / (1+p5689);
138
            model->mb_type[ctx][type][8] = 1 + 255 * p[5] / (1+p56);
139
            model->mb_type[ctx][type][9] = 1 + 255 * p[8] / (1+p89);
140 5ce117c3 Aurelien Jacobs
141
            /* restore initial value */
142 247df384 Aurelien Jacobs
            p[type] = 100 * model->mb_types_stats[ctx][type][1];
143 5ce117c3 Aurelien Jacobs
        }
144
    }
145
}
146
147 3d52bca6 Aurelien Jacobs
static VP56mb vp56_parse_mb_type(VP56Context *s,
148 76025d91 Aurelien Jacobs
                                 VP56mb prev_type, int ctx)
149 5ce117c3 Aurelien Jacobs
{
150 247df384 Aurelien Jacobs
    uint8_t *mb_type_model = s->modelp->mb_type[ctx][prev_type];
151 3d52bca6 Aurelien Jacobs
    VP56RangeCoder *c = &s->c;
152 5ce117c3 Aurelien Jacobs
153
    if (vp56_rac_get_prob(c, mb_type_model[0]))
154
        return prev_type;
155
    else
156
        return vp56_rac_get_tree(c, vp56_pmbt_tree, mb_type_model);
157
}
158
159 3d52bca6 Aurelien Jacobs
static void vp56_decode_4mv(VP56Context *s, int row, int col)
160 5ce117c3 Aurelien Jacobs
{
161 3d52bca6 Aurelien Jacobs
    VP56mv mv = {0,0};
162 5ce117c3 Aurelien Jacobs
    int type[4];
163
    int b;
164
165
    /* parse each block type */
166
    for (b=0; b<4; b++) {
167
        type[b] = vp56_rac_gets(&s->c, 2);
168
        if (type[b])
169
            type[b]++;  /* only returns 0, 2, 3 or 4 (all INTER_PF) */
170
    }
171
172
    /* get vectors */
173
    for (b=0; b<4; b++) {
174
        switch (type[b]) {
175
            case VP56_MB_INTER_NOVEC_PF:
176 3d52bca6 Aurelien Jacobs
                s->mv[b] = (VP56mv) {0,0};
177 5ce117c3 Aurelien Jacobs
                break;
178
            case VP56_MB_INTER_DELTA_PF:
179
                s->parse_vector_adjustment(s, &s->mv[b]);
180
                break;
181
            case VP56_MB_INTER_V1_PF:
182
                s->mv[b] = s->vector_candidate[0];
183
                break;
184
            case VP56_MB_INTER_V2_PF:
185
                s->mv[b] = s->vector_candidate[1];
186
                break;
187
        }
188
        mv.x += s->mv[b].x;
189
        mv.y += s->mv[b].y;
190
    }
191
192
    /* this is the one selected for the whole MB for prediction */
193
    s->macroblocks[row * s->mb_width + col].mv = s->mv[3];
194
195
    /* chroma vectors are average luma vectors */
196
    if (s->avctx->codec->id == CODEC_ID_VP5) {
197
        s->mv[4].x = s->mv[5].x = RSHIFT(mv.x,2);
198
        s->mv[4].y = s->mv[5].y = RSHIFT(mv.y,2);
199
    } else {
200 3d52bca6 Aurelien Jacobs
        s->mv[4] = s->mv[5] = (VP56mv) {mv.x/4, mv.y/4};
201 5ce117c3 Aurelien Jacobs
    }
202
}
203
204 3d52bca6 Aurelien Jacobs
static VP56mb vp56_decode_mv(VP56Context *s, int row, int col)
205 5ce117c3 Aurelien Jacobs
{
206 3d52bca6 Aurelien Jacobs
    VP56mv *mv, vect = {0,0};
207 5ce117c3 Aurelien Jacobs
    int ctx, b;
208
209
    ctx = vp56_get_vectors_predictors(s, row, col, VP56_FRAME_PREVIOUS);
210
    s->mb_type = vp56_parse_mb_type(s, s->mb_type, ctx);
211
    s->macroblocks[row * s->mb_width + col].type = s->mb_type;
212
213
    switch (s->mb_type) {
214
        case VP56_MB_INTER_V1_PF:
215
            mv = &s->vector_candidate[0];
216
            break;
217
218
        case VP56_MB_INTER_V2_PF:
219
            mv = &s->vector_candidate[1];
220
            break;
221
222
        case VP56_MB_INTER_V1_GF:
223
            vp56_get_vectors_predictors(s, row, col, VP56_FRAME_GOLDEN);
224
            mv = &s->vector_candidate[0];
225
            break;
226
227
        case VP56_MB_INTER_V2_GF:
228
            vp56_get_vectors_predictors(s, row, col, VP56_FRAME_GOLDEN);
229
            mv = &s->vector_candidate[1];
230
            break;
231
232
        case VP56_MB_INTER_DELTA_PF:
233 d120e402 Aurelien Jacobs
            s->parse_vector_adjustment(s, &vect);
234
            mv = &vect;
235 5ce117c3 Aurelien Jacobs
            break;
236
237
        case VP56_MB_INTER_DELTA_GF:
238
            vp56_get_vectors_predictors(s, row, col, VP56_FRAME_GOLDEN);
239 d120e402 Aurelien Jacobs
            s->parse_vector_adjustment(s, &vect);
240
            mv = &vect;
241 5ce117c3 Aurelien Jacobs
            break;
242
243
        case VP56_MB_INTER_4V:
244
            vp56_decode_4mv(s, row, col);
245
            return s->mb_type;
246
247
        default:
248 d120e402 Aurelien Jacobs
            mv = &vect;
249 5ce117c3 Aurelien Jacobs
            break;
250
    }
251
252
    s->macroblocks[row*s->mb_width + col].mv = *mv;
253
254
    /* same vector for all blocks */
255
    for (b=0; b<6; b++)
256
        s->mv[b] = *mv;
257
258
    return s->mb_type;
259
}
260
261 3d52bca6 Aurelien Jacobs
static void vp56_add_predictors_dc(VP56Context *s, VP56Frame ref_frame)
262 5ce117c3 Aurelien Jacobs
{
263
    int idx = s->scantable.permutated[0];
264 d7af6a9d Aurelien Jacobs
    int b;
265 5ce117c3 Aurelien Jacobs
266 d7af6a9d Aurelien Jacobs
    for (b=0; b<6; b++) {
267 3d52bca6 Aurelien Jacobs
        VP56RefDc *ab = &s->above_blocks[s->above_block_idx[b]];
268
        VP56RefDc *lb = &s->left_block[vp56_b6to4[b]];
269 5ce117c3 Aurelien Jacobs
        int count = 0;
270
        int dc = 0;
271 7ecae905 Aurelien Jacobs
        int i;
272 5ce117c3 Aurelien Jacobs
273
        if (ref_frame == lb->ref_frame) {
274
            dc += lb->dc_coeff;
275
            count++;
276
        }
277
        if (ref_frame == ab->ref_frame) {
278
            dc += ab->dc_coeff;
279
            count++;
280
        }
281 7ecae905 Aurelien Jacobs
        if (s->avctx->codec->id == CODEC_ID_VP5)
282
            for (i=0; i<2; i++)
283
                if (count < 2 && ref_frame == ab[-1+2*i].ref_frame) {
284
                    dc += ab[-1+2*i].dc_coeff;
285
                    count++;
286
                }
287 5ce117c3 Aurelien Jacobs
        if (count == 0)
288 d7af6a9d Aurelien Jacobs
            dc = s->prev_dc[vp56_b2p[b]][ref_frame];
289 5ce117c3 Aurelien Jacobs
        else if (count == 2)
290
            dc /= 2;
291
292 d7af6a9d Aurelien Jacobs
        s->block_coeff[b][idx] += dc;
293
        s->prev_dc[vp56_b2p[b]][ref_frame] = s->block_coeff[b][idx];
294
        ab->dc_coeff = s->block_coeff[b][idx];
295 5ce117c3 Aurelien Jacobs
        ab->ref_frame = ref_frame;
296 d7af6a9d Aurelien Jacobs
        lb->dc_coeff = s->block_coeff[b][idx];
297 5ce117c3 Aurelien Jacobs
        lb->ref_frame = ref_frame;
298 d7af6a9d Aurelien Jacobs
        s->block_coeff[b][idx] *= s->dequant_dc;
299 5ce117c3 Aurelien Jacobs
    }
300
}
301
302 3d52bca6 Aurelien Jacobs
static void vp56_edge_filter(VP56Context *s, uint8_t *yuv,
303 5ce117c3 Aurelien Jacobs
                             int pix_inc, int line_inc, int t)
304
{
305
    int pix2_inc = 2 * pix_inc;
306
    int i, v;
307
308
    for (i=0; i<12; i++) {
309
        v = (yuv[-pix2_inc] + 3*(yuv[0]-yuv[-pix_inc]) - yuv[pix_inc] + 4) >>3;
310
        v = s->adjust(v, t);
311 f66e4f5f Reimar Döffinger
        yuv[-pix_inc] = av_clip_uint8(yuv[-pix_inc] + v);
312
        yuv[0] = av_clip_uint8(yuv[0] - v);
313 5ce117c3 Aurelien Jacobs
        yuv += line_inc;
314
    }
315
}
316
317 3d52bca6 Aurelien Jacobs
static void vp56_deblock_filter(VP56Context *s, uint8_t *yuv,
318 5ce117c3 Aurelien Jacobs
                                int stride, int dx, int dy)
319
{
320
    int t = vp56_filter_threshold[s->quantizer];
321
    if (dx)  vp56_edge_filter(s, yuv +         10-dx ,      1, stride, t);
322
    if (dy)  vp56_edge_filter(s, yuv + stride*(10-dy), stride,      1, t);
323
}
324
325 3d52bca6 Aurelien Jacobs
static void vp56_mc(VP56Context *s, int b, int plane, uint8_t *src,
326 5ce117c3 Aurelien Jacobs
                    int stride, int x, int y)
327
{
328 704a2881 Aurelien Jacobs
    uint8_t *dst=s->framep[VP56_FRAME_CURRENT]->data[plane]+s->block_offset[b];
329 5ce117c3 Aurelien Jacobs
    uint8_t *src_block;
330
    int src_offset;
331
    int overlap_offset = 0;
332
    int mask = s->vp56_coord_div[b] - 1;
333
    int deblock_filtering = s->deblock_filtering;
334
    int dx;
335
    int dy;
336
337
    if (s->avctx->skip_loop_filter >= AVDISCARD_ALL ||
338
        (s->avctx->skip_loop_filter >= AVDISCARD_NONKEY
339 704a2881 Aurelien Jacobs
         && !s->framep[VP56_FRAME_CURRENT]->key_frame))
340 5ce117c3 Aurelien Jacobs
        deblock_filtering = 0;
341
342
    dx = s->mv[b].x / s->vp56_coord_div[b];
343
    dy = s->mv[b].y / s->vp56_coord_div[b];
344
345
    if (b >= 4) {
346
        x /= 2;
347
        y /= 2;
348
    }
349
    x += dx - 2;
350
    y += dy - 2;
351
352
    if (x<0 || x+12>=s->plane_width[plane] ||
353
        y<0 || y+12>=s->plane_height[plane]) {
354
        ff_emulated_edge_mc(s->edge_emu_buffer,
355
                            src + s->block_offset[b] + (dy-2)*stride + (dx-2),
356
                            stride, 12, 12, x, y,
357
                            s->plane_width[plane],
358
                            s->plane_height[plane]);
359
        src_block = s->edge_emu_buffer;
360
        src_offset = 2 + 2*stride;
361
    } else if (deblock_filtering) {
362 a8678a3a Aurelien Jacobs
        /* only need a 12x12 block, but there is no such dsp function, */
363
        /* so copy a 16x12 block */
364
        s->dsp.put_pixels_tab[0][0](s->edge_emu_buffer,
365
                                    src + s->block_offset[b] + (dy-2)*stride + (dx-2),
366
                                    stride, 12);
367 5ce117c3 Aurelien Jacobs
        src_block = s->edge_emu_buffer;
368
        src_offset = 2 + 2*stride;
369
    } else {
370
        src_block = src;
371
        src_offset = s->block_offset[b] + dy*stride + dx;
372
    }
373
374
    if (deblock_filtering)
375
        vp56_deblock_filter(s, src_block, stride, dx&7, dy&7);
376
377
    if (s->mv[b].x & mask)
378
        overlap_offset += (s->mv[b].x > 0) ? 1 : -1;
379
    if (s->mv[b].y & mask)
380
        overlap_offset += (s->mv[b].y > 0) ? stride : -stride;
381
382
    if (overlap_offset) {
383
        if (s->filter)
384
            s->filter(s, dst, src_block, src_offset, src_offset+overlap_offset,
385
                      stride, s->mv[b], mask, s->filter_selection, b<4);
386
        else
387
            s->dsp.put_no_rnd_pixels_l2[1](dst, src_block+src_offset,
388
                                           src_block+src_offset+overlap_offset,
389
                                           stride, 8);
390
    } else {
391
        s->dsp.put_pixels_tab[1][0](dst, src_block+src_offset, stride, 8);
392
    }
393
}
394
395 3d52bca6 Aurelien Jacobs
static void vp56_decode_mb(VP56Context *s, int row, int col, int is_alpha)
396 5ce117c3 Aurelien Jacobs
{
397
    AVFrame *frame_current, *frame_ref;
398 3d52bca6 Aurelien Jacobs
    VP56mb mb_type;
399
    VP56Frame ref_frame;
400 442b145a Aurelien Jacobs
    int b, ab, b_max, plane, off;
401 5ce117c3 Aurelien Jacobs
402 704a2881 Aurelien Jacobs
    if (s->framep[VP56_FRAME_CURRENT]->key_frame)
403 5ce117c3 Aurelien Jacobs
        mb_type = VP56_MB_INTRA;
404
    else
405
        mb_type = vp56_decode_mv(s, row, col);
406
    ref_frame = vp56_reference_frame[mb_type];
407
408 5fecfb7d Loren Merritt
    s->dsp.clear_blocks(*s->block_coeff);
409 5ce117c3 Aurelien Jacobs
410
    s->parse_coeff(s);
411
412
    vp56_add_predictors_dc(s, ref_frame);
413
414 704a2881 Aurelien Jacobs
    frame_current = s->framep[VP56_FRAME_CURRENT];
415
    frame_ref = s->framep[ref_frame];
416 5ce117c3 Aurelien Jacobs
417 91fc2cf1 Aurelien Jacobs
    ab = 6*is_alpha;
418
    b_max = 6 - 2*is_alpha;
419
420 5ce117c3 Aurelien Jacobs
    switch (mb_type) {
421
        case VP56_MB_INTRA:
422 91fc2cf1 Aurelien Jacobs
            for (b=0; b<b_max; b++) {
423 442b145a Aurelien Jacobs
                plane = vp56_b2p[b+ab];
424
                s->dsp.idct_put(frame_current->data[plane] + s->block_offset[b],
425
                                s->stride[plane], s->block_coeff[b]);
426 5ce117c3 Aurelien Jacobs
            }
427
            break;
428
429
        case VP56_MB_INTER_NOVEC_PF:
430
        case VP56_MB_INTER_NOVEC_GF:
431 91fc2cf1 Aurelien Jacobs
            for (b=0; b<b_max; b++) {
432 442b145a Aurelien Jacobs
                plane = vp56_b2p[b+ab];
433 5ce117c3 Aurelien Jacobs
                off = s->block_offset[b];
434 442b145a Aurelien Jacobs
                s->dsp.put_pixels_tab[1][0](frame_current->data[plane] + off,
435
                                            frame_ref->data[plane] + off,
436
                                            s->stride[plane], 8);
437
                s->dsp.idct_add(frame_current->data[plane] + off,
438
                                s->stride[plane], s->block_coeff[b]);
439 5ce117c3 Aurelien Jacobs
            }
440
            break;
441
442
        case VP56_MB_INTER_DELTA_PF:
443
        case VP56_MB_INTER_V1_PF:
444
        case VP56_MB_INTER_V2_PF:
445
        case VP56_MB_INTER_DELTA_GF:
446
        case VP56_MB_INTER_4V:
447
        case VP56_MB_INTER_V1_GF:
448
        case VP56_MB_INTER_V2_GF:
449 91fc2cf1 Aurelien Jacobs
            for (b=0; b<b_max; b++) {
450 5ce117c3 Aurelien Jacobs
                int x_off = b==1 || b==3 ? 8 : 0;
451
                int y_off = b==2 || b==3 ? 8 : 0;
452 442b145a Aurelien Jacobs
                plane = vp56_b2p[b+ab];
453
                vp56_mc(s, b, plane, frame_ref->data[plane], s->stride[plane],
454 5ce117c3 Aurelien Jacobs
                        16*col+x_off, 16*row+y_off);
455 442b145a Aurelien Jacobs
                s->dsp.idct_add(frame_current->data[plane] + s->block_offset[b],
456
                                s->stride[plane], s->block_coeff[b]);
457 5ce117c3 Aurelien Jacobs
            }
458
            break;
459
    }
460
}
461
462 1457516f Aurelien Jacobs
static int vp56_size_changed(AVCodecContext *avctx)
463 5ce117c3 Aurelien Jacobs
{
464 3d52bca6 Aurelien Jacobs
    VP56Context *s = avctx->priv_data;
465 704a2881 Aurelien Jacobs
    int stride = s->framep[VP56_FRAME_CURRENT]->linesize[0];
466 5ce117c3 Aurelien Jacobs
    int i;
467
468 91fc2cf1 Aurelien Jacobs
    s->plane_width[0]  = s->plane_width[3]  = avctx->coded_width;
469 1457516f Aurelien Jacobs
    s->plane_width[1]  = s->plane_width[2]  = avctx->coded_width/2;
470 91fc2cf1 Aurelien Jacobs
    s->plane_height[0] = s->plane_height[3] = avctx->coded_height;
471 1457516f Aurelien Jacobs
    s->plane_height[1] = s->plane_height[2] = avctx->coded_height/2;
472 5ce117c3 Aurelien Jacobs
473 91fc2cf1 Aurelien Jacobs
    for (i=0; i<4; i++)
474 704a2881 Aurelien Jacobs
        s->stride[i] = s->flip * s->framep[VP56_FRAME_CURRENT]->linesize[i];
475 5ce117c3 Aurelien Jacobs
476 1457516f Aurelien Jacobs
    s->mb_width  = (avctx->coded_width +15) / 16;
477
    s->mb_height = (avctx->coded_height+15) / 16;
478 5ce117c3 Aurelien Jacobs
479
    if (s->mb_width > 1000 || s->mb_height > 1000) {
480
        av_log(avctx, AV_LOG_ERROR, "picture too big\n");
481
        return -1;
482
    }
483
484
    s->above_blocks = av_realloc(s->above_blocks,
485
                                 (4*s->mb_width+6) * sizeof(*s->above_blocks));
486
    s->macroblocks = av_realloc(s->macroblocks,
487
                                s->mb_width*s->mb_height*sizeof(*s->macroblocks));
488 a8678a3a Aurelien Jacobs
    av_free(s->edge_emu_buffer_alloc);
489
    s->edge_emu_buffer_alloc = av_malloc(16*stride);
490 5ce117c3 Aurelien Jacobs
    s->edge_emu_buffer = s->edge_emu_buffer_alloc;
491
    if (s->flip < 0)
492
        s->edge_emu_buffer += 15 * stride;
493
494
    return 0;
495
}
496
497
int vp56_decode_frame(AVCodecContext *avctx, void *data, int *data_size,
498 012f9308 Aurelien Jacobs
                      AVPacket *avpkt)
499 5ce117c3 Aurelien Jacobs
{
500 012f9308 Aurelien Jacobs
    const uint8_t *buf = avpkt->data;
501 3d52bca6 Aurelien Jacobs
    VP56Context *s = avctx->priv_data;
502 704a2881 Aurelien Jacobs
    AVFrame *const p = s->framep[VP56_FRAME_CURRENT];
503 012f9308 Aurelien Jacobs
    int remaining_buf_size = avpkt->size;
504 ed761067 Aurelien Jacobs
    int is_alpha, av_uninit(alpha_offset);
505 91fc2cf1 Aurelien Jacobs
506
    if (s->has_alpha) {
507 68a4d349 Laurent Aimar
        if (remaining_buf_size < 3)
508
            return -1;
509 91fc2cf1 Aurelien Jacobs
        alpha_offset = bytestream_get_be24(&buf);
510 5210529e Aurelien Jacobs
        remaining_buf_size -= 3;
511 68a4d349 Laurent Aimar
        if (remaining_buf_size < alpha_offset)
512
            return -1;
513 91fc2cf1 Aurelien Jacobs
    }
514
515
    for (is_alpha=0; is_alpha < 1+s->has_alpha; is_alpha++) {
516 f62a2b61 Aurelien Jacobs
        int mb_row, mb_col, mb_row_flip, mb_offset = 0;
517
        int block, y, uv, stride_y, stride_uv;
518
        int golden_frame = 0;
519
        int res;
520 5ce117c3 Aurelien Jacobs
521 f62a2b61 Aurelien Jacobs
        s->modelp = &s->models[is_alpha];
522 247df384 Aurelien Jacobs
523 5210529e Aurelien Jacobs
        res = s->parse_header(s, buf, remaining_buf_size, &golden_frame);
524 f62a2b61 Aurelien Jacobs
        if (!res)
525
            return -1;
526 5ce117c3 Aurelien Jacobs
527 f62a2b61 Aurelien Jacobs
        if (!is_alpha) {
528
            p->reference = 1;
529
            if (avctx->get_buffer(avctx, p) < 0) {
530
                av_log(avctx, AV_LOG_ERROR, "get_buffer() failed\n");
531
                return -1;
532
            }
533 5ce117c3 Aurelien Jacobs
534 f62a2b61 Aurelien Jacobs
            if (res == 2)
535
                if (vp56_size_changed(avctx)) {
536
                    avctx->release_buffer(avctx, p);
537
                    return -1;
538
                }
539 5ce117c3 Aurelien Jacobs
        }
540
541 f62a2b61 Aurelien Jacobs
        if (p->key_frame) {
542
            p->pict_type = FF_I_TYPE;
543
            s->default_models_init(s);
544
            for (block=0; block<s->mb_height*s->mb_width; block++)
545
                s->macroblocks[block].type = VP56_MB_INTRA;
546
        } else {
547
            p->pict_type = FF_P_TYPE;
548
            vp56_parse_mb_type_models(s);
549
            s->parse_vector_models(s);
550
            s->mb_type = VP56_MB_INTER_NOVEC_PF;
551
        }
552 5ce117c3 Aurelien Jacobs
553 f62a2b61 Aurelien Jacobs
        s->parse_coeff_models(s);
554 5ce117c3 Aurelien Jacobs
555 f62a2b61 Aurelien Jacobs
        memset(s->prev_dc, 0, sizeof(s->prev_dc));
556
        s->prev_dc[1][VP56_FRAME_CURRENT] = 128;
557
        s->prev_dc[2][VP56_FRAME_CURRENT] = 128;
558 5ce117c3 Aurelien Jacobs
559 f62a2b61 Aurelien Jacobs
        for (block=0; block < 4*s->mb_width+6; block++) {
560 738a89b9 Carl Eugen Hoyos
            s->above_blocks[block].ref_frame = VP56_FRAME_NONE;
561 f62a2b61 Aurelien Jacobs
            s->above_blocks[block].dc_coeff = 0;
562
            s->above_blocks[block].not_null_dc = 0;
563
        }
564 738a89b9 Carl Eugen Hoyos
        s->above_blocks[2*s->mb_width + 2].ref_frame = VP56_FRAME_CURRENT;
565
        s->above_blocks[3*s->mb_width + 4].ref_frame = VP56_FRAME_CURRENT;
566 5ce117c3 Aurelien Jacobs
567 f62a2b61 Aurelien Jacobs
        stride_y  = p->linesize[0];
568
        stride_uv = p->linesize[1];
569 5ce117c3 Aurelien Jacobs
570
        if (s->flip < 0)
571 f62a2b61 Aurelien Jacobs
            mb_offset = 7;
572
573
        /* main macroblocks loop */
574
        for (mb_row=0; mb_row<s->mb_height; mb_row++) {
575
            if (s->flip < 0)
576
                mb_row_flip = s->mb_height - mb_row - 1;
577
            else
578
                mb_row_flip = mb_row;
579
580
            for (block=0; block<4; block++) {
581 738a89b9 Carl Eugen Hoyos
                s->left_block[block].ref_frame = VP56_FRAME_NONE;
582 f62a2b61 Aurelien Jacobs
                s->left_block[block].dc_coeff = 0;
583
                s->left_block[block].not_null_dc = 0;
584 5ce117c3 Aurelien Jacobs
            }
585 d3f9edba Aurelien Jacobs
            memset(s->coeff_ctx, 0, sizeof(s->coeff_ctx));
586 f62a2b61 Aurelien Jacobs
            memset(s->coeff_ctx_last, 24, sizeof(s->coeff_ctx_last));
587
588
            s->above_block_idx[0] = 1;
589
            s->above_block_idx[1] = 2;
590
            s->above_block_idx[2] = 1;
591
            s->above_block_idx[3] = 2;
592
            s->above_block_idx[4] = 2*s->mb_width + 2 + 1;
593
            s->above_block_idx[5] = 3*s->mb_width + 4 + 1;
594
595
            s->block_offset[s->frbi] = (mb_row_flip*16 + mb_offset) * stride_y;
596
            s->block_offset[s->srbi] = s->block_offset[s->frbi] + 8*stride_y;
597
            s->block_offset[1] = s->block_offset[0] + 8;
598
            s->block_offset[3] = s->block_offset[2] + 8;
599
            s->block_offset[4] = (mb_row_flip*8 + mb_offset) * stride_uv;
600
            s->block_offset[5] = s->block_offset[4];
601
602
            for (mb_col=0; mb_col<s->mb_width; mb_col++) {
603
                vp56_decode_mb(s, mb_row, mb_col, is_alpha);
604
605
                for (y=0; y<4; y++) {
606
                    s->above_block_idx[y] += 2;
607
                    s->block_offset[y] += 16;
608
                }
609 5ce117c3 Aurelien Jacobs
610 f62a2b61 Aurelien Jacobs
                for (uv=4; uv<6; uv++) {
611
                    s->above_block_idx[uv] += 1;
612
                    s->block_offset[uv] += 8;
613
                }
614 5ce117c3 Aurelien Jacobs
            }
615
        }
616
617 f62a2b61 Aurelien Jacobs
        if (p->key_frame || golden_frame) {
618
            if (s->framep[VP56_FRAME_GOLDEN]->data[0] &&
619
                s->framep[VP56_FRAME_GOLDEN] != s->framep[VP56_FRAME_GOLDEN2])
620
                avctx->release_buffer(avctx, s->framep[VP56_FRAME_GOLDEN]);
621
            s->framep[VP56_FRAME_GOLDEN] = p;
622
        }
623 91fc2cf1 Aurelien Jacobs
624 f62a2b61 Aurelien Jacobs
        if (s->has_alpha) {
625
            FFSWAP(AVFrame *, s->framep[VP56_FRAME_GOLDEN],
626
                              s->framep[VP56_FRAME_GOLDEN2]);
627
            buf += alpha_offset;
628 5210529e Aurelien Jacobs
            remaining_buf_size -= alpha_offset;
629 f62a2b61 Aurelien Jacobs
        }
630 91fc2cf1 Aurelien Jacobs
    }
631
632
    if (s->framep[VP56_FRAME_PREVIOUS] == s->framep[VP56_FRAME_GOLDEN] ||
633
        s->framep[VP56_FRAME_PREVIOUS] == s->framep[VP56_FRAME_GOLDEN2]) {
634
        if (s->framep[VP56_FRAME_UNUSED] != s->framep[VP56_FRAME_GOLDEN] &&
635
            s->framep[VP56_FRAME_UNUSED] != s->framep[VP56_FRAME_GOLDEN2])
636
            FFSWAP(AVFrame *, s->framep[VP56_FRAME_PREVIOUS],
637
                              s->framep[VP56_FRAME_UNUSED]);
638
        else
639
            FFSWAP(AVFrame *, s->framep[VP56_FRAME_PREVIOUS],
640
                              s->framep[VP56_FRAME_UNUSED2]);
641
    } else if (s->framep[VP56_FRAME_PREVIOUS]->data[0])
642
        avctx->release_buffer(avctx, s->framep[VP56_FRAME_PREVIOUS]);
643 704a2881 Aurelien Jacobs
    FFSWAP(AVFrame *, s->framep[VP56_FRAME_CURRENT],
644
                      s->framep[VP56_FRAME_PREVIOUS]);
645 5ce117c3 Aurelien Jacobs
646 eba0fcad Aurelien Jacobs
    *(AVFrame*)data = *p;
647
    *data_size = sizeof(AVFrame);
648 5ce117c3 Aurelien Jacobs
649 012f9308 Aurelien Jacobs
    return avpkt->size;
650 5ce117c3 Aurelien Jacobs
}
651
652 98a6fff9 Zuxy Meng
av_cold void vp56_init(AVCodecContext *avctx, int flip, int has_alpha)
653 5ce117c3 Aurelien Jacobs
{
654 3d52bca6 Aurelien Jacobs
    VP56Context *s = avctx->priv_data;
655 5ce117c3 Aurelien Jacobs
    int i;
656
657
    s->avctx = avctx;
658 91fc2cf1 Aurelien Jacobs
    avctx->pix_fmt = has_alpha ? PIX_FMT_YUVA420P : PIX_FMT_YUV420P;
659 5ce117c3 Aurelien Jacobs
660 1457516f Aurelien Jacobs
    if (avctx->idct_algo == FF_IDCT_AUTO)
661
        avctx->idct_algo = FF_IDCT_VP3;
662
    dsputil_init(&s->dsp, avctx);
663 5ce117c3 Aurelien Jacobs
    ff_init_scantable(s->dsp.idct_permutation, &s->scantable,ff_zigzag_direct);
664
665 91fc2cf1 Aurelien Jacobs
    for (i=0; i<4; i++)
666 704a2881 Aurelien Jacobs
        s->framep[i] = &s->frames[i];
667
    s->framep[VP56_FRAME_UNUSED] = s->framep[VP56_FRAME_GOLDEN];
668 91fc2cf1 Aurelien Jacobs
    s->framep[VP56_FRAME_UNUSED2] = s->framep[VP56_FRAME_GOLDEN2];
669 5ce117c3 Aurelien Jacobs
    s->edge_emu_buffer_alloc = NULL;
670
671
    s->above_blocks = NULL;
672
    s->macroblocks = NULL;
673
    s->quantizer = -1;
674
    s->deblock_filtering = 1;
675
676
    s->filter = NULL;
677
678 91fc2cf1 Aurelien Jacobs
    s->has_alpha = has_alpha;
679 5ce117c3 Aurelien Jacobs
    if (flip) {
680
        s->flip = -1;
681
        s->frbi = 2;
682
        s->srbi = 0;
683
    } else {
684
        s->flip = 1;
685
        s->frbi = 0;
686
        s->srbi = 2;
687
    }
688
}
689
690 98a6fff9 Zuxy Meng
av_cold int vp56_free(AVCodecContext *avctx)
691 5ce117c3 Aurelien Jacobs
{
692 3d52bca6 Aurelien Jacobs
    VP56Context *s = avctx->priv_data;
693 5ce117c3 Aurelien Jacobs
694 6242b1c4 Reimar Döffinger
    av_freep(&s->above_blocks);
695
    av_freep(&s->macroblocks);
696
    av_freep(&s->edge_emu_buffer_alloc);
697 1457516f Aurelien Jacobs
    if (s->framep[VP56_FRAME_GOLDEN]->data[0])
698 704a2881 Aurelien Jacobs
        avctx->release_buffer(avctx, s->framep[VP56_FRAME_GOLDEN]);
699 91fc2cf1 Aurelien Jacobs
    if (s->framep[VP56_FRAME_GOLDEN2]->data[0])
700
        avctx->release_buffer(avctx, s->framep[VP56_FRAME_GOLDEN2]);
701 704a2881 Aurelien Jacobs
    if (s->framep[VP56_FRAME_PREVIOUS]->data[0])
702
        avctx->release_buffer(avctx, s->framep[VP56_FRAME_PREVIOUS]);
703 5ce117c3 Aurelien Jacobs
    return 0;
704
}