annotate imgresample.c @ 3436:3ccec186532f libavcodec

moving utf8 reading function to libavutil
author michael
date Fri, 07 Jul 2006 22:43:32 +0000
parents ad3d12d04645
children f60df50f50aa
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
0
986e461dc072 Initial revision
glantau
parents:
diff changeset
1 /*
2967
ef2149182f1c COSMETICS: Remove all trailing whitespace.
diego
parents: 2904
diff changeset
2 * High quality image resampling with polyphase filters
429
718a22dc121f license/copyright change
glantau
parents: 412
diff changeset
3 * Copyright (c) 2001 Fabrice Bellard.
0
986e461dc072 Initial revision
glantau
parents:
diff changeset
4 *
429
718a22dc121f license/copyright change
glantau
parents: 412
diff changeset
5 * This library is free software; you can redistribute it and/or
718a22dc121f license/copyright change
glantau
parents: 412
diff changeset
6 * modify it under the terms of the GNU Lesser General Public
718a22dc121f license/copyright change
glantau
parents: 412
diff changeset
7 * License as published by the Free Software Foundation; either
718a22dc121f license/copyright change
glantau
parents: 412
diff changeset
8 * version 2 of the License, or (at your option) any later version.
0
986e461dc072 Initial revision
glantau
parents:
diff changeset
9 *
429
718a22dc121f license/copyright change
glantau
parents: 412
diff changeset
10 * This library is distributed in the hope that it will be useful,
0
986e461dc072 Initial revision
glantau
parents:
diff changeset
11 * but WITHOUT ANY WARRANTY; without even the implied warranty of
429
718a22dc121f license/copyright change
glantau
parents: 412
diff changeset
12 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
718a22dc121f license/copyright change
glantau
parents: 412
diff changeset
13 * Lesser General Public License for more details.
0
986e461dc072 Initial revision
glantau
parents:
diff changeset
14 *
429
718a22dc121f license/copyright change
glantau
parents: 412
diff changeset
15 * You should have received a copy of the GNU Lesser General Public
718a22dc121f license/copyright change
glantau
parents: 412
diff changeset
16 * License along with this library; if not, write to the Free Software
3036
0b546eab515d Update licensing information: The FSF changed postal address.
diego
parents: 2979
diff changeset
17 * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
0
986e461dc072 Initial revision
glantau
parents:
diff changeset
18 */
2967
ef2149182f1c COSMETICS: Remove all trailing whitespace.
diego
parents: 2904
diff changeset
19
1106
1e39f273ecd6 per file doxy
michaelni
parents: 1064
diff changeset
20 /**
1e39f273ecd6 per file doxy
michaelni
parents: 1064
diff changeset
21 * @file imgresample.c
1e39f273ecd6 per file doxy
michaelni
parents: 1064
diff changeset
22 * High quality image resampling with polyphase filters .
1e39f273ecd6 per file doxy
michaelni
parents: 1064
diff changeset
23 */
2967
ef2149182f1c COSMETICS: Remove all trailing whitespace.
diego
parents: 2904
diff changeset
24
396
fce0a2520551 removed useless header includes - use av memory functions
glantau
parents: 18
diff changeset
25 #include "avcodec.h"
3249
ad3d12d04645 Introduce swscale interface in libavcodec
lucabe
parents: 3089
diff changeset
26 #include "swscale.h"
0
986e461dc072 Initial revision
glantau
parents:
diff changeset
27 #include "dsputil.h"
986e461dc072 Initial revision
glantau
parents:
diff changeset
28
17
b69fe46fd708 Adding fastmemcpy stuff to speedup mplayer project
nickols_k
parents: 2
diff changeset
29 #ifdef USE_FASTMEMCPY
b69fe46fd708 Adding fastmemcpy stuff to speedup mplayer project
nickols_k
parents: 2
diff changeset
30 #include "fastmemcpy.h"
b69fe46fd708 Adding fastmemcpy stuff to speedup mplayer project
nickols_k
parents: 2
diff changeset
31 #endif
b69fe46fd708 Adding fastmemcpy stuff to speedup mplayer project
nickols_k
parents: 2
diff changeset
32
0
986e461dc072 Initial revision
glantau
parents:
diff changeset
33 #define NB_COMPONENTS 3
986e461dc072 Initial revision
glantau
parents:
diff changeset
34
986e461dc072 Initial revision
glantau
parents:
diff changeset
35 #define PHASE_BITS 4
986e461dc072 Initial revision
glantau
parents:
diff changeset
36 #define NB_PHASES (1 << PHASE_BITS)
986e461dc072 Initial revision
glantau
parents:
diff changeset
37 #define NB_TAPS 4
986e461dc072 Initial revision
glantau
parents:
diff changeset
38 #define FCENTER 1 /* index of the center of the filter */
630
b4ee42142ad1 croping patch by (talus25 at speakeasy dot net) with fixes from atmos & me
michaelni
parents: 429
diff changeset
39 //#define TEST 1 /* Test it */
0
986e461dc072 Initial revision
glantau
parents:
diff changeset
40
986e461dc072 Initial revision
glantau
parents:
diff changeset
41 #define POS_FRAC_BITS 16
986e461dc072 Initial revision
glantau
parents:
diff changeset
42 #define POS_FRAC (1 << POS_FRAC_BITS)
986e461dc072 Initial revision
glantau
parents:
diff changeset
43 /* 6 bits precision is needed for MMX */
986e461dc072 Initial revision
glantau
parents:
diff changeset
44 #define FILTER_BITS 8
986e461dc072 Initial revision
glantau
parents:
diff changeset
45
986e461dc072 Initial revision
glantau
parents:
diff changeset
46 #define LINE_BUF_HEIGHT (NB_TAPS * 4)
986e461dc072 Initial revision
glantau
parents:
diff changeset
47
986e461dc072 Initial revision
glantau
parents:
diff changeset
48 struct ImgReSampleContext {
1928
0c23a5564489 padding support in ffmpeg patch by (Todd Kirby <doubleshot at pacbell dot net>)
michael
parents: 1488
diff changeset
49 int iwidth, iheight, owidth, oheight;
0c23a5564489 padding support in ffmpeg patch by (Todd Kirby <doubleshot at pacbell dot net>)
michael
parents: 1488
diff changeset
50 int topBand, bottomBand, leftBand, rightBand;
0c23a5564489 padding support in ffmpeg patch by (Todd Kirby <doubleshot at pacbell dot net>)
michael
parents: 1488
diff changeset
51 int padtop, padbottom, padleft, padright;
0c23a5564489 padding support in ffmpeg patch by (Todd Kirby <doubleshot at pacbell dot net>)
michael
parents: 1488
diff changeset
52 int pad_owidth, pad_oheight;
0
986e461dc072 Initial revision
glantau
parents:
diff changeset
53 int h_incr, v_incr;
3089
072dbc669253 MSVC-compatible __align8/__align16 declaration
diego
parents: 3036
diff changeset
54 DECLARE_ALIGNED_8(int16_t, h_filters[NB_PHASES][NB_TAPS]); /* horizontal filters */
072dbc669253 MSVC-compatible __align8/__align16 declaration
diego
parents: 3036
diff changeset
55 DECLARE_ALIGNED_8(int16_t, v_filters[NB_PHASES][NB_TAPS]); /* vertical filters */
1064
b32afefe7d33 * UINTX -> uintx_t INTX -> intx_t
kabi
parents: 1057
diff changeset
56 uint8_t *line_buf;
0
986e461dc072 Initial revision
glantau
parents:
diff changeset
57 };
986e461dc072 Initial revision
glantau
parents:
diff changeset
58
2082
3dc9bbe1b152 polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents: 2064
diff changeset
59 void av_build_filter(int16_t *filter, double factor, int tap_count, int phase_count, int scale, int type);
3dc9bbe1b152 polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents: 2064
diff changeset
60
0
986e461dc072 Initial revision
glantau
parents:
diff changeset
61 static inline int get_phase(int pos)
986e461dc072 Initial revision
glantau
parents:
diff changeset
62 {
986e461dc072 Initial revision
glantau
parents:
diff changeset
63 return ((pos) >> (POS_FRAC_BITS - PHASE_BITS)) & ((1 << PHASE_BITS) - 1);
986e461dc072 Initial revision
glantau
parents:
diff changeset
64 }
986e461dc072 Initial revision
glantau
parents:
diff changeset
65
986e461dc072 Initial revision
glantau
parents:
diff changeset
66 /* This function must be optimized */
1488
766a2f4edbea avcodec const correctness patch by (Drew Hess <dhess at ilm dot com>)
michaelni
parents: 1106
diff changeset
67 static void h_resample_fast(uint8_t *dst, int dst_width, const uint8_t *src,
2979
bfabfdf9ce55 COSMETICS: tabs --> spaces, some prettyprinting
diego
parents: 2967
diff changeset
68 int src_width, int src_start, int src_incr,
bfabfdf9ce55 COSMETICS: tabs --> spaces, some prettyprinting
diego
parents: 2967
diff changeset
69 int16_t *filters)
0
986e461dc072 Initial revision
glantau
parents:
diff changeset
70 {
986e461dc072 Initial revision
glantau
parents:
diff changeset
71 int src_pos, phase, sum, i;
1488
766a2f4edbea avcodec const correctness patch by (Drew Hess <dhess at ilm dot com>)
michaelni
parents: 1106
diff changeset
72 const uint8_t *s;
1064
b32afefe7d33 * UINTX -> uintx_t INTX -> intx_t
kabi
parents: 1057
diff changeset
73 int16_t *filter;
0
986e461dc072 Initial revision
glantau
parents:
diff changeset
74
986e461dc072 Initial revision
glantau
parents:
diff changeset
75 src_pos = src_start;
986e461dc072 Initial revision
glantau
parents:
diff changeset
76 for(i=0;i<dst_width;i++) {
986e461dc072 Initial revision
glantau
parents:
diff changeset
77 #ifdef TEST
986e461dc072 Initial revision
glantau
parents:
diff changeset
78 /* test */
986e461dc072 Initial revision
glantau
parents:
diff changeset
79 if ((src_pos >> POS_FRAC_BITS) < 0 ||
986e461dc072 Initial revision
glantau
parents:
diff changeset
80 (src_pos >> POS_FRAC_BITS) > (src_width - NB_TAPS))
653
714795876872 Change abort() calls to av_abort() calls.
philipjsg
parents: 644
diff changeset
81 av_abort();
0
986e461dc072 Initial revision
glantau
parents:
diff changeset
82 #endif
986e461dc072 Initial revision
glantau
parents:
diff changeset
83 s = src + (src_pos >> POS_FRAC_BITS);
986e461dc072 Initial revision
glantau
parents:
diff changeset
84 phase = get_phase(src_pos);
986e461dc072 Initial revision
glantau
parents:
diff changeset
85 filter = filters + phase * NB_TAPS;
986e461dc072 Initial revision
glantau
parents:
diff changeset
86 #if NB_TAPS == 4
986e461dc072 Initial revision
glantau
parents:
diff changeset
87 sum = s[0] * filter[0] +
986e461dc072 Initial revision
glantau
parents:
diff changeset
88 s[1] * filter[1] +
986e461dc072 Initial revision
glantau
parents:
diff changeset
89 s[2] * filter[2] +
986e461dc072 Initial revision
glantau
parents:
diff changeset
90 s[3] * filter[3];
986e461dc072 Initial revision
glantau
parents:
diff changeset
91 #else
986e461dc072 Initial revision
glantau
parents:
diff changeset
92 {
986e461dc072 Initial revision
glantau
parents:
diff changeset
93 int j;
986e461dc072 Initial revision
glantau
parents:
diff changeset
94 sum = 0;
986e461dc072 Initial revision
glantau
parents:
diff changeset
95 for(j=0;j<NB_TAPS;j++)
986e461dc072 Initial revision
glantau
parents:
diff changeset
96 sum += s[j] * filter[j];
986e461dc072 Initial revision
glantau
parents:
diff changeset
97 }
986e461dc072 Initial revision
glantau
parents:
diff changeset
98 #endif
986e461dc072 Initial revision
glantau
parents:
diff changeset
99 sum = sum >> FILTER_BITS;
986e461dc072 Initial revision
glantau
parents:
diff changeset
100 if (sum < 0)
986e461dc072 Initial revision
glantau
parents:
diff changeset
101 sum = 0;
986e461dc072 Initial revision
glantau
parents:
diff changeset
102 else if (sum > 255)
986e461dc072 Initial revision
glantau
parents:
diff changeset
103 sum = 255;
986e461dc072 Initial revision
glantau
parents:
diff changeset
104 dst[0] = sum;
986e461dc072 Initial revision
glantau
parents:
diff changeset
105 src_pos += src_incr;
986e461dc072 Initial revision
glantau
parents:
diff changeset
106 dst++;
986e461dc072 Initial revision
glantau
parents:
diff changeset
107 }
986e461dc072 Initial revision
glantau
parents:
diff changeset
108 }
986e461dc072 Initial revision
glantau
parents:
diff changeset
109
986e461dc072 Initial revision
glantau
parents:
diff changeset
110 /* This function must be optimized */
1488
766a2f4edbea avcodec const correctness patch by (Drew Hess <dhess at ilm dot com>)
michaelni
parents: 1106
diff changeset
111 static void v_resample(uint8_t *dst, int dst_width, const uint8_t *src,
2979
bfabfdf9ce55 COSMETICS: tabs --> spaces, some prettyprinting
diego
parents: 2967
diff changeset
112 int wrap, int16_t *filter)
0
986e461dc072 Initial revision
glantau
parents:
diff changeset
113 {
986e461dc072 Initial revision
glantau
parents:
diff changeset
114 int sum, i;
1488
766a2f4edbea avcodec const correctness patch by (Drew Hess <dhess at ilm dot com>)
michaelni
parents: 1106
diff changeset
115 const uint8_t *s;
0
986e461dc072 Initial revision
glantau
parents:
diff changeset
116
986e461dc072 Initial revision
glantau
parents:
diff changeset
117 s = src;
986e461dc072 Initial revision
glantau
parents:
diff changeset
118 for(i=0;i<dst_width;i++) {
986e461dc072 Initial revision
glantau
parents:
diff changeset
119 #if NB_TAPS == 4
986e461dc072 Initial revision
glantau
parents:
diff changeset
120 sum = s[0 * wrap] * filter[0] +
986e461dc072 Initial revision
glantau
parents:
diff changeset
121 s[1 * wrap] * filter[1] +
986e461dc072 Initial revision
glantau
parents:
diff changeset
122 s[2 * wrap] * filter[2] +
986e461dc072 Initial revision
glantau
parents:
diff changeset
123 s[3 * wrap] * filter[3];
986e461dc072 Initial revision
glantau
parents:
diff changeset
124 #else
986e461dc072 Initial revision
glantau
parents:
diff changeset
125 {
986e461dc072 Initial revision
glantau
parents:
diff changeset
126 int j;
1064
b32afefe7d33 * UINTX -> uintx_t INTX -> intx_t
kabi
parents: 1057
diff changeset
127 uint8_t *s1 = s;
0
986e461dc072 Initial revision
glantau
parents:
diff changeset
128
986e461dc072 Initial revision
glantau
parents:
diff changeset
129 sum = 0;
986e461dc072 Initial revision
glantau
parents:
diff changeset
130 for(j=0;j<NB_TAPS;j++) {
986e461dc072 Initial revision
glantau
parents:
diff changeset
131 sum += s1[0] * filter[j];
986e461dc072 Initial revision
glantau
parents:
diff changeset
132 s1 += wrap;
986e461dc072 Initial revision
glantau
parents:
diff changeset
133 }
986e461dc072 Initial revision
glantau
parents:
diff changeset
134 }
986e461dc072 Initial revision
glantau
parents:
diff changeset
135 #endif
986e461dc072 Initial revision
glantau
parents:
diff changeset
136 sum = sum >> FILTER_BITS;
986e461dc072 Initial revision
glantau
parents:
diff changeset
137 if (sum < 0)
986e461dc072 Initial revision
glantau
parents:
diff changeset
138 sum = 0;
986e461dc072 Initial revision
glantau
parents:
diff changeset
139 else if (sum > 255)
986e461dc072 Initial revision
glantau
parents:
diff changeset
140 sum = 255;
986e461dc072 Initial revision
glantau
parents:
diff changeset
141 dst[0] = sum;
986e461dc072 Initial revision
glantau
parents:
diff changeset
142 dst++;
986e461dc072 Initial revision
glantau
parents:
diff changeset
143 s++;
986e461dc072 Initial revision
glantau
parents:
diff changeset
144 }
986e461dc072 Initial revision
glantau
parents:
diff changeset
145 }
986e461dc072 Initial revision
glantau
parents:
diff changeset
146
2
2e2c46c87460 fixed config for direct mplayer build compatibility
glantau
parents: 0
diff changeset
147 #ifdef HAVE_MMX
0
986e461dc072 Initial revision
glantau
parents:
diff changeset
148
986e461dc072 Initial revision
glantau
parents:
diff changeset
149 #include "i386/mmx.h"
986e461dc072 Initial revision
glantau
parents:
diff changeset
150
986e461dc072 Initial revision
glantau
parents:
diff changeset
151 #define FILTER4(reg) \
986e461dc072 Initial revision
glantau
parents:
diff changeset
152 {\
986e461dc072 Initial revision
glantau
parents:
diff changeset
153 s = src + (src_pos >> POS_FRAC_BITS);\
986e461dc072 Initial revision
glantau
parents:
diff changeset
154 phase = get_phase(src_pos);\
986e461dc072 Initial revision
glantau
parents:
diff changeset
155 filter = filters + phase * NB_TAPS;\
986e461dc072 Initial revision
glantau
parents:
diff changeset
156 movq_m2r(*s, reg);\
986e461dc072 Initial revision
glantau
parents:
diff changeset
157 punpcklbw_r2r(mm7, reg);\
986e461dc072 Initial revision
glantau
parents:
diff changeset
158 movq_m2r(*filter, mm6);\
986e461dc072 Initial revision
glantau
parents:
diff changeset
159 pmaddwd_r2r(reg, mm6);\
986e461dc072 Initial revision
glantau
parents:
diff changeset
160 movq_r2r(mm6, reg);\
986e461dc072 Initial revision
glantau
parents:
diff changeset
161 psrlq_i2r(32, reg);\
986e461dc072 Initial revision
glantau
parents:
diff changeset
162 paddd_r2r(mm6, reg);\
986e461dc072 Initial revision
glantau
parents:
diff changeset
163 psrad_i2r(FILTER_BITS, reg);\
986e461dc072 Initial revision
glantau
parents:
diff changeset
164 src_pos += src_incr;\
986e461dc072 Initial revision
glantau
parents:
diff changeset
165 }
986e461dc072 Initial revision
glantau
parents:
diff changeset
166
986e461dc072 Initial revision
glantau
parents:
diff changeset
167 #define DUMP(reg) movq_r2m(reg, tmp); printf(#reg "=%016Lx\n", tmp.uq);
986e461dc072 Initial revision
glantau
parents:
diff changeset
168
986e461dc072 Initial revision
glantau
parents:
diff changeset
169 /* XXX: do four pixels at a time */
1488
766a2f4edbea avcodec const correctness patch by (Drew Hess <dhess at ilm dot com>)
michaelni
parents: 1106
diff changeset
170 static void h_resample_fast4_mmx(uint8_t *dst, int dst_width,
2979
bfabfdf9ce55 COSMETICS: tabs --> spaces, some prettyprinting
diego
parents: 2967
diff changeset
171 const uint8_t *src, int src_width,
1064
b32afefe7d33 * UINTX -> uintx_t INTX -> intx_t
kabi
parents: 1057
diff changeset
172 int src_start, int src_incr, int16_t *filters)
0
986e461dc072 Initial revision
glantau
parents:
diff changeset
173 {
986e461dc072 Initial revision
glantau
parents:
diff changeset
174 int src_pos, phase;
1488
766a2f4edbea avcodec const correctness patch by (Drew Hess <dhess at ilm dot com>)
michaelni
parents: 1106
diff changeset
175 const uint8_t *s;
1064
b32afefe7d33 * UINTX -> uintx_t INTX -> intx_t
kabi
parents: 1057
diff changeset
176 int16_t *filter;
0
986e461dc072 Initial revision
glantau
parents:
diff changeset
177 mmx_t tmp;
2967
ef2149182f1c COSMETICS: Remove all trailing whitespace.
diego
parents: 2904
diff changeset
178
0
986e461dc072 Initial revision
glantau
parents:
diff changeset
179 src_pos = src_start;
986e461dc072 Initial revision
glantau
parents:
diff changeset
180 pxor_r2r(mm7, mm7);
986e461dc072 Initial revision
glantau
parents:
diff changeset
181
986e461dc072 Initial revision
glantau
parents:
diff changeset
182 while (dst_width >= 4) {
986e461dc072 Initial revision
glantau
parents:
diff changeset
183
986e461dc072 Initial revision
glantau
parents:
diff changeset
184 FILTER4(mm0);
986e461dc072 Initial revision
glantau
parents:
diff changeset
185 FILTER4(mm1);
986e461dc072 Initial revision
glantau
parents:
diff changeset
186 FILTER4(mm2);
986e461dc072 Initial revision
glantau
parents:
diff changeset
187 FILTER4(mm3);
986e461dc072 Initial revision
glantau
parents:
diff changeset
188
986e461dc072 Initial revision
glantau
parents:
diff changeset
189 packuswb_r2r(mm7, mm0);
986e461dc072 Initial revision
glantau
parents:
diff changeset
190 packuswb_r2r(mm7, mm1);
986e461dc072 Initial revision
glantau
parents:
diff changeset
191 packuswb_r2r(mm7, mm3);
986e461dc072 Initial revision
glantau
parents:
diff changeset
192 packuswb_r2r(mm7, mm2);
986e461dc072 Initial revision
glantau
parents:
diff changeset
193 movq_r2m(mm0, tmp);
986e461dc072 Initial revision
glantau
parents:
diff changeset
194 dst[0] = tmp.ub[0];
986e461dc072 Initial revision
glantau
parents:
diff changeset
195 movq_r2m(mm1, tmp);
986e461dc072 Initial revision
glantau
parents:
diff changeset
196 dst[1] = tmp.ub[0];
986e461dc072 Initial revision
glantau
parents:
diff changeset
197 movq_r2m(mm2, tmp);
986e461dc072 Initial revision
glantau
parents:
diff changeset
198 dst[2] = tmp.ub[0];
986e461dc072 Initial revision
glantau
parents:
diff changeset
199 movq_r2m(mm3, tmp);
986e461dc072 Initial revision
glantau
parents:
diff changeset
200 dst[3] = tmp.ub[0];
986e461dc072 Initial revision
glantau
parents:
diff changeset
201 dst += 4;
986e461dc072 Initial revision
glantau
parents:
diff changeset
202 dst_width -= 4;
986e461dc072 Initial revision
glantau
parents:
diff changeset
203 }
986e461dc072 Initial revision
glantau
parents:
diff changeset
204 while (dst_width > 0) {
986e461dc072 Initial revision
glantau
parents:
diff changeset
205 FILTER4(mm0);
986e461dc072 Initial revision
glantau
parents:
diff changeset
206 packuswb_r2r(mm7, mm0);
986e461dc072 Initial revision
glantau
parents:
diff changeset
207 movq_r2m(mm0, tmp);
986e461dc072 Initial revision
glantau
parents:
diff changeset
208 dst[0] = tmp.ub[0];
986e461dc072 Initial revision
glantau
parents:
diff changeset
209 dst++;
986e461dc072 Initial revision
glantau
parents:
diff changeset
210 dst_width--;
986e461dc072 Initial revision
glantau
parents:
diff changeset
211 }
986e461dc072 Initial revision
glantau
parents:
diff changeset
212 emms();
986e461dc072 Initial revision
glantau
parents:
diff changeset
213 }
986e461dc072 Initial revision
glantau
parents:
diff changeset
214
1488
766a2f4edbea avcodec const correctness patch by (Drew Hess <dhess at ilm dot com>)
michaelni
parents: 1106
diff changeset
215 static void v_resample4_mmx(uint8_t *dst, int dst_width, const uint8_t *src,
2979
bfabfdf9ce55 COSMETICS: tabs --> spaces, some prettyprinting
diego
parents: 2967
diff changeset
216 int wrap, int16_t *filter)
0
986e461dc072 Initial revision
glantau
parents:
diff changeset
217 {
986e461dc072 Initial revision
glantau
parents:
diff changeset
218 int sum, i, v;
1488
766a2f4edbea avcodec const correctness patch by (Drew Hess <dhess at ilm dot com>)
michaelni
parents: 1106
diff changeset
219 const uint8_t *s;
0
986e461dc072 Initial revision
glantau
parents:
diff changeset
220 mmx_t tmp;
986e461dc072 Initial revision
glantau
parents:
diff changeset
221 mmx_t coefs[4];
2967
ef2149182f1c COSMETICS: Remove all trailing whitespace.
diego
parents: 2904
diff changeset
222
0
986e461dc072 Initial revision
glantau
parents:
diff changeset
223 for(i=0;i<4;i++) {
986e461dc072 Initial revision
glantau
parents:
diff changeset
224 v = filter[i];
986e461dc072 Initial revision
glantau
parents:
diff changeset
225 coefs[i].uw[0] = v;
986e461dc072 Initial revision
glantau
parents:
diff changeset
226 coefs[i].uw[1] = v;
986e461dc072 Initial revision
glantau
parents:
diff changeset
227 coefs[i].uw[2] = v;
986e461dc072 Initial revision
glantau
parents:
diff changeset
228 coefs[i].uw[3] = v;
986e461dc072 Initial revision
glantau
parents:
diff changeset
229 }
2967
ef2149182f1c COSMETICS: Remove all trailing whitespace.
diego
parents: 2904
diff changeset
230
0
986e461dc072 Initial revision
glantau
parents:
diff changeset
231 pxor_r2r(mm7, mm7);
986e461dc072 Initial revision
glantau
parents:
diff changeset
232 s = src;
986e461dc072 Initial revision
glantau
parents:
diff changeset
233 while (dst_width >= 4) {
986e461dc072 Initial revision
glantau
parents:
diff changeset
234 movq_m2r(s[0 * wrap], mm0);
986e461dc072 Initial revision
glantau
parents:
diff changeset
235 punpcklbw_r2r(mm7, mm0);
986e461dc072 Initial revision
glantau
parents:
diff changeset
236 movq_m2r(s[1 * wrap], mm1);
986e461dc072 Initial revision
glantau
parents:
diff changeset
237 punpcklbw_r2r(mm7, mm1);
986e461dc072 Initial revision
glantau
parents:
diff changeset
238 movq_m2r(s[2 * wrap], mm2);
986e461dc072 Initial revision
glantau
parents:
diff changeset
239 punpcklbw_r2r(mm7, mm2);
986e461dc072 Initial revision
glantau
parents:
diff changeset
240 movq_m2r(s[3 * wrap], mm3);
986e461dc072 Initial revision
glantau
parents:
diff changeset
241 punpcklbw_r2r(mm7, mm3);
986e461dc072 Initial revision
glantau
parents:
diff changeset
242
986e461dc072 Initial revision
glantau
parents:
diff changeset
243 pmullw_m2r(coefs[0], mm0);
986e461dc072 Initial revision
glantau
parents:
diff changeset
244 pmullw_m2r(coefs[1], mm1);
986e461dc072 Initial revision
glantau
parents:
diff changeset
245 pmullw_m2r(coefs[2], mm2);
986e461dc072 Initial revision
glantau
parents:
diff changeset
246 pmullw_m2r(coefs[3], mm3);
986e461dc072 Initial revision
glantau
parents:
diff changeset
247
986e461dc072 Initial revision
glantau
parents:
diff changeset
248 paddw_r2r(mm1, mm0);
986e461dc072 Initial revision
glantau
parents:
diff changeset
249 paddw_r2r(mm3, mm2);
986e461dc072 Initial revision
glantau
parents:
diff changeset
250 paddw_r2r(mm2, mm0);
986e461dc072 Initial revision
glantau
parents:
diff changeset
251 psraw_i2r(FILTER_BITS, mm0);
2967
ef2149182f1c COSMETICS: Remove all trailing whitespace.
diego
parents: 2904
diff changeset
252
0
986e461dc072 Initial revision
glantau
parents:
diff changeset
253 packuswb_r2r(mm7, mm0);
986e461dc072 Initial revision
glantau
parents:
diff changeset
254 movq_r2m(mm0, tmp);
986e461dc072 Initial revision
glantau
parents:
diff changeset
255
1064
b32afefe7d33 * UINTX -> uintx_t INTX -> intx_t
kabi
parents: 1057
diff changeset
256 *(uint32_t *)dst = tmp.ud[0];
0
986e461dc072 Initial revision
glantau
parents:
diff changeset
257 dst += 4;
986e461dc072 Initial revision
glantau
parents:
diff changeset
258 s += 4;
986e461dc072 Initial revision
glantau
parents:
diff changeset
259 dst_width -= 4;
986e461dc072 Initial revision
glantau
parents:
diff changeset
260 }
986e461dc072 Initial revision
glantau
parents:
diff changeset
261 while (dst_width > 0) {
986e461dc072 Initial revision
glantau
parents:
diff changeset
262 sum = s[0 * wrap] * filter[0] +
986e461dc072 Initial revision
glantau
parents:
diff changeset
263 s[1 * wrap] * filter[1] +
986e461dc072 Initial revision
glantau
parents:
diff changeset
264 s[2 * wrap] * filter[2] +
986e461dc072 Initial revision
glantau
parents:
diff changeset
265 s[3 * wrap] * filter[3];
986e461dc072 Initial revision
glantau
parents:
diff changeset
266 sum = sum >> FILTER_BITS;
986e461dc072 Initial revision
glantau
parents:
diff changeset
267 if (sum < 0)
986e461dc072 Initial revision
glantau
parents:
diff changeset
268 sum = 0;
986e461dc072 Initial revision
glantau
parents:
diff changeset
269 else if (sum > 255)
986e461dc072 Initial revision
glantau
parents:
diff changeset
270 sum = 255;
986e461dc072 Initial revision
glantau
parents:
diff changeset
271 dst[0] = sum;
986e461dc072 Initial revision
glantau
parents:
diff changeset
272 dst++;
986e461dc072 Initial revision
glantau
parents:
diff changeset
273 s++;
986e461dc072 Initial revision
glantau
parents:
diff changeset
274 dst_width--;
986e461dc072 Initial revision
glantau
parents:
diff changeset
275 }
986e461dc072 Initial revision
glantau
parents:
diff changeset
276 emms();
986e461dc072 Initial revision
glantau
parents:
diff changeset
277 }
986e461dc072 Initial revision
glantau
parents:
diff changeset
278 #endif
986e461dc072 Initial revision
glantau
parents:
diff changeset
279
894
a408778eff87 altivec accelerated v-resample patch by (Brian Foley <bfoley at compsoc dot nuigalway dot ie>)
michaelni
parents: 653
diff changeset
280 #ifdef HAVE_ALTIVEC
2979
bfabfdf9ce55 COSMETICS: tabs --> spaces, some prettyprinting
diego
parents: 2967
diff changeset
281 typedef union {
894
a408778eff87 altivec accelerated v-resample patch by (Brian Foley <bfoley at compsoc dot nuigalway dot ie>)
michaelni
parents: 653
diff changeset
282 vector unsigned char v;
a408778eff87 altivec accelerated v-resample patch by (Brian Foley <bfoley at compsoc dot nuigalway dot ie>)
michaelni
parents: 653
diff changeset
283 unsigned char c[16];
a408778eff87 altivec accelerated v-resample patch by (Brian Foley <bfoley at compsoc dot nuigalway dot ie>)
michaelni
parents: 653
diff changeset
284 } vec_uc_t;
a408778eff87 altivec accelerated v-resample patch by (Brian Foley <bfoley at compsoc dot nuigalway dot ie>)
michaelni
parents: 653
diff changeset
285
2979
bfabfdf9ce55 COSMETICS: tabs --> spaces, some prettyprinting
diego
parents: 2967
diff changeset
286 typedef union {
894
a408778eff87 altivec accelerated v-resample patch by (Brian Foley <bfoley at compsoc dot nuigalway dot ie>)
michaelni
parents: 653
diff changeset
287 vector signed short v;
a408778eff87 altivec accelerated v-resample patch by (Brian Foley <bfoley at compsoc dot nuigalway dot ie>)
michaelni
parents: 653
diff changeset
288 signed short s[8];
a408778eff87 altivec accelerated v-resample patch by (Brian Foley <bfoley at compsoc dot nuigalway dot ie>)
michaelni
parents: 653
diff changeset
289 } vec_ss_t;
a408778eff87 altivec accelerated v-resample patch by (Brian Foley <bfoley at compsoc dot nuigalway dot ie>)
michaelni
parents: 653
diff changeset
290
1488
766a2f4edbea avcodec const correctness patch by (Drew Hess <dhess at ilm dot com>)
michaelni
parents: 1106
diff changeset
291 void v_resample16_altivec(uint8_t *dst, int dst_width, const uint8_t *src,
2979
bfabfdf9ce55 COSMETICS: tabs --> spaces, some prettyprinting
diego
parents: 2967
diff changeset
292 int wrap, int16_t *filter)
894
a408778eff87 altivec accelerated v-resample patch by (Brian Foley <bfoley at compsoc dot nuigalway dot ie>)
michaelni
parents: 653
diff changeset
293 {
a408778eff87 altivec accelerated v-resample patch by (Brian Foley <bfoley at compsoc dot nuigalway dot ie>)
michaelni
parents: 653
diff changeset
294 int sum, i;
1488
766a2f4edbea avcodec const correctness patch by (Drew Hess <dhess at ilm dot com>)
michaelni
parents: 1106
diff changeset
295 const uint8_t *s;
894
a408778eff87 altivec accelerated v-resample patch by (Brian Foley <bfoley at compsoc dot nuigalway dot ie>)
michaelni
parents: 653
diff changeset
296 vector unsigned char *tv, tmp, dstv, zero;
a408778eff87 altivec accelerated v-resample patch by (Brian Foley <bfoley at compsoc dot nuigalway dot ie>)
michaelni
parents: 653
diff changeset
297 vec_ss_t srchv[4], srclv[4], fv[4];
2967
ef2149182f1c COSMETICS: Remove all trailing whitespace.
diego
parents: 2904
diff changeset
298 vector signed short zeros, sumhv, sumlv;
894
a408778eff87 altivec accelerated v-resample patch by (Brian Foley <bfoley at compsoc dot nuigalway dot ie>)
michaelni
parents: 653
diff changeset
299 s = src;
a408778eff87 altivec accelerated v-resample patch by (Brian Foley <bfoley at compsoc dot nuigalway dot ie>)
michaelni
parents: 653
diff changeset
300
a408778eff87 altivec accelerated v-resample patch by (Brian Foley <bfoley at compsoc dot nuigalway dot ie>)
michaelni
parents: 653
diff changeset
301 for(i=0;i<4;i++)
a408778eff87 altivec accelerated v-resample patch by (Brian Foley <bfoley at compsoc dot nuigalway dot ie>)
michaelni
parents: 653
diff changeset
302 {
a408778eff87 altivec accelerated v-resample patch by (Brian Foley <bfoley at compsoc dot nuigalway dot ie>)
michaelni
parents: 653
diff changeset
303 /*
a408778eff87 altivec accelerated v-resample patch by (Brian Foley <bfoley at compsoc dot nuigalway dot ie>)
michaelni
parents: 653
diff changeset
304 The vec_madds later on does an implicit >>15 on the result.
a408778eff87 altivec accelerated v-resample patch by (Brian Foley <bfoley at compsoc dot nuigalway dot ie>)
michaelni
parents: 653
diff changeset
305 Since FILTER_BITS is 8, and we have 15 bits of magnitude in
a408778eff87 altivec accelerated v-resample patch by (Brian Foley <bfoley at compsoc dot nuigalway dot ie>)
michaelni
parents: 653
diff changeset
306 a signed short, we have just enough bits to pre-shift our
a408778eff87 altivec accelerated v-resample patch by (Brian Foley <bfoley at compsoc dot nuigalway dot ie>)
michaelni
parents: 653
diff changeset
307 filter constants <<7 to compensate for vec_madds.
a408778eff87 altivec accelerated v-resample patch by (Brian Foley <bfoley at compsoc dot nuigalway dot ie>)
michaelni
parents: 653
diff changeset
308 */
a408778eff87 altivec accelerated v-resample patch by (Brian Foley <bfoley at compsoc dot nuigalway dot ie>)
michaelni
parents: 653
diff changeset
309 fv[i].s[0] = filter[i] << (15-FILTER_BITS);
a408778eff87 altivec accelerated v-resample patch by (Brian Foley <bfoley at compsoc dot nuigalway dot ie>)
michaelni
parents: 653
diff changeset
310 fv[i].v = vec_splat(fv[i].v, 0);
a408778eff87 altivec accelerated v-resample patch by (Brian Foley <bfoley at compsoc dot nuigalway dot ie>)
michaelni
parents: 653
diff changeset
311 }
2967
ef2149182f1c COSMETICS: Remove all trailing whitespace.
diego
parents: 2904
diff changeset
312
894
a408778eff87 altivec accelerated v-resample patch by (Brian Foley <bfoley at compsoc dot nuigalway dot ie>)
michaelni
parents: 653
diff changeset
313 zero = vec_splat_u8(0);
a408778eff87 altivec accelerated v-resample patch by (Brian Foley <bfoley at compsoc dot nuigalway dot ie>)
michaelni
parents: 653
diff changeset
314 zeros = vec_splat_s16(0);
a408778eff87 altivec accelerated v-resample patch by (Brian Foley <bfoley at compsoc dot nuigalway dot ie>)
michaelni
parents: 653
diff changeset
315
a408778eff87 altivec accelerated v-resample patch by (Brian Foley <bfoley at compsoc dot nuigalway dot ie>)
michaelni
parents: 653
diff changeset
316
a408778eff87 altivec accelerated v-resample patch by (Brian Foley <bfoley at compsoc dot nuigalway dot ie>)
michaelni
parents: 653
diff changeset
317 /*
a408778eff87 altivec accelerated v-resample patch by (Brian Foley <bfoley at compsoc dot nuigalway dot ie>)
michaelni
parents: 653
diff changeset
318 When we're resampling, we'd ideally like both our input buffers,
a408778eff87 altivec accelerated v-resample patch by (Brian Foley <bfoley at compsoc dot nuigalway dot ie>)
michaelni
parents: 653
diff changeset
319 and output buffers to be 16-byte aligned, so we can do both aligned
a408778eff87 altivec accelerated v-resample patch by (Brian Foley <bfoley at compsoc dot nuigalway dot ie>)
michaelni
parents: 653
diff changeset
320 reads and writes. Sadly we can't always have this at the moment, so
a408778eff87 altivec accelerated v-resample patch by (Brian Foley <bfoley at compsoc dot nuigalway dot ie>)
michaelni
parents: 653
diff changeset
321 we opt for aligned writes, as unaligned writes have a huge overhead.
a408778eff87 altivec accelerated v-resample patch by (Brian Foley <bfoley at compsoc dot nuigalway dot ie>)
michaelni
parents: 653
diff changeset
322 To do this, do enough scalar resamples to get dst 16-byte aligned.
a408778eff87 altivec accelerated v-resample patch by (Brian Foley <bfoley at compsoc dot nuigalway dot ie>)
michaelni
parents: 653
diff changeset
323 */
898
6d5e3fe7aea1 Simplify an expression and eliminate a compile warning
philipjsg
parents: 894
diff changeset
324 i = (-(int)dst) & 0xf;
894
a408778eff87 altivec accelerated v-resample patch by (Brian Foley <bfoley at compsoc dot nuigalway dot ie>)
michaelni
parents: 653
diff changeset
325 while(i>0) {
a408778eff87 altivec accelerated v-resample patch by (Brian Foley <bfoley at compsoc dot nuigalway dot ie>)
michaelni
parents: 653
diff changeset
326 sum = s[0 * wrap] * filter[0] +
a408778eff87 altivec accelerated v-resample patch by (Brian Foley <bfoley at compsoc dot nuigalway dot ie>)
michaelni
parents: 653
diff changeset
327 s[1 * wrap] * filter[1] +
a408778eff87 altivec accelerated v-resample patch by (Brian Foley <bfoley at compsoc dot nuigalway dot ie>)
michaelni
parents: 653
diff changeset
328 s[2 * wrap] * filter[2] +
a408778eff87 altivec accelerated v-resample patch by (Brian Foley <bfoley at compsoc dot nuigalway dot ie>)
michaelni
parents: 653
diff changeset
329 s[3 * wrap] * filter[3];
a408778eff87 altivec accelerated v-resample patch by (Brian Foley <bfoley at compsoc dot nuigalway dot ie>)
michaelni
parents: 653
diff changeset
330 sum = sum >> FILTER_BITS;
a408778eff87 altivec accelerated v-resample patch by (Brian Foley <bfoley at compsoc dot nuigalway dot ie>)
michaelni
parents: 653
diff changeset
331 if (sum<0) sum = 0; else if (sum>255) sum=255;
a408778eff87 altivec accelerated v-resample patch by (Brian Foley <bfoley at compsoc dot nuigalway dot ie>)
michaelni
parents: 653
diff changeset
332 dst[0] = sum;
a408778eff87 altivec accelerated v-resample patch by (Brian Foley <bfoley at compsoc dot nuigalway dot ie>)
michaelni
parents: 653
diff changeset
333 dst++;
a408778eff87 altivec accelerated v-resample patch by (Brian Foley <bfoley at compsoc dot nuigalway dot ie>)
michaelni
parents: 653
diff changeset
334 s++;
a408778eff87 altivec accelerated v-resample patch by (Brian Foley <bfoley at compsoc dot nuigalway dot ie>)
michaelni
parents: 653
diff changeset
335 dst_width--;
a408778eff87 altivec accelerated v-resample patch by (Brian Foley <bfoley at compsoc dot nuigalway dot ie>)
michaelni
parents: 653
diff changeset
336 i--;
a408778eff87 altivec accelerated v-resample patch by (Brian Foley <bfoley at compsoc dot nuigalway dot ie>)
michaelni
parents: 653
diff changeset
337 }
2967
ef2149182f1c COSMETICS: Remove all trailing whitespace.
diego
parents: 2904
diff changeset
338
894
a408778eff87 altivec accelerated v-resample patch by (Brian Foley <bfoley at compsoc dot nuigalway dot ie>)
michaelni
parents: 653
diff changeset
339 /* Do our altivec resampling on 16 pixels at once. */
a408778eff87 altivec accelerated v-resample patch by (Brian Foley <bfoley at compsoc dot nuigalway dot ie>)
michaelni
parents: 653
diff changeset
340 while(dst_width>=16) {
a408778eff87 altivec accelerated v-resample patch by (Brian Foley <bfoley at compsoc dot nuigalway dot ie>)
michaelni
parents: 653
diff changeset
341 /*
a408778eff87 altivec accelerated v-resample patch by (Brian Foley <bfoley at compsoc dot nuigalway dot ie>)
michaelni
parents: 653
diff changeset
342 Read 16 (potentially unaligned) bytes from each of
a408778eff87 altivec accelerated v-resample patch by (Brian Foley <bfoley at compsoc dot nuigalway dot ie>)
michaelni
parents: 653
diff changeset
343 4 lines into 4 vectors, and split them into shorts.
a408778eff87 altivec accelerated v-resample patch by (Brian Foley <bfoley at compsoc dot nuigalway dot ie>)
michaelni
parents: 653
diff changeset
344 Interleave the multipy/accumulate for the resample
a408778eff87 altivec accelerated v-resample patch by (Brian Foley <bfoley at compsoc dot nuigalway dot ie>)
michaelni
parents: 653
diff changeset
345 filter with the loads to hide the 3 cycle latency
a408778eff87 altivec accelerated v-resample patch by (Brian Foley <bfoley at compsoc dot nuigalway dot ie>)
michaelni
parents: 653
diff changeset
346 the vec_madds have.
a408778eff87 altivec accelerated v-resample patch by (Brian Foley <bfoley at compsoc dot nuigalway dot ie>)
michaelni
parents: 653
diff changeset
347 */
a408778eff87 altivec accelerated v-resample patch by (Brian Foley <bfoley at compsoc dot nuigalway dot ie>)
michaelni
parents: 653
diff changeset
348 tv = (vector unsigned char *) &s[0 * wrap];
a408778eff87 altivec accelerated v-resample patch by (Brian Foley <bfoley at compsoc dot nuigalway dot ie>)
michaelni
parents: 653
diff changeset
349 tmp = vec_perm(tv[0], tv[1], vec_lvsl(0, &s[i * wrap]));
a408778eff87 altivec accelerated v-resample patch by (Brian Foley <bfoley at compsoc dot nuigalway dot ie>)
michaelni
parents: 653
diff changeset
350 srchv[0].v = (vector signed short) vec_mergeh(zero, tmp);
a408778eff87 altivec accelerated v-resample patch by (Brian Foley <bfoley at compsoc dot nuigalway dot ie>)
michaelni
parents: 653
diff changeset
351 srclv[0].v = (vector signed short) vec_mergel(zero, tmp);
a408778eff87 altivec accelerated v-resample patch by (Brian Foley <bfoley at compsoc dot nuigalway dot ie>)
michaelni
parents: 653
diff changeset
352 sumhv = vec_madds(srchv[0].v, fv[0].v, zeros);
a408778eff87 altivec accelerated v-resample patch by (Brian Foley <bfoley at compsoc dot nuigalway dot ie>)
michaelni
parents: 653
diff changeset
353 sumlv = vec_madds(srclv[0].v, fv[0].v, zeros);
a408778eff87 altivec accelerated v-resample patch by (Brian Foley <bfoley at compsoc dot nuigalway dot ie>)
michaelni
parents: 653
diff changeset
354
a408778eff87 altivec accelerated v-resample patch by (Brian Foley <bfoley at compsoc dot nuigalway dot ie>)
michaelni
parents: 653
diff changeset
355 tv = (vector unsigned char *) &s[1 * wrap];
a408778eff87 altivec accelerated v-resample patch by (Brian Foley <bfoley at compsoc dot nuigalway dot ie>)
michaelni
parents: 653
diff changeset
356 tmp = vec_perm(tv[0], tv[1], vec_lvsl(0, &s[1 * wrap]));
a408778eff87 altivec accelerated v-resample patch by (Brian Foley <bfoley at compsoc dot nuigalway dot ie>)
michaelni
parents: 653
diff changeset
357 srchv[1].v = (vector signed short) vec_mergeh(zero, tmp);
a408778eff87 altivec accelerated v-resample patch by (Brian Foley <bfoley at compsoc dot nuigalway dot ie>)
michaelni
parents: 653
diff changeset
358 srclv[1].v = (vector signed short) vec_mergel(zero, tmp);
a408778eff87 altivec accelerated v-resample patch by (Brian Foley <bfoley at compsoc dot nuigalway dot ie>)
michaelni
parents: 653
diff changeset
359 sumhv = vec_madds(srchv[1].v, fv[1].v, sumhv);
a408778eff87 altivec accelerated v-resample patch by (Brian Foley <bfoley at compsoc dot nuigalway dot ie>)
michaelni
parents: 653
diff changeset
360 sumlv = vec_madds(srclv[1].v, fv[1].v, sumlv);
a408778eff87 altivec accelerated v-resample patch by (Brian Foley <bfoley at compsoc dot nuigalway dot ie>)
michaelni
parents: 653
diff changeset
361
a408778eff87 altivec accelerated v-resample patch by (Brian Foley <bfoley at compsoc dot nuigalway dot ie>)
michaelni
parents: 653
diff changeset
362 tv = (vector unsigned char *) &s[2 * wrap];
a408778eff87 altivec accelerated v-resample patch by (Brian Foley <bfoley at compsoc dot nuigalway dot ie>)
michaelni
parents: 653
diff changeset
363 tmp = vec_perm(tv[0], tv[1], vec_lvsl(0, &s[2 * wrap]));
a408778eff87 altivec accelerated v-resample patch by (Brian Foley <bfoley at compsoc dot nuigalway dot ie>)
michaelni
parents: 653
diff changeset
364 srchv[2].v = (vector signed short) vec_mergeh(zero, tmp);
a408778eff87 altivec accelerated v-resample patch by (Brian Foley <bfoley at compsoc dot nuigalway dot ie>)
michaelni
parents: 653
diff changeset
365 srclv[2].v = (vector signed short) vec_mergel(zero, tmp);
a408778eff87 altivec accelerated v-resample patch by (Brian Foley <bfoley at compsoc dot nuigalway dot ie>)
michaelni
parents: 653
diff changeset
366 sumhv = vec_madds(srchv[2].v, fv[2].v, sumhv);
a408778eff87 altivec accelerated v-resample patch by (Brian Foley <bfoley at compsoc dot nuigalway dot ie>)
michaelni
parents: 653
diff changeset
367 sumlv = vec_madds(srclv[2].v, fv[2].v, sumlv);
a408778eff87 altivec accelerated v-resample patch by (Brian Foley <bfoley at compsoc dot nuigalway dot ie>)
michaelni
parents: 653
diff changeset
368
a408778eff87 altivec accelerated v-resample patch by (Brian Foley <bfoley at compsoc dot nuigalway dot ie>)
michaelni
parents: 653
diff changeset
369 tv = (vector unsigned char *) &s[3 * wrap];
a408778eff87 altivec accelerated v-resample patch by (Brian Foley <bfoley at compsoc dot nuigalway dot ie>)
michaelni
parents: 653
diff changeset
370 tmp = vec_perm(tv[0], tv[1], vec_lvsl(0, &s[3 * wrap]));
a408778eff87 altivec accelerated v-resample patch by (Brian Foley <bfoley at compsoc dot nuigalway dot ie>)
michaelni
parents: 653
diff changeset
371 srchv[3].v = (vector signed short) vec_mergeh(zero, tmp);
a408778eff87 altivec accelerated v-resample patch by (Brian Foley <bfoley at compsoc dot nuigalway dot ie>)
michaelni
parents: 653
diff changeset
372 srclv[3].v = (vector signed short) vec_mergel(zero, tmp);
a408778eff87 altivec accelerated v-resample patch by (Brian Foley <bfoley at compsoc dot nuigalway dot ie>)
michaelni
parents: 653
diff changeset
373 sumhv = vec_madds(srchv[3].v, fv[3].v, sumhv);
a408778eff87 altivec accelerated v-resample patch by (Brian Foley <bfoley at compsoc dot nuigalway dot ie>)
michaelni
parents: 653
diff changeset
374 sumlv = vec_madds(srclv[3].v, fv[3].v, sumlv);
2967
ef2149182f1c COSMETICS: Remove all trailing whitespace.
diego
parents: 2904
diff changeset
375
894
a408778eff87 altivec accelerated v-resample patch by (Brian Foley <bfoley at compsoc dot nuigalway dot ie>)
michaelni
parents: 653
diff changeset
376 /*
a408778eff87 altivec accelerated v-resample patch by (Brian Foley <bfoley at compsoc dot nuigalway dot ie>)
michaelni
parents: 653
diff changeset
377 Pack the results into our destination vector,
a408778eff87 altivec accelerated v-resample patch by (Brian Foley <bfoley at compsoc dot nuigalway dot ie>)
michaelni
parents: 653
diff changeset
378 and do an aligned write of that back to memory.
a408778eff87 altivec accelerated v-resample patch by (Brian Foley <bfoley at compsoc dot nuigalway dot ie>)
michaelni
parents: 653
diff changeset
379 */
a408778eff87 altivec accelerated v-resample patch by (Brian Foley <bfoley at compsoc dot nuigalway dot ie>)
michaelni
parents: 653
diff changeset
380 dstv = vec_packsu(sumhv, sumlv) ;
a408778eff87 altivec accelerated v-resample patch by (Brian Foley <bfoley at compsoc dot nuigalway dot ie>)
michaelni
parents: 653
diff changeset
381 vec_st(dstv, 0, (vector unsigned char *) dst);
2967
ef2149182f1c COSMETICS: Remove all trailing whitespace.
diego
parents: 2904
diff changeset
382
894
a408778eff87 altivec accelerated v-resample patch by (Brian Foley <bfoley at compsoc dot nuigalway dot ie>)
michaelni
parents: 653
diff changeset
383 dst+=16;
a408778eff87 altivec accelerated v-resample patch by (Brian Foley <bfoley at compsoc dot nuigalway dot ie>)
michaelni
parents: 653
diff changeset
384 s+=16;
a408778eff87 altivec accelerated v-resample patch by (Brian Foley <bfoley at compsoc dot nuigalway dot ie>)
michaelni
parents: 653
diff changeset
385 dst_width-=16;
a408778eff87 altivec accelerated v-resample patch by (Brian Foley <bfoley at compsoc dot nuigalway dot ie>)
michaelni
parents: 653
diff changeset
386 }
a408778eff87 altivec accelerated v-resample patch by (Brian Foley <bfoley at compsoc dot nuigalway dot ie>)
michaelni
parents: 653
diff changeset
387
a408778eff87 altivec accelerated v-resample patch by (Brian Foley <bfoley at compsoc dot nuigalway dot ie>)
michaelni
parents: 653
diff changeset
388 /*
a408778eff87 altivec accelerated v-resample patch by (Brian Foley <bfoley at compsoc dot nuigalway dot ie>)
michaelni
parents: 653
diff changeset
389 If there are any leftover pixels, resample them
a408778eff87 altivec accelerated v-resample patch by (Brian Foley <bfoley at compsoc dot nuigalway dot ie>)
michaelni
parents: 653
diff changeset
390 with the slow scalar method.
a408778eff87 altivec accelerated v-resample patch by (Brian Foley <bfoley at compsoc dot nuigalway dot ie>)
michaelni
parents: 653
diff changeset
391 */
a408778eff87 altivec accelerated v-resample patch by (Brian Foley <bfoley at compsoc dot nuigalway dot ie>)
michaelni
parents: 653
diff changeset
392 while(dst_width>0) {
a408778eff87 altivec accelerated v-resample patch by (Brian Foley <bfoley at compsoc dot nuigalway dot ie>)
michaelni
parents: 653
diff changeset
393 sum = s[0 * wrap] * filter[0] +
a408778eff87 altivec accelerated v-resample patch by (Brian Foley <bfoley at compsoc dot nuigalway dot ie>)
michaelni
parents: 653
diff changeset
394 s[1 * wrap] * filter[1] +
a408778eff87 altivec accelerated v-resample patch by (Brian Foley <bfoley at compsoc dot nuigalway dot ie>)
michaelni
parents: 653
diff changeset
395 s[2 * wrap] * filter[2] +
a408778eff87 altivec accelerated v-resample patch by (Brian Foley <bfoley at compsoc dot nuigalway dot ie>)
michaelni
parents: 653
diff changeset
396 s[3 * wrap] * filter[3];
a408778eff87 altivec accelerated v-resample patch by (Brian Foley <bfoley at compsoc dot nuigalway dot ie>)
michaelni
parents: 653
diff changeset
397 sum = sum >> FILTER_BITS;
a408778eff87 altivec accelerated v-resample patch by (Brian Foley <bfoley at compsoc dot nuigalway dot ie>)
michaelni
parents: 653
diff changeset
398 if (sum<0) sum = 0; else if (sum>255) sum=255;
a408778eff87 altivec accelerated v-resample patch by (Brian Foley <bfoley at compsoc dot nuigalway dot ie>)
michaelni
parents: 653
diff changeset
399 dst[0] = sum;
a408778eff87 altivec accelerated v-resample patch by (Brian Foley <bfoley at compsoc dot nuigalway dot ie>)
michaelni
parents: 653
diff changeset
400 dst++;
a408778eff87 altivec accelerated v-resample patch by (Brian Foley <bfoley at compsoc dot nuigalway dot ie>)
michaelni
parents: 653
diff changeset
401 s++;
a408778eff87 altivec accelerated v-resample patch by (Brian Foley <bfoley at compsoc dot nuigalway dot ie>)
michaelni
parents: 653
diff changeset
402 dst_width--;
a408778eff87 altivec accelerated v-resample patch by (Brian Foley <bfoley at compsoc dot nuigalway dot ie>)
michaelni
parents: 653
diff changeset
403 }
a408778eff87 altivec accelerated v-resample patch by (Brian Foley <bfoley at compsoc dot nuigalway dot ie>)
michaelni
parents: 653
diff changeset
404 }
a408778eff87 altivec accelerated v-resample patch by (Brian Foley <bfoley at compsoc dot nuigalway dot ie>)
michaelni
parents: 653
diff changeset
405 #endif
a408778eff87 altivec accelerated v-resample patch by (Brian Foley <bfoley at compsoc dot nuigalway dot ie>)
michaelni
parents: 653
diff changeset
406
0
986e461dc072 Initial revision
glantau
parents:
diff changeset
407 /* slow version to handle limit cases. Does not need optimisation */
1488
766a2f4edbea avcodec const correctness patch by (Drew Hess <dhess at ilm dot com>)
michaelni
parents: 1106
diff changeset
408 static void h_resample_slow(uint8_t *dst, int dst_width,
2979
bfabfdf9ce55 COSMETICS: tabs --> spaces, some prettyprinting
diego
parents: 2967
diff changeset
409 const uint8_t *src, int src_width,
1064
b32afefe7d33 * UINTX -> uintx_t INTX -> intx_t
kabi
parents: 1057
diff changeset
410 int src_start, int src_incr, int16_t *filters)
0
986e461dc072 Initial revision
glantau
parents:
diff changeset
411 {
986e461dc072 Initial revision
glantau
parents:
diff changeset
412 int src_pos, phase, sum, j, v, i;
1488
766a2f4edbea avcodec const correctness patch by (Drew Hess <dhess at ilm dot com>)
michaelni
parents: 1106
diff changeset
413 const uint8_t *s, *src_end;
1064
b32afefe7d33 * UINTX -> uintx_t INTX -> intx_t
kabi
parents: 1057
diff changeset
414 int16_t *filter;
0
986e461dc072 Initial revision
glantau
parents:
diff changeset
415
986e461dc072 Initial revision
glantau
parents:
diff changeset
416 src_end = src + src_width;
986e461dc072 Initial revision
glantau
parents:
diff changeset
417 src_pos = src_start;
986e461dc072 Initial revision
glantau
parents:
diff changeset
418 for(i=0;i<dst_width;i++) {
986e461dc072 Initial revision
glantau
parents:
diff changeset
419 s = src + (src_pos >> POS_FRAC_BITS);
986e461dc072 Initial revision
glantau
parents:
diff changeset
420 phase = get_phase(src_pos);
986e461dc072 Initial revision
glantau
parents:
diff changeset
421 filter = filters + phase * NB_TAPS;
986e461dc072 Initial revision
glantau
parents:
diff changeset
422 sum = 0;
986e461dc072 Initial revision
glantau
parents:
diff changeset
423 for(j=0;j<NB_TAPS;j++) {
986e461dc072 Initial revision
glantau
parents:
diff changeset
424 if (s < src)
986e461dc072 Initial revision
glantau
parents:
diff changeset
425 v = src[0];
986e461dc072 Initial revision
glantau
parents:
diff changeset
426 else if (s >= src_end)
986e461dc072 Initial revision
glantau
parents:
diff changeset
427 v = src_end[-1];
986e461dc072 Initial revision
glantau
parents:
diff changeset
428 else
986e461dc072 Initial revision
glantau
parents:
diff changeset
429 v = s[0];
986e461dc072 Initial revision
glantau
parents:
diff changeset
430 sum += v * filter[j];
986e461dc072 Initial revision
glantau
parents:
diff changeset
431 s++;
986e461dc072 Initial revision
glantau
parents:
diff changeset
432 }
986e461dc072 Initial revision
glantau
parents:
diff changeset
433 sum = sum >> FILTER_BITS;
986e461dc072 Initial revision
glantau
parents:
diff changeset
434 if (sum < 0)
986e461dc072 Initial revision
glantau
parents:
diff changeset
435 sum = 0;
986e461dc072 Initial revision
glantau
parents:
diff changeset
436 else if (sum > 255)
986e461dc072 Initial revision
glantau
parents:
diff changeset
437 sum = 255;
986e461dc072 Initial revision
glantau
parents:
diff changeset
438 dst[0] = sum;
986e461dc072 Initial revision
glantau
parents:
diff changeset
439 src_pos += src_incr;
986e461dc072 Initial revision
glantau
parents:
diff changeset
440 dst++;
986e461dc072 Initial revision
glantau
parents:
diff changeset
441 }
986e461dc072 Initial revision
glantau
parents:
diff changeset
442 }
986e461dc072 Initial revision
glantau
parents:
diff changeset
443
1488
766a2f4edbea avcodec const correctness patch by (Drew Hess <dhess at ilm dot com>)
michaelni
parents: 1106
diff changeset
444 static void h_resample(uint8_t *dst, int dst_width, const uint8_t *src,
2979
bfabfdf9ce55 COSMETICS: tabs --> spaces, some prettyprinting
diego
parents: 2967
diff changeset
445 int src_width, int src_start, int src_incr,
bfabfdf9ce55 COSMETICS: tabs --> spaces, some prettyprinting
diego
parents: 2967
diff changeset
446 int16_t *filters)
0
986e461dc072 Initial revision
glantau
parents:
diff changeset
447 {
986e461dc072 Initial revision
glantau
parents:
diff changeset
448 int n, src_end;
986e461dc072 Initial revision
glantau
parents:
diff changeset
449
986e461dc072 Initial revision
glantau
parents:
diff changeset
450 if (src_start < 0) {
986e461dc072 Initial revision
glantau
parents:
diff changeset
451 n = (0 - src_start + src_incr - 1) / src_incr;
986e461dc072 Initial revision
glantau
parents:
diff changeset
452 h_resample_slow(dst, n, src, src_width, src_start, src_incr, filters);
986e461dc072 Initial revision
glantau
parents:
diff changeset
453 dst += n;
986e461dc072 Initial revision
glantau
parents:
diff changeset
454 dst_width -= n;
986e461dc072 Initial revision
glantau
parents:
diff changeset
455 src_start += n * src_incr;
986e461dc072 Initial revision
glantau
parents:
diff changeset
456 }
986e461dc072 Initial revision
glantau
parents:
diff changeset
457 src_end = src_start + dst_width * src_incr;
986e461dc072 Initial revision
glantau
parents:
diff changeset
458 if (src_end > ((src_width - NB_TAPS) << POS_FRAC_BITS)) {
2967
ef2149182f1c COSMETICS: Remove all trailing whitespace.
diego
parents: 2904
diff changeset
459 n = (((src_width - NB_TAPS + 1) << POS_FRAC_BITS) - 1 - src_start) /
0
986e461dc072 Initial revision
glantau
parents:
diff changeset
460 src_incr;
986e461dc072 Initial revision
glantau
parents:
diff changeset
461 } else {
986e461dc072 Initial revision
glantau
parents:
diff changeset
462 n = dst_width;
986e461dc072 Initial revision
glantau
parents:
diff changeset
463 }
2
2e2c46c87460 fixed config for direct mplayer build compatibility
glantau
parents: 0
diff changeset
464 #ifdef HAVE_MMX
0
986e461dc072 Initial revision
glantau
parents:
diff changeset
465 if ((mm_flags & MM_MMX) && NB_TAPS == 4)
2967
ef2149182f1c COSMETICS: Remove all trailing whitespace.
diego
parents: 2904
diff changeset
466 h_resample_fast4_mmx(dst, n,
0
986e461dc072 Initial revision
glantau
parents:
diff changeset
467 src, src_width, src_start, src_incr, filters);
986e461dc072 Initial revision
glantau
parents:
diff changeset
468 else
986e461dc072 Initial revision
glantau
parents:
diff changeset
469 #endif
2967
ef2149182f1c COSMETICS: Remove all trailing whitespace.
diego
parents: 2904
diff changeset
470 h_resample_fast(dst, n,
0
986e461dc072 Initial revision
glantau
parents:
diff changeset
471 src, src_width, src_start, src_incr, filters);
986e461dc072 Initial revision
glantau
parents:
diff changeset
472 if (n < dst_width) {
986e461dc072 Initial revision
glantau
parents:
diff changeset
473 dst += n;
986e461dc072 Initial revision
glantau
parents:
diff changeset
474 dst_width -= n;
986e461dc072 Initial revision
glantau
parents:
diff changeset
475 src_start += n * src_incr;
2967
ef2149182f1c COSMETICS: Remove all trailing whitespace.
diego
parents: 2904
diff changeset
476 h_resample_slow(dst, dst_width,
0
986e461dc072 Initial revision
glantau
parents:
diff changeset
477 src, src_width, src_start, src_incr, filters);
986e461dc072 Initial revision
glantau
parents:
diff changeset
478 }
986e461dc072 Initial revision
glantau
parents:
diff changeset
479 }
986e461dc072 Initial revision
glantau
parents:
diff changeset
480
2967
ef2149182f1c COSMETICS: Remove all trailing whitespace.
diego
parents: 2904
diff changeset
481 static void component_resample(ImgReSampleContext *s,
1064
b32afefe7d33 * UINTX -> uintx_t INTX -> intx_t
kabi
parents: 1057
diff changeset
482 uint8_t *output, int owrap, int owidth, int oheight,
b32afefe7d33 * UINTX -> uintx_t INTX -> intx_t
kabi
parents: 1057
diff changeset
483 uint8_t *input, int iwrap, int iwidth, int iheight)
0
986e461dc072 Initial revision
glantau
parents:
diff changeset
484 {
986e461dc072 Initial revision
glantau
parents:
diff changeset
485 int src_y, src_y1, last_src_y, ring_y, phase_y, y1, y;
1064
b32afefe7d33 * UINTX -> uintx_t INTX -> intx_t
kabi
parents: 1057
diff changeset
486 uint8_t *new_line, *src_line;
0
986e461dc072 Initial revision
glantau
parents:
diff changeset
487
986e461dc072 Initial revision
glantau
parents:
diff changeset
488 last_src_y = - FCENTER - 1;
986e461dc072 Initial revision
glantau
parents:
diff changeset
489 /* position of the bottom of the filter in the source image */
2967
ef2149182f1c COSMETICS: Remove all trailing whitespace.
diego
parents: 2904
diff changeset
490 src_y = (last_src_y + NB_TAPS) * POS_FRAC;
0
986e461dc072 Initial revision
glantau
parents:
diff changeset
491 ring_y = NB_TAPS; /* position in ring buffer */
986e461dc072 Initial revision
glantau
parents:
diff changeset
492 for(y=0;y<oheight;y++) {
986e461dc072 Initial revision
glantau
parents:
diff changeset
493 /* apply horizontal filter on new lines from input if needed */
986e461dc072 Initial revision
glantau
parents:
diff changeset
494 src_y1 = src_y >> POS_FRAC_BITS;
986e461dc072 Initial revision
glantau
parents:
diff changeset
495 while (last_src_y < src_y1) {
986e461dc072 Initial revision
glantau
parents:
diff changeset
496 if (++ring_y >= LINE_BUF_HEIGHT + NB_TAPS)
986e461dc072 Initial revision
glantau
parents:
diff changeset
497 ring_y = NB_TAPS;
986e461dc072 Initial revision
glantau
parents:
diff changeset
498 last_src_y++;
630
b4ee42142ad1 croping patch by (talus25 at speakeasy dot net) with fixes from atmos & me
michaelni
parents: 429
diff changeset
499 /* handle limit conditions : replicate line (slightly
b4ee42142ad1 croping patch by (talus25 at speakeasy dot net) with fixes from atmos & me
michaelni
parents: 429
diff changeset
500 inefficient because we filter multiple times) */
0
986e461dc072 Initial revision
glantau
parents:
diff changeset
501 y1 = last_src_y;
986e461dc072 Initial revision
glantau
parents:
diff changeset
502 if (y1 < 0) {
986e461dc072 Initial revision
glantau
parents:
diff changeset
503 y1 = 0;
986e461dc072 Initial revision
glantau
parents:
diff changeset
504 } else if (y1 >= iheight) {
986e461dc072 Initial revision
glantau
parents:
diff changeset
505 y1 = iheight - 1;
986e461dc072 Initial revision
glantau
parents:
diff changeset
506 }
986e461dc072 Initial revision
glantau
parents:
diff changeset
507 src_line = input + y1 * iwrap;
986e461dc072 Initial revision
glantau
parents:
diff changeset
508 new_line = s->line_buf + ring_y * owidth;
986e461dc072 Initial revision
glantau
parents:
diff changeset
509 /* apply filter and handle limit cases correctly */
2967
ef2149182f1c COSMETICS: Remove all trailing whitespace.
diego
parents: 2904
diff changeset
510 h_resample(new_line, owidth,
ef2149182f1c COSMETICS: Remove all trailing whitespace.
diego
parents: 2904
diff changeset
511 src_line, iwidth, - FCENTER * POS_FRAC, s->h_incr,
0
986e461dc072 Initial revision
glantau
parents:
diff changeset
512 &s->h_filters[0][0]);
986e461dc072 Initial revision
glantau
parents:
diff changeset
513 /* handle ring buffer wraping */
986e461dc072 Initial revision
glantau
parents:
diff changeset
514 if (ring_y >= LINE_BUF_HEIGHT) {
986e461dc072 Initial revision
glantau
parents:
diff changeset
515 memcpy(s->line_buf + (ring_y - LINE_BUF_HEIGHT) * owidth,
986e461dc072 Initial revision
glantau
parents:
diff changeset
516 new_line, owidth);
986e461dc072 Initial revision
glantau
parents:
diff changeset
517 }
986e461dc072 Initial revision
glantau
parents:
diff changeset
518 }
986e461dc072 Initial revision
glantau
parents:
diff changeset
519 /* apply vertical filter */
986e461dc072 Initial revision
glantau
parents:
diff changeset
520 phase_y = get_phase(src_y);
2
2e2c46c87460 fixed config for direct mplayer build compatibility
glantau
parents: 0
diff changeset
521 #ifdef HAVE_MMX
0
986e461dc072 Initial revision
glantau
parents:
diff changeset
522 /* desactivated MMX because loss of precision */
986e461dc072 Initial revision
glantau
parents:
diff changeset
523 if ((mm_flags & MM_MMX) && NB_TAPS == 4 && 0)
2967
ef2149182f1c COSMETICS: Remove all trailing whitespace.
diego
parents: 2904
diff changeset
524 v_resample4_mmx(output, owidth,
ef2149182f1c COSMETICS: Remove all trailing whitespace.
diego
parents: 2904
diff changeset
525 s->line_buf + (ring_y - NB_TAPS + 1) * owidth, owidth,
0
986e461dc072 Initial revision
glantau
parents:
diff changeset
526 &s->v_filters[phase_y][0]);
986e461dc072 Initial revision
glantau
parents:
diff changeset
527 else
986e461dc072 Initial revision
glantau
parents:
diff changeset
528 #endif
894
a408778eff87 altivec accelerated v-resample patch by (Brian Foley <bfoley at compsoc dot nuigalway dot ie>)
michaelni
parents: 653
diff changeset
529 #ifdef HAVE_ALTIVEC
920
a0ad8e3452f2 practically disabling altivec resampling code (some ppl said its broken) patch by (Dieter Shirley <dieters at schemasoft dot com>)
michaelni
parents: 898
diff changeset
530 if ((mm_flags & MM_ALTIVEC) && NB_TAPS == 4 && FILTER_BITS <= 6)
894
a408778eff87 altivec accelerated v-resample patch by (Brian Foley <bfoley at compsoc dot nuigalway dot ie>)
michaelni
parents: 653
diff changeset
531 v_resample16_altivec(output, owidth,
a408778eff87 altivec accelerated v-resample patch by (Brian Foley <bfoley at compsoc dot nuigalway dot ie>)
michaelni
parents: 653
diff changeset
532 s->line_buf + (ring_y - NB_TAPS + 1) * owidth, owidth,
a408778eff87 altivec accelerated v-resample patch by (Brian Foley <bfoley at compsoc dot nuigalway dot ie>)
michaelni
parents: 653
diff changeset
533 &s->v_filters[phase_y][0]);
a408778eff87 altivec accelerated v-resample patch by (Brian Foley <bfoley at compsoc dot nuigalway dot ie>)
michaelni
parents: 653
diff changeset
534 else
a408778eff87 altivec accelerated v-resample patch by (Brian Foley <bfoley at compsoc dot nuigalway dot ie>)
michaelni
parents: 653
diff changeset
535 #endif
2967
ef2149182f1c COSMETICS: Remove all trailing whitespace.
diego
parents: 2904
diff changeset
536 v_resample(output, owidth,
ef2149182f1c COSMETICS: Remove all trailing whitespace.
diego
parents: 2904
diff changeset
537 s->line_buf + (ring_y - NB_TAPS + 1) * owidth, owidth,
0
986e461dc072 Initial revision
glantau
parents:
diff changeset
538 &s->v_filters[phase_y][0]);
2967
ef2149182f1c COSMETICS: Remove all trailing whitespace.
diego
parents: 2904
diff changeset
539
0
986e461dc072 Initial revision
glantau
parents:
diff changeset
540 src_y += s->v_incr;
2967
ef2149182f1c COSMETICS: Remove all trailing whitespace.
diego
parents: 2904
diff changeset
541
0
986e461dc072 Initial revision
glantau
parents:
diff changeset
542 output += owrap;
986e461dc072 Initial revision
glantau
parents:
diff changeset
543 }
986e461dc072 Initial revision
glantau
parents:
diff changeset
544 }
986e461dc072 Initial revision
glantau
parents:
diff changeset
545
986e461dc072 Initial revision
glantau
parents:
diff changeset
546 ImgReSampleContext *img_resample_init(int owidth, int oheight,
986e461dc072 Initial revision
glantau
parents:
diff changeset
547 int iwidth, int iheight)
986e461dc072 Initial revision
glantau
parents:
diff changeset
548 {
2967
ef2149182f1c COSMETICS: Remove all trailing whitespace.
diego
parents: 2904
diff changeset
549 return img_resample_full_init(owidth, oheight, iwidth, iheight,
1928
0c23a5564489 padding support in ffmpeg patch by (Todd Kirby <doubleshot at pacbell dot net>)
michael
parents: 1488
diff changeset
550 0, 0, 0, 0, 0, 0, 0, 0);
630
b4ee42142ad1 croping patch by (talus25 at speakeasy dot net) with fixes from atmos & me
michaelni
parents: 429
diff changeset
551 }
b4ee42142ad1 croping patch by (talus25 at speakeasy dot net) with fixes from atmos & me
michaelni
parents: 429
diff changeset
552
b4ee42142ad1 croping patch by (talus25 at speakeasy dot net) with fixes from atmos & me
michaelni
parents: 429
diff changeset
553 ImgReSampleContext *img_resample_full_init(int owidth, int oheight,
b4ee42142ad1 croping patch by (talus25 at speakeasy dot net) with fixes from atmos & me
michaelni
parents: 429
diff changeset
554 int iwidth, int iheight,
b4ee42142ad1 croping patch by (talus25 at speakeasy dot net) with fixes from atmos & me
michaelni
parents: 429
diff changeset
555 int topBand, int bottomBand,
1928
0c23a5564489 padding support in ffmpeg patch by (Todd Kirby <doubleshot at pacbell dot net>)
michael
parents: 1488
diff changeset
556 int leftBand, int rightBand,
0c23a5564489 padding support in ffmpeg patch by (Todd Kirby <doubleshot at pacbell dot net>)
michael
parents: 1488
diff changeset
557 int padtop, int padbottom,
0c23a5564489 padding support in ffmpeg patch by (Todd Kirby <doubleshot at pacbell dot net>)
michael
parents: 1488
diff changeset
558 int padleft, int padright)
630
b4ee42142ad1 croping patch by (talus25 at speakeasy dot net) with fixes from atmos & me
michaelni
parents: 429
diff changeset
559 {
0
986e461dc072 Initial revision
glantau
parents:
diff changeset
560 ImgReSampleContext *s;
986e461dc072 Initial revision
glantau
parents:
diff changeset
561
2904
84404858c529 sanity check whether dimensions are non-null
alex
parents: 2846
diff changeset
562 if (!owidth || !oheight || !iwidth || !iheight)
2979
bfabfdf9ce55 COSMETICS: tabs --> spaces, some prettyprinting
diego
parents: 2967
diff changeset
563 return NULL;
2904
84404858c529 sanity check whether dimensions are non-null
alex
parents: 2846
diff changeset
564
0
986e461dc072 Initial revision
glantau
parents:
diff changeset
565 s = av_mallocz(sizeof(ImgReSampleContext));
986e461dc072 Initial revision
glantau
parents:
diff changeset
566 if (!s)
986e461dc072 Initial revision
glantau
parents:
diff changeset
567 return NULL;
2422
18b8b2dcc037 various security fixes and precautionary checks
michael
parents: 2400
diff changeset
568 if((unsigned)owidth >= UINT_MAX / (LINE_BUF_HEIGHT + NB_TAPS))
18b8b2dcc037 various security fixes and precautionary checks
michael
parents: 2400
diff changeset
569 return NULL;
0
986e461dc072 Initial revision
glantau
parents:
diff changeset
570 s->line_buf = av_mallocz(owidth * (LINE_BUF_HEIGHT + NB_TAPS));
2967
ef2149182f1c COSMETICS: Remove all trailing whitespace.
diego
parents: 2904
diff changeset
571 if (!s->line_buf)
0
986e461dc072 Initial revision
glantau
parents:
diff changeset
572 goto fail;
2967
ef2149182f1c COSMETICS: Remove all trailing whitespace.
diego
parents: 2904
diff changeset
573
0
986e461dc072 Initial revision
glantau
parents:
diff changeset
574 s->owidth = owidth;
986e461dc072 Initial revision
glantau
parents:
diff changeset
575 s->oheight = oheight;
986e461dc072 Initial revision
glantau
parents:
diff changeset
576 s->iwidth = iwidth;
986e461dc072 Initial revision
glantau
parents:
diff changeset
577 s->iheight = iheight;
2967
ef2149182f1c COSMETICS: Remove all trailing whitespace.
diego
parents: 2904
diff changeset
578
630
b4ee42142ad1 croping patch by (talus25 at speakeasy dot net) with fixes from atmos & me
michaelni
parents: 429
diff changeset
579 s->topBand = topBand;
b4ee42142ad1 croping patch by (talus25 at speakeasy dot net) with fixes from atmos & me
michaelni
parents: 429
diff changeset
580 s->bottomBand = bottomBand;
b4ee42142ad1 croping patch by (talus25 at speakeasy dot net) with fixes from atmos & me
michaelni
parents: 429
diff changeset
581 s->leftBand = leftBand;
b4ee42142ad1 croping patch by (talus25 at speakeasy dot net) with fixes from atmos & me
michaelni
parents: 429
diff changeset
582 s->rightBand = rightBand;
2967
ef2149182f1c COSMETICS: Remove all trailing whitespace.
diego
parents: 2904
diff changeset
583
1928
0c23a5564489 padding support in ffmpeg patch by (Todd Kirby <doubleshot at pacbell dot net>)
michael
parents: 1488
diff changeset
584 s->padtop = padtop;
0c23a5564489 padding support in ffmpeg patch by (Todd Kirby <doubleshot at pacbell dot net>)
michael
parents: 1488
diff changeset
585 s->padbottom = padbottom;
0c23a5564489 padding support in ffmpeg patch by (Todd Kirby <doubleshot at pacbell dot net>)
michael
parents: 1488
diff changeset
586 s->padleft = padleft;
0c23a5564489 padding support in ffmpeg patch by (Todd Kirby <doubleshot at pacbell dot net>)
michael
parents: 1488
diff changeset
587 s->padright = padright;
0c23a5564489 padding support in ffmpeg patch by (Todd Kirby <doubleshot at pacbell dot net>)
michael
parents: 1488
diff changeset
588
0c23a5564489 padding support in ffmpeg patch by (Todd Kirby <doubleshot at pacbell dot net>)
michael
parents: 1488
diff changeset
589 s->pad_owidth = owidth - (padleft + padright);
0c23a5564489 padding support in ffmpeg patch by (Todd Kirby <doubleshot at pacbell dot net>)
michael
parents: 1488
diff changeset
590 s->pad_oheight = oheight - (padtop + padbottom);
0c23a5564489 padding support in ffmpeg patch by (Todd Kirby <doubleshot at pacbell dot net>)
michael
parents: 1488
diff changeset
591
0c23a5564489 padding support in ffmpeg patch by (Todd Kirby <doubleshot at pacbell dot net>)
michael
parents: 1488
diff changeset
592 s->h_incr = ((iwidth - leftBand - rightBand) * POS_FRAC) / s->pad_owidth;
2967
ef2149182f1c COSMETICS: Remove all trailing whitespace.
diego
parents: 2904
diff changeset
593 s->v_incr = ((iheight - topBand - bottomBand) * POS_FRAC) / s->pad_oheight;
1928
0c23a5564489 padding support in ffmpeg patch by (Todd Kirby <doubleshot at pacbell dot net>)
michael
parents: 1488
diff changeset
594
2967
ef2149182f1c COSMETICS: Remove all trailing whitespace.
diego
parents: 2904
diff changeset
595 av_build_filter(&s->h_filters[0][0], (float) s->pad_owidth /
2082
3dc9bbe1b152 polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents: 2064
diff changeset
596 (float) (iwidth - leftBand - rightBand), NB_TAPS, NB_PHASES, 1<<FILTER_BITS, 0);
2967
ef2149182f1c COSMETICS: Remove all trailing whitespace.
diego
parents: 2904
diff changeset
597 av_build_filter(&s->v_filters[0][0], (float) s->pad_oheight /
2082
3dc9bbe1b152 polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents: 2064
diff changeset
598 (float) (iheight - topBand - bottomBand), NB_TAPS, NB_PHASES, 1<<FILTER_BITS, 0);
0
986e461dc072 Initial revision
glantau
parents:
diff changeset
599
986e461dc072 Initial revision
glantau
parents:
diff changeset
600 return s;
1928
0c23a5564489 padding support in ffmpeg patch by (Todd Kirby <doubleshot at pacbell dot net>)
michael
parents: 1488
diff changeset
601 fail:
396
fce0a2520551 removed useless header includes - use av memory functions
glantau
parents: 18
diff changeset
602 av_free(s);
0
986e461dc072 Initial revision
glantau
parents:
diff changeset
603 return NULL;
986e461dc072 Initial revision
glantau
parents:
diff changeset
604 }
986e461dc072 Initial revision
glantau
parents:
diff changeset
605
2967
ef2149182f1c COSMETICS: Remove all trailing whitespace.
diego
parents: 2904
diff changeset
606 void img_resample(ImgReSampleContext *s,
1488
766a2f4edbea avcodec const correctness patch by (Drew Hess <dhess at ilm dot com>)
michaelni
parents: 1106
diff changeset
607 AVPicture *output, const AVPicture *input)
0
986e461dc072 Initial revision
glantau
parents:
diff changeset
608 {
986e461dc072 Initial revision
glantau
parents:
diff changeset
609 int i, shift;
1928
0c23a5564489 padding support in ffmpeg patch by (Todd Kirby <doubleshot at pacbell dot net>)
michael
parents: 1488
diff changeset
610 uint8_t* optr;
0
986e461dc072 Initial revision
glantau
parents:
diff changeset
611
1928
0c23a5564489 padding support in ffmpeg patch by (Todd Kirby <doubleshot at pacbell dot net>)
michael
parents: 1488
diff changeset
612 for (i=0;i<3;i++) {
0
986e461dc072 Initial revision
glantau
parents:
diff changeset
613 shift = (i == 0) ? 0 : 1;
1928
0c23a5564489 padding support in ffmpeg patch by (Todd Kirby <doubleshot at pacbell dot net>)
michael
parents: 1488
diff changeset
614
2967
ef2149182f1c COSMETICS: Remove all trailing whitespace.
diego
parents: 2904
diff changeset
615 optr = output->data[i] + (((output->linesize[i] *
1928
0c23a5564489 padding support in ffmpeg patch by (Todd Kirby <doubleshot at pacbell dot net>)
michael
parents: 1488
diff changeset
616 s->padtop) + s->padleft) >> shift);
0c23a5564489 padding support in ffmpeg patch by (Todd Kirby <doubleshot at pacbell dot net>)
michael
parents: 1488
diff changeset
617
2967
ef2149182f1c COSMETICS: Remove all trailing whitespace.
diego
parents: 2904
diff changeset
618 component_resample(s, optr, output->linesize[i],
1928
0c23a5564489 padding support in ffmpeg patch by (Todd Kirby <doubleshot at pacbell dot net>)
michael
parents: 1488
diff changeset
619 s->pad_owidth >> shift, s->pad_oheight >> shift,
2967
ef2149182f1c COSMETICS: Remove all trailing whitespace.
diego
parents: 2904
diff changeset
620 input->data[i] + (input->linesize[i] *
1928
0c23a5564489 padding support in ffmpeg patch by (Todd Kirby <doubleshot at pacbell dot net>)
michael
parents: 1488
diff changeset
621 (s->topBand >> shift)) + (s->leftBand >> shift),
2967
ef2149182f1c COSMETICS: Remove all trailing whitespace.
diego
parents: 2904
diff changeset
622 input->linesize[i], ((s->iwidth - s->leftBand -
1928
0c23a5564489 padding support in ffmpeg patch by (Todd Kirby <doubleshot at pacbell dot net>)
michael
parents: 1488
diff changeset
623 s->rightBand) >> shift),
630
b4ee42142ad1 croping patch by (talus25 at speakeasy dot net) with fixes from atmos & me
michaelni
parents: 429
diff changeset
624 (s->iheight - s->topBand - s->bottomBand) >> shift);
0
986e461dc072 Initial revision
glantau
parents:
diff changeset
625 }
986e461dc072 Initial revision
glantau
parents:
diff changeset
626 }
986e461dc072 Initial revision
glantau
parents:
diff changeset
627
986e461dc072 Initial revision
glantau
parents:
diff changeset
628 void img_resample_close(ImgReSampleContext *s)
986e461dc072 Initial revision
glantau
parents:
diff changeset
629 {
396
fce0a2520551 removed useless header includes - use av memory functions
glantau
parents: 18
diff changeset
630 av_free(s->line_buf);
fce0a2520551 removed useless header includes - use av memory functions
glantau
parents: 18
diff changeset
631 av_free(s);
0
986e461dc072 Initial revision
glantau
parents:
diff changeset
632 }
986e461dc072 Initial revision
glantau
parents:
diff changeset
633
3249
ad3d12d04645 Introduce swscale interface in libavcodec
lucabe
parents: 3089
diff changeset
634 struct SwsContext *sws_getContext(int srcW, int srcH, int srcFormat,
ad3d12d04645 Introduce swscale interface in libavcodec
lucabe
parents: 3089
diff changeset
635 int dstW, int dstH, int dstFormat,
ad3d12d04645 Introduce swscale interface in libavcodec
lucabe
parents: 3089
diff changeset
636 int flags, SwsFilter *srcFilter,
ad3d12d04645 Introduce swscale interface in libavcodec
lucabe
parents: 3089
diff changeset
637 SwsFilter *dstFilter, double *param)
ad3d12d04645 Introduce swscale interface in libavcodec
lucabe
parents: 3089
diff changeset
638 {
ad3d12d04645 Introduce swscale interface in libavcodec
lucabe
parents: 3089
diff changeset
639 struct SwsContext *ctx;
ad3d12d04645 Introduce swscale interface in libavcodec
lucabe
parents: 3089
diff changeset
640
ad3d12d04645 Introduce swscale interface in libavcodec
lucabe
parents: 3089
diff changeset
641 ctx = av_malloc(sizeof(struct SwsContext));
ad3d12d04645 Introduce swscale interface in libavcodec
lucabe
parents: 3089
diff changeset
642 if (ctx == NULL) {
ad3d12d04645 Introduce swscale interface in libavcodec
lucabe
parents: 3089
diff changeset
643 av_log(NULL, AV_LOG_ERROR, "Cannot allocate a resampling context!\n");
ad3d12d04645 Introduce swscale interface in libavcodec
lucabe
parents: 3089
diff changeset
644
ad3d12d04645 Introduce swscale interface in libavcodec
lucabe
parents: 3089
diff changeset
645 return NULL;
ad3d12d04645 Introduce swscale interface in libavcodec
lucabe
parents: 3089
diff changeset
646 }
ad3d12d04645 Introduce swscale interface in libavcodec
lucabe
parents: 3089
diff changeset
647
ad3d12d04645 Introduce swscale interface in libavcodec
lucabe
parents: 3089
diff changeset
648 if ((srcH != dstH) || (srcW != dstW)) {
ad3d12d04645 Introduce swscale interface in libavcodec
lucabe
parents: 3089
diff changeset
649 if ((srcFormat != PIX_FMT_YUV420P) || (dstFormat != PIX_FMT_YUV420P)) {
ad3d12d04645 Introduce swscale interface in libavcodec
lucabe
parents: 3089
diff changeset
650 av_log(NULL, AV_LOG_INFO, "PIX_FMT_YUV420P will be used as an intermediate format for rescaling\n");
ad3d12d04645 Introduce swscale interface in libavcodec
lucabe
parents: 3089
diff changeset
651 }
ad3d12d04645 Introduce swscale interface in libavcodec
lucabe
parents: 3089
diff changeset
652 ctx->resampling_ctx = img_resample_init(dstW, dstH, srcW, srcH);
ad3d12d04645 Introduce swscale interface in libavcodec
lucabe
parents: 3089
diff changeset
653 } else {
ad3d12d04645 Introduce swscale interface in libavcodec
lucabe
parents: 3089
diff changeset
654 ctx->resampling_ctx = av_malloc(sizeof(ImgReSampleContext));
ad3d12d04645 Introduce swscale interface in libavcodec
lucabe
parents: 3089
diff changeset
655 ctx->resampling_ctx->iheight = srcH;
ad3d12d04645 Introduce swscale interface in libavcodec
lucabe
parents: 3089
diff changeset
656 ctx->resampling_ctx->iwidth = srcW;
ad3d12d04645 Introduce swscale interface in libavcodec
lucabe
parents: 3089
diff changeset
657 ctx->resampling_ctx->oheight = dstH;
ad3d12d04645 Introduce swscale interface in libavcodec
lucabe
parents: 3089
diff changeset
658 ctx->resampling_ctx->owidth = dstW;
ad3d12d04645 Introduce swscale interface in libavcodec
lucabe
parents: 3089
diff changeset
659 }
ad3d12d04645 Introduce swscale interface in libavcodec
lucabe
parents: 3089
diff changeset
660 ctx->src_pix_fmt = srcFormat;
ad3d12d04645 Introduce swscale interface in libavcodec
lucabe
parents: 3089
diff changeset
661 ctx->dst_pix_fmt = dstFormat;
ad3d12d04645 Introduce swscale interface in libavcodec
lucabe
parents: 3089
diff changeset
662
ad3d12d04645 Introduce swscale interface in libavcodec
lucabe
parents: 3089
diff changeset
663 return ctx;
ad3d12d04645 Introduce swscale interface in libavcodec
lucabe
parents: 3089
diff changeset
664 }
ad3d12d04645 Introduce swscale interface in libavcodec
lucabe
parents: 3089
diff changeset
665
ad3d12d04645 Introduce swscale interface in libavcodec
lucabe
parents: 3089
diff changeset
666 void sws_freeContext(struct SwsContext *ctx)
ad3d12d04645 Introduce swscale interface in libavcodec
lucabe
parents: 3089
diff changeset
667 {
ad3d12d04645 Introduce swscale interface in libavcodec
lucabe
parents: 3089
diff changeset
668 if ((ctx->resampling_ctx->iwidth != ctx->resampling_ctx->owidth) ||
ad3d12d04645 Introduce swscale interface in libavcodec
lucabe
parents: 3089
diff changeset
669 (ctx->resampling_ctx->iheight != ctx->resampling_ctx->oheight)) {
ad3d12d04645 Introduce swscale interface in libavcodec
lucabe
parents: 3089
diff changeset
670 img_resample_close(ctx->resampling_ctx);
ad3d12d04645 Introduce swscale interface in libavcodec
lucabe
parents: 3089
diff changeset
671 } else {
ad3d12d04645 Introduce swscale interface in libavcodec
lucabe
parents: 3089
diff changeset
672 av_free(ctx->resampling_ctx);
ad3d12d04645 Introduce swscale interface in libavcodec
lucabe
parents: 3089
diff changeset
673 }
ad3d12d04645 Introduce swscale interface in libavcodec
lucabe
parents: 3089
diff changeset
674 av_free(ctx);
ad3d12d04645 Introduce swscale interface in libavcodec
lucabe
parents: 3089
diff changeset
675 }
ad3d12d04645 Introduce swscale interface in libavcodec
lucabe
parents: 3089
diff changeset
676
ad3d12d04645 Introduce swscale interface in libavcodec
lucabe
parents: 3089
diff changeset
677 int sws_scale(struct SwsContext *ctx, uint8_t* src[], int srcStride[],
ad3d12d04645 Introduce swscale interface in libavcodec
lucabe
parents: 3089
diff changeset
678 int srcSliceY, int srcSliceH, uint8_t* dst[], int dstStride[])
ad3d12d04645 Introduce swscale interface in libavcodec
lucabe
parents: 3089
diff changeset
679 {
ad3d12d04645 Introduce swscale interface in libavcodec
lucabe
parents: 3089
diff changeset
680 AVPicture src_pict, dst_pict;
ad3d12d04645 Introduce swscale interface in libavcodec
lucabe
parents: 3089
diff changeset
681 int i, res = 0;
ad3d12d04645 Introduce swscale interface in libavcodec
lucabe
parents: 3089
diff changeset
682 AVPicture picture_format_temp;
ad3d12d04645 Introduce swscale interface in libavcodec
lucabe
parents: 3089
diff changeset
683 AVPicture picture_resample_temp, *formatted_picture, *resampled_picture;
ad3d12d04645 Introduce swscale interface in libavcodec
lucabe
parents: 3089
diff changeset
684 uint8_t *buf1 = NULL, *buf2 = NULL;
ad3d12d04645 Introduce swscale interface in libavcodec
lucabe
parents: 3089
diff changeset
685 enum PixelFormat current_pix_fmt;
ad3d12d04645 Introduce swscale interface in libavcodec
lucabe
parents: 3089
diff changeset
686
ad3d12d04645 Introduce swscale interface in libavcodec
lucabe
parents: 3089
diff changeset
687 for (i = 0; i < 3; i++) {
ad3d12d04645 Introduce swscale interface in libavcodec
lucabe
parents: 3089
diff changeset
688 src_pict.data[i] = src[i];
ad3d12d04645 Introduce swscale interface in libavcodec
lucabe
parents: 3089
diff changeset
689 src_pict.linesize[i] = srcStride[i];
ad3d12d04645 Introduce swscale interface in libavcodec
lucabe
parents: 3089
diff changeset
690 dst_pict.data[i] = dst[i];
ad3d12d04645 Introduce swscale interface in libavcodec
lucabe
parents: 3089
diff changeset
691 dst_pict.linesize[i] = dstStride[i];
ad3d12d04645 Introduce swscale interface in libavcodec
lucabe
parents: 3089
diff changeset
692 }
ad3d12d04645 Introduce swscale interface in libavcodec
lucabe
parents: 3089
diff changeset
693 if ((ctx->resampling_ctx->iwidth != ctx->resampling_ctx->owidth) ||
ad3d12d04645 Introduce swscale interface in libavcodec
lucabe
parents: 3089
diff changeset
694 (ctx->resampling_ctx->iheight != ctx->resampling_ctx->oheight)) {
ad3d12d04645 Introduce swscale interface in libavcodec
lucabe
parents: 3089
diff changeset
695 /* We have to rescale the picture, but only YUV420P rescaling is supported... */
ad3d12d04645 Introduce swscale interface in libavcodec
lucabe
parents: 3089
diff changeset
696
ad3d12d04645 Introduce swscale interface in libavcodec
lucabe
parents: 3089
diff changeset
697 if (ctx->src_pix_fmt != PIX_FMT_YUV420P) {
ad3d12d04645 Introduce swscale interface in libavcodec
lucabe
parents: 3089
diff changeset
698 int size;
ad3d12d04645 Introduce swscale interface in libavcodec
lucabe
parents: 3089
diff changeset
699
ad3d12d04645 Introduce swscale interface in libavcodec
lucabe
parents: 3089
diff changeset
700 /* create temporary picture for rescaling input*/
ad3d12d04645 Introduce swscale interface in libavcodec
lucabe
parents: 3089
diff changeset
701 size = avpicture_get_size(PIX_FMT_YUV420P, ctx->resampling_ctx->iwidth, ctx->resampling_ctx->iheight);
ad3d12d04645 Introduce swscale interface in libavcodec
lucabe
parents: 3089
diff changeset
702 buf1 = av_malloc(size);
ad3d12d04645 Introduce swscale interface in libavcodec
lucabe
parents: 3089
diff changeset
703 if (!buf1) {
ad3d12d04645 Introduce swscale interface in libavcodec
lucabe
parents: 3089
diff changeset
704 res = -1;
ad3d12d04645 Introduce swscale interface in libavcodec
lucabe
parents: 3089
diff changeset
705 goto the_end;
ad3d12d04645 Introduce swscale interface in libavcodec
lucabe
parents: 3089
diff changeset
706 }
ad3d12d04645 Introduce swscale interface in libavcodec
lucabe
parents: 3089
diff changeset
707 formatted_picture = &picture_format_temp;
ad3d12d04645 Introduce swscale interface in libavcodec
lucabe
parents: 3089
diff changeset
708 avpicture_fill((AVPicture*)formatted_picture, buf1,
ad3d12d04645 Introduce swscale interface in libavcodec
lucabe
parents: 3089
diff changeset
709 PIX_FMT_YUV420P, ctx->resampling_ctx->iwidth, ctx->resampling_ctx->iheight);
ad3d12d04645 Introduce swscale interface in libavcodec
lucabe
parents: 3089
diff changeset
710
ad3d12d04645 Introduce swscale interface in libavcodec
lucabe
parents: 3089
diff changeset
711 if (img_convert((AVPicture*)formatted_picture, PIX_FMT_YUV420P,
ad3d12d04645 Introduce swscale interface in libavcodec
lucabe
parents: 3089
diff changeset
712 &src_pict, ctx->src_pix_fmt,
ad3d12d04645 Introduce swscale interface in libavcodec
lucabe
parents: 3089
diff changeset
713 ctx->resampling_ctx->iwidth, ctx->resampling_ctx->iheight) < 0) {
ad3d12d04645 Introduce swscale interface in libavcodec
lucabe
parents: 3089
diff changeset
714
ad3d12d04645 Introduce swscale interface in libavcodec
lucabe
parents: 3089
diff changeset
715 av_log(NULL, AV_LOG_ERROR, "pixel format conversion not handled\n");
ad3d12d04645 Introduce swscale interface in libavcodec
lucabe
parents: 3089
diff changeset
716 res = -1;
ad3d12d04645 Introduce swscale interface in libavcodec
lucabe
parents: 3089
diff changeset
717 goto the_end;
ad3d12d04645 Introduce swscale interface in libavcodec
lucabe
parents: 3089
diff changeset
718 }
ad3d12d04645 Introduce swscale interface in libavcodec
lucabe
parents: 3089
diff changeset
719 } else {
ad3d12d04645 Introduce swscale interface in libavcodec
lucabe
parents: 3089
diff changeset
720 formatted_picture = &src_pict;
ad3d12d04645 Introduce swscale interface in libavcodec
lucabe
parents: 3089
diff changeset
721 }
ad3d12d04645 Introduce swscale interface in libavcodec
lucabe
parents: 3089
diff changeset
722
ad3d12d04645 Introduce swscale interface in libavcodec
lucabe
parents: 3089
diff changeset
723 if (ctx->dst_pix_fmt != PIX_FMT_YUV420P) {
ad3d12d04645 Introduce swscale interface in libavcodec
lucabe
parents: 3089
diff changeset
724 int size;
ad3d12d04645 Introduce swscale interface in libavcodec
lucabe
parents: 3089
diff changeset
725
ad3d12d04645 Introduce swscale interface in libavcodec
lucabe
parents: 3089
diff changeset
726 /* create temporary picture for rescaling output*/
ad3d12d04645 Introduce swscale interface in libavcodec
lucabe
parents: 3089
diff changeset
727 size = avpicture_get_size(PIX_FMT_YUV420P, ctx->resampling_ctx->owidth, ctx->resampling_ctx->oheight);
ad3d12d04645 Introduce swscale interface in libavcodec
lucabe
parents: 3089
diff changeset
728 buf2 = av_malloc(size);
ad3d12d04645 Introduce swscale interface in libavcodec
lucabe
parents: 3089
diff changeset
729 if (!buf2) {
ad3d12d04645 Introduce swscale interface in libavcodec
lucabe
parents: 3089
diff changeset
730 res = -1;
ad3d12d04645 Introduce swscale interface in libavcodec
lucabe
parents: 3089
diff changeset
731 goto the_end;
ad3d12d04645 Introduce swscale interface in libavcodec
lucabe
parents: 3089
diff changeset
732 }
ad3d12d04645 Introduce swscale interface in libavcodec
lucabe
parents: 3089
diff changeset
733 resampled_picture = &picture_resample_temp;
ad3d12d04645 Introduce swscale interface in libavcodec
lucabe
parents: 3089
diff changeset
734 avpicture_fill((AVPicture*)resampled_picture, buf2,
ad3d12d04645 Introduce swscale interface in libavcodec
lucabe
parents: 3089
diff changeset
735 PIX_FMT_YUV420P, ctx->resampling_ctx->owidth, ctx->resampling_ctx->oheight);
ad3d12d04645 Introduce swscale interface in libavcodec
lucabe
parents: 3089
diff changeset
736
ad3d12d04645 Introduce swscale interface in libavcodec
lucabe
parents: 3089
diff changeset
737 } else {
ad3d12d04645 Introduce swscale interface in libavcodec
lucabe
parents: 3089
diff changeset
738 resampled_picture = &dst_pict;
ad3d12d04645 Introduce swscale interface in libavcodec
lucabe
parents: 3089
diff changeset
739 }
ad3d12d04645 Introduce swscale interface in libavcodec
lucabe
parents: 3089
diff changeset
740
ad3d12d04645 Introduce swscale interface in libavcodec
lucabe
parents: 3089
diff changeset
741 /* ...and finally rescale!!! */
ad3d12d04645 Introduce swscale interface in libavcodec
lucabe
parents: 3089
diff changeset
742 img_resample(ctx->resampling_ctx, resampled_picture, formatted_picture);
ad3d12d04645 Introduce swscale interface in libavcodec
lucabe
parents: 3089
diff changeset
743 current_pix_fmt = PIX_FMT_YUV420P;
ad3d12d04645 Introduce swscale interface in libavcodec
lucabe
parents: 3089
diff changeset
744 } else {
ad3d12d04645 Introduce swscale interface in libavcodec
lucabe
parents: 3089
diff changeset
745 resampled_picture = &src_pict;
ad3d12d04645 Introduce swscale interface in libavcodec
lucabe
parents: 3089
diff changeset
746 current_pix_fmt = ctx->src_pix_fmt;
ad3d12d04645 Introduce swscale interface in libavcodec
lucabe
parents: 3089
diff changeset
747 }
ad3d12d04645 Introduce swscale interface in libavcodec
lucabe
parents: 3089
diff changeset
748
ad3d12d04645 Introduce swscale interface in libavcodec
lucabe
parents: 3089
diff changeset
749 if (current_pix_fmt != ctx->dst_pix_fmt) {
ad3d12d04645 Introduce swscale interface in libavcodec
lucabe
parents: 3089
diff changeset
750 if (img_convert(&dst_pict, ctx->dst_pix_fmt,
ad3d12d04645 Introduce swscale interface in libavcodec
lucabe
parents: 3089
diff changeset
751 resampled_picture, current_pix_fmt,
ad3d12d04645 Introduce swscale interface in libavcodec
lucabe
parents: 3089
diff changeset
752 ctx->resampling_ctx->owidth, ctx->resampling_ctx->oheight) < 0) {
ad3d12d04645 Introduce swscale interface in libavcodec
lucabe
parents: 3089
diff changeset
753
ad3d12d04645 Introduce swscale interface in libavcodec
lucabe
parents: 3089
diff changeset
754 av_log(NULL, AV_LOG_ERROR, "pixel format conversion not handled\n");
ad3d12d04645 Introduce swscale interface in libavcodec
lucabe
parents: 3089
diff changeset
755
ad3d12d04645 Introduce swscale interface in libavcodec
lucabe
parents: 3089
diff changeset
756 res = -1;
ad3d12d04645 Introduce swscale interface in libavcodec
lucabe
parents: 3089
diff changeset
757 goto the_end;
ad3d12d04645 Introduce swscale interface in libavcodec
lucabe
parents: 3089
diff changeset
758 }
ad3d12d04645 Introduce swscale interface in libavcodec
lucabe
parents: 3089
diff changeset
759 }
ad3d12d04645 Introduce swscale interface in libavcodec
lucabe
parents: 3089
diff changeset
760
ad3d12d04645 Introduce swscale interface in libavcodec
lucabe
parents: 3089
diff changeset
761 the_end:
ad3d12d04645 Introduce swscale interface in libavcodec
lucabe
parents: 3089
diff changeset
762 av_free(buf1);
ad3d12d04645 Introduce swscale interface in libavcodec
lucabe
parents: 3089
diff changeset
763 av_free(buf2);
ad3d12d04645 Introduce swscale interface in libavcodec
lucabe
parents: 3089
diff changeset
764 return res;
ad3d12d04645 Introduce swscale interface in libavcodec
lucabe
parents: 3089
diff changeset
765 }
ad3d12d04645 Introduce swscale interface in libavcodec
lucabe
parents: 3089
diff changeset
766
ad3d12d04645 Introduce swscale interface in libavcodec
lucabe
parents: 3089
diff changeset
767
0
986e461dc072 Initial revision
glantau
parents:
diff changeset
768 #ifdef TEST
2400
17ec73c65748 imgresample test cleanup patch by (Panagiotis Issaris <takis )( lumumba d0t luc d0t ac.be>)
michael
parents: 2082
diff changeset
769 #include <stdio.h>
630
b4ee42142ad1 croping patch by (talus25 at speakeasy dot net) with fixes from atmos & me
michaelni
parents: 429
diff changeset
770
0
986e461dc072 Initial revision
glantau
parents:
diff changeset
771 /* input */
986e461dc072 Initial revision
glantau
parents:
diff changeset
772 #define XSIZE 256
986e461dc072 Initial revision
glantau
parents:
diff changeset
773 #define YSIZE 256
1064
b32afefe7d33 * UINTX -> uintx_t INTX -> intx_t
kabi
parents: 1057
diff changeset
774 uint8_t img[XSIZE * YSIZE];
0
986e461dc072 Initial revision
glantau
parents:
diff changeset
775
986e461dc072 Initial revision
glantau
parents:
diff changeset
776 /* output */
986e461dc072 Initial revision
glantau
parents:
diff changeset
777 #define XSIZE1 512
986e461dc072 Initial revision
glantau
parents:
diff changeset
778 #define YSIZE1 512
1064
b32afefe7d33 * UINTX -> uintx_t INTX -> intx_t
kabi
parents: 1057
diff changeset
779 uint8_t img1[XSIZE1 * YSIZE1];
b32afefe7d33 * UINTX -> uintx_t INTX -> intx_t
kabi
parents: 1057
diff changeset
780 uint8_t img2[XSIZE1 * YSIZE1];
0
986e461dc072 Initial revision
glantau
parents:
diff changeset
781
1064
b32afefe7d33 * UINTX -> uintx_t INTX -> intx_t
kabi
parents: 1057
diff changeset
782 void save_pgm(const char *filename, uint8_t *img, int xsize, int ysize)
0
986e461dc072 Initial revision
glantau
parents:
diff changeset
783 {
2846
40765c51a7a9 Compilation fixes part 1 patch by (Arvind R. and Burkhard Plaum, plaum, ipf uni-stuttgart de)
michael
parents: 2423
diff changeset
784 #undef fprintf
0
986e461dc072 Initial revision
glantau
parents:
diff changeset
785 FILE *f;
986e461dc072 Initial revision
glantau
parents:
diff changeset
786 f=fopen(filename,"w");
986e461dc072 Initial revision
glantau
parents:
diff changeset
787 fprintf(f,"P5\n%d %d\n%d\n", xsize, ysize, 255);
986e461dc072 Initial revision
glantau
parents:
diff changeset
788 fwrite(img,1, xsize * ysize,f);
986e461dc072 Initial revision
glantau
parents:
diff changeset
789 fclose(f);
2846
40765c51a7a9 Compilation fixes part 1 patch by (Arvind R. and Burkhard Plaum, plaum, ipf uni-stuttgart de)
michael
parents: 2423
diff changeset
790 #define fprintf please_use_av_log
0
986e461dc072 Initial revision
glantau
parents:
diff changeset
791 }
986e461dc072 Initial revision
glantau
parents:
diff changeset
792
1064
b32afefe7d33 * UINTX -> uintx_t INTX -> intx_t
kabi
parents: 1057
diff changeset
793 static void dump_filter(int16_t *filter)
0
986e461dc072 Initial revision
glantau
parents:
diff changeset
794 {
986e461dc072 Initial revision
glantau
parents:
diff changeset
795 int i, ph;
986e461dc072 Initial revision
glantau
parents:
diff changeset
796
986e461dc072 Initial revision
glantau
parents:
diff changeset
797 for(ph=0;ph<NB_PHASES;ph++) {
2400
17ec73c65748 imgresample test cleanup patch by (Panagiotis Issaris <takis )( lumumba d0t luc d0t ac.be>)
michael
parents: 2082
diff changeset
798 av_log(NULL, AV_LOG_INFO, "%2d: ", ph);
0
986e461dc072 Initial revision
glantau
parents:
diff changeset
799 for(i=0;i<NB_TAPS;i++) {
2400
17ec73c65748 imgresample test cleanup patch by (Panagiotis Issaris <takis )( lumumba d0t luc d0t ac.be>)
michael
parents: 2082
diff changeset
800 av_log(NULL, AV_LOG_INFO, " %5.2f", filter[ph * NB_TAPS + i] / 256.0);
0
986e461dc072 Initial revision
glantau
parents:
diff changeset
801 }
2400
17ec73c65748 imgresample test cleanup patch by (Panagiotis Issaris <takis )( lumumba d0t luc d0t ac.be>)
michael
parents: 2082
diff changeset
802 av_log(NULL, AV_LOG_INFO, "\n");
0
986e461dc072 Initial revision
glantau
parents:
diff changeset
803 }
986e461dc072 Initial revision
glantau
parents:
diff changeset
804 }
986e461dc072 Initial revision
glantau
parents:
diff changeset
805
2
2e2c46c87460 fixed config for direct mplayer build compatibility
glantau
parents: 0
diff changeset
806 #ifdef HAVE_MMX
644
641708b5c33c fixing tests
michaelni
parents: 630
diff changeset
807 int mm_flags;
0
986e461dc072 Initial revision
glantau
parents:
diff changeset
808 #endif
986e461dc072 Initial revision
glantau
parents:
diff changeset
809
986e461dc072 Initial revision
glantau
parents:
diff changeset
810 int main(int argc, char **argv)
986e461dc072 Initial revision
glantau
parents:
diff changeset
811 {
986e461dc072 Initial revision
glantau
parents:
diff changeset
812 int x, y, v, i, xsize, ysize;
986e461dc072 Initial revision
glantau
parents:
diff changeset
813 ImgReSampleContext *s;
986e461dc072 Initial revision
glantau
parents:
diff changeset
814 float fact, factors[] = { 1/2.0, 3.0/4.0, 1.0, 4.0/3.0, 16.0/9.0, 2.0 };
986e461dc072 Initial revision
glantau
parents:
diff changeset
815 char buf[256];
986e461dc072 Initial revision
glantau
parents:
diff changeset
816
986e461dc072 Initial revision
glantau
parents:
diff changeset
817 /* build test image */
986e461dc072 Initial revision
glantau
parents:
diff changeset
818 for(y=0;y<YSIZE;y++) {
986e461dc072 Initial revision
glantau
parents:
diff changeset
819 for(x=0;x<XSIZE;x++) {
986e461dc072 Initial revision
glantau
parents:
diff changeset
820 if (x < XSIZE/2 && y < YSIZE/2) {
986e461dc072 Initial revision
glantau
parents:
diff changeset
821 if (x < XSIZE/4 && y < YSIZE/4) {
986e461dc072 Initial revision
glantau
parents:
diff changeset
822 if ((x % 10) <= 6 &&
986e461dc072 Initial revision
glantau
parents:
diff changeset
823 (y % 10) <= 6)
986e461dc072 Initial revision
glantau
parents:
diff changeset
824 v = 0xff;
986e461dc072 Initial revision
glantau
parents:
diff changeset
825 else
986e461dc072 Initial revision
glantau
parents:
diff changeset
826 v = 0x00;
986e461dc072 Initial revision
glantau
parents:
diff changeset
827 } else if (x < XSIZE/4) {
2967
ef2149182f1c COSMETICS: Remove all trailing whitespace.
diego
parents: 2904
diff changeset
828 if (x & 1)
0
986e461dc072 Initial revision
glantau
parents:
diff changeset
829 v = 0xff;
2967
ef2149182f1c COSMETICS: Remove all trailing whitespace.
diego
parents: 2904
diff changeset
830 else
0
986e461dc072 Initial revision
glantau
parents:
diff changeset
831 v = 0;
986e461dc072 Initial revision
glantau
parents:
diff changeset
832 } else if (y < XSIZE/4) {
2967
ef2149182f1c COSMETICS: Remove all trailing whitespace.
diego
parents: 2904
diff changeset
833 if (y & 1)
0
986e461dc072 Initial revision
glantau
parents:
diff changeset
834 v = 0xff;
2967
ef2149182f1c COSMETICS: Remove all trailing whitespace.
diego
parents: 2904
diff changeset
835 else
0
986e461dc072 Initial revision
glantau
parents:
diff changeset
836 v = 0;
986e461dc072 Initial revision
glantau
parents:
diff changeset
837 } else {
986e461dc072 Initial revision
glantau
parents:
diff changeset
838 if (y < YSIZE*3/8) {
2967
ef2149182f1c COSMETICS: Remove all trailing whitespace.
diego
parents: 2904
diff changeset
839 if ((y+x) & 1)
0
986e461dc072 Initial revision
glantau
parents:
diff changeset
840 v = 0xff;
2967
ef2149182f1c COSMETICS: Remove all trailing whitespace.
diego
parents: 2904
diff changeset
841 else
0
986e461dc072 Initial revision
glantau
parents:
diff changeset
842 v = 0;
986e461dc072 Initial revision
glantau
parents:
diff changeset
843 } else {
986e461dc072 Initial revision
glantau
parents:
diff changeset
844 if (((x+3) % 4) <= 1 &&
986e461dc072 Initial revision
glantau
parents:
diff changeset
845 ((y+3) % 4) <= 1)
986e461dc072 Initial revision
glantau
parents:
diff changeset
846 v = 0xff;
986e461dc072 Initial revision
glantau
parents:
diff changeset
847 else
986e461dc072 Initial revision
glantau
parents:
diff changeset
848 v = 0x00;
986e461dc072 Initial revision
glantau
parents:
diff changeset
849 }
986e461dc072 Initial revision
glantau
parents:
diff changeset
850 }
986e461dc072 Initial revision
glantau
parents:
diff changeset
851 } else if (x < XSIZE/2) {
986e461dc072 Initial revision
glantau
parents:
diff changeset
852 v = ((x - (XSIZE/2)) * 255) / (XSIZE/2);
986e461dc072 Initial revision
glantau
parents:
diff changeset
853 } else if (y < XSIZE/2) {
986e461dc072 Initial revision
glantau
parents:
diff changeset
854 v = ((y - (XSIZE/2)) * 255) / (XSIZE/2);
986e461dc072 Initial revision
glantau
parents:
diff changeset
855 } else {
986e461dc072 Initial revision
glantau
parents:
diff changeset
856 v = ((x + y - XSIZE) * 255) / XSIZE;
986e461dc072 Initial revision
glantau
parents:
diff changeset
857 }
630
b4ee42142ad1 croping patch by (talus25 at speakeasy dot net) with fixes from atmos & me
michaelni
parents: 429
diff changeset
858 img[(YSIZE - y) * XSIZE + (XSIZE - x)] = v;
0
986e461dc072 Initial revision
glantau
parents:
diff changeset
859 }
986e461dc072 Initial revision
glantau
parents:
diff changeset
860 }
986e461dc072 Initial revision
glantau
parents:
diff changeset
861 save_pgm("/tmp/in.pgm", img, XSIZE, YSIZE);
986e461dc072 Initial revision
glantau
parents:
diff changeset
862 for(i=0;i<sizeof(factors)/sizeof(float);i++) {
986e461dc072 Initial revision
glantau
parents:
diff changeset
863 fact = factors[i];
986e461dc072 Initial revision
glantau
parents:
diff changeset
864 xsize = (int)(XSIZE * fact);
630
b4ee42142ad1 croping patch by (talus25 at speakeasy dot net) with fixes from atmos & me
michaelni
parents: 429
diff changeset
865 ysize = (int)((YSIZE - 100) * fact);
2400
17ec73c65748 imgresample test cleanup patch by (Panagiotis Issaris <takis )( lumumba d0t luc d0t ac.be>)
michael
parents: 2082
diff changeset
866 s = img_resample_full_init(xsize, ysize, XSIZE, YSIZE, 50 ,50, 0, 0, 0, 0, 0, 0);
17ec73c65748 imgresample test cleanup patch by (Panagiotis Issaris <takis )( lumumba d0t luc d0t ac.be>)
michael
parents: 2082
diff changeset
867 av_log(NULL, AV_LOG_INFO, "Factor=%0.2f\n", fact);
0
986e461dc072 Initial revision
glantau
parents:
diff changeset
868 dump_filter(&s->h_filters[0][0]);
986e461dc072 Initial revision
glantau
parents:
diff changeset
869 component_resample(s, img1, xsize, xsize, ysize,
630
b4ee42142ad1 croping patch by (talus25 at speakeasy dot net) with fixes from atmos & me
michaelni
parents: 429
diff changeset
870 img + 50 * XSIZE, XSIZE, XSIZE, YSIZE - 100);
0
986e461dc072 Initial revision
glantau
parents:
diff changeset
871 img_resample_close(s);
986e461dc072 Initial revision
glantau
parents:
diff changeset
872
2423
87b7fbed8609 dissallow sprintf
michael
parents: 2422
diff changeset
873 snprintf(buf, sizeof(buf), "/tmp/out%d.pgm", i);
0
986e461dc072 Initial revision
glantau
parents:
diff changeset
874 save_pgm(buf, img1, xsize, ysize);
986e461dc072 Initial revision
glantau
parents:
diff changeset
875 }
986e461dc072 Initial revision
glantau
parents:
diff changeset
876
986e461dc072 Initial revision
glantau
parents:
diff changeset
877 /* mmx test */
2
2e2c46c87460 fixed config for direct mplayer build compatibility
glantau
parents: 0
diff changeset
878 #ifdef HAVE_MMX
2400
17ec73c65748 imgresample test cleanup patch by (Panagiotis Issaris <takis )( lumumba d0t luc d0t ac.be>)
michael
parents: 2082
diff changeset
879 av_log(NULL, AV_LOG_INFO, "MMX test\n");
0
986e461dc072 Initial revision
glantau
parents:
diff changeset
880 fact = 0.72;
986e461dc072 Initial revision
glantau
parents:
diff changeset
881 xsize = (int)(XSIZE * fact);
986e461dc072 Initial revision
glantau
parents:
diff changeset
882 ysize = (int)(YSIZE * fact);
986e461dc072 Initial revision
glantau
parents:
diff changeset
883 mm_flags = MM_MMX;
986e461dc072 Initial revision
glantau
parents:
diff changeset
884 s = img_resample_init(xsize, ysize, XSIZE, YSIZE);
986e461dc072 Initial revision
glantau
parents:
diff changeset
885 component_resample(s, img1, xsize, xsize, ysize,
986e461dc072 Initial revision
glantau
parents:
diff changeset
886 img, XSIZE, XSIZE, YSIZE);
986e461dc072 Initial revision
glantau
parents:
diff changeset
887
986e461dc072 Initial revision
glantau
parents:
diff changeset
888 mm_flags = 0;
986e461dc072 Initial revision
glantau
parents:
diff changeset
889 s = img_resample_init(xsize, ysize, XSIZE, YSIZE);
986e461dc072 Initial revision
glantau
parents:
diff changeset
890 component_resample(s, img2, xsize, xsize, ysize,
986e461dc072 Initial revision
glantau
parents:
diff changeset
891 img, XSIZE, XSIZE, YSIZE);
986e461dc072 Initial revision
glantau
parents:
diff changeset
892 if (memcmp(img1, img2, xsize * ysize) != 0) {
2400
17ec73c65748 imgresample test cleanup patch by (Panagiotis Issaris <takis )( lumumba d0t luc d0t ac.be>)
michael
parents: 2082
diff changeset
893 av_log(NULL, AV_LOG_ERROR, "mmx error\n");
0
986e461dc072 Initial revision
glantau
parents:
diff changeset
894 exit(1);
986e461dc072 Initial revision
glantau
parents:
diff changeset
895 }
2400
17ec73c65748 imgresample test cleanup patch by (Panagiotis Issaris <takis )( lumumba d0t luc d0t ac.be>)
michael
parents: 2082
diff changeset
896 av_log(NULL, AV_LOG_INFO, "MMX OK\n");
0
986e461dc072 Initial revision
glantau
parents:
diff changeset
897 #endif
986e461dc072 Initial revision
glantau
parents:
diff changeset
898 return 0;
986e461dc072 Initial revision
glantau
parents:
diff changeset
899 }
986e461dc072 Initial revision
glantau
parents:
diff changeset
900
986e461dc072 Initial revision
glantau
parents:
diff changeset
901 #endif