repositories / deb_ffmpeg.git / blame
? search:
summary | shortlog | log | commit | commitdiff | tree
blob | blame (incremental) | history | HEAD
Imported Debian version 2.4.3~trusty1
[deb_ffmpeg.git] / ffmpeg / libavcodec / qpeldsp.c
CommitLineData
2ba45a60
DM		 1/*
2 * quarterpel DSP functions
3 * Copyright (c) 2000, 2001 Fabrice Bellard
4 * Copyright (c) 2002-2004 Michael Niedermayer <michaelni@gmx.at>
5 *
6 * This file is part of FFmpeg.
7 *
8 * FFmpeg 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 * FFmpeg 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 FFmpeg; if not, write to the Free Software
20 * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
21 */
22
23/**
24 * @file
25 * quarterpel DSP functions
26 */
27
28#include <stddef.h>
29#include <stdint.h>
30
31#include "config.h"
32#include "libavutil/attributes.h"
33#include "copy_block.h"
34#include "qpeldsp.h"
35#include "diracdsp.h"
36
37#define BIT_DEPTH 8
38#include "hpel_template.c"
39#include "pel_template.c"
40#include "qpel_template.c"
41
42#define QPEL_MC(r, OPNAME, RND, OP) \
43static void OPNAME ## mpeg4_qpel8_h_lowpass(uint8_t *dst, const uint8_t *src, \
44 int dstStride, int srcStride, \
45 int h) \
46{ \
47 const uint8_t *cm = ff_crop_tab + MAX_NEG_CROP; \
48 int i; \
49 \
50 for (i = 0; i < h; i++) { \
51 OP(dst[0], (src[0] + src[1]) * 20 - (src[0] + src[2]) * 6 + (src[1] + src[3]) * 3 - (src[2] + src[4])); \
52 OP(dst[1], (src[1] + src[2]) * 20 - (src[0] + src[3]) * 6 + (src[0] + src[4]) * 3 - (src[1] + src[5])); \
53 OP(dst[2], (src[2] + src[3]) * 20 - (src[1] + src[4]) * 6 + (src[0] + src[5]) * 3 - (src[0] + src[6])); \
54 OP(dst[3], (src[3] + src[4]) * 20 - (src[2] + src[5]) * 6 + (src[1] + src[6]) * 3 - (src[0] + src[7])); \
55 OP(dst[4], (src[4] + src[5]) * 20 - (src[3] + src[6]) * 6 + (src[2] + src[7]) * 3 - (src[1] + src[8])); \
56 OP(dst[5], (src[5] + src[6]) * 20 - (src[4] + src[7]) * 6 + (src[3] + src[8]) * 3 - (src[2] + src[8])); \
57 OP(dst[6], (src[6] + src[7]) * 20 - (src[5] + src[8]) * 6 + (src[4] + src[8]) * 3 - (src[3] + src[7])); \
58 OP(dst[7], (src[7] + src[8]) * 20 - (src[6] + src[8]) * 6 + (src[5] + src[7]) * 3 - (src[4] + src[6])); \
59 dst += dstStride; \
60 src += srcStride; \
61 } \
62} \
63 \
64static void OPNAME ## mpeg4_qpel8_v_lowpass(uint8_t *dst, const uint8_t *src, \
65 int dstStride, int srcStride) \
66{ \
67 const uint8_t *cm = ff_crop_tab + MAX_NEG_CROP; \
68 const int w = 8; \
69 int i; \
70 \
71 for (i = 0; i < w; i++) { \
72 const int src0 = src[0 * srcStride]; \
73 const int src1 = src[1 * srcStride]; \
74 const int src2 = src[2 * srcStride]; \
75 const int src3 = src[3 * srcStride]; \
76 const int src4 = src[4 * srcStride]; \
77 const int src5 = src[5 * srcStride]; \
78 const int src6 = src[6 * srcStride]; \
79 const int src7 = src[7 * srcStride]; \
80 const int src8 = src[8 * srcStride]; \
81 OP(dst[0 * dstStride], (src0 + src1) * 20 - (src0 + src2) * 6 + (src1 + src3) * 3 - (src2 + src4)); \
82 OP(dst[1 * dstStride], (src1 + src2) * 20 - (src0 + src3) * 6 + (src0 + src4) * 3 - (src1 + src5)); \
83 OP(dst[2 * dstStride], (src2 + src3) * 20 - (src1 + src4) * 6 + (src0 + src5) * 3 - (src0 + src6)); \
84 OP(dst[3 * dstStride], (src3 + src4) * 20 - (src2 + src5) * 6 + (src1 + src6) * 3 - (src0 + src7)); \
85 OP(dst[4 * dstStride], (src4 + src5) * 20 - (src3 + src6) * 6 + (src2 + src7) * 3 - (src1 + src8)); \
86 OP(dst[5 * dstStride], (src5 + src6) * 20 - (src4 + src7) * 6 + (src3 + src8) * 3 - (src2 + src8)); \
87 OP(dst[6 * dstStride], (src6 + src7) * 20 - (src5 + src8) * 6 + (src4 + src8) * 3 - (src3 + src7)); \
88 OP(dst[7 * dstStride], (src7 + src8) * 20 - (src6 + src8) * 6 + (src5 + src7) * 3 - (src4 + src6)); \
89 dst++; \
90 src++; \
91 } \
92} \
93 \
94static void OPNAME ## mpeg4_qpel16_h_lowpass(uint8_t *dst, \
95 const uint8_t *src, \
96 int dstStride, int srcStride, \
97 int h) \
98{ \
99 const uint8_t *cm = ff_crop_tab + MAX_NEG_CROP; \
100 int i; \
101 \
102 for (i = 0; i < h; i++) { \
103 OP(dst[0], (src[0] + src[1]) * 20 - (src[0] + src[2]) * 6 + (src[1] + src[3]) * 3 - (src[2] + src[4])); \
104 OP(dst[1], (src[1] + src[2]) * 20 - (src[0] + src[3]) * 6 + (src[0] + src[4]) * 3 - (src[1] + src[5])); \
105 OP(dst[2], (src[2] + src[3]) * 20 - (src[1] + src[4]) * 6 + (src[0] + src[5]) * 3 - (src[0] + src[6])); \
106 OP(dst[3], (src[3] + src[4]) * 20 - (src[2] + src[5]) * 6 + (src[1] + src[6]) * 3 - (src[0] + src[7])); \
107 OP(dst[4], (src[4] + src[5]) * 20 - (src[3] + src[6]) * 6 + (src[2] + src[7]) * 3 - (src[1] + src[8])); \
108 OP(dst[5], (src[5] + src[6]) * 20 - (src[4] + src[7]) * 6 + (src[3] + src[8]) * 3 - (src[2] + src[9])); \
109 OP(dst[6], (src[6] + src[7]) * 20 - (src[5] + src[8]) * 6 + (src[4] + src[9]) * 3 - (src[3] + src[10])); \
110 OP(dst[7], (src[7] + src[8]) * 20 - (src[6] + src[9]) * 6 + (src[5] + src[10]) * 3 - (src[4] + src[11])); \
111 OP(dst[8], (src[8] + src[9]) * 20 - (src[7] + src[10]) * 6 + (src[6] + src[11]) * 3 - (src[5] + src[12])); \
112 OP(dst[9], (src[9] + src[10]) * 20 - (src[8] + src[11]) * 6 + (src[7] + src[12]) * 3 - (src[6] + src[13])); \
113 OP(dst[10], (src[10] + src[11]) * 20 - (src[9] + src[12]) * 6 + (src[8] + src[13]) * 3 - (src[7] + src[14])); \
114 OP(dst[11], (src[11] + src[12]) * 20 - (src[10] + src[13]) * 6 + (src[9] + src[14]) * 3 - (src[8] + src[15])); \
115 OP(dst[12], (src[12] + src[13]) * 20 - (src[11] + src[14]) * 6 + (src[10] + src[15]) * 3 - (src[9] + src[16])); \
116 OP(dst[13], (src[13] + src[14]) * 20 - (src[12] + src[15]) * 6 + (src[11] + src[16]) * 3 - (src[10] + src[16])); \
117 OP(dst[14], (src[14] + src[15]) * 20 - (src[13] + src[16]) * 6 + (src[12] + src[16]) * 3 - (src[11] + src[15])); \
118 OP(dst[15], (src[15] + src[16]) * 20 - (src[14] + src[16]) * 6 + (src[13] + src[15]) * 3 - (src[12] + src[14])); \
119 dst += dstStride; \
120 src += srcStride; \
121 } \
122} \
123 \
124static void OPNAME ## mpeg4_qpel16_v_lowpass(uint8_t *dst, \
125 const uint8_t *src, \
126 int dstStride, int srcStride) \
127{ \
128 const uint8_t *cm = ff_crop_tab + MAX_NEG_CROP; \
129 const int w = 16; \
130 int i; \
131 \
132 for (i = 0; i < w; i++) { \
133 const int src0 = src[0 * srcStride]; \
134 const int src1 = src[1 * srcStride]; \
135 const int src2 = src[2 * srcStride]; \
136 const int src3 = src[3 * srcStride]; \
137 const int src4 = src[4 * srcStride]; \
138 const int src5 = src[5 * srcStride]; \
139 const int src6 = src[6 * srcStride]; \
140 const int src7 = src[7 * srcStride]; \
141 const int src8 = src[8 * srcStride]; \
142 const int src9 = src[9 * srcStride]; \
143 const int src10 = src[10 * srcStride]; \
144 const int src11 = src[11 * srcStride]; \
145 const int src12 = src[12 * srcStride]; \
146 const int src13 = src[13 * srcStride]; \
147 const int src14 = src[14 * srcStride]; \
148 const int src15 = src[15 * srcStride]; \
149 const int src16 = src[16 * srcStride]; \
150 OP(dst[0 * dstStride], (src0 + src1) * 20 - (src0 + src2) * 6 + (src1 + src3) * 3 - (src2 + src4)); \
151 OP(dst[1 * dstStride], (src1 + src2) * 20 - (src0 + src3) * 6 + (src0 + src4) * 3 - (src1 + src5)); \
152 OP(dst[2 * dstStride], (src2 + src3) * 20 - (src1 + src4) * 6 + (src0 + src5) * 3 - (src0 + src6)); \
153 OP(dst[3 * dstStride], (src3 + src4) * 20 - (src2 + src5) * 6 + (src1 + src6) * 3 - (src0 + src7)); \
154 OP(dst[4 * dstStride], (src4 + src5) * 20 - (src3 + src6) * 6 + (src2 + src7) * 3 - (src1 + src8)); \
155 OP(dst[5 * dstStride], (src5 + src6) * 20 - (src4 + src7) * 6 + (src3 + src8) * 3 - (src2 + src9)); \
156 OP(dst[6 * dstStride], (src6 + src7) * 20 - (src5 + src8) * 6 + (src4 + src9) * 3 - (src3 + src10)); \
157 OP(dst[7 * dstStride], (src7 + src8) * 20 - (src6 + src9) * 6 + (src5 + src10) * 3 - (src4 + src11)); \
158 OP(dst[8 * dstStride], (src8 + src9) * 20 - (src7 + src10) * 6 + (src6 + src11) * 3 - (src5 + src12)); \
159 OP(dst[9 * dstStride], (src9 + src10) * 20 - (src8 + src11) * 6 + (src7 + src12) * 3 - (src6 + src13)); \
160 OP(dst[10 * dstStride], (src10 + src11) * 20 - (src9 + src12) * 6 + (src8 + src13) * 3 - (src7 + src14)); \
161 OP(dst[11 * dstStride], (src11 + src12) * 20 - (src10 + src13) * 6 + (src9 + src14) * 3 - (src8 + src15)); \
162 OP(dst[12 * dstStride], (src12 + src13) * 20 - (src11 + src14) * 6 + (src10 + src15) * 3 - (src9 + src16)); \
163 OP(dst[13 * dstStride], (src13 + src14) * 20 - (src12 + src15) * 6 + (src11 + src16) * 3 - (src10 + src16)); \
164 OP(dst[14 * dstStride], (src14 + src15) * 20 - (src13 + src16) * 6 + (src12 + src16) * 3 - (src11 + src15)); \
165 OP(dst[15 * dstStride], (src15 + src16) * 20 - (src14 + src16) * 6 + (src13 + src15) * 3 - (src12 + src14)); \
166 dst++; \
167 src++; \
168 } \
169} \
170 \
171static void OPNAME ## qpel8_mc10_c(uint8_t *dst, const uint8_t *src, \
172 ptrdiff_t stride) \
173{ \
174 uint8_t half[64]; \
175 \
176 put ## RND ## mpeg4_qpel8_h_lowpass(half, src, 8, stride, 8); \
177 OPNAME ## pixels8_l2_8(dst, src, half, stride, stride, 8, 8); \
178} \
179 \
180static void OPNAME ## qpel8_mc20_c(uint8_t *dst, const uint8_t *src, \
181 ptrdiff_t stride) \
182{ \
183 OPNAME ## mpeg4_qpel8_h_lowpass(dst, src, stride, stride, 8); \
184} \
185 \
186static void OPNAME ## qpel8_mc30_c(uint8_t *dst, const uint8_t *src, \
187 ptrdiff_t stride) \
188{ \
189 uint8_t half[64]; \
190 \
191 put ## RND ## mpeg4_qpel8_h_lowpass(half, src, 8, stride, 8); \
192 OPNAME ## pixels8_l2_8(dst, src + 1, half, stride, stride, 8, 8); \
193} \
194 \
195static void OPNAME ## qpel8_mc01_c(uint8_t *dst, const uint8_t *src, \
196 ptrdiff_t stride) \
197{ \
198 uint8_t full[16 * 9]; \
199 uint8_t half[64]; \
200 \
201 copy_block9(full, src, 16, stride, 9); \
202 put ## RND ## mpeg4_qpel8_v_lowpass(half, full, 8, 16); \
203 OPNAME ## pixels8_l2_8(dst, full, half, stride, 16, 8, 8); \
204} \
205 \
206static void OPNAME ## qpel8_mc02_c(uint8_t *dst, const uint8_t *src, \
207 ptrdiff_t stride) \
208{ \
209 uint8_t full[16 * 9]; \
210 \
211 copy_block9(full, src, 16, stride, 9); \
212 OPNAME ## mpeg4_qpel8_v_lowpass(dst, full, stride, 16); \
213} \
214 \
215static void OPNAME ## qpel8_mc03_c(uint8_t *dst, const uint8_t *src, \
216 ptrdiff_t stride) \
217{ \
218 uint8_t full[16 * 9]; \
219 uint8_t half[64]; \
220 \
221 copy_block9(full, src, 16, stride, 9); \
222 put ## RND ## mpeg4_qpel8_v_lowpass(half, full, 8, 16); \
223 OPNAME ## pixels8_l2_8(dst, full + 16, half, stride, 16, 8, 8); \
224} \
225 \
226void ff_ ## OPNAME ## qpel8_mc11_old_c(uint8_t *dst, const uint8_t *src, \
227 ptrdiff_t stride) \
228{ \
229 uint8_t full[16 * 9]; \
230 uint8_t halfH[72]; \
231 uint8_t halfV[64]; \
232 uint8_t halfHV[64]; \
233 \
234 copy_block9(full, src, 16, stride, 9); \
235 put ## RND ## mpeg4_qpel8_h_lowpass(halfH, full, 8, 16, 9); \
236 put ## RND ## mpeg4_qpel8_v_lowpass(halfV, full, 8, 16); \
237 put ## RND ## mpeg4_qpel8_v_lowpass(halfHV, halfH, 8, 8); \
238 OPNAME ## pixels8_l4_8(dst, full, halfH, halfV, halfHV, \
239 stride, 16, 8, 8, 8, 8); \
240} \
241 \
242static void OPNAME ## qpel8_mc11_c(uint8_t *dst, const uint8_t *src, \
243 ptrdiff_t stride) \
244{ \
245 uint8_t full[16 * 9]; \
246 uint8_t halfH[72]; \
247 uint8_t halfHV[64]; \
248 \
249 copy_block9(full, src, 16, stride, 9); \
250 put ## RND ## mpeg4_qpel8_h_lowpass(halfH, full, 8, 16, 9); \
251 put ## RND ## pixels8_l2_8(halfH, halfH, full, 8, 8, 16, 9); \
252 put ## RND ## mpeg4_qpel8_v_lowpass(halfHV, halfH, 8, 8); \
253 OPNAME ## pixels8_l2_8(dst, halfH, halfHV, stride, 8, 8, 8); \
254} \
255 \
256void ff_ ## OPNAME ## qpel8_mc31_old_c(uint8_t *dst, const uint8_t *src, \
257 ptrdiff_t stride) \
258{ \
259 uint8_t full[16 * 9]; \
260 uint8_t halfH[72]; \
261 uint8_t halfV[64]; \
262 uint8_t halfHV[64]; \
263 \
264 copy_block9(full, src, 16, stride, 9); \
265 put ## RND ## mpeg4_qpel8_h_lowpass(halfH, full, 8, 16, 9); \
266 put ## RND ## mpeg4_qpel8_v_lowpass(halfV, full + 1, 8, 16); \
267 put ## RND ## mpeg4_qpel8_v_lowpass(halfHV, halfH, 8, 8); \
268 OPNAME ## pixels8_l4_8(dst, full + 1, halfH, halfV, halfHV, \
269 stride, 16, 8, 8, 8, 8); \
270} \
271 \
272static void OPNAME ## qpel8_mc31_c(uint8_t *dst, const uint8_t *src, \
273 ptrdiff_t stride) \
274{ \
275 uint8_t full[16 * 9]; \
276 uint8_t halfH[72]; \
277 uint8_t halfHV[64]; \
278 \
279 copy_block9(full, src, 16, stride, 9); \
280 put ## RND ## mpeg4_qpel8_h_lowpass(halfH, full, 8, 16, 9); \
281 put ## RND ## pixels8_l2_8(halfH, halfH, full + 1, 8, 8, 16, 9); \
282 put ## RND ## mpeg4_qpel8_v_lowpass(halfHV, halfH, 8, 8); \
283 OPNAME ## pixels8_l2_8(dst, halfH, halfHV, stride, 8, 8, 8); \
284} \
285 \
286void ff_ ## OPNAME ## qpel8_mc13_old_c(uint8_t *dst, const uint8_t *src, \
287 ptrdiff_t stride) \
288{ \
289 uint8_t full[16 * 9]; \
290 uint8_t halfH[72]; \
291 uint8_t halfV[64]; \
292 uint8_t halfHV[64]; \
293 \
294 copy_block9(full, src, 16, stride, 9); \
295 put ## RND ## mpeg4_qpel8_h_lowpass(halfH, full, 8, 16, 9); \
296 put ## RND ## mpeg4_qpel8_v_lowpass(halfV, full, 8, 16); \
297 put ## RND ## mpeg4_qpel8_v_lowpass(halfHV, halfH, 8, 8); \
298 OPNAME ## pixels8_l4_8(dst, full + 16, halfH + 8, halfV, halfHV, \
299 stride, 16, 8, 8, 8, 8); \
300} \
301 \
302static void OPNAME ## qpel8_mc13_c(uint8_t *dst, const uint8_t *src, \
303 ptrdiff_t stride) \
304{ \
305 uint8_t full[16 * 9]; \
306 uint8_t halfH[72]; \
307 uint8_t halfHV[64]; \
308 \
309 copy_block9(full, src, 16, stride, 9); \
310 put ## RND ## mpeg4_qpel8_h_lowpass(halfH, full, 8, 16, 9); \
311 put ## RND ## pixels8_l2_8(halfH, halfH, full, 8, 8, 16, 9); \
312 put ## RND ## mpeg4_qpel8_v_lowpass(halfHV, halfH, 8, 8); \
313 OPNAME ## pixels8_l2_8(dst, halfH + 8, halfHV, stride, 8, 8, 8); \
314} \
315 \
316void ff_ ## OPNAME ## qpel8_mc33_old_c(uint8_t *dst, const uint8_t *src, \
317 ptrdiff_t stride) \
318{ \
319 uint8_t full[16 * 9]; \
320 uint8_t halfH[72]; \
321 uint8_t halfV[64]; \
322 uint8_t halfHV[64]; \
323 \
324 copy_block9(full, src, 16, stride, 9); \
325 put ## RND ## mpeg4_qpel8_h_lowpass(halfH, full, 8, 16, 9); \
326 put ## RND ## mpeg4_qpel8_v_lowpass(halfV, full + 1, 8, 16); \
327 put ## RND ## mpeg4_qpel8_v_lowpass(halfHV, halfH, 8, 8); \
328 OPNAME ## pixels8_l4_8(dst, full + 17, halfH + 8, halfV, halfHV, \
329 stride, 16, 8, 8, 8, 8); \
330} \
331 \
332static void OPNAME ## qpel8_mc33_c(uint8_t *dst, const uint8_t *src, \
333 ptrdiff_t stride) \
334{ \
335 uint8_t full[16 * 9]; \
336 uint8_t halfH[72]; \
337 uint8_t halfHV[64]; \
338 \
339 copy_block9(full, src, 16, stride, 9); \
340 put ## RND ## mpeg4_qpel8_h_lowpass(halfH, full, 8, 16, 9); \
341 put ## RND ## pixels8_l2_8(halfH, halfH, full + 1, 8, 8, 16, 9); \
342 put ## RND ## mpeg4_qpel8_v_lowpass(halfHV, halfH, 8, 8); \
343 OPNAME ## pixels8_l2_8(dst, halfH + 8, halfHV, stride, 8, 8, 8); \
344} \
345 \
346static void OPNAME ## qpel8_mc21_c(uint8_t *dst, const uint8_t *src, \
347 ptrdiff_t stride) \
348{ \
349 uint8_t halfH[72]; \
350 uint8_t halfHV[64]; \
351 \
352 put ## RND ## mpeg4_qpel8_h_lowpass(halfH, src, 8, stride, 9); \
353 put ## RND ## mpeg4_qpel8_v_lowpass(halfHV, halfH, 8, 8); \
354 OPNAME ## pixels8_l2_8(dst, halfH, halfHV, stride, 8, 8, 8); \
355} \
356 \
357static void OPNAME ## qpel8_mc23_c(uint8_t *dst, const uint8_t *src, \
358 ptrdiff_t stride) \
359{ \
360 uint8_t halfH[72]; \
361 uint8_t halfHV[64]; \
362 \
363 put ## RND ## mpeg4_qpel8_h_lowpass(halfH, src, 8, stride, 9); \
364 put ## RND ## mpeg4_qpel8_v_lowpass(halfHV, halfH, 8, 8); \
365 OPNAME ## pixels8_l2_8(dst, halfH + 8, halfHV, stride, 8, 8, 8); \
366} \
367 \
368void ff_ ## OPNAME ## qpel8_mc12_old_c(uint8_t *dst, const uint8_t *src, \
369 ptrdiff_t stride) \
370{ \
371 uint8_t full[16 * 9]; \
372 uint8_t halfH[72]; \
373 uint8_t halfV[64]; \
374 uint8_t halfHV[64]; \
375 \
376 copy_block9(full, src, 16, stride, 9); \
377 put ## RND ## mpeg4_qpel8_h_lowpass(halfH, full, 8, 16, 9); \
378 put ## RND ## mpeg4_qpel8_v_lowpass(halfV, full, 8, 16); \
379 put ## RND ## mpeg4_qpel8_v_lowpass(halfHV, halfH, 8, 8); \
380 OPNAME ## pixels8_l2_8(dst, halfV, halfHV, stride, 8, 8, 8); \
381} \
382 \
383static void OPNAME ## qpel8_mc12_c(uint8_t *dst, const uint8_t *src, \
384 ptrdiff_t stride) \
385{ \
386 uint8_t full[16 * 9]; \
387 uint8_t halfH[72]; \
388 \
389 copy_block9(full, src, 16, stride, 9); \
390 put ## RND ## mpeg4_qpel8_h_lowpass(halfH, full, 8, 16, 9); \
391 put ## RND ## pixels8_l2_8(halfH, halfH, full, 8, 8, 16, 9); \
392 OPNAME ## mpeg4_qpel8_v_lowpass(dst, halfH, stride, 8); \
393} \
394 \
395void ff_ ## OPNAME ## qpel8_mc32_old_c(uint8_t *dst, const uint8_t *src, \
396 ptrdiff_t stride) \
397{ \
398 uint8_t full[16 * 9]; \
399 uint8_t halfH[72]; \
400 uint8_t halfV[64]; \
401 uint8_t halfHV[64]; \
402 \
403 copy_block9(full, src, 16, stride, 9); \
404 put ## RND ## mpeg4_qpel8_h_lowpass(halfH, full, 8, 16, 9); \
405 put ## RND ## mpeg4_qpel8_v_lowpass(halfV, full + 1, 8, 16); \
406 put ## RND ## mpeg4_qpel8_v_lowpass(halfHV, halfH, 8, 8); \
407 OPNAME ## pixels8_l2_8(dst, halfV, halfHV, stride, 8, 8, 8); \
408} \
409 \
410static void OPNAME ## qpel8_mc32_c(uint8_t *dst, const uint8_t *src, \
411 ptrdiff_t stride) \
412{ \
413 uint8_t full[16 * 9]; \
414 uint8_t halfH[72]; \
415 \
416 copy_block9(full, src, 16, stride, 9); \
417 put ## RND ## mpeg4_qpel8_h_lowpass(halfH, full, 8, 16, 9); \
418 put ## RND ## pixels8_l2_8(halfH, halfH, full + 1, 8, 8, 16, 9); \
419 OPNAME ## mpeg4_qpel8_v_lowpass(dst, halfH, stride, 8); \
420} \
421 \
422static void OPNAME ## qpel8_mc22_c(uint8_t *dst, const uint8_t *src, \
423 ptrdiff_t stride) \
424{ \
425 uint8_t halfH[72]; \
426 \
427 put ## RND ## mpeg4_qpel8_h_lowpass(halfH, src, 8, stride, 9); \
428 OPNAME ## mpeg4_qpel8_v_lowpass(dst, halfH, stride, 8); \
429} \
430 \
431static void OPNAME ## qpel16_mc10_c(uint8_t *dst, const uint8_t *src, \
432 ptrdiff_t stride) \
433{ \
434 uint8_t half[256]; \
435 \
436 put ## RND ## mpeg4_qpel16_h_lowpass(half, src, 16, stride, 16); \
437 OPNAME ## pixels16_l2_8(dst, src, half, stride, stride, 16, 16); \
438} \
439 \
440static void OPNAME ## qpel16_mc20_c(uint8_t *dst, const uint8_t *src, \
441 ptrdiff_t stride) \
442{ \
443 OPNAME ## mpeg4_qpel16_h_lowpass(dst, src, stride, stride, 16); \
444} \
445 \
446static void OPNAME ## qpel16_mc30_c(uint8_t *dst, const uint8_t *src, \
447 ptrdiff_t stride) \
448{ \
449 uint8_t half[256]; \
450 \
451 put ## RND ## mpeg4_qpel16_h_lowpass(half, src, 16, stride, 16); \
452 OPNAME ## pixels16_l2_8(dst, src + 1, half, stride, stride, 16, 16); \
453} \
454 \
455static void OPNAME ## qpel16_mc01_c(uint8_t *dst, const uint8_t *src, \
456 ptrdiff_t stride) \
457{ \
458 uint8_t full[24 * 17]; \
459 uint8_t half[256]; \
460 \
461 copy_block17(full, src, 24, stride, 17); \
462 put ## RND ## mpeg4_qpel16_v_lowpass(half, full, 16, 24); \
463 OPNAME ## pixels16_l2_8(dst, full, half, stride, 24, 16, 16); \
464} \
465 \
466static void OPNAME ## qpel16_mc02_c(uint8_t *dst, const uint8_t *src, \
467 ptrdiff_t stride) \
468{ \
469 uint8_t full[24 * 17]; \
470 \
471 copy_block17(full, src, 24, stride, 17); \
472 OPNAME ## mpeg4_qpel16_v_lowpass(dst, full, stride, 24); \
473} \
474 \
475static void OPNAME ## qpel16_mc03_c(uint8_t *dst, const uint8_t *src, \
476 ptrdiff_t stride) \
477{ \
478 uint8_t full[24 * 17]; \
479 uint8_t half[256]; \
480 \
481 copy_block17(full, src, 24, stride, 17); \
482 put ## RND ## mpeg4_qpel16_v_lowpass(half, full, 16, 24); \
483 OPNAME ## pixels16_l2_8(dst, full + 24, half, stride, 24, 16, 16); \
484} \
485 \
486void ff_ ## OPNAME ## qpel16_mc11_old_c(uint8_t *dst, const uint8_t *src, \
487 ptrdiff_t stride) \
488{ \
489 uint8_t full[24 * 17]; \
490 uint8_t halfH[272]; \
491 uint8_t halfV[256]; \
492 uint8_t halfHV[256]; \
493 \
494 copy_block17(full, src, 24, stride, 17); \
495 put ## RND ## mpeg4_qpel16_h_lowpass(halfH, full, 16, 24, 17); \
496 put ## RND ## mpeg4_qpel16_v_lowpass(halfV, full, 16, 24); \
497 put ## RND ## mpeg4_qpel16_v_lowpass(halfHV, halfH, 16, 16); \
498 OPNAME ## pixels16_l4_8(dst, full, halfH, halfV, halfHV, \
499 stride, 24, 16, 16, 16, 16); \
500} \
501 \
502static void OPNAME ## qpel16_mc11_c(uint8_t *dst, const uint8_t *src, \
503 ptrdiff_t stride) \
504{ \
505 uint8_t full[24 * 17]; \
506 uint8_t halfH[272]; \
507 uint8_t halfHV[256]; \
508 \
509 copy_block17(full, src, 24, stride, 17); \
510 put ## RND ## mpeg4_qpel16_h_lowpass(halfH, full, 16, 24, 17); \
511 put ## RND ## pixels16_l2_8(halfH, halfH, full, 16, 16, 24, 17); \
512 put ## RND ## mpeg4_qpel16_v_lowpass(halfHV, halfH, 16, 16); \
513 OPNAME ## pixels16_l2_8(dst, halfH, halfHV, stride, 16, 16, 16); \
514} \
515 \
516void ff_ ## OPNAME ## qpel16_mc31_old_c(uint8_t *dst, const uint8_t *src, \
517 ptrdiff_t stride) \
518{ \
519 uint8_t full[24 * 17]; \
520 uint8_t halfH[272]; \
521 uint8_t halfV[256]; \
522 uint8_t halfHV[256]; \
523 \
524 copy_block17(full, src, 24, stride, 17); \
525 put ## RND ## mpeg4_qpel16_h_lowpass(halfH, full, 16, 24, 17); \
526 put ## RND ## mpeg4_qpel16_v_lowpass(halfV, full + 1, 16, 24); \
527 put ## RND ## mpeg4_qpel16_v_lowpass(halfHV, halfH, 16, 16); \
528 OPNAME ## pixels16_l4_8(dst, full + 1, halfH, halfV, halfHV, \
529 stride, 24, 16, 16, 16, 16); \
530} \
531 \
532static void OPNAME ## qpel16_mc31_c(uint8_t *dst, const uint8_t *src, \
533 ptrdiff_t stride) \
534{ \
535 uint8_t full[24 * 17]; \
536 uint8_t halfH[272]; \
537 uint8_t halfHV[256]; \
538 \
539 copy_block17(full, src, 24, stride, 17); \
540 put ## RND ## mpeg4_qpel16_h_lowpass(halfH, full, 16, 24, 17); \
541 put ## RND ## pixels16_l2_8(halfH, halfH, full + 1, 16, 16, 24, 17); \
542 put ## RND ## mpeg4_qpel16_v_lowpass(halfHV, halfH, 16, 16); \
543 OPNAME ## pixels16_l2_8(dst, halfH, halfHV, stride, 16, 16, 16); \
544} \
545 \
546void ff_ ## OPNAME ## qpel16_mc13_old_c(uint8_t *dst, const uint8_t *src, \
547 ptrdiff_t stride) \
548{ \
549 uint8_t full[24 * 17]; \
550 uint8_t halfH[272]; \
551 uint8_t halfV[256]; \
552 uint8_t halfHV[256]; \
553 \
554 copy_block17(full, src, 24, stride, 17); \
555 put ## RND ## mpeg4_qpel16_h_lowpass(halfH, full, 16, 24, 17); \
556 put ## RND ## mpeg4_qpel16_v_lowpass(halfV, full, 16, 24); \
557 put ## RND ## mpeg4_qpel16_v_lowpass(halfHV, halfH, 16, 16); \
558 OPNAME ## pixels16_l4_8(dst, full + 24, halfH + 16, halfV, halfHV, \
559 stride, 24, 16, 16, 16, 16); \
560} \
561 \
562static void OPNAME ## qpel16_mc13_c(uint8_t *dst, const uint8_t *src, \
563 ptrdiff_t stride) \
564{ \
565 uint8_t full[24 * 17]; \
566 uint8_t halfH[272]; \
567 uint8_t halfHV[256]; \
568 \
569 copy_block17(full, src, 24, stride, 17); \
570 put ## RND ## mpeg4_qpel16_h_lowpass(halfH, full, 16, 24, 17); \
571 put ## RND ## pixels16_l2_8(halfH, halfH, full, 16, 16, 24, 17); \
572 put ## RND ## mpeg4_qpel16_v_lowpass(halfHV, halfH, 16, 16); \
573 OPNAME ## pixels16_l2_8(dst, halfH + 16, halfHV, stride, 16, 16, 16); \
574} \
575 \
576void ff_ ## OPNAME ## qpel16_mc33_old_c(uint8_t *dst, const uint8_t *src, \
577 ptrdiff_t stride) \
578{ \
579 uint8_t full[24 * 17]; \
580 uint8_t halfH[272]; \
581 uint8_t halfV[256]; \
582 uint8_t halfHV[256]; \
583 \
584 copy_block17(full, src, 24, stride, 17); \
585 put ## RND ## mpeg4_qpel16_h_lowpass(halfH, full, 16, 24, 17); \
586 put ## RND ## mpeg4_qpel16_v_lowpass(halfV, full + 1, 16, 24); \
587 put ## RND ## mpeg4_qpel16_v_lowpass(halfHV, halfH, 16, 16); \
588 OPNAME ## pixels16_l4_8(dst, full + 25, halfH + 16, halfV, halfHV, \
589 stride, 24, 16, 16, 16, 16); \
590} \
591 \
592static void OPNAME ## qpel16_mc33_c(uint8_t *dst, const uint8_t *src, \
593 ptrdiff_t stride) \
594{ \
595 uint8_t full[24 * 17]; \
596 uint8_t halfH[272]; \
597 uint8_t halfHV[256]; \
598 \
599 copy_block17(full, src, 24, stride, 17); \
600 put ## RND ## mpeg4_qpel16_h_lowpass(halfH, full, 16, 24, 17); \
601 put ## RND ## pixels16_l2_8(halfH, halfH, full + 1, 16, 16, 24, 17); \
602 put ## RND ## mpeg4_qpel16_v_lowpass(halfHV, halfH, 16, 16); \
603 OPNAME ## pixels16_l2_8(dst, halfH + 16, halfHV, stride, 16, 16, 16); \
604} \
605 \
606static void OPNAME ## qpel16_mc21_c(uint8_t *dst, const uint8_t *src, \
607 ptrdiff_t stride) \
608{ \
609 uint8_t halfH[272]; \
610 uint8_t halfHV[256]; \
611 \
612 put ## RND ## mpeg4_qpel16_h_lowpass(halfH, src, 16, stride, 17); \
613 put ## RND ## mpeg4_qpel16_v_lowpass(halfHV, halfH, 16, 16); \
614 OPNAME ## pixels16_l2_8(dst, halfH, halfHV, stride, 16, 16, 16); \
615} \
616 \
617static void OPNAME ## qpel16_mc23_c(uint8_t *dst, const uint8_t *src, \
618 ptrdiff_t stride) \
619{ \
620 uint8_t halfH[272]; \
621 uint8_t halfHV[256]; \
622 \
623 put ## RND ## mpeg4_qpel16_h_lowpass(halfH, src, 16, stride, 17); \
624 put ## RND ## mpeg4_qpel16_v_lowpass(halfHV, halfH, 16, 16); \
625 OPNAME ## pixels16_l2_8(dst, halfH + 16, halfHV, stride, 16, 16, 16); \
626} \
627 \
628void ff_ ## OPNAME ## qpel16_mc12_old_c(uint8_t *dst, const uint8_t *src, \
629 ptrdiff_t stride) \
630{ \
631 uint8_t full[24 * 17]; \
632 uint8_t halfH[272]; \
633 uint8_t halfV[256]; \
634 uint8_t halfHV[256]; \
635 \
636 copy_block17(full, src, 24, stride, 17); \
637 put ## RND ## mpeg4_qpel16_h_lowpass(halfH, full, 16, 24, 17); \
638 put ## RND ## mpeg4_qpel16_v_lowpass(halfV, full, 16, 24); \
639 put ## RND ## mpeg4_qpel16_v_lowpass(halfHV, halfH, 16, 16); \
640 OPNAME ## pixels16_l2_8(dst, halfV, halfHV, stride, 16, 16, 16); \
641} \
642 \
643static void OPNAME ## qpel16_mc12_c(uint8_t *dst, const uint8_t *src, \
644 ptrdiff_t stride) \
645{ \
646 uint8_t full[24 * 17]; \
647 uint8_t halfH[272]; \
648 \
649 copy_block17(full, src, 24, stride, 17); \
650 put ## RND ## mpeg4_qpel16_h_lowpass(halfH, full, 16, 24, 17); \
651 put ## RND ## pixels16_l2_8(halfH, halfH, full, 16, 16, 24, 17); \
652 OPNAME ## mpeg4_qpel16_v_lowpass(dst, halfH, stride, 16); \
653} \
654 \
655void ff_ ## OPNAME ## qpel16_mc32_old_c(uint8_t *dst, const uint8_t *src, \
656 ptrdiff_t stride) \
657{ \
658 uint8_t full[24 * 17]; \
659 uint8_t halfH[272]; \
660 uint8_t halfV[256]; \
661 uint8_t halfHV[256]; \
662 \
663 copy_block17(full, src, 24, stride, 17); \
664 put ## RND ## mpeg4_qpel16_h_lowpass(halfH, full, 16, 24, 17); \
665 put ## RND ## mpeg4_qpel16_v_lowpass(halfV, full + 1, 16, 24); \
666 put ## RND ## mpeg4_qpel16_v_lowpass(halfHV, halfH, 16, 16); \
667 OPNAME ## pixels16_l2_8(dst, halfV, halfHV, stride, 16, 16, 16); \
668} \
669 \
670static void OPNAME ## qpel16_mc32_c(uint8_t *dst, const uint8_t *src, \
671 ptrdiff_t stride) \
672{ \
673 uint8_t full[24 * 17]; \
674 uint8_t halfH[272]; \
675 \
676 copy_block17(full, src, 24, stride, 17); \
677 put ## RND ## mpeg4_qpel16_h_lowpass(halfH, full, 16, 24, 17); \
678 put ## RND ## pixels16_l2_8(halfH, halfH, full + 1, 16, 16, 24, 17); \
679 OPNAME ## mpeg4_qpel16_v_lowpass(dst, halfH, stride, 16); \
680} \
681 \
682static void OPNAME ## qpel16_mc22_c(uint8_t *dst, const uint8_t *src, \
683 ptrdiff_t stride) \
684{ \
685 uint8_t halfH[272]; \
686 \
687 put ## RND ## mpeg4_qpel16_h_lowpass(halfH, src, 16, stride, 17); \
688 OPNAME ## mpeg4_qpel16_v_lowpass(dst, halfH, stride, 16); \
689}
690
691#define op_avg(a, b) a = (((a) + cm[((b) + 16) >> 5] + 1) >> 1)
692#define op_put(a, b) a = cm[((b) + 16) >> 5]
693#define op_put_no_rnd(a, b) a = cm[((b) + 15) >> 5]
694
695QPEL_MC(0, put_, _, op_put)
696QPEL_MC(1, put_no_rnd_, _no_rnd_, op_put_no_rnd)
697QPEL_MC(0, avg_, _, op_avg)
698
699#undef op_avg
700#undef op_put
701#undef op_put_no_rnd
702
703void ff_put_pixels8x8_c(uint8_t *dst, const uint8_t *src, ptrdiff_t stride)
704{
705 put_pixels8_8_c(dst, src, stride, 8);
706}
707
708void ff_avg_pixels8x8_c(uint8_t *dst, const uint8_t *src, ptrdiff_t stride)
709{
710 avg_pixels8_8_c(dst, src, stride, 8);
711}
712
713void ff_put_pixels16x16_c(uint8_t *dst, const uint8_t *src, ptrdiff_t stride)
714{
715 put_pixels16_8_c(dst, src, stride, 16);
716}
717
718void ff_avg_pixels16x16_c(uint8_t *dst, const uint8_t *src, ptrdiff_t stride)
719{
720 avg_pixels16_8_c(dst, src, stride, 16);
721}
722
723#define put_qpel8_mc00_c ff_put_pixels8x8_c
724#define avg_qpel8_mc00_c ff_avg_pixels8x8_c
725#define put_qpel16_mc00_c ff_put_pixels16x16_c
726#define avg_qpel16_mc00_c ff_avg_pixels16x16_c
727#define put_no_rnd_qpel8_mc00_c ff_put_pixels8x8_c
728#define put_no_rnd_qpel16_mc00_c ff_put_pixels16x16_c
729
730void ff_put_pixels8_l2_8(uint8_t *dst, const uint8_t *src1, const uint8_t *src2,
731 int dst_stride, int src_stride1, int src_stride2,
732 int h)
733{
734 put_pixels8_l2_8(dst, src1, src2, dst_stride, src_stride1, src_stride2, h);
735
736}
737
738#if CONFIG_DIRAC_DECODER
739#define DIRAC_MC(OPNAME)\
740void ff_ ## OPNAME ## _dirac_pixels8_c(uint8_t *dst, const uint8_t *src[5], int stride, int h)\
741{\
742 OPNAME ## _pixels8_8_c(dst, src[0], stride, h);\
743}\
744void ff_ ## OPNAME ## _dirac_pixels16_c(uint8_t *dst, const uint8_t *src[5], int stride, int h)\
745{\
746 OPNAME ## _pixels16_8_c(dst, src[0], stride, h);\
747}\
748void ff_ ## OPNAME ## _dirac_pixels32_c(uint8_t *dst, const uint8_t *src[5], int stride, int h)\
749{\
750 OPNAME ## _pixels16_8_c(dst , src[0] , stride, h);\
751 OPNAME ## _pixels16_8_c(dst+16, src[0]+16, stride, h);\
752}\
753void ff_ ## OPNAME ## _dirac_pixels8_l2_c(uint8_t *dst, const uint8_t *src[5], int stride, int h)\
754{\
755 OPNAME ## _pixels8_l2_8(dst, src[0], src[1], stride, stride, stride, h);\
756}\
757void ff_ ## OPNAME ## _dirac_pixels16_l2_c(uint8_t *dst, const uint8_t *src[5], int stride, int h)\
758{\
759 OPNAME ## _pixels16_l2_8(dst, src[0], src[1], stride, stride, stride, h);\
760}\
761void ff_ ## OPNAME ## _dirac_pixels32_l2_c(uint8_t *dst, const uint8_t *src[5], int stride, int h)\
762{\
763 OPNAME ## _pixels16_l2_8(dst , src[0] , src[1] , stride, stride, stride, h);\
764 OPNAME ## _pixels16_l2_8(dst+16, src[0]+16, src[1]+16, stride, stride, stride, h);\
765}\
766void ff_ ## OPNAME ## _dirac_pixels8_l4_c(uint8_t *dst, const uint8_t *src[5], int stride, int h)\
767{\
768 OPNAME ## _pixels8_l4_8(dst, src[0], src[1], src[2], src[3], stride, stride, stride, stride, stride, h);\
769}\
770void ff_ ## OPNAME ## _dirac_pixels16_l4_c(uint8_t *dst, const uint8_t *src[5], int stride, int h)\
771{\
772 OPNAME ## _pixels16_l4_8(dst, src[0], src[1], src[2], src[3], stride, stride, stride, stride, stride, h);\
773}\
774void ff_ ## OPNAME ## _dirac_pixels32_l4_c(uint8_t *dst, const uint8_t *src[5], int stride, int h)\
775{\
776 OPNAME ## _pixels16_l4_8(dst , src[0] , src[1] , src[2] , src[3] , stride, stride, stride, stride, stride, h);\
777 OPNAME ## _pixels16_l4_8(dst+16, src[0]+16, src[1]+16, src[2]+16, src[3]+16, stride, stride, stride, stride, stride, h);\
778}
779DIRAC_MC(put)
780DIRAC_MC(avg)
781#endif
782
783av_cold void ff_qpeldsp_init(QpelDSPContext *c)
784{
785#define dspfunc(PFX, IDX, NUM) \
786 c->PFX ## _pixels_tab[IDX][0] = PFX ## NUM ## _mc00_c; \
787 c->PFX ## _pixels_tab[IDX][1] = PFX ## NUM ## _mc10_c; \
788 c->PFX ## _pixels_tab[IDX][2] = PFX ## NUM ## _mc20_c; \
789 c->PFX ## _pixels_tab[IDX][3] = PFX ## NUM ## _mc30_c; \
790 c->PFX ## _pixels_tab[IDX][4] = PFX ## NUM ## _mc01_c; \
791 c->PFX ## _pixels_tab[IDX][5] = PFX ## NUM ## _mc11_c; \
792 c->PFX ## _pixels_tab[IDX][6] = PFX ## NUM ## _mc21_c; \
793 c->PFX ## _pixels_tab[IDX][7] = PFX ## NUM ## _mc31_c; \
794 c->PFX ## _pixels_tab[IDX][8] = PFX ## NUM ## _mc02_c; \
795 c->PFX ## _pixels_tab[IDX][9] = PFX ## NUM ## _mc12_c; \
796 c->PFX ## _pixels_tab[IDX][10] = PFX ## NUM ## _mc22_c; \
797 c->PFX ## _pixels_tab[IDX][11] = PFX ## NUM ## _mc32_c; \
798 c->PFX ## _pixels_tab[IDX][12] = PFX ## NUM ## _mc03_c; \
799 c->PFX ## _pixels_tab[IDX][13] = PFX ## NUM ## _mc13_c; \
800 c->PFX ## _pixels_tab[IDX][14] = PFX ## NUM ## _mc23_c; \
801 c->PFX ## _pixels_tab[IDX][15] = PFX ## NUM ## _mc33_c
802
803 dspfunc(put_qpel, 0, 16);
804 dspfunc(put_qpel, 1, 8);
805
806 dspfunc(put_no_rnd_qpel, 0, 16);
807 dspfunc(put_no_rnd_qpel, 1, 8);
808
809 dspfunc(avg_qpel, 0, 16);
810 dspfunc(avg_qpel, 1, 8);
811
812 if (ARCH_X86)
813 ff_qpeldsp_init_x86(c);
814}
ffmpeg Piment Noir Debian package repository