Merge remote-tracking branch 'qatar/master'
[ffmpeg.git] / libavcodec / proresdsp.c
1 /*
2  * Apple ProRes compatible decoder
3  *
4  * Copyright (c) 2010-2011 Maxim Poliakovski
5  *
6  * This file is part of Libav.
7  *
8  * Libav is free software; you can redistribute it and/or
9  * modify it under the terms of the GNU Lesser General Public
10  * License as published by the Free Software Foundation; either
11  * version 2.1 of the License, or (at your option) any later version.
12  *
13  * Libav is distributed in the hope that it will be useful,
14  * but WITHOUT ANY WARRANTY; without even the implied warranty of
15  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
16  * Lesser General Public License for more details.
17  *
18  * You should have received a copy of the GNU Lesser General Public
19  * License along with Libav; if not, write to the Free Software
20  * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
21  */
22
23 #include "proresdsp.h"
24 #include "simple_idct.h"
25
26 #define BIAS     (1 << (PRORES_BITS_PER_SAMPLE - 1))           ///< bias value for converting signed pixels into unsigned ones
27 #define CLIP_MIN (1 << (PRORES_BITS_PER_SAMPLE - 8))           ///< minimum value for clipping resulting pixels
28 #define CLIP_MAX (1 << PRORES_BITS_PER_SAMPLE) - CLIP_MIN - 1  ///< maximum value for clipping resulting pixels
29
30 #define CLIP_AND_BIAS(x) (av_clip((x) + BIAS, CLIP_MIN, CLIP_MAX))
31
32 /**
33  * Add bias value, clamp and output pixels of a slice
34  */
35 static void put_pixels(uint16_t *dst, int stride, const DCTELEM *in)
36 {
37     int x, y, src_offset, dst_offset;
38
39     for (y = 0, dst_offset = 0; y < 8; y++, dst_offset += stride) {
40         for (x = 0; x < 8; x++) {
41             src_offset = (y << 3) + x;
42
43             dst[dst_offset + x] = CLIP_AND_BIAS(in[src_offset]);
44         }
45     }
46 }
47
48 static void prores_idct_put_c(uint16_t *out, int linesize, DCTELEM *block, const int16_t *qmat)
49 {
50     ff_prores_idct(block, qmat);
51     put_pixels(out, linesize >> 1, block);
52 }
53
54 static void prores_fdct_c(const uint16_t *src, int linesize, DCTELEM *block)
55 {
56     int x, y;
57     const uint16_t *tsrc = src;
58
59     for (y = 0; y < 8; y++) {
60         for (x = 0; x < 8; x++)
61             block[y * 8 + x] = tsrc[x];
62         tsrc += linesize >> 1;
63     }
64     ff_jpeg_fdct_islow_10(block);
65 }
66
67 void ff_proresdsp_init(ProresDSPContext *dsp, AVCodecContext *avctx)
68 {
69     dsp->idct_put = prores_idct_put_c;
70     dsp->idct_permutation_type = FF_NO_IDCT_PERM;
71     dsp->fdct     = prores_fdct_c;
72     dsp->dct_permutation_type  = FF_NO_IDCT_PERM;
73
74     if (HAVE_MMX) ff_proresdsp_x86_init(dsp, avctx);
75
76     ff_init_scantable_permutation(dsp->idct_permutation,
77                                   dsp->idct_permutation_type);
78     ff_init_scantable_permutation(dsp->dct_permutation,
79                                   dsp->dct_permutation_type);
80 }