comparison x86/vp6dsp_mmx.c @ 8817:492f8911992c libavcodec

add MMX version of vp6_filter_diag original patch by Sebastien Lucas sebastien.lucas _at_ gmail _dot_ com
author aurel
date Thu, 12 Feb 2009 23:52:52 +0000
parents
children 7dd2a45249a9
comparison
equal deleted inserted replaced
8816:53f9f3994ec8 8817:492f8911992c
1 /**
2 * @file libavcodec/x86/vp6dsp_mmx.c
3 * MMX-optimized functions for the VP6 decoder
4 *
5 * Copyright (C) 2009 Sebastien Lucas <sebastien.lucas@gmail.com>
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 #include "libavutil/x86_cpu.h"
25 #include "libavcodec/dsputil.h"
26 #include "dsputil_mmx.h"
27 #include "vp6dsp_mmx.h"
28
29
30 #define DIAG4_MMX(in1,in2,in3,in4) \
31 "movq "#in1"(%0), %%mm0 \n\t" \
32 "movq "#in2"(%0), %%mm1 \n\t" \
33 "movq %%mm0, %%mm3 \n\t" \
34 "movq %%mm1, %%mm4 \n\t" \
35 "punpcklbw %%mm7, %%mm0 \n\t" \
36 "punpcklbw %%mm7, %%mm1 \n\t" \
37 "punpckhbw %%mm7, %%mm3 \n\t" \
38 "punpckhbw %%mm7, %%mm4 \n\t" \
39 "pmullw 0(%2), %%mm0 \n\t" /* src[x-8 ] * biweight [0] */ \
40 "pmullw 8(%2), %%mm1 \n\t" /* src[x ] * biweight [1] */ \
41 "pmullw 0(%2), %%mm3 \n\t" /* src[x-8 ] * biweight [0] */ \
42 "pmullw 8(%2), %%mm4 \n\t" /* src[x ] * biweight [1] */ \
43 "paddw %%mm1, %%mm0 \n\t" \
44 "paddw %%mm4, %%mm3 \n\t" \
45 "movq "#in3"(%0), %%mm1 \n\t" \
46 "movq "#in4"(%0), %%mm2 \n\t" \
47 "movq %%mm1, %%mm4 \n\t" \
48 "movq %%mm2, %%mm5 \n\t" \
49 "punpcklbw %%mm7, %%mm1 \n\t" \
50 "punpcklbw %%mm7, %%mm2 \n\t" \
51 "punpckhbw %%mm7, %%mm4 \n\t" \
52 "punpckhbw %%mm7, %%mm5 \n\t" \
53 "pmullw 16(%2), %%mm1 \n\t" /* src[x+8 ] * biweight [2] */ \
54 "pmullw 24(%2), %%mm2 \n\t" /* src[x+16] * biweight [3] */ \
55 "pmullw 16(%2), %%mm4 \n\t" /* src[x+8 ] * biweight [2] */ \
56 "pmullw 24(%2), %%mm5 \n\t" /* src[x+16] * biweight [3] */ \
57 "paddw %%mm2, %%mm1 \n\t" \
58 "paddw %%mm5, %%mm4 \n\t" \
59 "paddsw %%mm1, %%mm0 \n\t" \
60 "paddsw %%mm4, %%mm3 \n\t" \
61 "paddsw %%mm6, %%mm0 \n\t" /* Add 64 */ \
62 "paddsw %%mm6, %%mm3 \n\t" /* Add 64 */ \
63 "psraw $7, %%mm0 \n\t" \
64 "psraw $7, %%mm3 \n\t" \
65 "packuswb %%mm3, %%mm0 \n\t" \
66 "movq %%mm0, (%1) \n\t"
67
68 void ff_vp6_filter_diag4_mmx(uint8_t *dst, uint8_t *src, int stride,
69 const int16_t *h_weights, const int16_t *v_weights)
70 {
71 uint8_t tmp[8*11], *t = tmp;
72 int16_t weights[4*4];
73 int i;
74 src -= stride;
75
76 for (i=0; i<4*4; i++)
77 weights[i] = h_weights[i>>2];
78
79 __asm__ volatile(
80 "pxor %%mm7, %%mm7 \n\t"
81 "movq "MANGLE(ff_pw_64)", %%mm6 \n\t"
82 "1: \n\t"
83 DIAG4_MMX(-1,0,1,2)
84 "add $8, %1 \n\t"
85 "add %3, %0 \n\t"
86 "decl %4 \n\t"
87 "jnz 1b \n\t"
88 : "+r"(src), "+r"(t)
89 : "r"(weights), "r"((x86_reg)stride), "r"(11)
90 : "memory");
91
92 t = tmp + 8;
93 for (i=0; i<4*4; i++)
94 weights[i] = v_weights[i>>2];
95
96 __asm__ volatile(
97 "pxor %%mm7, %%mm7 \n\t"
98 "movq "MANGLE(ff_pw_64)", %%mm6 \n\t"
99 "1: \n\t"
100 DIAG4_MMX(-8,0,8,16)
101 "add $8, %0 \n\t"
102 "add %3, %1 \n\t"
103 "decl %4 \n\t"
104 "jnz 1b \n\t"
105 : "+r"(t), "+r"(dst)
106 : "r"(weights), "r"((x86_reg)stride), "r"(8)
107 : "memory");
108 }