?  search:
summary | log | graphiclog1 | graphiclog2 | commit | commitdiff | tree | refs | edit | fork
blame | history | raw | HEAD
Merge commit 'origin/master' into ordered_chapters
[FFMpeg-mirror/ordered_chapters.git] / libavcodec / sh4 / qpel.c
blob2069bd38c1b0c77b4202d812ac1b1c34b7479d7f
1 /*
2 * This is optimized for sh, which have post increment addressing (*p++).
3 * Some CPU may be index (p[n]) faster than post increment (*p++).
4 *
5 * copyright (c) 2001-2003 BERO <bero@geocities.co.jp>
6 *
7 * This file is part of FFmpeg.
8 *
9 * FFmpeg is free software; you can redistribute it and/or
10 * modify it under the terms of the GNU Lesser General Public
11 * License as published by the Free Software Foundation; either
12 * version 2.1 of the License, or (at your option) any later version.
13 *
14 * FFmpeg is distributed in the hope that it will be useful,
15 * but WITHOUT ANY WARRANTY; without even the implied warranty of
16 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
17 * Lesser General Public License for more details.
18 *
19 * You should have received a copy of the GNU Lesser General Public
20 * License along with FFmpeg; if not, write to the Free Software
21 * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
22 */
23
24 #define PIXOP2(OPNAME, OP) \
25 \
26 static inline void OPNAME ## _pixels4_l2_aligned(uint8_t *dst, const uint8_t *src1, const uint8_t *src2, int dst_stride, int src_stride1, int src_stride2, int h) \
27 {\
28 do {\
29 OP(LP(dst ),rnd_avg32(LP(src1 ),LP(src2 )) ); \
30 src1+=src_stride1; \
31 src2+=src_stride2; \
32 dst+=dst_stride; \
33 } while(--h); \
34 }\
35 \
36 static inline void OPNAME ## _pixels4_l2_aligned2(uint8_t *dst, const uint8_t *src1, const uint8_t *src2, int dst_stride, int src_stride1, int src_stride2, int h) \
37 {\
38 do {\
39 OP(LP(dst ),rnd_avg32(AV_RN32(src1 ),LP(src2 )) ); \
40 src1+=src_stride1; \
41 src2+=src_stride2; \
42 dst+=dst_stride; \
43 } while(--h); \
44 }\
45 \
46 static inline void OPNAME ## _no_rnd_pixels16_l2_aligned2(uint8_t *dst, const uint8_t *src1, const uint8_t *src2, int dst_stride, int src_stride1, int src_stride2, int h) \
47 {\
48 do {\
49 OP(LP(dst ),no_rnd_avg32(AV_RN32(src1 ),LP(src2 )) ); \
50 OP(LP(dst+4),no_rnd_avg32(AV_RN32(src1+4),LP(src2+4)) ); \
51 OP(LP(dst+8),no_rnd_avg32(AV_RN32(src1+8),LP(src2+8)) ); \
52 OP(LP(dst+12),no_rnd_avg32(AV_RN32(src1+12),LP(src2+12)) ); \
53 src1+=src_stride1; \
54 src2+=src_stride2; \
55 dst+=dst_stride; \
56 } while(--h); \
57 }\
58 \
59 static inline void OPNAME ## _pixels16_l2_aligned2(uint8_t *dst, const uint8_t *src1, const uint8_t *src2, int dst_stride, int src_stride1, int src_stride2, int h) \
60 {\
61 do {\
62 OP(LP(dst ),rnd_avg32(AV_RN32(src1 ),LP(src2 )) ); \
63 OP(LP(dst+4),rnd_avg32(AV_RN32(src1+4),LP(src2+4)) ); \
64 OP(LP(dst+8),rnd_avg32(AV_RN32(src1+8),LP(src2+8)) ); \
65 OP(LP(dst+12),rnd_avg32(AV_RN32(src1+12),LP(src2+12)) ); \
66 src1+=src_stride1; \
67 src2+=src_stride2; \
68 dst+=dst_stride; \
69 } while(--h); \
70 }\
71 \
72 static inline void OPNAME ## _no_rnd_pixels8_l2_aligned2(uint8_t *dst, const uint8_t *src1, const uint8_t *src2, int dst_stride, int src_stride1, int src_stride2, int h) \
73 {\
74 do { /* onlye src2 aligned */\
75 OP(LP(dst ),no_rnd_avg32(AV_RN32(src1 ),LP(src2 )) ); \
76 OP(LP(dst+4),no_rnd_avg32(AV_RN32(src1+4),LP(src2+4)) ); \
77 src1+=src_stride1; \
78 src2+=src_stride2; \
79 dst+=dst_stride; \
80 } while(--h); \
81 }\
82 \
83 static inline void OPNAME ## _pixels8_l2_aligned2(uint8_t *dst, const uint8_t *src1, const uint8_t *src2, int dst_stride, int src_stride1, int src_stride2, int h) \
84 {\
85 do {\
86 OP(LP(dst ),rnd_avg32(AV_RN32(src1 ),LP(src2 )) ); \
87 OP(LP(dst+4),rnd_avg32(AV_RN32(src1+4),LP(src2+4)) ); \
88 src1+=src_stride1; \
89 src2+=src_stride2; \
90 dst+=dst_stride; \
91 } while(--h); \
92 }\
93 \
94 static inline void OPNAME ## _no_rnd_pixels8_l2_aligned(uint8_t *dst, const uint8_t *src1, const uint8_t *src2, int dst_stride, int src_stride1, int src_stride2, int h) \
95 {\
96 do {\
97 OP(LP(dst ),no_rnd_avg32(LP(src1 ),LP(src2 )) ); \
98 OP(LP(dst+4),no_rnd_avg32(LP(src1+4),LP(src2+4)) ); \
99 src1+=src_stride1; \
100 src2+=src_stride2; \
101 dst+=dst_stride; \
102 } while(--h); \
103 }\
104 \
105 static inline void OPNAME ## _pixels8_l2_aligned(uint8_t *dst, const uint8_t *src1, const uint8_t *src2, int dst_stride, int src_stride1, int src_stride2, int h) \
106 {\
107 do {\
108 OP(LP(dst ),rnd_avg32(LP(src1 ),LP(src2 )) ); \
109 OP(LP(dst+4),rnd_avg32(LP(src1+4),LP(src2+4)) ); \
110 src1+=src_stride1; \
111 src2+=src_stride2; \
112 dst+=dst_stride; \
113 } while(--h); \
114 }\
115 \
116 static inline void OPNAME ## _no_rnd_pixels16_l2_aligned(uint8_t *dst, const uint8_t *src1, const uint8_t *src2, int dst_stride, int src_stride1, int src_stride2, int h) \
117 {\
118 do {\
119 OP(LP(dst ),no_rnd_avg32(LP(src1 ),LP(src2 )) ); \
120 OP(LP(dst+4),no_rnd_avg32(LP(src1+4),LP(src2+4)) ); \
121 OP(LP(dst+8),no_rnd_avg32(LP(src1+8),LP(src2+8)) ); \
122 OP(LP(dst+12),no_rnd_avg32(LP(src1+12),LP(src2+12)) ); \
123 src1+=src_stride1; \
124 src2+=src_stride2; \
125 dst+=dst_stride; \
126 } while(--h); \
127 }\
128 \
129 static inline void OPNAME ## _pixels16_l2_aligned(uint8_t *dst, const uint8_t *src1, const uint8_t *src2, int dst_stride, int src_stride1, int src_stride2, int h) \
130 {\
131 do {\
132 OP(LP(dst ),rnd_avg32(LP(src1 ),LP(src2 )) ); \
133 OP(LP(dst+4),rnd_avg32(LP(src1+4),LP(src2+4)) ); \
134 OP(LP(dst+8),rnd_avg32(LP(src1+8),LP(src2+8)) ); \
135 OP(LP(dst+12),rnd_avg32(LP(src1+12),LP(src2+12)) ); \
136 src1+=src_stride1; \
137 src2+=src_stride2; \
138 dst+=dst_stride; \
139 } while(--h); \
140 }\
141 \
142 static inline void OPNAME ## _no_rnd_pixels16_l2_aligned1(uint8_t *dst, const uint8_t *src1, const uint8_t *src2, int dst_stride, int src_stride1, int src_stride2, int h) \
143 { OPNAME ## _no_rnd_pixels16_l2_aligned2(dst,src2,src1,dst_stride,src_stride2,src_stride1,h); } \
144 \
145 static inline void OPNAME ## _pixels16_l2_aligned1(uint8_t *dst, const uint8_t *src1, const uint8_t *src2, int dst_stride, int src_stride1, int src_stride2, int h) \
146 { OPNAME ## _pixels16_l2_aligned2(dst,src2,src1,dst_stride,src_stride2,src_stride1,h); } \
147 \
148 static inline void OPNAME ## _no_rnd_pixels8_l2_aligned1(uint8_t *dst, const uint8_t *src1, const uint8_t *src2, int dst_stride, int src_stride1, int src_stride2, int h) \
149 { OPNAME ## _no_rnd_pixels8_l2_aligned2(dst,src2,src1,dst_stride,src_stride2,src_stride1,h); } \
150 \
151 static inline void OPNAME ## _pixels8_l2_aligned1(uint8_t *dst, const uint8_t *src1, const uint8_t *src2, int dst_stride, int src_stride1, int src_stride2, int h) \
152 { OPNAME ## _pixels8_l2_aligned2(dst,src2,src1,dst_stride,src_stride2,src_stride1,h); } \
153 \
154 static inline void OPNAME ## _pixels8_l4_aligned(uint8_t *dst, const uint8_t *src1, uint8_t *src2, uint8_t *src3, uint8_t *src4,int dst_stride, int src_stride1, int src_stride2,int src_stride3,int src_stride4, int h){\
155 do { \
156 uint32_t a0,a1,a2,a3; \
157 UNPACK(a0,a1,LP(src1),LP(src2)); \
158 UNPACK(a2,a3,LP(src3),LP(src4)); \
159 OP(LP(dst),rnd_PACK(a0,a1,a2,a3)); \
160 UNPACK(a0,a1,LP(src1+4),LP(src2+4)); \
161 UNPACK(a2,a3,LP(src3+4),LP(src4+4)); \
162 OP(LP(dst+4),rnd_PACK(a0,a1,a2,a3)); \
163 src1+=src_stride1;\
164 src2+=src_stride2;\
165 src3+=src_stride3;\
166 src4+=src_stride4;\
167 dst+=dst_stride;\
168 } while(--h); \
169 } \
170 \
171 static inline void OPNAME ## _no_rnd_pixels8_l4_aligned(uint8_t *dst, const uint8_t *src1, uint8_t *src2, uint8_t *src3, uint8_t *src4,int dst_stride, int src_stride1, int src_stride2,int src_stride3,int src_stride4, int h){\
172 do { \
173 uint32_t a0,a1,a2,a3; \
174 UNPACK(a0,a1,LP(src1),LP(src2)); \
175 UNPACK(a2,a3,LP(src3),LP(src4)); \
176 OP(LP(dst),no_rnd_PACK(a0,a1,a2,a3)); \
177 UNPACK(a0,a1,LP(src1+4),LP(src2+4)); \
178 UNPACK(a2,a3,LP(src3+4),LP(src4+4)); \
179 OP(LP(dst+4),no_rnd_PACK(a0,a1,a2,a3)); \
180 src1+=src_stride1;\
181 src2+=src_stride2;\
182 src3+=src_stride3;\
183 src4+=src_stride4;\
184 dst+=dst_stride;\
185 } while(--h); \
186 } \
187 \
188 static inline void OPNAME ## _pixels8_l4_aligned0(uint8_t *dst, const uint8_t *src1, uint8_t *src2, uint8_t *src3, uint8_t *src4,int dst_stride, int src_stride1, int src_stride2,int src_stride3,int src_stride4, int h){\
189 do { \
190 uint32_t a0,a1,a2,a3; /* src1 only not aligned */\
191 UNPACK(a0,a1,AV_RN32(src1),LP(src2)); \
192 UNPACK(a2,a3,LP(src3),LP(src4)); \
193 OP(LP(dst),rnd_PACK(a0,a1,a2,a3)); \
194 UNPACK(a0,a1,AV_RN32(src1+4),LP(src2+4)); \
195 UNPACK(a2,a3,LP(src3+4),LP(src4+4)); \
196 OP(LP(dst+4),rnd_PACK(a0,a1,a2,a3)); \
197 src1+=src_stride1;\
198 src2+=src_stride2;\
199 src3+=src_stride3;\
200 src4+=src_stride4;\
201 dst+=dst_stride;\
202 } while(--h); \
203 } \
204 \
205 static inline void OPNAME ## _no_rnd_pixels8_l4_aligned0(uint8_t *dst, const uint8_t *src1, uint8_t *src2, uint8_t *src3, uint8_t *src4,int dst_stride, int src_stride1, int src_stride2,int src_stride3,int src_stride4, int h){\
206 do { \
207 uint32_t a0,a1,a2,a3; \
208 UNPACK(a0,a1,AV_RN32(src1),LP(src2)); \
209 UNPACK(a2,a3,LP(src3),LP(src4)); \
210 OP(LP(dst),no_rnd_PACK(a0,a1,a2,a3)); \
211 UNPACK(a0,a1,AV_RN32(src1+4),LP(src2+4)); \
212 UNPACK(a2,a3,LP(src3+4),LP(src4+4)); \
213 OP(LP(dst+4),no_rnd_PACK(a0,a1,a2,a3)); \
214 src1+=src_stride1;\
215 src2+=src_stride2;\
216 src3+=src_stride3;\
217 src4+=src_stride4;\
218 dst+=dst_stride;\
219 } while(--h); \
220 } \
221 \
222 static inline void OPNAME ## _pixels16_l4_aligned(uint8_t *dst, const uint8_t *src1, uint8_t *src2, uint8_t *src3, uint8_t *src4,int dst_stride, int src_stride1, int src_stride2,int src_stride3,int src_stride4, int h){\
223 do { \
224 uint32_t a0,a1,a2,a3; \
225 UNPACK(a0,a1,LP(src1),LP(src2)); \
226 UNPACK(a2,a3,LP(src3),LP(src4)); \
227 OP(LP(dst),rnd_PACK(a0,a1,a2,a3)); \
228 UNPACK(a0,a1,LP(src1+4),LP(src2+4)); \
229 UNPACK(a2,a3,LP(src3+4),LP(src4+4)); \
230 OP(LP(dst+8),rnd_PACK(a0,a1,a2,a3)); \
231 UNPACK(a0,a1,LP(src1+8),LP(src2+8)); \
232 UNPACK(a2,a3,LP(src3+8),LP(src4+8)); \
233 OP(LP(dst+8),rnd_PACK(a0,a1,a2,a3)); \
234 UNPACK(a0,a1,LP(src1+12),LP(src2+12)); \
235 UNPACK(a2,a3,LP(src3+12),LP(src4+12)); \
236 OP(LP(dst+12),rnd_PACK(a0,a1,a2,a3)); \
237 src1+=src_stride1;\
238 src2+=src_stride2;\
239 src3+=src_stride3;\
240 src4+=src_stride4;\
241 dst+=dst_stride;\
242 } while(--h); \
243 } \
244 \
245 static inline void OPNAME ## _no_rnd_pixels16_l4_aligned(uint8_t *dst, const uint8_t *src1, uint8_t *src2, uint8_t *src3, uint8_t *src4,int dst_stride, int src_stride1, int src_stride2,int src_stride3,int src_stride4, int h){\
246 do { \
247 uint32_t a0,a1,a2,a3; \
248 UNPACK(a0,a1,LP(src1),LP(src2)); \
249 UNPACK(a2,a3,LP(src3),LP(src4)); \
250 OP(LP(dst),no_rnd_PACK(a0,a1,a2,a3)); \
251 UNPACK(a0,a1,LP(src1+4),LP(src2+4)); \
252 UNPACK(a2,a3,LP(src3+4),LP(src4+4)); \
253 OP(LP(dst+4),no_rnd_PACK(a0,a1,a2,a3)); \
254 UNPACK(a0,a1,LP(src1+8),LP(src2+8)); \
255 UNPACK(a2,a3,LP(src3+8),LP(src4+8)); \
256 OP(LP(dst+8),no_rnd_PACK(a0,a1,a2,a3)); \
257 UNPACK(a0,a1,LP(src1+12),LP(src2+12)); \
258 UNPACK(a2,a3,LP(src3+12),LP(src4+12)); \
259 OP(LP(dst+12),no_rnd_PACK(a0,a1,a2,a3)); \
260 src1+=src_stride1;\
261 src2+=src_stride2;\
262 src3+=src_stride3;\
263 src4+=src_stride4;\
264 dst+=dst_stride;\
265 } while(--h); \
266 } \
267 \
268 static inline void OPNAME ## _pixels16_l4_aligned0(uint8_t *dst, const uint8_t *src1, uint8_t *src2, uint8_t *src3, uint8_t *src4,int dst_stride, int src_stride1, int src_stride2,int src_stride3,int src_stride4, int h){\
269 do { /* src1 is unaligned */\
270 uint32_t a0,a1,a2,a3; \
271 UNPACK(a0,a1,AV_RN32(src1),LP(src2)); \
272 UNPACK(a2,a3,LP(src3),LP(src4)); \
273 OP(LP(dst),rnd_PACK(a0,a1,a2,a3)); \
274 UNPACK(a0,a1,AV_RN32(src1+4),LP(src2+4)); \
275 UNPACK(a2,a3,LP(src3+4),LP(src4+4)); \
276 OP(LP(dst+8),rnd_PACK(a0,a1,a2,a3)); \
277 UNPACK(a0,a1,AV_RN32(src1+8),LP(src2+8)); \
278 UNPACK(a2,a3,LP(src3+8),LP(src4+8)); \
279 OP(LP(dst+8),rnd_PACK(a0,a1,a2,a3)); \
280 UNPACK(a0,a1,AV_RN32(src1+12),LP(src2+12)); \
281 UNPACK(a2,a3,LP(src3+12),LP(src4+12)); \
282 OP(LP(dst+12),rnd_PACK(a0,a1,a2,a3)); \
283 src1+=src_stride1;\
284 src2+=src_stride2;\
285 src3+=src_stride3;\
286 src4+=src_stride4;\
287 dst+=dst_stride;\
288 } while(--h); \
289 } \
290 \
291 static inline void OPNAME ## _no_rnd_pixels16_l4_aligned0(uint8_t *dst, const uint8_t *src1, uint8_t *src2, uint8_t *src3, uint8_t *src4,int dst_stride, int src_stride1, int src_stride2,int src_stride3,int src_stride4, int h){\
292 do { \
293 uint32_t a0,a1,a2,a3; \
294 UNPACK(a0,a1,AV_RN32(src1),LP(src2)); \
295 UNPACK(a2,a3,LP(src3),LP(src4)); \
296 OP(LP(dst),no_rnd_PACK(a0,a1,a2,a3)); \
297 UNPACK(a0,a1,AV_RN32(src1+4),LP(src2+4)); \
298 UNPACK(a2,a3,LP(src3+4),LP(src4+4)); \
299 OP(LP(dst+4),no_rnd_PACK(a0,a1,a2,a3)); \
300 UNPACK(a0,a1,AV_RN32(src1+8),LP(src2+8)); \
301 UNPACK(a2,a3,LP(src3+8),LP(src4+8)); \
302 OP(LP(dst+8),no_rnd_PACK(a0,a1,a2,a3)); \
303 UNPACK(a0,a1,AV_RN32(src1+12),LP(src2+12)); \
304 UNPACK(a2,a3,LP(src3+12),LP(src4+12)); \
305 OP(LP(dst+12),no_rnd_PACK(a0,a1,a2,a3)); \
306 src1+=src_stride1;\
307 src2+=src_stride2;\
308 src3+=src_stride3;\
309 src4+=src_stride4;\
310 dst+=dst_stride;\
311 } while(--h); \
312 } \
313 \
314
315 #define op_avg(a, b) a = rnd_avg32(a,b)
316 #define op_put(a, b) a = b
317
318 PIXOP2(avg, op_avg)
319 PIXOP2(put, op_put)
320 #undef op_avg
321 #undef op_put
322
323 #define avg2(a,b) ((a+b+1)>>1)
324 #define avg4(a,b,c,d) ((a+b+c+d+2)>>2)
325
326
327 static void gmc1_c(uint8_t *dst, uint8_t *src, int stride, int h, int x16, int y16, int rounder)
328 {
329 const int A=(16-x16)*(16-y16);
330 const int B=( x16)*(16-y16);
331 const int C=(16-x16)*( y16);
332 const int D=( x16)*( y16);
333
334 do {
335 int t0,t1,t2,t3;
336 uint8_t *s0 = src;
337 uint8_t *s1 = src+stride;
338 t0 = *s0++; t2 = *s1++;
339 t1 = *s0++; t3 = *s1++;
340 dst[0]= (A*t0 + B*t1 + C*t2 + D*t3 + rounder)>>8;
341 t0 = *s0++; t2 = *s1++;
342 dst[1]= (A*t1 + B*t0 + C*t3 + D*t2 + rounder)>>8;
343 t1 = *s0++; t3 = *s1++;
344 dst[2]= (A*t0 + B*t1 + C*t2 + D*t3 + rounder)>>8;
345 t0 = *s0++; t2 = *s1++;
346 dst[3]= (A*t1 + B*t0 + C*t3 + D*t2 + rounder)>>8;
347 t1 = *s0++; t3 = *s1++;
348 dst[4]= (A*t0 + B*t1 + C*t2 + D*t3 + rounder)>>8;
349 t0 = *s0++; t2 = *s1++;
350 dst[5]= (A*t1 + B*t0 + C*t3 + D*t2 + rounder)>>8;
351 t1 = *s0++; t3 = *s1++;
352 dst[6]= (A*t0 + B*t1 + C*t2 + D*t3 + rounder)>>8;
353 t0 = *s0++; t2 = *s1++;
354 dst[7]= (A*t1 + B*t0 + C*t3 + D*t2 + rounder)>>8;
355 dst+= stride;
356 src+= stride;
357 }while(--h);
358 }
359
360 static void gmc_c(uint8_t *dst, uint8_t *src, int stride, int h, int ox, int oy,
361 int dxx, int dxy, int dyx, int dyy, int shift, int r, int width, int height)
362 {
363 int y, vx, vy;
364 const int s= 1<<shift;
365
366 width--;
367 height--;
368
369 for(y=0; y<h; y++){
370 int x;
371
372 vx= ox;
373 vy= oy;
374 for(x=0; x<8; x++){ //XXX FIXME optimize
375 int src_x, src_y, frac_x, frac_y, index;
376
377 src_x= vx>>16;
378 src_y= vy>>16;
379 frac_x= src_x&(s-1);
380 frac_y= src_y&(s-1);
381 src_x>>=shift;
382 src_y>>=shift;
383
384 if((unsigned)src_x < width){
385 if((unsigned)src_y < height){
386 index= src_x + src_y*stride;
387 dst[y*stride + x]= ( ( src[index ]*(s-frac_x)
388 + src[index +1]* frac_x )*(s-frac_y)
389 + ( src[index+stride ]*(s-frac_x)
390 + src[index+stride+1]* frac_x )* frac_y
391 + r)>>(shift*2);
392 }else{
393 index= src_x + av_clip(src_y, 0, height)*stride;
394 dst[y*stride + x]= ( ( src[index ]*(s-frac_x)
395 + src[index +1]* frac_x )*s
396 + r)>>(shift*2);
397 }
398 }else{
399 if((unsigned)src_y < height){
400 index= av_clip(src_x, 0, width) + src_y*stride;
401 dst[y*stride + x]= ( ( src[index ]*(s-frac_y)
402 + src[index+stride ]* frac_y )*s
403 + r)>>(shift*2);
404 }else{
405 index= av_clip(src_x, 0, width) + av_clip(src_y, 0, height)*stride;
406 dst[y*stride + x]= src[index ];
407 }
408 }
409
410 vx+= dxx;
411 vy+= dyx;
412 }
413 ox += dxy;
414 oy += dyy;
415 }
416 }
417 #define H264_CHROMA_MC(OPNAME, OP)\
418 static void OPNAME ## h264_chroma_mc2_sh4(uint8_t *dst/*align 8*/, uint8_t *src/*align 1*/, int stride, int h, int x, int y){\
419 const int A=(8-x)*(8-y);\
420 const int B=( x)*(8-y);\
421 const int C=(8-x)*( y);\
422 const int D=( x)*( y);\
423 \
424 assert(x<8 && y<8 && x>=0 && y>=0);\
425 \
426 do {\
427 int t0,t1,t2,t3; \
428 uint8_t *s0 = src; \
429 uint8_t *s1 = src+stride; \
430 t0 = *s0++; t2 = *s1++; \
431 t1 = *s0++; t3 = *s1++; \
432 OP(dst[0], (A*t0 + B*t1 + C*t2 + D*t3));\
433 t0 = *s0++; t2 = *s1++; \
434 OP(dst[1], (A*t1 + B*t0 + C*t3 + D*t2));\
435 dst+= stride;\
436 src+= stride;\
437 }while(--h);\
438 }\
439 \
440 static void OPNAME ## h264_chroma_mc4_sh4(uint8_t *dst/*align 8*/, uint8_t *src/*align 1*/, int stride, int h, int x, int y){\
441 const int A=(8-x)*(8-y);\
442 const int B=( x)*(8-y);\
443 const int C=(8-x)*( y);\
444 const int D=( x)*( y);\
445 \
446 assert(x<8 && y<8 && x>=0 && y>=0);\
447 \
448 do {\
449 int t0,t1,t2,t3; \
450 uint8_t *s0 = src; \
451 uint8_t *s1 = src+stride; \
452 t0 = *s0++; t2 = *s1++; \
453 t1 = *s0++; t3 = *s1++; \
454 OP(dst[0], (A*t0 + B*t1 + C*t2 + D*t3));\
455 t0 = *s0++; t2 = *s1++; \
456 OP(dst[1], (A*t1 + B*t0 + C*t3 + D*t2));\
457 t1 = *s0++; t3 = *s1++; \
458 OP(dst[2], (A*t0 + B*t1 + C*t2 + D*t3));\
459 t0 = *s0++; t2 = *s1++; \
460 OP(dst[3], (A*t1 + B*t0 + C*t3 + D*t2));\
461 dst+= stride;\
462 src+= stride;\
463 }while(--h);\
464 }\
465 \
466 static void OPNAME ## h264_chroma_mc8_sh4(uint8_t *dst/*align 8*/, uint8_t *src/*align 1*/, int stride, int h, int x, int y){\
467 const int A=(8-x)*(8-y);\
468 const int B=( x)*(8-y);\
469 const int C=(8-x)*( y);\
470 const int D=( x)*( y);\
471 \
472 assert(x<8 && y<8 && x>=0 && y>=0);\
473 \
474 do {\
475 int t0,t1,t2,t3; \
476 uint8_t *s0 = src; \
477 uint8_t *s1 = src+stride; \
478 t0 = *s0++; t2 = *s1++; \
479 t1 = *s0++; t3 = *s1++; \
480 OP(dst[0], (A*t0 + B*t1 + C*t2 + D*t3));\
481 t0 = *s0++; t2 = *s1++; \
482 OP(dst[1], (A*t1 + B*t0 + C*t3 + D*t2));\
483 t1 = *s0++; t3 = *s1++; \
484 OP(dst[2], (A*t0 + B*t1 + C*t2 + D*t3));\
485 t0 = *s0++; t2 = *s1++; \
486 OP(dst[3], (A*t1 + B*t0 + C*t3 + D*t2));\
487 t1 = *s0++; t3 = *s1++; \
488 OP(dst[4], (A*t0 + B*t1 + C*t2 + D*t3));\
489 t0 = *s0++; t2 = *s1++; \
490 OP(dst[5], (A*t1 + B*t0 + C*t3 + D*t2));\
491 t1 = *s0++; t3 = *s1++; \
492 OP(dst[6], (A*t0 + B*t1 + C*t2 + D*t3));\
493 t0 = *s0++; t2 = *s1++; \
494 OP(dst[7], (A*t1 + B*t0 + C*t3 + D*t2));\
495 dst+= stride;\
496 src+= stride;\
497 }while(--h);\
498 }
499
500 #define op_avg(a, b) a = (((a)+(((b) + 32)>>6)+1)>>1)
501 #define op_put(a, b) a = (((b) + 32)>>6)
502
503 H264_CHROMA_MC(put_ , op_put)
504 H264_CHROMA_MC(avg_ , op_avg)
505 #undef op_avg
506 #undef op_put
507
508 #define QPEL_MC(r, OPNAME, RND, OP) \
509 static void OPNAME ## mpeg4_qpel8_h_lowpass(uint8_t *dst, uint8_t *src, int dstStride, int srcStride, int h){\
510 uint8_t *cm = ff_cropTbl + MAX_NEG_CROP;\
511 do {\
512 uint8_t *s = src; \
513 int src0,src1,src2,src3,src4,src5,src6,src7,src8;\
514 src0= *s++;\
515 src1= *s++;\
516 src2= *s++;\
517 src3= *s++;\
518 src4= *s++;\
519 OP(dst[0], (src0+src1)*20 - (src0+src2)*6 + (src1+src3)*3 - (src2+src4));\
520 src5= *s++;\
521 OP(dst[1], (src1+src2)*20 - (src0+src3)*6 + (src0+src4)*3 - (src1+src5));\
522 src6= *s++;\
523 OP(dst[2], (src2+src3)*20 - (src1+src4)*6 + (src0+src5)*3 - (src0+src6));\
524 src7= *s++;\
525 OP(dst[3], (src3+src4)*20 - (src2+src5)*6 + (src1+src6)*3 - (src0+src7));\
526 src8= *s++;\
527 OP(dst[4], (src4+src5)*20 - (src3+src6)*6 + (src2+src7)*3 - (src1+src8));\
528 OP(dst[5], (src5+src6)*20 - (src4+src7)*6 + (src3+src8)*3 - (src2+src8));\
529 OP(dst[6], (src6+src7)*20 - (src5+src8)*6 + (src4+src8)*3 - (src3+src7));\
530 OP(dst[7], (src7+src8)*20 - (src6+src8)*6 + (src5+src7)*3 - (src4+src6));\
531 dst+=dstStride;\
532 src+=srcStride;\
533 }while(--h);\
534 }\
535 \
536 static void OPNAME ## mpeg4_qpel8_v_lowpass(uint8_t *dst, uint8_t *src, int dstStride, int srcStride){\
537 uint8_t *cm = ff_cropTbl + MAX_NEG_CROP;\
538 int w=8;\
539 do{\
540 uint8_t *s = src, *d=dst;\
541 int src0,src1,src2,src3,src4,src5,src6,src7,src8;\
542 src0 = *s; s+=srcStride; \
543 src1 = *s; s+=srcStride; \
544 src2 = *s; s+=srcStride; \
545 src3 = *s; s+=srcStride; \
546 src4 = *s; s+=srcStride; \
547 OP(*d, (src0+src1)*20 - (src0+src2)*6 + (src1+src3)*3 - (src2+src4));d+=dstStride;\
548 src5 = *s; s+=srcStride; \
549 OP(*d, (src1+src2)*20 - (src0+src3)*6 + (src0+src4)*3 - (src1+src5));d+=dstStride;\
550 src6 = *s; s+=srcStride; \
551 OP(*d, (src2+src3)*20 - (src1+src4)*6 + (src0+src5)*3 - (src0+src6));d+=dstStride;\
552 src7 = *s; s+=srcStride; \
553 OP(*d, (src3+src4)*20 - (src2+src5)*6 + (src1+src6)*3 - (src0+src7));d+=dstStride;\
554 src8 = *s; \
555 OP(*d, (src4+src5)*20 - (src3+src6)*6 + (src2+src7)*3 - (src1+src8));d+=dstStride;\
556 OP(*d, (src5+src6)*20 - (src4+src7)*6 + (src3+src8)*3 - (src2+src8));d+=dstStride;\
557 OP(*d, (src6+src7)*20 - (src5+src8)*6 + (src4+src8)*3 - (src3+src7));d+=dstStride;\
558 OP(*d, (src7+src8)*20 - (src6+src8)*6 + (src5+src7)*3 - (src4+src6));\
559 dst++;\
560 src++;\
561 }while(--w);\
562 }\
563 \
564 static void OPNAME ## mpeg4_qpel16_h_lowpass(uint8_t *dst, uint8_t *src, int dstStride, int srcStride, int h){\
565 uint8_t *cm = ff_cropTbl + MAX_NEG_CROP;\
566 do {\
567 uint8_t *s = src;\
568 int src0,src1,src2,src3,src4,src5,src6,src7,src8;\
569 int src9,src10,src11,src12,src13,src14,src15,src16;\
570 src0= *s++;\
571 src1= *s++;\
572 src2= *s++;\
573 src3= *s++;\
574 src4= *s++;\
575 OP(dst[ 0], (src0 +src1 )*20 - (src0 +src2 )*6 + (src1 +src3 )*3 - (src2 +src4 ));\
576 src5= *s++;\
577 OP(dst[ 1], (src1 +src2 )*20 - (src0 +src3 )*6 + (src0 +src4 )*3 - (src1 +src5 ));\
578 src6= *s++;\
579 OP(dst[ 2], (src2 +src3 )*20 - (src1 +src4 )*6 + (src0 +src5 )*3 - (src0 +src6 ));\
580 src7= *s++;\
581 OP(dst[ 3], (src3 +src4 )*20 - (src2 +src5 )*6 + (src1 +src6 )*3 - (src0 +src7 ));\
582 src8= *s++;\
583 OP(dst[ 4], (src4 +src5 )*20 - (src3 +src6 )*6 + (src2 +src7 )*3 - (src1 +src8 ));\
584 src9= *s++;\
585 OP(dst[ 5], (src5 +src6 )*20 - (src4 +src7 )*6 + (src3 +src8 )*3 - (src2 +src9 ));\
586 src10= *s++;\
587 OP(dst[ 6], (src6 +src7 )*20 - (src5 +src8 )*6 + (src4 +src9 )*3 - (src3 +src10));\
588 src11= *s++;\
589 OP(dst[ 7], (src7 +src8 )*20 - (src6 +src9 )*6 + (src5 +src10)*3 - (src4 +src11));\
590 src12= *s++;\
591 OP(dst[ 8], (src8 +src9 )*20 - (src7 +src10)*6 + (src6 +src11)*3 - (src5 +src12));\
592 src13= *s++;\
593 OP(dst[ 9], (src9 +src10)*20 - (src8 +src11)*6 + (src7 +src12)*3 - (src6 +src13));\
594 src14= *s++;\
595 OP(dst[10], (src10+src11)*20 - (src9 +src12)*6 + (src8 +src13)*3 - (src7 +src14));\
596 src15= *s++;\
597 OP(dst[11], (src11+src12)*20 - (src10+src13)*6 + (src9 +src14)*3 - (src8 +src15));\
598 src16= *s++;\
599 OP(dst[12], (src12+src13)*20 - (src11+src14)*6 + (src10+src15)*3 - (src9 +src16));\
600 OP(dst[13], (src13+src14)*20 - (src12+src15)*6 + (src11+src16)*3 - (src10+src16));\
601 OP(dst[14], (src14+src15)*20 - (src13+src16)*6 + (src12+src16)*3 - (src11+src15));\
602 OP(dst[15], (src15+src16)*20 - (src14+src16)*6 + (src13+src15)*3 - (src12+src14));\
603 dst+=dstStride;\
604 src+=srcStride;\
605 }while(--h);\
606 }\
607 \
608 static void OPNAME ## mpeg4_qpel16_v_lowpass(uint8_t *dst, uint8_t *src, int dstStride, int srcStride){\
609 uint8_t *cm = ff_cropTbl + MAX_NEG_CROP;\
610 int w=16;\
611 do {\
612 uint8_t *s = src, *d=dst;\
613 int src0,src1,src2,src3,src4,src5,src6,src7,src8;\
614 int src9,src10,src11,src12,src13,src14,src15,src16;\
615 src0 = *s; s+=srcStride; \
616 src1 = *s; s+=srcStride; \
617 src2 = *s; s+=srcStride; \
618 src3 = *s; s+=srcStride; \
619 src4 = *s; s+=srcStride; \
620 OP(*d, (src0 +src1 )*20 - (src0 +src2 )*6 + (src1 +src3 )*3 - (src2 +src4 ));d+=dstStride;\
621 src5 = *s; s+=srcStride; \
622 OP(*d, (src1 +src2 )*20 - (src0 +src3 )*6 + (src0 +src4 )*3 - (src1 +src5 ));d+=dstStride;\
623 src6 = *s; s+=srcStride; \
624 OP(*d, (src2 +src3 )*20 - (src1 +src4 )*6 + (src0 +src5 )*3 - (src0 +src6 ));d+=dstStride;\
625 src7 = *s; s+=srcStride; \
626 OP(*d, (src3 +src4 )*20 - (src2 +src5 )*6 + (src1 +src6 )*3 - (src0 +src7 ));d+=dstStride;\
627 src8 = *s; s+=srcStride; \
628 OP(*d, (src4 +src5 )*20 - (src3 +src6 )*6 + (src2 +src7 )*3 - (src1 +src8 ));d+=dstStride;\
629 src9 = *s; s+=srcStride; \
630 OP(*d, (src5 +src6 )*20 - (src4 +src7 )*6 + (src3 +src8 )*3 - (src2 +src9 ));d+=dstStride;\
631 src10 = *s; s+=srcStride; \
632 OP(*d, (src6 +src7 )*20 - (src5 +src8 )*6 + (src4 +src9 )*3 - (src3 +src10));d+=dstStride;\
633 src11 = *s; s+=srcStride; \
634 OP(*d, (src7 +src8 )*20 - (src6 +src9 )*6 + (src5 +src10)*3 - (src4 +src11));d+=dstStride;\
635 src12 = *s; s+=srcStride; \
636 OP(*d, (src8 +src9 )*20 - (src7 +src10)*6 + (src6 +src11)*3 - (src5 +src12));d+=dstStride;\
637 src13 = *s; s+=srcStride; \
638 OP(*d, (src9 +src10)*20 - (src8 +src11)*6 + (src7 +src12)*3 - (src6 +src13));d+=dstStride;\
639 src14 = *s; s+=srcStride; \
640 OP(*d, (src10+src11)*20 - (src9 +src12)*6 + (src8 +src13)*3 - (src7 +src14));d+=dstStride;\
641 src15 = *s; s+=srcStride; \
642 OP(*d, (src11+src12)*20 - (src10+src13)*6 + (src9 +src14)*3 - (src8 +src15));d+=dstStride;\
643 src16 = *s; \
644 OP(*d, (src12+src13)*20 - (src11+src14)*6 + (src10+src15)*3 - (src9 +src16));d+=dstStride;\
645 OP(*d, (src13+src14)*20 - (src12+src15)*6 + (src11+src16)*3 - (src10+src16));d+=dstStride;\
646 OP(*d, (src14+src15)*20 - (src13+src16)*6 + (src12+src16)*3 - (src11+src15));d+=dstStride;\
647 OP(*d, (src15+src16)*20 - (src14+src16)*6 + (src13+src15)*3 - (src12+src14));\
648 dst++;\
649 src++;\
650 }while(--w);\
651 }\
652 \
653 static void OPNAME ## qpel8_mc00_sh4 (uint8_t *dst, uint8_t *src, int stride){\
654 OPNAME ## pixels8_c(dst, src, stride, 8);\
655 }\
656 \
657 static void OPNAME ## qpel8_mc10_sh4(uint8_t *dst, uint8_t *src, int stride){\
658 uint8_t half[64];\
659 put ## RND ## mpeg4_qpel8_h_lowpass(half, src, 8, stride, 8);\
660 OPNAME ## pixels8_l2_aligned2(dst, src, half, stride, stride, 8, 8);\
661 }\
662 \
663 static void OPNAME ## qpel8_mc20_sh4(uint8_t *dst, uint8_t *src, int stride){\
664 OPNAME ## mpeg4_qpel8_h_lowpass(dst, src, stride, stride, 8);\
665 }\
666 \
667 static void OPNAME ## qpel8_mc30_sh4(uint8_t *dst, uint8_t *src, int stride){\
668 uint8_t half[64];\
669 put ## RND ## mpeg4_qpel8_h_lowpass(half, src, 8, stride, 8);\
670 OPNAME ## pixels8_l2_aligned2(dst, src+1, half, stride, stride, 8, 8);\
671 }\
672 \
673 static void OPNAME ## qpel8_mc01_sh4(uint8_t *dst, uint8_t *src, int stride){\
674 uint8_t full[16*9];\
675 uint8_t half[64];\
676 copy_block9(full, src, 16, stride, 9);\
677 put ## RND ## mpeg4_qpel8_v_lowpass(half, full, 8, 16);\
678 OPNAME ## pixels8_l2_aligned(dst, full, half, stride, 16, 8, 8);\
679 }\
680 \
681 static void OPNAME ## qpel8_mc02_sh4(uint8_t *dst, uint8_t *src, int stride){\
682 uint8_t full[16*9];\
683 copy_block9(full, src, 16, stride, 9);\
684 OPNAME ## mpeg4_qpel8_v_lowpass(dst, full, stride, 16);\
685 }\
686 \
687 static void OPNAME ## qpel8_mc03_sh4(uint8_t *dst, uint8_t *src, int stride){\
688 uint8_t full[16*9];\
689 uint8_t half[64];\
690 copy_block9(full, src, 16, stride, 9);\
691 put ## RND ## mpeg4_qpel8_v_lowpass(half, full, 8, 16);\
692 OPNAME ## pixels8_l2_aligned(dst, full+16, half, stride, 16, 8, 8);\
693 }\
694 static void OPNAME ## qpel8_mc11_sh4(uint8_t *dst, uint8_t *src, int stride){\
695 uint8_t full[16*9];\
696 uint8_t halfH[72];\
697 uint8_t halfHV[64];\
698 copy_block9(full, src, 16, stride, 9);\
699 put ## RND ## mpeg4_qpel8_h_lowpass(halfH, full, 8, 16, 9);\
700 put ## RND ## pixels8_l2_aligned(halfH, halfH, full, 8, 8, 16, 9);\
701 put ## RND ## mpeg4_qpel8_v_lowpass(halfHV, halfH, 8, 8);\
702 OPNAME ## pixels8_l2_aligned(dst, halfH, halfHV, stride, 8, 8, 8);\
703 }\
704 static void OPNAME ## qpel8_mc31_sh4(uint8_t *dst, uint8_t *src, int stride){\
705 uint8_t full[16*9];\
706 uint8_t halfH[72];\
707 uint8_t halfHV[64];\
708 copy_block9(full, src, 16, stride, 9);\
709 put ## RND ## mpeg4_qpel8_h_lowpass(halfH, full, 8, 16, 9);\
710 put ## RND ## pixels8_l2_aligned1(halfH, halfH, full+1, 8, 8, 16, 9);\
711 put ## RND ## mpeg4_qpel8_v_lowpass(halfHV, halfH, 8, 8);\
712 OPNAME ## pixels8_l2_aligned(dst, halfH, halfHV, stride, 8, 8, 8);\
713 }\
714 static void OPNAME ## qpel8_mc13_sh4(uint8_t *dst, uint8_t *src, int stride){\
715 uint8_t full[16*9];\
716 uint8_t halfH[72];\
717 uint8_t halfHV[64];\
718 copy_block9(full, src, 16, stride, 9);\
719 put ## RND ## mpeg4_qpel8_h_lowpass(halfH, full, 8, 16, 9);\
720 put ## RND ## pixels8_l2_aligned(halfH, halfH, full, 8, 8, 16, 9);\
721 put ## RND ## mpeg4_qpel8_v_lowpass(halfHV, halfH, 8, 8);\
722 OPNAME ## pixels8_l2_aligned(dst, halfH+8, halfHV, stride, 8, 8, 8);\
723 }\
724 static void OPNAME ## qpel8_mc33_sh4(uint8_t *dst, uint8_t *src, int stride){\
725 uint8_t full[16*9];\
726 uint8_t halfH[72];\
727 uint8_t halfHV[64];\
728 copy_block9(full, src, 16, stride, 9);\
729 put ## RND ## mpeg4_qpel8_h_lowpass(halfH, full, 8, 16, 9);\
730 put ## RND ## pixels8_l2_aligned1(halfH, halfH, full+1, 8, 8, 16, 9);\
731 put ## RND ## mpeg4_qpel8_v_lowpass(halfHV, halfH, 8, 8);\
732 OPNAME ## pixels8_l2_aligned(dst, halfH+8, halfHV, stride, 8, 8, 8);\
733 }\
734 static void OPNAME ## qpel8_mc21_sh4(uint8_t *dst, uint8_t *src, int stride){\
735 uint8_t halfH[72];\
736 uint8_t halfHV[64];\
737 put ## RND ## mpeg4_qpel8_h_lowpass(halfH, src, 8, stride, 9);\
738 put ## RND ## mpeg4_qpel8_v_lowpass(halfHV, halfH, 8, 8);\
739 OPNAME ## pixels8_l2_aligned(dst, halfH, halfHV, stride, 8, 8, 8);\
740 }\
741 static void OPNAME ## qpel8_mc23_sh4(uint8_t *dst, uint8_t *src, int stride){\
742 uint8_t halfH[72];\
743 uint8_t halfHV[64];\
744 put ## RND ## mpeg4_qpel8_h_lowpass(halfH, src, 8, stride, 9);\
745 put ## RND ## mpeg4_qpel8_v_lowpass(halfHV, halfH, 8, 8);\
746 OPNAME ## pixels8_l2_aligned(dst, halfH+8, halfHV, stride, 8, 8, 8);\
747 }\
748 static void OPNAME ## qpel8_mc12_sh4(uint8_t *dst, uint8_t *src, int stride){\
749 uint8_t full[16*9];\
750 uint8_t halfH[72];\
751 copy_block9(full, src, 16, stride, 9);\
752 put ## RND ## mpeg4_qpel8_h_lowpass(halfH, full, 8, 16, 9);\
753 put ## RND ## pixels8_l2_aligned(halfH, halfH, full, 8, 8, 16, 9);\
754 OPNAME ## mpeg4_qpel8_v_lowpass(dst, halfH, stride, 8);\
755 }\
756 static void OPNAME ## qpel8_mc32_sh4(uint8_t *dst, uint8_t *src, int stride){\
757 uint8_t full[16*9];\
758 uint8_t halfH[72];\
759 copy_block9(full, src, 16, stride, 9);\
760 put ## RND ## mpeg4_qpel8_h_lowpass(halfH, full, 8, 16, 9);\
761 put ## RND ## pixels8_l2_aligned1(halfH, halfH, full+1, 8, 8, 16, 9);\
762 OPNAME ## mpeg4_qpel8_v_lowpass(dst, halfH, stride, 8);\
763 }\
764 static void OPNAME ## qpel8_mc22_sh4(uint8_t *dst, uint8_t *src, int stride){\
765 uint8_t halfH[72];\
766 put ## RND ## mpeg4_qpel8_h_lowpass(halfH, src, 8, stride, 9);\
767 OPNAME ## mpeg4_qpel8_v_lowpass(dst, halfH, stride, 8);\
768 }\
769 static void OPNAME ## qpel16_mc00_sh4 (uint8_t *dst, uint8_t *src, int stride){\
770 OPNAME ## pixels16_c(dst, src, stride, 16);\
771 }\
772 \
773 static void OPNAME ## qpel16_mc10_sh4(uint8_t *dst, uint8_t *src, int stride){\
774 uint8_t half[256];\
775 put ## RND ## mpeg4_qpel16_h_lowpass(half, src, 16, stride, 16);\
776 OPNAME ## pixels16_l2_aligned2(dst, src, half, stride, stride, 16, 16);\
777 }\
778 \
779 static void OPNAME ## qpel16_mc20_sh4(uint8_t *dst, uint8_t *src, int stride){\
780 OPNAME ## mpeg4_qpel16_h_lowpass(dst, src, stride, stride, 16);\
781 }\
782 \
783 static void OPNAME ## qpel16_mc30_sh4(uint8_t *dst, uint8_t *src, int stride){\
784 uint8_t half[256];\
785 put ## RND ## mpeg4_qpel16_h_lowpass(half, src, 16, stride, 16);\
786 OPNAME ## pixels16_l2_aligned2(dst, src+1, half, stride, stride, 16, 16);\
787 }\
788 \
789 static void OPNAME ## qpel16_mc01_sh4(uint8_t *dst, uint8_t *src, int stride){\
790 uint8_t full[24*17];\
791 uint8_t half[256];\
792 copy_block17(full, src, 24, stride, 17);\
793 put ## RND ## mpeg4_qpel16_v_lowpass(half, full, 16, 24);\
794 OPNAME ## pixels16_l2_aligned(dst, full, half, stride, 24, 16, 16);\
795 }\
796 \
797 static void OPNAME ## qpel16_mc02_sh4(uint8_t *dst, uint8_t *src, int stride){\
798 uint8_t full[24*17];\
799 copy_block17(full, src, 24, stride, 17);\
800 OPNAME ## mpeg4_qpel16_v_lowpass(dst, full, stride, 24);\
801 }\
802 \
803 static void OPNAME ## qpel16_mc03_sh4(uint8_t *dst, uint8_t *src, int stride){\
804 uint8_t full[24*17];\
805 uint8_t half[256];\
806 copy_block17(full, src, 24, stride, 17);\
807 put ## RND ## mpeg4_qpel16_v_lowpass(half, full, 16, 24);\
808 OPNAME ## pixels16_l2_aligned(dst, full+24, half, stride, 24, 16, 16);\
809 }\
810 static void OPNAME ## qpel16_mc11_sh4(uint8_t *dst, uint8_t *src, int stride){\
811 uint8_t full[24*17];\
812 uint8_t halfH[272];\
813 uint8_t halfHV[256];\
814 copy_block17(full, src, 24, stride, 17);\
815 put ## RND ## mpeg4_qpel16_h_lowpass(halfH, full, 16, 24, 17);\
816 put ## RND ## pixels16_l2_aligned(halfH, halfH, full, 16, 16, 24, 17);\
817 put ## RND ## mpeg4_qpel16_v_lowpass(halfHV, halfH, 16, 16);\
818 OPNAME ## pixels16_l2_aligned(dst, halfH, halfHV, stride, 16, 16, 16);\
819 }\
820 static void OPNAME ## qpel16_mc31_sh4(uint8_t *dst, uint8_t *src, int stride){\
821 uint8_t full[24*17];\
822 uint8_t halfH[272];\
823 uint8_t halfHV[256];\
824 copy_block17(full, src, 24, stride, 17);\
825 put ## RND ## mpeg4_qpel16_h_lowpass(halfH, full, 16, 24, 17);\
826 put ## RND ## pixels16_l2_aligned1(halfH, halfH, full+1, 16, 16, 24, 17);\
827 put ## RND ## mpeg4_qpel16_v_lowpass(halfHV, halfH, 16, 16);\
828 OPNAME ## pixels16_l2_aligned(dst, halfH, halfHV, stride, 16, 16, 16);\
829 }\
830 static void OPNAME ## qpel16_mc13_sh4(uint8_t *dst, uint8_t *src, int stride){\
831 uint8_t full[24*17];\
832 uint8_t halfH[272];\
833 uint8_t halfHV[256];\
834 copy_block17(full, src, 24, stride, 17);\
835 put ## RND ## mpeg4_qpel16_h_lowpass(halfH, full, 16, 24, 17);\
836 put ## RND ## pixels16_l2_aligned(halfH, halfH, full, 16, 16, 24, 17);\
837 put ## RND ## mpeg4_qpel16_v_lowpass(halfHV, halfH, 16, 16);\
838 OPNAME ## pixels16_l2_aligned(dst, halfH+16, halfHV, stride, 16, 16, 16);\
839 }\
840 static void OPNAME ## qpel16_mc33_sh4(uint8_t *dst, uint8_t *src, int stride){\
841 uint8_t full[24*17];\
842 uint8_t halfH[272];\
843 uint8_t halfHV[256];\
844 copy_block17(full, src, 24, stride, 17);\
845 put ## RND ## mpeg4_qpel16_h_lowpass(halfH, full, 16, 24, 17);\
846 put ## RND ## pixels16_l2_aligned1(halfH, halfH, full+1, 16, 16, 24, 17);\
847 put ## RND ## mpeg4_qpel16_v_lowpass(halfHV, halfH, 16, 16);\
848 OPNAME ## pixels16_l2_aligned(dst, halfH+16, halfHV, stride, 16, 16, 16);\
849 }\
850 static void OPNAME ## qpel16_mc21_sh4(uint8_t *dst, uint8_t *src, int stride){\
851 uint8_t halfH[272];\
852 uint8_t halfHV[256];\
853 put ## RND ## mpeg4_qpel16_h_lowpass(halfH, src, 16, stride, 17);\
854 put ## RND ## mpeg4_qpel16_v_lowpass(halfHV, halfH, 16, 16);\
855 OPNAME ## pixels16_l2_aligned(dst, halfH, halfHV, stride, 16, 16, 16);\
856 }\
857 static void OPNAME ## qpel16_mc23_sh4(uint8_t *dst, uint8_t *src, int stride){\
858 uint8_t halfH[272];\
859 uint8_t halfHV[256];\
860 put ## RND ## mpeg4_qpel16_h_lowpass(halfH, src, 16, stride, 17);\
861 put ## RND ## mpeg4_qpel16_v_lowpass(halfHV, halfH, 16, 16);\
862 OPNAME ## pixels16_l2_aligned(dst, halfH+16, halfHV, stride, 16, 16, 16);\
863 }\
864 static void OPNAME ## qpel16_mc12_sh4(uint8_t *dst, uint8_t *src, int stride){\
865 uint8_t full[24*17];\
866 uint8_t halfH[272];\
867 copy_block17(full, src, 24, stride, 17);\
868 put ## RND ## mpeg4_qpel16_h_lowpass(halfH, full, 16, 24, 17);\
869 put ## RND ## pixels16_l2_aligned(halfH, halfH, full, 16, 16, 24, 17);\
870 OPNAME ## mpeg4_qpel16_v_lowpass(dst, halfH, stride, 16);\
871 }\
872 static void OPNAME ## qpel16_mc32_sh4(uint8_t *dst, uint8_t *src, int stride){\
873 uint8_t full[24*17];\
874 uint8_t halfH[272];\
875 copy_block17(full, src, 24, stride, 17);\
876 put ## RND ## mpeg4_qpel16_h_lowpass(halfH, full, 16, 24, 17);\
877 put ## RND ## pixels16_l2_aligned1(halfH, halfH, full+1, 16, 16, 24, 17);\
878 OPNAME ## mpeg4_qpel16_v_lowpass(dst, halfH, stride, 16);\
879 }\
880 static void OPNAME ## qpel16_mc22_sh4(uint8_t *dst, uint8_t *src, int stride){\
881 uint8_t halfH[272];\
882 put ## RND ## mpeg4_qpel16_h_lowpass(halfH, src, 16, stride, 17);\
883 OPNAME ## mpeg4_qpel16_v_lowpass(dst, halfH, stride, 16);\
884 }
885
886 #define op_avg(a, b) a = (((a)+cm[((b) + 16)>>5]+1)>>1)
887 #define op_avg_no_rnd(a, b) a = (((a)+cm[((b) + 15)>>5])>>1)
888 #define op_put(a, b) a = cm[((b) + 16)>>5]
889 #define op_put_no_rnd(a, b) a = cm[((b) + 15)>>5]
890
891 QPEL_MC(0, put_ , _ , op_put)
892 QPEL_MC(1, put_no_rnd_, _no_rnd_, op_put_no_rnd)
893 QPEL_MC(0, avg_ , _ , op_avg)
894 //QPEL_MC(1, avg_no_rnd , _ , op_avg)
895 #undef op_avg
896 #undef op_avg_no_rnd
897 #undef op_put
898 #undef op_put_no_rnd
899
900 #if 1
901 #define H264_LOWPASS(OPNAME, OP, OP2) \
902 static inline void OPNAME ## h264_qpel_h_lowpass(uint8_t *dst, uint8_t *src, int dstStride, int srcStride,int w,int h){\
903 uint8_t *cm = ff_cropTbl + MAX_NEG_CROP;\
904 do {\
905 int srcB,srcA,src0,src1,src2,src3,src4,src5,src6;\
906 uint8_t *s = src-2;\
907 srcB = *s++;\
908 srcA = *s++;\
909 src0 = *s++;\
910 src1 = *s++;\
911 src2 = *s++;\
912 src3 = *s++;\
913 OP(dst[0], (src0+src1)*20 - (srcA+src2)*5 + (srcB+src3));\
914 src4 = *s++;\
915 OP(dst[1], (src1+src2)*20 - (src0+src3)*5 + (srcA+src4));\
916 src5 = *s++;\
917 OP(dst[2], (src2+src3)*20 - (src1+src4)*5 + (src0+src5));\
918 src6 = *s++;\
919 OP(dst[3], (src3+src4)*20 - (src2+src5)*5 + (src1+src6));\
920 if (w>4) { /* it optimized */ \
921 int src7,src8,src9,src10; \
922 src7 = *s++;\
923 OP(dst[4], (src4+src5)*20 - (src3+src6)*5 + (src2+src7));\
924 src8 = *s++;\
925 OP(dst[5], (src5+src6)*20 - (src4+src7)*5 + (src3+src8));\
926 src9 = *s++;\
927 OP(dst[6], (src6+src7)*20 - (src5+src8)*5 + (src4+src9));\
928 src10 = *s++;\
929 OP(dst[7], (src7+src8)*20 - (src6+src9)*5 + (src5+src10));\
930 if (w>8) { \
931 int src11,src12,src13,src14,src15,src16,src17,src18; \
932 src11 = *s++;\
933 OP(dst[8] , (src8 +src9 )*20 - (src7 +src10)*5 + (src6 +src11));\
934 src12 = *s++;\
935 OP(dst[9] , (src9 +src10)*20 - (src8 +src11)*5 + (src7 +src12));\
936 src13 = *s++;\
937 OP(dst[10], (src10+src11)*20 - (src9 +src12)*5 + (src8 +src13));\
938 src14 = *s++;\
939 OP(dst[11], (src11+src12)*20 - (src10+src13)*5 + (src9 +src14));\
940 src15 = *s++;\
941 OP(dst[12], (src12+src13)*20 - (src11+src14)*5 + (src10+src15));\
942 src16 = *s++;\
943 OP(dst[13], (src13+src14)*20 - (src12+src15)*5 + (src11+src16));\
944 src17 = *s++;\
945 OP(dst[14], (src14+src15)*20 - (src13+src16)*5 + (src12+src17));\
946 src18 = *s++;\
947 OP(dst[15], (src15+src16)*20 - (src14+src17)*5 + (src13+src18));\
948 } \
949 } \
950 dst+=dstStride;\
951 src+=srcStride;\
952 }while(--h);\
953 }\
954 \
955 static inline void OPNAME ## h264_qpel_v_lowpass(uint8_t *dst, uint8_t *src, int dstStride, int srcStride,int w,int h){\
956 uint8_t *cm = ff_cropTbl + MAX_NEG_CROP;\
957 do{\
958 int srcB,srcA,src0,src1,src2,src3,src4,src5,src6;\
959 uint8_t *s = src-2*srcStride,*d=dst;\
960 srcB = *s; s+=srcStride;\
961 srcA = *s; s+=srcStride;\
962 src0 = *s; s+=srcStride;\
963 src1 = *s; s+=srcStride;\
964 src2 = *s; s+=srcStride;\
965 src3 = *s; s+=srcStride;\
966 OP(*d, (src0+src1)*20 - (srcA+src2)*5 + (srcB+src3));d+=dstStride;\
967 src4 = *s; s+=srcStride;\
968 OP(*d, (src1+src2)*20 - (src0+src3)*5 + (srcA+src4));d+=dstStride;\
969 src5 = *s; s+=srcStride;\
970 OP(*d, (src2+src3)*20 - (src1+src4)*5 + (src0+src5));d+=dstStride;\
971 src6 = *s; s+=srcStride;\
972 OP(*d, (src3+src4)*20 - (src2+src5)*5 + (src1+src6));d+=dstStride;\
973 if (h>4) { \
974 int src7,src8,src9,src10; \
975 src7 = *s; s+=srcStride;\
976 OP(*d, (src4+src5)*20 - (src3+src6)*5 + (src2+src7));d+=dstStride;\
977 src8 = *s; s+=srcStride;\
978 OP(*d, (src5+src6)*20 - (src4+src7)*5 + (src3+src8));d+=dstStride;\
979 src9 = *s; s+=srcStride;\
980 OP(*d, (src6+src7)*20 - (src5+src8)*5 + (src4+src9));d+=dstStride;\
981 src10 = *s; s+=srcStride;\
982 OP(*d, (src7+src8)*20 - (src6+src9)*5 + (src5+src10));d+=dstStride;\
983 if (h>8) { \
984 int src11,src12,src13,src14,src15,src16,src17,src18; \
985 src11 = *s; s+=srcStride;\
986 OP(*d , (src8 +src9 )*20 - (src7 +src10)*5 + (src6 +src11));d+=dstStride;\
987 src12 = *s; s+=srcStride;\
988 OP(*d , (src9 +src10)*20 - (src8 +src11)*5 + (src7 +src12));d+=dstStride;\
989 src13 = *s; s+=srcStride;\
990 OP(*d, (src10+src11)*20 - (src9 +src12)*5 + (src8 +src13));d+=dstStride;\
991 src14 = *s; s+=srcStride;\
992 OP(*d, (src11+src12)*20 - (src10+src13)*5 + (src9 +src14));d+=dstStride;\
993 src15 = *s; s+=srcStride;\
994 OP(*d, (src12+src13)*20 - (src11+src14)*5 + (src10+src15));d+=dstStride;\
995 src16 = *s; s+=srcStride;\
996 OP(*d, (src13+src14)*20 - (src12+src15)*5 + (src11+src16));d+=dstStride;\
997 src17 = *s; s+=srcStride;\
998 OP(*d, (src14+src15)*20 - (src13+src16)*5 + (src12+src17));d+=dstStride;\
999 src18 = *s; s+=srcStride;\
1000 OP(*d, (src15+src16)*20 - (src14+src17)*5 + (src13+src18));d+=dstStride;\
1001 } \
1002 } \
1003 dst++;\
1004 src++;\
1005 }while(--w);\
1006 }\
1007 \
1008 static inline void OPNAME ## h264_qpel_hv_lowpass(uint8_t *dst, int16_t *tmp, uint8_t *src, int dstStride, int tmpStride, int srcStride,int w,int h){\
1009 uint8_t *cm = ff_cropTbl + MAX_NEG_CROP;\
1010 int i;\
1011 src -= 2*srcStride;\
1012 i= h+5; \
1013 do {\
1014 int srcB,srcA,src0,src1,src2,src3,src4,src5,src6;\
1015 uint8_t *s = src-2;\
1016 srcB = *s++;\
1017 srcA = *s++;\
1018 src0 = *s++;\
1019 src1 = *s++;\
1020 src2 = *s++;\
1021 src3 = *s++;\
1022 tmp[0] = ((src0+src1)*20 - (srcA+src2)*5 + (srcB+src3));\
1023 src4 = *s++;\
1024 tmp[1] = ((src1+src2)*20 - (src0+src3)*5 + (srcA+src4));\
1025 src5 = *s++;\
1026 tmp[2] = ((src2+src3)*20 - (src1+src4)*5 + (src0+src5));\
1027 src6 = *s++;\
1028 tmp[3] = ((src3+src4)*20 - (src2+src5)*5 + (src1+src6));\
1029 if (w>4) { /* it optimized */ \
1030 int src7,src8,src9,src10; \
1031 src7 = *s++;\
1032 tmp[4] = ((src4+src5)*20 - (src3+src6)*5 + (src2+src7));\
1033 src8 = *s++;\
1034 tmp[5] = ((src5+src6)*20 - (src4+src7)*5 + (src3+src8));\
1035 src9 = *s++;\
1036 tmp[6] = ((src6+src7)*20 - (src5+src8)*5 + (src4+src9));\
1037 src10 = *s++;\
1038 tmp[7] = ((src7+src8)*20 - (src6+src9)*5 + (src5+src10));\
1039 if (w>8) { \
1040 int src11,src12,src13,src14,src15,src16,src17,src18; \
1041 src11 = *s++;\
1042 tmp[8] = ((src8 +src9 )*20 - (src7 +src10)*5 + (src6 +src11));\
1043 src12 = *s++;\
1044 tmp[9] = ((src9 +src10)*20 - (src8 +src11)*5 + (src7 +src12));\
1045 src13 = *s++;\
1046 tmp[10] = ((src10+src11)*20 - (src9 +src12)*5 + (src8 +src13));\
1047 src14 = *s++;\
1048 tmp[11] = ((src11+src12)*20 - (src10+src13)*5 + (src9 +src14));\
1049 src15 = *s++;\
1050 tmp[12] = ((src12+src13)*20 - (src11+src14)*5 + (src10+src15));\
1051 src16 = *s++;\
1052 tmp[13] = ((src13+src14)*20 - (src12+src15)*5 + (src11+src16));\
1053 src17 = *s++;\
1054 tmp[14] = ((src14+src15)*20 - (src13+src16)*5 + (src12+src17));\
1055 src18 = *s++;\
1056 tmp[15] = ((src15+src16)*20 - (src14+src17)*5 + (src13+src18));\
1057 } \
1058 } \
1059 tmp+=tmpStride;\
1060 src+=srcStride;\
1061 }while(--i);\
1062 tmp -= tmpStride*(h+5-2);\
1063 i = w; \
1064 do {\
1065 int tmpB,tmpA,tmp0,tmp1,tmp2,tmp3,tmp4,tmp5,tmp6;\
1066 int16_t *s = tmp-2*tmpStride; \
1067 uint8_t *d=dst;\
1068 tmpB = *s; s+=tmpStride;\
1069 tmpA = *s; s+=tmpStride;\
1070 tmp0 = *s; s+=tmpStride;\
1071 tmp1 = *s; s+=tmpStride;\
1072 tmp2 = *s; s+=tmpStride;\
1073 tmp3 = *s; s+=tmpStride;\
1074 OP2(*d, (tmp0+tmp1)*20 - (tmpA+tmp2)*5 + (tmpB+tmp3));d+=dstStride;\
1075 tmp4 = *s; s+=tmpStride;\
1076 OP2(*d, (tmp1+tmp2)*20 - (tmp0+tmp3)*5 + (tmpA+tmp4));d+=dstStride;\
1077 tmp5 = *s; s+=tmpStride;\
1078 OP2(*d, (tmp2+tmp3)*20 - (tmp1+tmp4)*5 + (tmp0+tmp5));d+=dstStride;\
1079 tmp6 = *s; s+=tmpStride;\
1080 OP2(*d, (tmp3+tmp4)*20 - (tmp2+tmp5)*5 + (tmp1+tmp6));d+=dstStride;\
1081 if (h>4) { \
1082 int tmp7,tmp8,tmp9,tmp10; \
1083 tmp7 = *s; s+=tmpStride;\
1084 OP2(*d, (tmp4+tmp5)*20 - (tmp3+tmp6)*5 + (tmp2+tmp7));d+=dstStride;\
1085 tmp8 = *s; s+=tmpStride;\
1086 OP2(*d, (tmp5+tmp6)*20 - (tmp4+tmp7)*5 + (tmp3+tmp8));d+=dstStride;\
1087 tmp9 = *s; s+=tmpStride;\
1088 OP2(*d, (tmp6+tmp7)*20 - (tmp5+tmp8)*5 + (tmp4+tmp9));d+=dstStride;\
1089 tmp10 = *s; s+=tmpStride;\
1090 OP2(*d, (tmp7+tmp8)*20 - (tmp6+tmp9)*5 + (tmp5+tmp10));d+=dstStride;\
1091 if (h>8) { \
1092 int tmp11,tmp12,tmp13,tmp14,tmp15,tmp16,tmp17,tmp18; \
1093 tmp11 = *s; s+=tmpStride;\
1094 OP2(*d , (tmp8 +tmp9 )*20 - (tmp7 +tmp10)*5 + (tmp6 +tmp11));d+=dstStride;\
1095 tmp12 = *s; s+=tmpStride;\
1096 OP2(*d , (tmp9 +tmp10)*20 - (tmp8 +tmp11)*5 + (tmp7 +tmp12));d+=dstStride;\
1097 tmp13 = *s; s+=tmpStride;\
1098 OP2(*d, (tmp10+tmp11)*20 - (tmp9 +tmp12)*5 + (tmp8 +tmp13));d+=dstStride;\
1099 tmp14 = *s; s+=tmpStride;\
1100 OP2(*d, (tmp11+tmp12)*20 - (tmp10+tmp13)*5 + (tmp9 +tmp14));d+=dstStride;\
1101 tmp15 = *s; s+=tmpStride;\
1102 OP2(*d, (tmp12+tmp13)*20 - (tmp11+tmp14)*5 + (tmp10+tmp15));d+=dstStride;\
1103 tmp16 = *s; s+=tmpStride;\
1104 OP2(*d, (tmp13+tmp14)*20 - (tmp12+tmp15)*5 + (tmp11+tmp16));d+=dstStride;\
1105 tmp17 = *s; s+=tmpStride;\
1106 OP2(*d, (tmp14+tmp15)*20 - (tmp13+tmp16)*5 + (tmp12+tmp17));d+=dstStride;\
1107 tmp18 = *s; s+=tmpStride;\
1108 OP2(*d, (tmp15+tmp16)*20 - (tmp14+tmp17)*5 + (tmp13+tmp18));d+=dstStride;\
1109 } \
1110 } \
1111 dst++;\
1112 tmp++;\
1113 }while(--i);\
1114 }\
1115 \
1116 static void OPNAME ## h264_qpel4_h_lowpass(uint8_t *dst, uint8_t *src, int dstStride, int srcStride){\
1117 OPNAME ## h264_qpel_h_lowpass(dst,src,dstStride,srcStride,4,4); \
1118 }\
1119 static void OPNAME ## h264_qpel8_h_lowpass(uint8_t *dst, uint8_t *src, int dstStride, int srcStride){\
1120 OPNAME ## h264_qpel_h_lowpass(dst,src,dstStride,srcStride,8,8); \
1121 }\
1122 static void OPNAME ## h264_qpel16_h_lowpass(uint8_t *dst, uint8_t *src, int dstStride, int srcStride){\
1123 OPNAME ## h264_qpel_h_lowpass(dst,src,dstStride,srcStride,16,16); \
1124 }\
1125 \
1126 static void OPNAME ## h264_qpel4_v_lowpass(uint8_t *dst, uint8_t *src, int dstStride, int srcStride){\
1127 OPNAME ## h264_qpel_v_lowpass(dst,src,dstStride,srcStride,4,4); \
1128 }\
1129 static void OPNAME ## h264_qpel8_v_lowpass(uint8_t *dst, uint8_t *src, int dstStride, int srcStride){\
1130 OPNAME ## h264_qpel_v_lowpass(dst,src,dstStride,srcStride,8,8); \
1131 }\
1132 static void OPNAME ## h264_qpel16_v_lowpass(uint8_t *dst, uint8_t *src, int dstStride, int srcStride){\
1133 OPNAME ## h264_qpel_v_lowpass(dst,src,dstStride,srcStride,16,16); \
1134 }\
1135 static void OPNAME ## h264_qpel4_hv_lowpass(uint8_t *dst, int16_t *tmp, uint8_t *src, int dstStride, int tmpStride, int srcStride){\
1136 OPNAME ## h264_qpel_hv_lowpass(dst,tmp,src,dstStride,tmpStride,srcStride,4,4); \
1137 }\
1138 static void OPNAME ## h264_qpel8_hv_lowpass(uint8_t *dst, int16_t *tmp, uint8_t *src, int dstStride, int tmpStride, int srcStride){\
1139 OPNAME ## h264_qpel_hv_lowpass(dst,tmp,src,dstStride,tmpStride,srcStride,8,8); \
1140 }\
1141 static void OPNAME ## h264_qpel16_hv_lowpass(uint8_t *dst, int16_t *tmp, uint8_t *src, int dstStride, int tmpStride, int srcStride){\
1142 OPNAME ## h264_qpel_hv_lowpass(dst,tmp,src,dstStride,tmpStride,srcStride,16,16); \
1143 }\
1144
1145 #define H264_MC(OPNAME, SIZE) \
1146 static void OPNAME ## h264_qpel ## SIZE ## _mc00_sh4 (uint8_t *dst, uint8_t *src, int stride){\
1147 OPNAME ## pixels ## SIZE ## _c(dst, src, stride, SIZE);\
1148 }\
1149 \
1150 static void OPNAME ## h264_qpel ## SIZE ## _mc10_sh4(uint8_t *dst, uint8_t *src, int stride){\
1151 uint8_t half[SIZE*SIZE];\
1152 put_h264_qpel ## SIZE ## _h_lowpass(half, src, SIZE, stride);\
1153 OPNAME ## pixels ## SIZE ## _l2_aligned2(dst, src, half, stride, stride, SIZE, SIZE);\
1154 }\
1155 \
1156 static void OPNAME ## h264_qpel ## SIZE ## _mc20_sh4(uint8_t *dst, uint8_t *src, int stride){\
1157 OPNAME ## h264_qpel ## SIZE ## _h_lowpass(dst, src, stride, stride);\
1158 }\
1159 \
1160 static void OPNAME ## h264_qpel ## SIZE ## _mc30_sh4(uint8_t *dst, uint8_t *src, int stride){\
1161 uint8_t half[SIZE*SIZE];\
1162 put_h264_qpel ## SIZE ## _h_lowpass(half, src, SIZE, stride);\
1163 OPNAME ## pixels ## SIZE ## _l2_aligned2(dst, src+1, half, stride, stride, SIZE, SIZE);\
1164 }\
1165 \
1166 static void OPNAME ## h264_qpel ## SIZE ## _mc01_sh4(uint8_t *dst, uint8_t *src, int stride){\
1167 uint8_t full[SIZE*(SIZE+5)];\
1168 uint8_t * const full_mid= full + SIZE*2;\
1169 uint8_t half[SIZE*SIZE];\
1170 copy_block ## SIZE (full, src - stride*2, SIZE, stride, SIZE + 5);\
1171 put_h264_qpel ## SIZE ## _v_lowpass(half, full_mid, SIZE, SIZE);\
1172 OPNAME ## pixels ## SIZE ## _l2_aligned(dst, full_mid, half, stride, SIZE, SIZE, SIZE);\
1173 }\
1174 \
1175 static void OPNAME ## h264_qpel ## SIZE ## _mc02_sh4(uint8_t *dst, uint8_t *src, int stride){\
1176 uint8_t full[SIZE*(SIZE+5)];\
1177 uint8_t * const full_mid= full + SIZE*2;\
1178 copy_block ## SIZE (full, src - stride*2, SIZE, stride, SIZE + 5);\
1179 OPNAME ## h264_qpel ## SIZE ## _v_lowpass(dst, full_mid, stride, SIZE);\
1180 }\
1181 \
1182 static void OPNAME ## h264_qpel ## SIZE ## _mc03_sh4(uint8_t *dst, uint8_t *src, int stride){\
1183 uint8_t full[SIZE*(SIZE+5)];\
1184 uint8_t * const full_mid= full + SIZE*2;\
1185 uint8_t half[SIZE*SIZE];\
1186 copy_block ## SIZE (full, src - stride*2, SIZE, stride, SIZE + 5);\
1187 put_h264_qpel ## SIZE ## _v_lowpass(half, full_mid, SIZE, SIZE);\
1188 OPNAME ## pixels ## SIZE ## _l2_aligned(dst, full_mid+SIZE, half, stride, SIZE, SIZE, SIZE);\
1189 }\
1190 \
1191 static void OPNAME ## h264_qpel ## SIZE ## _mc11_sh4(uint8_t *dst, uint8_t *src, int stride){\
1192 uint8_t full[SIZE*(SIZE+5)];\
1193 uint8_t * const full_mid= full + SIZE*2;\
1194 uint8_t halfH[SIZE*SIZE];\
1195 uint8_t halfV[SIZE*SIZE];\
1196 put_h264_qpel ## SIZE ## _h_lowpass(halfH, src, SIZE, stride);\
1197 copy_block ## SIZE (full, src - stride*2, SIZE, stride, SIZE + 5);\
1198 put_h264_qpel ## SIZE ## _v_lowpass(halfV, full_mid, SIZE, SIZE);\
1199 OPNAME ## pixels ## SIZE ## _l2_aligned(dst, halfH, halfV, stride, SIZE, SIZE, SIZE);\
1200 }\
1201 \
1202 static void OPNAME ## h264_qpel ## SIZE ## _mc31_sh4(uint8_t *dst, uint8_t *src, int stride){\
1203 uint8_t full[SIZE*(SIZE+5)];\
1204 uint8_t * const full_mid= full + SIZE*2;\
1205 uint8_t halfH[SIZE*SIZE];\
1206 uint8_t halfV[SIZE*SIZE];\
1207 put_h264_qpel ## SIZE ## _h_lowpass(halfH, src, SIZE, stride);\
1208 copy_block ## SIZE (full, src - stride*2 + 1, SIZE, stride, SIZE + 5);\
1209 put_h264_qpel ## SIZE ## _v_lowpass(halfV, full_mid, SIZE, SIZE);\
1210 OPNAME ## pixels ## SIZE ## _l2_aligned(dst, halfH, halfV, stride, SIZE, SIZE, SIZE);\
1211 }\
1212 \
1213 static void OPNAME ## h264_qpel ## SIZE ## _mc13_sh4(uint8_t *dst, uint8_t *src, int stride){\
1214 uint8_t full[SIZE*(SIZE+5)];\
1215 uint8_t * const full_mid= full + SIZE*2;\
1216 uint8_t halfH[SIZE*SIZE];\
1217 uint8_t halfV[SIZE*SIZE];\
1218 put_h264_qpel ## SIZE ## _h_lowpass(halfH, src + stride, SIZE, stride);\
1219 copy_block ## SIZE (full, src - stride*2, SIZE, stride, SIZE + 5);\
1220 put_h264_qpel ## SIZE ## _v_lowpass(halfV, full_mid, SIZE, SIZE);\
1221 OPNAME ## pixels ## SIZE ## _l2_aligned(dst, halfH, halfV, stride, SIZE, SIZE, SIZE);\
1222 }\
1223 \
1224 static void OPNAME ## h264_qpel ## SIZE ## _mc33_sh4(uint8_t *dst, uint8_t *src, int stride){\
1225 uint8_t full[SIZE*(SIZE+5)];\
1226 uint8_t * const full_mid= full + SIZE*2;\
1227 uint8_t halfH[SIZE*SIZE];\
1228 uint8_t halfV[SIZE*SIZE];\
1229 put_h264_qpel ## SIZE ## _h_lowpass(halfH, src + stride, SIZE, stride);\
1230 copy_block ## SIZE (full, src - stride*2 + 1, SIZE, stride, SIZE + 5);\
1231 put_h264_qpel ## SIZE ## _v_lowpass(halfV, full_mid, SIZE, SIZE);\
1232 OPNAME ## pixels ## SIZE ## _l2_aligned(dst, halfH, halfV, stride, SIZE, SIZE, SIZE);\
1233 }\
1234 \
1235 static void OPNAME ## h264_qpel ## SIZE ## _mc22_sh4(uint8_t *dst, uint8_t *src, int stride){\
1236 int16_t tmp[SIZE*(SIZE+5)];\
1237 OPNAME ## h264_qpel ## SIZE ## _hv_lowpass(dst, tmp, src, stride, SIZE, stride);\
1238 }\
1239 \
1240 static void OPNAME ## h264_qpel ## SIZE ## _mc21_sh4(uint8_t *dst, uint8_t *src, int stride){\
1241 int16_t tmp[SIZE*(SIZE+5)];\
1242 uint8_t halfH[SIZE*SIZE];\
1243 uint8_t halfHV[SIZE*SIZE];\
1244 put_h264_qpel ## SIZE ## _h_lowpass(halfH, src, SIZE, stride);\
1245 put_h264_qpel ## SIZE ## _hv_lowpass(halfHV, tmp, src, SIZE, SIZE, stride);\
1246 OPNAME ## pixels ## SIZE ## _l2_aligned(dst, halfH, halfHV, stride, SIZE, SIZE, SIZE);\
1247 }\
1248 \
1249 static void OPNAME ## h264_qpel ## SIZE ## _mc23_sh4(uint8_t *dst, uint8_t *src, int stride){\
1250 int16_t tmp[SIZE*(SIZE+5)];\
1251 uint8_t halfH[SIZE*SIZE];\
1252 uint8_t halfHV[SIZE*SIZE];\
1253 put_h264_qpel ## SIZE ## _h_lowpass(halfH, src + stride, SIZE, stride);\
1254 put_h264_qpel ## SIZE ## _hv_lowpass(halfHV, tmp, src, SIZE, SIZE, stride);\
1255 OPNAME ## pixels ## SIZE ## _l2_aligned(dst, halfH, halfHV, stride, SIZE, SIZE, SIZE);\
1256 }\
1257 \
1258 static void OPNAME ## h264_qpel ## SIZE ## _mc12_sh4(uint8_t *dst, uint8_t *src, int stride){\
1259 uint8_t full[SIZE*(SIZE+5)];\
1260 uint8_t * const full_mid= full + SIZE*2;\
1261 int16_t tmp[SIZE*(SIZE+5)];\
1262 uint8_t halfV[SIZE*SIZE];\
1263 uint8_t halfHV[SIZE*SIZE];\
1264 copy_block ## SIZE (full, src - stride*2, SIZE, stride, SIZE + 5);\
1265 put_h264_qpel ## SIZE ## _v_lowpass(halfV, full_mid, SIZE, SIZE);\
1266 put_h264_qpel ## SIZE ## _hv_lowpass(halfHV, tmp, src, SIZE, SIZE, stride);\
1267 OPNAME ## pixels ## SIZE ## _l2_aligned(dst, halfV, halfHV, stride, SIZE, SIZE, SIZE);\
1268 }\
1269 \
1270 static void OPNAME ## h264_qpel ## SIZE ## _mc32_sh4(uint8_t *dst, uint8_t *src, int stride){\
1271 uint8_t full[SIZE*(SIZE+5)];\
1272 uint8_t * const full_mid= full + SIZE*2;\
1273 int16_t tmp[SIZE*(SIZE+5)];\
1274 uint8_t halfV[SIZE*SIZE];\
1275 uint8_t halfHV[SIZE*SIZE];\
1276 copy_block ## SIZE (full, src - stride*2 + 1, SIZE, stride, SIZE + 5);\
1277 put_h264_qpel ## SIZE ## _v_lowpass(halfV, full_mid, SIZE, SIZE);\
1278 put_h264_qpel ## SIZE ## _hv_lowpass(halfHV, tmp, src, SIZE, SIZE, stride);\
1279 OPNAME ## pixels ## SIZE ## _l2_aligned(dst, halfV, halfHV, stride, SIZE, SIZE, SIZE);\
1280 }\
1281
1282 #define op_avg(a, b) a = (((a)+cm[((b) + 16)>>5]+1)>>1)
1283 //#define op_avg2(a, b) a = (((a)*w1+cm[((b) + 16)>>5]*w2 + o + 64)>>7)
1284 #define op_put(a, b) a = cm[((b) + 16)>>5]
1285 #define op2_avg(a, b) a = (((a)+cm[((b) + 512)>>10]+1)>>1)
1286 #define op2_put(a, b) a = cm[((b) + 512)>>10]
1287
1288 H264_LOWPASS(put_ , op_put, op2_put)
1289 H264_LOWPASS(avg_ , op_avg, op2_avg)
1290 H264_MC(put_, 4)
1291 H264_MC(put_, 8)
1292 H264_MC(put_, 16)
1293 H264_MC(avg_, 4)
1294 H264_MC(avg_, 8)
1295 H264_MC(avg_, 16)
1296
1297 #undef op_avg
1298 #undef op_put
1299 #undef op2_avg
1300 #undef op2_put
1301 #endif
1302
1303 static void wmv2_mspel8_h_lowpass(uint8_t *dst, uint8_t *src, int dstStride, int srcStride, int h){
1304 uint8_t *cm = ff_cropTbl + MAX_NEG_CROP;
1305
1306 do{
1307 int src_1,src0,src1,src2,src3,src4,src5,src6,src7,src8,src9;
1308 uint8_t *s = src;
1309 src_1 = s[-1];
1310 src0 = *s++;
1311 src1 = *s++;
1312 src2 = *s++;
1313 dst[0]= cm[(9*(src0 + src1) - (src_1 + src2) + 8)>>4];
1314 src3 = *s++;
1315 dst[1]= cm[(9*(src1 + src2) - (src0 + src3) + 8)>>4];
1316 src4 = *s++;
1317 dst[2]= cm[(9*(src2 + src3) - (src1 + src4) + 8)>>4];
1318 src5 = *s++;
1319 dst[3]= cm[(9*(src3 + src4) - (src2 + src5) + 8)>>4];
1320 src6 = *s++;
1321 dst[4]= cm[(9*(src4 + src5) - (src3 + src6) + 8)>>4];
1322 src7 = *s++;
1323 dst[5]= cm[(9*(src5 + src6) - (src4 + src7) + 8)>>4];
1324 src8 = *s++;
1325 dst[6]= cm[(9*(src6 + src7) - (src5 + src8) + 8)>>4];
1326 src9 = *s++;
1327 dst[7]= cm[(9*(src7 + src8) - (src6 + src9) + 8)>>4];
1328 dst+=dstStride;
1329 src+=srcStride;
1330 }while(--h);
1331 }
1332
1333 static void wmv2_mspel8_v_lowpass(uint8_t *dst, uint8_t *src, int dstStride, int srcStride, int w){
1334 uint8_t *cm = ff_cropTbl + MAX_NEG_CROP;
1335
1336 do{
1337 int src_1,src0,src1,src2,src3,src4,src5,src6,src7,src8,src9;
1338 uint8_t *s = src,*d = dst;
1339 src_1 = *(s-srcStride);
1340 src0 = *s; s+=srcStride;
1341 src1 = *s; s+=srcStride;
1342 src2 = *s; s+=srcStride;
1343 *d= cm[(9*(src0 + src1) - (src_1 + src2) + 8)>>4]; d+=dstStride;
1344 src3 = *s; s+=srcStride;
1345 *d= cm[(9*(src1 + src2) - (src0 + src3) + 8)>>4]; d+=dstStride;
1346 src4 = *s; s+=srcStride;
1347 *d= cm[(9*(src2 + src3) - (src1 + src4) + 8)>>4]; d+=dstStride;
1348 src5 = *s; s+=srcStride;
1349 *d= cm[(9*(src3 + src4) - (src2 + src5) + 8)>>4]; d+=dstStride;
1350 src6 = *s; s+=srcStride;
1351 *d= cm[(9*(src4 + src5) - (src3 + src6) + 8)>>4]; d+=dstStride;
1352 src7 = *s; s+=srcStride;
1353 *d= cm[(9*(src5 + src6) - (src4 + src7) + 8)>>4]; d+=dstStride;
1354 src8 = *s; s+=srcStride;
1355 *d= cm[(9*(src6 + src7) - (src5 + src8) + 8)>>4]; d+=dstStride;
1356 src9 = *s;
1357 *d= cm[(9*(src7 + src8) - (src6 + src9) + 8)>>4]; d+=dstStride;
1358 src++;
1359 dst++;
1360 }while(--w);
1361 }
1362
1363 static void put_mspel8_mc00_sh4 (uint8_t *dst, uint8_t *src, int stride){
1364 put_pixels8_c(dst, src, stride, 8);
1365 }
1366
1367 static void put_mspel8_mc10_sh4(uint8_t *dst, uint8_t *src, int stride){
1368 uint8_t half[64];
1369 wmv2_mspel8_h_lowpass(half, src, 8, stride, 8);
1370 put_pixels8_l2_aligned2(dst, src, half, stride, stride, 8, 8);
1371 }
1372
1373 static void put_mspel8_mc20_sh4(uint8_t *dst, uint8_t *src, int stride){
1374 wmv2_mspel8_h_lowpass(dst, src, stride, stride, 8);
1375 }
1376
1377 static void put_mspel8_mc30_sh4(uint8_t *dst, uint8_t *src, int stride){
1378 uint8_t half[64];
1379 wmv2_mspel8_h_lowpass(half, src, 8, stride, 8);
1380 put_pixels8_l2_aligned2(dst, src+1, half, stride, stride, 8, 8);
1381 }
1382
1383 static void put_mspel8_mc02_sh4(uint8_t *dst, uint8_t *src, int stride){
1384 wmv2_mspel8_v_lowpass(dst, src, stride, stride, 8);
1385 }
1386
1387 static void put_mspel8_mc12_sh4(uint8_t *dst, uint8_t *src, int stride){
1388 uint8_t halfH[88];
1389 uint8_t halfV[64];
1390 uint8_t halfHV[64];
1391 wmv2_mspel8_h_lowpass(halfH, src-stride, 8, stride, 11);
1392 wmv2_mspel8_v_lowpass(halfV, src, 8, stride, 8);
1393 wmv2_mspel8_v_lowpass(halfHV, halfH+8, 8, 8, 8);
1394 put_pixels8_l2_aligned(dst, halfV, halfHV, stride, 8, 8, 8);
1395 }
1396 static void put_mspel8_mc32_sh4(uint8_t *dst, uint8_t *src, int stride){
1397 uint8_t halfH[88];
1398 uint8_t halfV[64];
1399 uint8_t halfHV[64];
1400 wmv2_mspel8_h_lowpass(halfH, src-stride, 8, stride, 11);
1401 wmv2_mspel8_v_lowpass(halfV, src+1, 8, stride, 8);
1402 wmv2_mspel8_v_lowpass(halfHV, halfH+8, 8, 8, 8);
1403 put_pixels8_l2_aligned(dst, halfV, halfHV, stride, 8, 8, 8);
1404 }
1405 static void put_mspel8_mc22_sh4(uint8_t *dst, uint8_t *src, int stride){
1406 uint8_t halfH[88];
1407 wmv2_mspel8_h_lowpass(halfH, src-stride, 8, stride, 11);
1408 wmv2_mspel8_v_lowpass(dst, halfH+8, stride, 8, 8);
1409 }
ffmpeg with support for matroska ordered chapters