annotate libmpcodecs/vf_yadif.c @ 20911:5ab1e8c8f749

bump date of NUT spec for previous commit
author ods15
date Tue, 14 Nov 2006 17:05:20 +0000
parents 5ae0a24fbed8
children b37a53da3b8a
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
18608
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
1 /*
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
2 Copyright (C) 2006 Michael Niedermayer <michaelni@gmx.at>
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
3
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
4 This program is free software; you can redistribute it and/or modify
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
5 it under the terms of the GNU General Public License as published by
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
6 the Free Software Foundation; either version 2 of the License, or
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
7 (at your option) any later version.
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
8
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
9 This program is distributed in the hope that it will be useful,
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
10 but WITHOUT ANY WARRANTY; without even the implied warranty of
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
11 MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
12 GNU General Public License for more details.
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
13
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
14 You should have received a copy of the GNU General Public License
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
15 along with this program; if not, write to the Free Software
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
16 Foundation, Inc., 51 Franklin St, Fifth Floor, Boston, MA 02110-1301 USA
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
17 */
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
18
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
19 #include <stdio.h>
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
20 #include <stdlib.h>
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
21 #include <string.h>
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
22 #include <inttypes.h>
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
23 #include <math.h>
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
24
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
25 #include "config.h"
19830
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
26 #include "cpudetect.h"
18608
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
27
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
28 #include "mp_msg.h"
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
29
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
30 #ifdef HAVE_MALLOC_H
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
31 #include <malloc.h>
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
32 #endif
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
33
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
34 #include "img_format.h"
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
35 #include "mp_image.h"
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
36 #include "vf.h"
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
37 #include "libvo/fastmemcpy.h"
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
38
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
39 #define MIN(a,b) ((a) > (b) ? (b) : (a))
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
40 #define MAX(a,b) ((a) < (b) ? (b) : (a))
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
41 #define ABS(a) ((a) > 0 ? (a) : (-(a)))
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
42
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
43 #define MIN3(a,b,c) MIN(MIN(a,b),c)
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
44 #define MAX3(a,b,c) MAX(MAX(a,b),c)
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
45
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
46 //===========================================================================//
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
47
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
48 struct vf_priv_s {
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
49 int mode;
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
50 int parity;
18917
d9a75b26da6c Add a new video pts tracking mode, enabled by option -correct-pts.
uau
parents: 18657
diff changeset
51 int buffered_i;
d9a75b26da6c Add a new video pts tracking mode, enabled by option -correct-pts.
uau
parents: 18657
diff changeset
52 int buffered_tff;
d9a75b26da6c Add a new video pts tracking mode, enabled by option -correct-pts.
uau
parents: 18657
diff changeset
53 double buffered_pts;
d9a75b26da6c Add a new video pts tracking mode, enabled by option -correct-pts.
uau
parents: 18657
diff changeset
54 mp_image_t *buffered_mpi;
18608
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
55 int stride[3];
18625
dea3a68e2589 simplify
michael
parents: 18623
diff changeset
56 uint8_t *ref[4][3];
18608
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
57 };
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
58
19830
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
59 static void (*filter_line)(struct vf_priv_s *p, uint8_t *dst, uint8_t *prev, uint8_t *cur, uint8_t *next, int w, int refs, int parity);
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
60
18608
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
61 static void store_ref(struct vf_priv_s *p, uint8_t *src[3], int src_stride[3], int width, int height){
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
62 int i;
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
63
18625
dea3a68e2589 simplify
michael
parents: 18623
diff changeset
64 memcpy (p->ref[3], p->ref[0], sizeof(uint8_t *)*3);
dea3a68e2589 simplify
michael
parents: 18623
diff changeset
65 memmove(p->ref[0], p->ref[1], sizeof(uint8_t *)*3*3);
18608
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
66
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
67 for(i=0; i<3; i++){
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
68 int is_chroma= !!i;
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
69
18625
dea3a68e2589 simplify
michael
parents: 18623
diff changeset
70 memcpy_pic(p->ref[2][i], src[i], width>>is_chroma, height>>is_chroma, p->stride[i], src_stride[i]);
18608
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
71 }
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
72 }
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
73
19888
2ce14efa8917 Replace preprocessor hacks to work around compilers not supporting named
diego
parents: 19830
diff changeset
74 #if defined(HAVE_MMX) && defined(NAMED_ASM_ARGS)
19830
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
75
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
76 #define LOAD4(mem,dst) \
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
77 "movd "mem", "#dst" \n\t"\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
78 "punpcklbw %%mm7, "#dst" \n\t"
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
79
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
80 #define PABS(tmp,dst) \
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
81 "pxor "#tmp", "#tmp" \n\t"\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
82 "psubw "#dst", "#tmp" \n\t"\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
83 "pmaxsw "#tmp", "#dst" \n\t"
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
84
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
85 #define CHECK(pj,mj) \
20010
5ae0a24fbed8 Compile fix for odd versions of binutils, namely the ones in Fedora.
gpoirier
parents: 19942
diff changeset
86 "movq "#pj"(%[cur],%[mrefs]), %%mm2 \n\t" /* cur[x-refs-1+j] */\
5ae0a24fbed8 Compile fix for odd versions of binutils, namely the ones in Fedora.
gpoirier
parents: 19942
diff changeset
87 "movq "#mj"(%[cur],%[prefs]), %%mm3 \n\t" /* cur[x+refs-1-j] */\
19830
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
88 "movq %%mm2, %%mm4 \n\t"\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
89 "movq %%mm2, %%mm5 \n\t"\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
90 "pxor %%mm3, %%mm4 \n\t"\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
91 "pavgb %%mm3, %%mm5 \n\t"\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
92 "pand %[pb1], %%mm4 \n\t"\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
93 "psubusb %%mm4, %%mm5 \n\t"\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
94 "psrlq $8, %%mm5 \n\t"\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
95 "punpcklbw %%mm7, %%mm5 \n\t" /* (cur[x-refs+j] + cur[x+refs-j])>>1 */\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
96 "movq %%mm2, %%mm4 \n\t"\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
97 "psubusb %%mm3, %%mm2 \n\t"\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
98 "psubusb %%mm4, %%mm3 \n\t"\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
99 "pmaxub %%mm3, %%mm2 \n\t"\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
100 "movq %%mm2, %%mm3 \n\t"\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
101 "movq %%mm2, %%mm4 \n\t" /* ABS(cur[x-refs-1+j] - cur[x+refs-1-j]) */\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
102 "psrlq $8, %%mm3 \n\t" /* ABS(cur[x-refs +j] - cur[x+refs -j]) */\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
103 "psrlq $16, %%mm4 \n\t" /* ABS(cur[x-refs+1+j] - cur[x+refs+1-j]) */\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
104 "punpcklbw %%mm7, %%mm2 \n\t"\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
105 "punpcklbw %%mm7, %%mm3 \n\t"\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
106 "punpcklbw %%mm7, %%mm4 \n\t"\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
107 "paddw %%mm3, %%mm2 \n\t"\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
108 "paddw %%mm4, %%mm2 \n\t" /* score */
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
109
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
110 #define CHECK1 \
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
111 "movq %%mm0, %%mm3 \n\t"\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
112 "pcmpgtw %%mm2, %%mm3 \n\t" /* if(score < spatial_score) */\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
113 "pminsw %%mm2, %%mm0 \n\t" /* spatial_score= score; */\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
114 "movq %%mm3, %%mm6 \n\t"\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
115 "pand %%mm3, %%mm5 \n\t"\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
116 "pandn %%mm1, %%mm3 \n\t"\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
117 "por %%mm5, %%mm3 \n\t"\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
118 "movq %%mm3, %%mm1 \n\t" /* spatial_pred= (cur[x-refs+j] + cur[x+refs-j])>>1; */
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
119
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
120 #define CHECK2 /* pretend not to have checked dir=2 if dir=1 was bad.\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
121 hurts both quality and speed, but matches the C version. */\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
122 "paddw %[pw1], %%mm6 \n\t"\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
123 "psllw $14, %%mm6 \n\t"\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
124 "paddsw %%mm6, %%mm2 \n\t"\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
125 "movq %%mm0, %%mm3 \n\t"\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
126 "pcmpgtw %%mm2, %%mm3 \n\t"\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
127 "pminsw %%mm2, %%mm0 \n\t"\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
128 "pand %%mm3, %%mm5 \n\t"\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
129 "pandn %%mm1, %%mm3 \n\t"\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
130 "por %%mm5, %%mm3 \n\t"\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
131 "movq %%mm3, %%mm1 \n\t"
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
132
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
133 static void filter_line_mmx2(struct vf_priv_s *p, uint8_t *dst, uint8_t *prev, uint8_t *cur, uint8_t *next, int w, int refs, int parity){
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
134 static const uint64_t pw_1 = 0x0001000100010001ULL;
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
135 static const uint64_t pb_1 = 0x0101010101010101ULL;
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
136 const int mode = p->mode;
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
137 uint64_t tmp0, tmp1, tmp2, tmp3;
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
138 int x;
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
139
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
140 #define FILTER\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
141 for(x=0; x<w; x+=4){\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
142 asm volatile(\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
143 "pxor %%mm7, %%mm7 \n\t"\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
144 LOAD4("(%[cur],%[mrefs])", %%mm0) /* c = cur[x-refs] */\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
145 LOAD4("(%[cur],%[prefs])", %%mm1) /* e = cur[x+refs] */\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
146 LOAD4("(%["prev2"])", %%mm2) /* prev2[x] */\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
147 LOAD4("(%["next2"])", %%mm3) /* next2[x] */\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
148 "movq %%mm3, %%mm4 \n\t"\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
149 "paddw %%mm2, %%mm3 \n\t"\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
150 "psraw $1, %%mm3 \n\t" /* d = (prev2[x] + next2[x])>>1 */\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
151 "movq %%mm0, %[tmp0] \n\t" /* c */\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
152 "movq %%mm3, %[tmp1] \n\t" /* d */\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
153 "movq %%mm1, %[tmp2] \n\t" /* e */\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
154 "psubw %%mm4, %%mm2 \n\t"\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
155 PABS( %%mm4, %%mm2) /* temporal_diff0 */\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
156 LOAD4("(%[prev],%[mrefs])", %%mm3) /* prev[x-refs] */\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
157 LOAD4("(%[prev],%[prefs])", %%mm4) /* prev[x+refs] */\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
158 "psubw %%mm0, %%mm3 \n\t"\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
159 "psubw %%mm1, %%mm4 \n\t"\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
160 PABS( %%mm5, %%mm3)\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
161 PABS( %%mm5, %%mm4)\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
162 "paddw %%mm4, %%mm3 \n\t" /* temporal_diff1 */\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
163 "psrlw $1, %%mm2 \n\t"\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
164 "psrlw $1, %%mm3 \n\t"\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
165 "pmaxsw %%mm3, %%mm2 \n\t"\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
166 LOAD4("(%[next],%[mrefs])", %%mm3) /* next[x-refs] */\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
167 LOAD4("(%[next],%[prefs])", %%mm4) /* next[x+refs] */\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
168 "psubw %%mm0, %%mm3 \n\t"\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
169 "psubw %%mm1, %%mm4 \n\t"\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
170 PABS( %%mm5, %%mm3)\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
171 PABS( %%mm5, %%mm4)\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
172 "paddw %%mm4, %%mm3 \n\t" /* temporal_diff2 */\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
173 "psrlw $1, %%mm3 \n\t"\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
174 "pmaxsw %%mm3, %%mm2 \n\t"\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
175 "movq %%mm2, %[tmp3] \n\t" /* diff */\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
176 \
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
177 "paddw %%mm0, %%mm1 \n\t"\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
178 "paddw %%mm0, %%mm0 \n\t"\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
179 "psubw %%mm1, %%mm0 \n\t"\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
180 "psrlw $1, %%mm1 \n\t" /* spatial_pred */\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
181 PABS( %%mm2, %%mm0) /* ABS(c-e) */\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
182 \
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
183 "movq -1(%[cur],%[mrefs]), %%mm2 \n\t" /* cur[x-refs-1] */\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
184 "movq -1(%[cur],%[prefs]), %%mm3 \n\t" /* cur[x+refs-1] */\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
185 "movq %%mm2, %%mm4 \n\t"\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
186 "psubusb %%mm3, %%mm2 \n\t"\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
187 "psubusb %%mm4, %%mm3 \n\t"\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
188 "pmaxub %%mm3, %%mm2 \n\t"\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
189 "pshufw $9,%%mm2, %%mm3 \n\t"\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
190 "punpcklbw %%mm7, %%mm2 \n\t" /* ABS(cur[x-refs-1] - cur[x+refs-1]) */\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
191 "punpcklbw %%mm7, %%mm3 \n\t" /* ABS(cur[x-refs+1] - cur[x+refs+1]) */\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
192 "paddw %%mm2, %%mm0 \n\t"\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
193 "paddw %%mm3, %%mm0 \n\t"\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
194 "psubw %[pw1], %%mm0 \n\t" /* spatial_score */\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
195 \
20010
5ae0a24fbed8 Compile fix for odd versions of binutils, namely the ones in Fedora.
gpoirier
parents: 19942
diff changeset
196 CHECK(-2,0)\
19830
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
197 CHECK1\
20010
5ae0a24fbed8 Compile fix for odd versions of binutils, namely the ones in Fedora.
gpoirier
parents: 19942
diff changeset
198 CHECK(-3,1)\
19830
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
199 CHECK2\
20010
5ae0a24fbed8 Compile fix for odd versions of binutils, namely the ones in Fedora.
gpoirier
parents: 19942
diff changeset
200 CHECK(0,-2)\
19830
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
201 CHECK1\
20010
5ae0a24fbed8 Compile fix for odd versions of binutils, namely the ones in Fedora.
gpoirier
parents: 19942
diff changeset
202 CHECK(1,-3)\
19830
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
203 CHECK2\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
204 \
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
205 /* if(p->mode<2) ... */\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
206 "movq %[tmp3], %%mm6 \n\t" /* diff */\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
207 "cmp $2, %[mode] \n\t"\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
208 "jge 1f \n\t"\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
209 LOAD4("(%["prev2"],%[mrefs],2)", %%mm2) /* prev2[x-2*refs] */\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
210 LOAD4("(%["next2"],%[mrefs],2)", %%mm4) /* next2[x-2*refs] */\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
211 LOAD4("(%["prev2"],%[prefs],2)", %%mm3) /* prev2[x+2*refs] */\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
212 LOAD4("(%["next2"],%[prefs],2)", %%mm5) /* next2[x+2*refs] */\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
213 "paddw %%mm4, %%mm2 \n\t"\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
214 "paddw %%mm5, %%mm3 \n\t"\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
215 "psrlw $1, %%mm2 \n\t" /* b */\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
216 "psrlw $1, %%mm3 \n\t" /* f */\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
217 "movq %[tmp0], %%mm4 \n\t" /* c */\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
218 "movq %[tmp1], %%mm5 \n\t" /* d */\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
219 "movq %[tmp2], %%mm7 \n\t" /* e */\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
220 "psubw %%mm4, %%mm2 \n\t" /* b-c */\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
221 "psubw %%mm7, %%mm3 \n\t" /* f-e */\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
222 "movq %%mm5, %%mm0 \n\t"\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
223 "psubw %%mm4, %%mm5 \n\t" /* d-c */\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
224 "psubw %%mm7, %%mm0 \n\t" /* d-e */\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
225 "movq %%mm2, %%mm4 \n\t"\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
226 "pminsw %%mm3, %%mm2 \n\t"\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
227 "pmaxsw %%mm4, %%mm3 \n\t"\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
228 "pmaxsw %%mm5, %%mm2 \n\t"\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
229 "pminsw %%mm5, %%mm3 \n\t"\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
230 "pmaxsw %%mm0, %%mm2 \n\t" /* max */\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
231 "pminsw %%mm0, %%mm3 \n\t" /* min */\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
232 "pxor %%mm4, %%mm4 \n\t"\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
233 "pmaxsw %%mm3, %%mm6 \n\t"\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
234 "psubw %%mm2, %%mm4 \n\t" /* -max */\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
235 "pmaxsw %%mm4, %%mm6 \n\t" /* diff= MAX3(diff, min, -max); */\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
236 "1: \n\t"\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
237 \
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
238 "movq %[tmp1], %%mm2 \n\t" /* d */\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
239 "movq %%mm2, %%mm3 \n\t"\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
240 "psubw %%mm6, %%mm2 \n\t" /* d-diff */\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
241 "paddw %%mm6, %%mm3 \n\t" /* d+diff */\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
242 "pmaxsw %%mm2, %%mm1 \n\t"\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
243 "pminsw %%mm3, %%mm1 \n\t" /* d = clip(spatial_pred, d-diff, d+diff); */\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
244 "packuswb %%mm1, %%mm1 \n\t"\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
245 \
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
246 :[tmp0]"=m"(tmp0),\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
247 [tmp1]"=m"(tmp1),\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
248 [tmp2]"=m"(tmp2),\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
249 [tmp3]"=m"(tmp3)\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
250 :[prev] "r"(prev),\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
251 [cur] "r"(cur),\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
252 [next] "r"(next),\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
253 [prefs]"r"((long)refs),\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
254 [mrefs]"r"((long)-refs),\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
255 [pw1] "m"(pw_1),\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
256 [pb1] "m"(pb_1),\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
257 [mode] "g"(mode)\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
258 );\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
259 asm volatile("movd %%mm1, %0" :"=m"(*dst));\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
260 dst += 4;\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
261 prev+= 4;\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
262 cur += 4;\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
263 next+= 4;\
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
264 }
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
265
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
266 if(parity){
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
267 #define prev2 "prev"
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
268 #define next2 "cur"
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
269 FILTER
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
270 #undef prev2
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
271 #undef next2
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
272 }else{
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
273 #define prev2 "cur"
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
274 #define next2 "next"
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
275 FILTER
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
276 #undef prev2
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
277 #undef next2
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
278 }
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
279 }
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
280 #undef LOAD4
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
281 #undef PABS
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
282 #undef CHECK
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
283 #undef CHECK1
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
284 #undef CHECK2
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
285 #undef FILTER
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
286
19888
2ce14efa8917 Replace preprocessor hacks to work around compilers not supporting named
diego
parents: 19830
diff changeset
287 #endif /* defined(HAVE_MMX) && defined(NAMED_ASM_ARGS) */
19830
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
288
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
289 static void filter_line_c(struct vf_priv_s *p, uint8_t *dst, uint8_t *prev, uint8_t *cur, uint8_t *next, int w, int refs, int parity){
19828
1d5832b4b204 cosmetics
lorenm
parents: 18917
diff changeset
290 int x;
1d5832b4b204 cosmetics
lorenm
parents: 18917
diff changeset
291 uint8_t *prev2= parity ? prev : cur ;
1d5832b4b204 cosmetics
lorenm
parents: 18917
diff changeset
292 uint8_t *next2= parity ? cur : next;
19829
4f4a58a3fba8 cosmetics (indentation)
lorenm
parents: 19828
diff changeset
293 for(x=0; x<w; x++){
4f4a58a3fba8 cosmetics (indentation)
lorenm
parents: 19828
diff changeset
294 int c= cur[-refs];
4f4a58a3fba8 cosmetics (indentation)
lorenm
parents: 19828
diff changeset
295 int d= (prev2[0] + next2[0])>>1;
4f4a58a3fba8 cosmetics (indentation)
lorenm
parents: 19828
diff changeset
296 int e= cur[+refs];
4f4a58a3fba8 cosmetics (indentation)
lorenm
parents: 19828
diff changeset
297 int temporal_diff0= ABS(prev2[0] - next2[0]);
4f4a58a3fba8 cosmetics (indentation)
lorenm
parents: 19828
diff changeset
298 int temporal_diff1=( ABS(prev[-refs] - c) + ABS(prev[+refs] - e) )>>1;
4f4a58a3fba8 cosmetics (indentation)
lorenm
parents: 19828
diff changeset
299 int temporal_diff2=( ABS(next[-refs] - c) + ABS(next[+refs] - e) )>>1;
4f4a58a3fba8 cosmetics (indentation)
lorenm
parents: 19828
diff changeset
300 int diff= MAX3(temporal_diff0>>1, temporal_diff1, temporal_diff2);
4f4a58a3fba8 cosmetics (indentation)
lorenm
parents: 19828
diff changeset
301 int spatial_pred= (c+e)>>1;
4f4a58a3fba8 cosmetics (indentation)
lorenm
parents: 19828
diff changeset
302 int spatial_score= ABS(cur[-refs-1] - cur[+refs-1]) + ABS(c-e)
4f4a58a3fba8 cosmetics (indentation)
lorenm
parents: 19828
diff changeset
303 + ABS(cur[-refs+1] - cur[+refs+1]) - 1;
18621
0001a9a9f1a9 cleanup
michael
parents: 18610
diff changeset
304
18657
26f54fad099c better spatial interpolation
michael
parents: 18626
diff changeset
305 #define CHECK(j)\
26f54fad099c better spatial interpolation
michael
parents: 18626
diff changeset
306 { int score= ABS(cur[-refs-1+j] - cur[+refs-1-j])\
26f54fad099c better spatial interpolation
michael
parents: 18626
diff changeset
307 + ABS(cur[-refs +j] - cur[+refs -j])\
26f54fad099c better spatial interpolation
michael
parents: 18626
diff changeset
308 + ABS(cur[-refs+1+j] - cur[+refs+1-j]);\
26f54fad099c better spatial interpolation
michael
parents: 18626
diff changeset
309 if(score < spatial_score){\
26f54fad099c better spatial interpolation
michael
parents: 18626
diff changeset
310 spatial_score= score;\
26f54fad099c better spatial interpolation
michael
parents: 18626
diff changeset
311 spatial_pred= (cur[-refs +j] + cur[+refs -j])>>1;\
18608
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
312
19829
4f4a58a3fba8 cosmetics (indentation)
lorenm
parents: 19828
diff changeset
313 CHECK(-1) CHECK(-2) }} }}
4f4a58a3fba8 cosmetics (indentation)
lorenm
parents: 19828
diff changeset
314 CHECK( 1) CHECK( 2) }} }}
18657
26f54fad099c better spatial interpolation
michael
parents: 18626
diff changeset
315
19829
4f4a58a3fba8 cosmetics (indentation)
lorenm
parents: 19828
diff changeset
316 if(p->mode<2){
4f4a58a3fba8 cosmetics (indentation)
lorenm
parents: 19828
diff changeset
317 int b= (prev2[-2*refs] + next2[-2*refs])>>1;
4f4a58a3fba8 cosmetics (indentation)
lorenm
parents: 19828
diff changeset
318 int f= (prev2[+2*refs] + next2[+2*refs])>>1;
18608
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
319 #if 0
19829
4f4a58a3fba8 cosmetics (indentation)
lorenm
parents: 19828
diff changeset
320 int a= cur[-3*refs];
4f4a58a3fba8 cosmetics (indentation)
lorenm
parents: 19828
diff changeset
321 int g= cur[+3*refs];
4f4a58a3fba8 cosmetics (indentation)
lorenm
parents: 19828
diff changeset
322 int max= MAX3(d-e, d-c, MIN3(MAX(b-c,f-e),MAX(b-c,b-a),MAX(f-g,f-e)) );
4f4a58a3fba8 cosmetics (indentation)
lorenm
parents: 19828
diff changeset
323 int min= MIN3(d-e, d-c, MAX3(MIN(b-c,f-e),MIN(b-c,b-a),MIN(f-g,f-e)) );
18608
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
324 #else
19829
4f4a58a3fba8 cosmetics (indentation)
lorenm
parents: 19828
diff changeset
325 int max= MAX3(d-e, d-c, MIN(b-c, f-e));
4f4a58a3fba8 cosmetics (indentation)
lorenm
parents: 19828
diff changeset
326 int min= MIN3(d-e, d-c, MAX(b-c, f-e));
18608
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
327 #endif
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
328
19829
4f4a58a3fba8 cosmetics (indentation)
lorenm
parents: 19828
diff changeset
329 diff= MAX3(diff, min, -max);
4f4a58a3fba8 cosmetics (indentation)
lorenm
parents: 19828
diff changeset
330 }
18608
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
331
19829
4f4a58a3fba8 cosmetics (indentation)
lorenm
parents: 19828
diff changeset
332 if(spatial_pred > d + diff)
4f4a58a3fba8 cosmetics (indentation)
lorenm
parents: 19828
diff changeset
333 spatial_pred = d + diff;
4f4a58a3fba8 cosmetics (indentation)
lorenm
parents: 19828
diff changeset
334 else if(spatial_pred < d - diff)
4f4a58a3fba8 cosmetics (indentation)
lorenm
parents: 19828
diff changeset
335 spatial_pred = d - diff;
19828
1d5832b4b204 cosmetics
lorenm
parents: 18917
diff changeset
336
19829
4f4a58a3fba8 cosmetics (indentation)
lorenm
parents: 19828
diff changeset
337 dst[0] = spatial_pred;
19828
1d5832b4b204 cosmetics
lorenm
parents: 18917
diff changeset
338
19829
4f4a58a3fba8 cosmetics (indentation)
lorenm
parents: 19828
diff changeset
339 dst++;
4f4a58a3fba8 cosmetics (indentation)
lorenm
parents: 19828
diff changeset
340 cur++;
4f4a58a3fba8 cosmetics (indentation)
lorenm
parents: 19828
diff changeset
341 prev++;
4f4a58a3fba8 cosmetics (indentation)
lorenm
parents: 19828
diff changeset
342 next++;
4f4a58a3fba8 cosmetics (indentation)
lorenm
parents: 19828
diff changeset
343 prev2++;
4f4a58a3fba8 cosmetics (indentation)
lorenm
parents: 19828
diff changeset
344 next2++;
4f4a58a3fba8 cosmetics (indentation)
lorenm
parents: 19828
diff changeset
345 }
19828
1d5832b4b204 cosmetics
lorenm
parents: 18917
diff changeset
346 }
18608
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
347
19828
1d5832b4b204 cosmetics
lorenm
parents: 18917
diff changeset
348 static void filter(struct vf_priv_s *p, uint8_t *dst[3], int dst_stride[3], int width, int height, int parity, int tff){
1d5832b4b204 cosmetics
lorenm
parents: 18917
diff changeset
349 int x, y, i;
1d5832b4b204 cosmetics
lorenm
parents: 18917
diff changeset
350
1d5832b4b204 cosmetics
lorenm
parents: 18917
diff changeset
351 for(i=0; i<3; i++){
1d5832b4b204 cosmetics
lorenm
parents: 18917
diff changeset
352 int is_chroma= !!i;
1d5832b4b204 cosmetics
lorenm
parents: 18917
diff changeset
353 int w= width >>is_chroma;
1d5832b4b204 cosmetics
lorenm
parents: 18917
diff changeset
354 int h= height>>is_chroma;
1d5832b4b204 cosmetics
lorenm
parents: 18917
diff changeset
355 int refs= p->stride[i];
1d5832b4b204 cosmetics
lorenm
parents: 18917
diff changeset
356
1d5832b4b204 cosmetics
lorenm
parents: 18917
diff changeset
357 for(y=0; y<h; y++){
1d5832b4b204 cosmetics
lorenm
parents: 18917
diff changeset
358 if((y ^ parity) & 1){
1d5832b4b204 cosmetics
lorenm
parents: 18917
diff changeset
359 uint8_t *prev= &p->ref[0][i][y*refs];
1d5832b4b204 cosmetics
lorenm
parents: 18917
diff changeset
360 uint8_t *cur = &p->ref[1][i][y*refs];
1d5832b4b204 cosmetics
lorenm
parents: 18917
diff changeset
361 uint8_t *next= &p->ref[2][i][y*refs];
1d5832b4b204 cosmetics
lorenm
parents: 18917
diff changeset
362 uint8_t *dst2= &dst[i][y*dst_stride[i]];
1d5832b4b204 cosmetics
lorenm
parents: 18917
diff changeset
363 filter_line(p, dst2, prev, cur, next, w, refs, parity ^ tff);
18608
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
364 }else{
18610
e0df832235c6 simplify
michael
parents: 18608
diff changeset
365 memcpy(&dst[i][y*dst_stride[i]], &p->ref[1][i][y*refs], w);
18608
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
366 }
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
367 }
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
368 }
19942
0983feab2dc2 Add forgotten emms which caused weird bugs like nan pts values.
reimar
parents: 19923
diff changeset
369 #if defined(HAVE_MMX) && defined(NAMED_ASM_ARGS)
0983feab2dc2 Add forgotten emms which caused weird bugs like nan pts values.
reimar
parents: 19923
diff changeset
370 if(gCpuCaps.hasMMX2) asm volatile("emms \n\t" : : : "memory");
0983feab2dc2 Add forgotten emms which caused weird bugs like nan pts values.
reimar
parents: 19923
diff changeset
371 #endif
18608
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
372 }
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
373
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
374 static int config(struct vf_instance_s* vf,
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
375 int width, int height, int d_width, int d_height,
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
376 unsigned int flags, unsigned int outfmt){
18625
dea3a68e2589 simplify
michael
parents: 18623
diff changeset
377 int i, j;
18608
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
378
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
379 for(i=0; i<3; i++){
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
380 int is_chroma= !!i;
18625
dea3a68e2589 simplify
michael
parents: 18623
diff changeset
381 int w= ((width + 31) & (~31))>>is_chroma;
dea3a68e2589 simplify
michael
parents: 18623
diff changeset
382 int h= ((height+6+ 31) & (~31))>>is_chroma;
18608
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
383
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
384 vf->priv->stride[i]= w;
18625
dea3a68e2589 simplify
michael
parents: 18623
diff changeset
385 for(j=0; j<3; j++)
dea3a68e2589 simplify
michael
parents: 18623
diff changeset
386 vf->priv->ref[j][i]= malloc(w*h*sizeof(uint8_t))+3*w;
18608
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
387 }
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
388
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
389 return vf_next_config(vf,width,height,d_width,d_height,flags,outfmt);
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
390 }
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
391
18917
d9a75b26da6c Add a new video pts tracking mode, enabled by option -correct-pts.
uau
parents: 18657
diff changeset
392 static int continue_buffered_image(struct vf_instance_s *vf);
d9a75b26da6c Add a new video pts tracking mode, enabled by option -correct-pts.
uau
parents: 18657
diff changeset
393 extern int correct_pts;
d9a75b26da6c Add a new video pts tracking mode, enabled by option -correct-pts.
uau
parents: 18657
diff changeset
394
18608
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
395 static int put_image(struct vf_instance_s* vf, mp_image_t *mpi, double pts){
18917
d9a75b26da6c Add a new video pts tracking mode, enabled by option -correct-pts.
uau
parents: 18657
diff changeset
396 int tff;
18608
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
397
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
398 if(vf->priv->parity < 0) {
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
399 if (mpi->fields & MP_IMGFIELD_ORDERED)
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
400 tff = !!(mpi->fields & MP_IMGFIELD_TOP_FIRST);
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
401 else
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
402 tff = 1;
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
403 }
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
404 else tff = (vf->priv->parity&1)^1;
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
405
18625
dea3a68e2589 simplify
michael
parents: 18623
diff changeset
406 store_ref(vf->priv, mpi->planes, mpi->stride, mpi->w, mpi->h);
dea3a68e2589 simplify
michael
parents: 18623
diff changeset
407
18917
d9a75b26da6c Add a new video pts tracking mode, enabled by option -correct-pts.
uau
parents: 18657
diff changeset
408 vf->priv->buffered_mpi = mpi;
d9a75b26da6c Add a new video pts tracking mode, enabled by option -correct-pts.
uau
parents: 18657
diff changeset
409 vf->priv->buffered_tff = tff;
d9a75b26da6c Add a new video pts tracking mode, enabled by option -correct-pts.
uau
parents: 18657
diff changeset
410 vf->priv->buffered_i = 0;
d9a75b26da6c Add a new video pts tracking mode, enabled by option -correct-pts.
uau
parents: 18657
diff changeset
411 vf->priv->buffered_pts = pts;
d9a75b26da6c Add a new video pts tracking mode, enabled by option -correct-pts.
uau
parents: 18657
diff changeset
412
d9a75b26da6c Add a new video pts tracking mode, enabled by option -correct-pts.
uau
parents: 18657
diff changeset
413 return continue_buffered_image(vf);
d9a75b26da6c Add a new video pts tracking mode, enabled by option -correct-pts.
uau
parents: 18657
diff changeset
414 }
d9a75b26da6c Add a new video pts tracking mode, enabled by option -correct-pts.
uau
parents: 18657
diff changeset
415
d9a75b26da6c Add a new video pts tracking mode, enabled by option -correct-pts.
uau
parents: 18657
diff changeset
416 static int continue_buffered_image(struct vf_instance_s *vf)
d9a75b26da6c Add a new video pts tracking mode, enabled by option -correct-pts.
uau
parents: 18657
diff changeset
417 {
d9a75b26da6c Add a new video pts tracking mode, enabled by option -correct-pts.
uau
parents: 18657
diff changeset
418 mp_image_t *mpi = vf->priv->buffered_mpi;
d9a75b26da6c Add a new video pts tracking mode, enabled by option -correct-pts.
uau
parents: 18657
diff changeset
419 int tff = vf->priv->buffered_tff;
d9a75b26da6c Add a new video pts tracking mode, enabled by option -correct-pts.
uau
parents: 18657
diff changeset
420 double pts = vf->priv->buffered_pts;
d9a75b26da6c Add a new video pts tracking mode, enabled by option -correct-pts.
uau
parents: 18657
diff changeset
421 int i;
d9a75b26da6c Add a new video pts tracking mode, enabled by option -correct-pts.
uau
parents: 18657
diff changeset
422 int ret=0;
d9a75b26da6c Add a new video pts tracking mode, enabled by option -correct-pts.
uau
parents: 18657
diff changeset
423 mp_image_t *dmpi;
d9a75b26da6c Add a new video pts tracking mode, enabled by option -correct-pts.
uau
parents: 18657
diff changeset
424
d9a75b26da6c Add a new video pts tracking mode, enabled by option -correct-pts.
uau
parents: 18657
diff changeset
425 pts += vf->priv->buffered_i * .02; // XXX not right
d9a75b26da6c Add a new video pts tracking mode, enabled by option -correct-pts.
uau
parents: 18657
diff changeset
426
d9a75b26da6c Add a new video pts tracking mode, enabled by option -correct-pts.
uau
parents: 18657
diff changeset
427 for(i = vf->priv->buffered_i; i<=(vf->priv->mode&1); i++){
18608
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
428 dmpi=vf_get_image(vf->next,mpi->imgfmt,
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
429 MP_IMGTYPE_TEMP,
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
430 MP_IMGFLAG_ACCEPT_STRIDE|MP_IMGFLAG_PREFER_ALIGNED_STRIDE,
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
431 mpi->width,mpi->height);
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
432 vf_clone_mpi_attributes(dmpi, mpi);
18625
dea3a68e2589 simplify
michael
parents: 18623
diff changeset
433 filter(vf->priv, dmpi->planes, dmpi->stride, mpi->w, mpi->h, i ^ tff ^ 1, tff);
18917
d9a75b26da6c Add a new video pts tracking mode, enabled by option -correct-pts.
uau
parents: 18657
diff changeset
434 if (correct_pts && i < (vf->priv->mode & 1))
d9a75b26da6c Add a new video pts tracking mode, enabled by option -correct-pts.
uau
parents: 18657
diff changeset
435 vf_queue_frame(vf, continue_buffered_image);
18608
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
436 ret |= vf_next_put_image(vf, dmpi, pts /*FIXME*/);
18917
d9a75b26da6c Add a new video pts tracking mode, enabled by option -correct-pts.
uau
parents: 18657
diff changeset
437 if (correct_pts)
d9a75b26da6c Add a new video pts tracking mode, enabled by option -correct-pts.
uau
parents: 18657
diff changeset
438 break;
18626
2e2b8f34a69e make spatial interlacing check optional
michael
parents: 18625
diff changeset
439 if(i<(vf->priv->mode&1))
18608
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
440 vf_next_control(vf, VFCTRL_FLIP_PAGE, NULL);
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
441 }
18917
d9a75b26da6c Add a new video pts tracking mode, enabled by option -correct-pts.
uau
parents: 18657
diff changeset
442 vf->priv->buffered_i = 1;
18608
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
443 return ret;
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
444 }
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
445
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
446 static void uninit(struct vf_instance_s* vf){
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
447 int i;
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
448 if(!vf->priv) return;
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
449
18625
dea3a68e2589 simplify
michael
parents: 18623
diff changeset
450 for(i=0; i<3*3; i++){
dea3a68e2589 simplify
michael
parents: 18623
diff changeset
451 uint8_t **p= &vf->priv->ref[i%3][i/3];
dea3a68e2589 simplify
michael
parents: 18623
diff changeset
452 if(*p) free(*p - 3*vf->priv->stride[i/3]);
dea3a68e2589 simplify
michael
parents: 18623
diff changeset
453 *p= NULL;
18608
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
454 }
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
455 free(vf->priv);
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
456 vf->priv=NULL;
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
457 }
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
458
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
459 //===========================================================================//
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
460 static int query_format(struct vf_instance_s* vf, unsigned int fmt){
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
461 switch(fmt){
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
462 case IMGFMT_YV12:
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
463 case IMGFMT_I420:
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
464 case IMGFMT_IYUV:
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
465 case IMGFMT_Y800:
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
466 case IMGFMT_Y8:
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
467 return vf_next_query_format(vf,fmt);
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
468 }
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
469 return 0;
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
470 }
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
471
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
472 static int open(vf_instance_t *vf, char* args){
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
473
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
474 vf->config=config;
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
475 vf->put_image=put_image;
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
476 vf->query_format=query_format;
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
477 vf->uninit=uninit;
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
478 vf->priv=malloc(sizeof(struct vf_priv_s));
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
479 memset(vf->priv, 0, sizeof(struct vf_priv_s));
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
480
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
481 vf->priv->mode=0;
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
482 vf->priv->parity= -1;
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
483
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
484 if (args) sscanf(args, "%d:%d", &vf->priv->mode, &vf->priv->parity);
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
485
19830
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
486 filter_line = filter_line_c;
19923
7dcd96645457 Fix compilation if HAVE_MMX is defined but NAMED_ASM_ARGS isn't
reimar
parents: 19888
diff changeset
487 #if defined(HAVE_MMX) && defined(NAMED_ASM_ARGS)
19830
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
488 if(gCpuCaps.hasMMX2) filter_line = filter_line_mmx2;
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
489 #endif
a1a997376658 mmx implementation of yadif. 3.5x faster.
lorenm
parents: 19829
diff changeset
490
18608
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
491 return 1;
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
492 }
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
493
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
494 vf_info_t vf_info_yadif = {
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
495 "Yet Another DeInterlacing Filter",
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
496 "yadif",
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
497 "Michael Niedermayer",
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
498 "",
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
499 open,
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
500 NULL
a80c7de8a4ba yet another deinterlacing filter
michael
parents:
diff changeset
501 };