Revert r31153. It failed to build on:
[ffmpeg.git] / libswscale / swscale.c
1 /*
2  * Copyright (C) 2001-2003 Michael Niedermayer <michaelni@gmx.at>
3  *
4  * This file is part of FFmpeg.
5  *
6  * FFmpeg is free software; you can redistribute it and/or
7  * modify it under the terms of the GNU Lesser General Public
8  * License as published by the Free Software Foundation; either
9  * version 2.1 of the License, or (at your option) any later version.
10  *
11  * FFmpeg is distributed in the hope that it will be useful,
12  * but WITHOUT ANY WARRANTY; without even the implied warranty of
13  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
14  * Lesser General Public License for more details.
15  *
16  * You should have received a copy of the GNU Lesser General Public
17  * License along with FFmpeg; if not, write to the Free Software
18  * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
19  */
20
21 /*
22   supported Input formats: YV12, I420/IYUV, YUY2, UYVY, BGR32, BGR32_1, BGR24, BGR16, BGR15, RGB32, RGB32_1, RGB24, Y8/Y800, YVU9/IF09, PAL8
23   supported output formats: YV12, I420/IYUV, YUY2, UYVY, {BGR,RGB}{1,4,8,15,16,24,32}, Y8/Y800, YVU9/IF09
24   {BGR,RGB}{1,4,8,15,16} support dithering
25
26   unscaled special converters (YV12=I420=IYUV, Y800=Y8)
27   YV12 -> {BGR,RGB}{1,4,8,12,15,16,24,32}
28   x -> x
29   YUV9 -> YV12
30   YUV9/YV12 -> Y800
31   Y800 -> YUV9/YV12
32   BGR24 -> BGR32 & RGB24 -> RGB32
33   BGR32 -> BGR24 & RGB32 -> RGB24
34   BGR15 -> BGR16
35 */
36
37 /*
38 tested special converters (most are tested actually, but I did not write it down ...)
39  YV12 -> BGR12/BGR16
40  YV12 -> YV12
41  BGR15 -> BGR16
42  BGR16 -> BGR16
43  YVU9 -> YV12
44
45 untested special converters
46   YV12/I420 -> BGR15/BGR24/BGR32 (it is the yuv2rgb stuff, so it should be OK)
47   YV12/I420 -> YV12/I420
48   YUY2/BGR15/BGR24/BGR32/RGB24/RGB32 -> same format
49   BGR24 -> BGR32 & RGB24 -> RGB32
50   BGR32 -> BGR24 & RGB32 -> RGB24
51   BGR24 -> YV12
52 */
53
54 #include <inttypes.h>
55 #include <string.h>
56 #include <math.h>
57 #include <stdio.h>
58 #include "config.h"
59 #include <assert.h>
60 #include "swscale.h"
61 #include "swscale_internal.h"
62 #include "rgb2rgb.h"
63 #include "libavutil/intreadwrite.h"
64 #include "libavutil/x86_cpu.h"
65 #include "libavutil/avutil.h"
66 #include "libavutil/bswap.h"
67 #include "libavutil/pixdesc.h"
68
69 #undef MOVNTQ
70 #undef PAVGB
71
72 //#undef HAVE_MMX2
73 //#define HAVE_AMD3DNOW
74 //#undef HAVE_MMX
75 //#undef ARCH_X86
76 #define DITHER1XBPP
77
78 #define FAST_BGR2YV12 // use 7 bit coefficients instead of 15 bit
79
80 #ifdef M_PI
81 #define PI M_PI
82 #else
83 #define PI 3.14159265358979323846
84 #endif
85
86 #define isPacked(x)         (       \
87            (x)==PIX_FMT_PAL8        \
88         || (x)==PIX_FMT_YUYV422     \
89         || (x)==PIX_FMT_UYVY422     \
90         || isAnyRGB(x)              \
91     )
92
93 #define RGB2YUV_SHIFT 15
94 #define BY ( (int)(0.114*219/255*(1<<RGB2YUV_SHIFT)+0.5))
95 #define BV (-(int)(0.081*224/255*(1<<RGB2YUV_SHIFT)+0.5))
96 #define BU ( (int)(0.500*224/255*(1<<RGB2YUV_SHIFT)+0.5))
97 #define GY ( (int)(0.587*219/255*(1<<RGB2YUV_SHIFT)+0.5))
98 #define GV (-(int)(0.419*224/255*(1<<RGB2YUV_SHIFT)+0.5))
99 #define GU (-(int)(0.331*224/255*(1<<RGB2YUV_SHIFT)+0.5))
100 #define RY ( (int)(0.299*219/255*(1<<RGB2YUV_SHIFT)+0.5))
101 #define RV ( (int)(0.500*224/255*(1<<RGB2YUV_SHIFT)+0.5))
102 #define RU (-(int)(0.169*224/255*(1<<RGB2YUV_SHIFT)+0.5))
103
104 static const double rgb2yuv_table[8][9]={
105     {0.7152, 0.0722, 0.2126, -0.386, 0.5, -0.115, -0.454, -0.046, 0.5},
106     {0.7152, 0.0722, 0.2126, -0.386, 0.5, -0.115, -0.454, -0.046, 0.5},
107     {0.587 , 0.114 , 0.299 , -0.331, 0.5, -0.169, -0.419, -0.081, 0.5},
108     {0.587 , 0.114 , 0.299 , -0.331, 0.5, -0.169, -0.419, -0.081, 0.5},
109     {0.59  , 0.11  , 0.30  , -0.331, 0.5, -0.169, -0.421, -0.079, 0.5}, //FCC
110     {0.587 , 0.114 , 0.299 , -0.331, 0.5, -0.169, -0.419, -0.081, 0.5},
111     {0.587 , 0.114 , 0.299 , -0.331, 0.5, -0.169, -0.419, -0.081, 0.5}, //SMPTE 170M
112     {0.701 , 0.087 , 0.212 , -0.384, 0.5  -0.116, -0.445, -0.055, 0.5}, //SMPTE 240M
113 };
114
115 /*
116 NOTES
117 Special versions: fast Y 1:1 scaling (no interpolation in y direction)
118
119 TODO
120 more intelligent misalignment avoidance for the horizontal scaler
121 write special vertical cubic upscale version
122 optimize C code (YV12 / minmax)
123 add support for packed pixel YUV input & output
124 add support for Y8 output
125 optimize BGR24 & BGR32
126 add BGR4 output support
127 write special BGR->BGR scaler
128 */
129
130 #if ARCH_X86
131 DECLARE_ASM_CONST(8, uint64_t, bF8)=       0xF8F8F8F8F8F8F8F8LL;
132 DECLARE_ASM_CONST(8, uint64_t, bFC)=       0xFCFCFCFCFCFCFCFCLL;
133 DECLARE_ASM_CONST(8, uint64_t, w10)=       0x0010001000100010LL;
134 DECLARE_ASM_CONST(8, uint64_t, w02)=       0x0002000200020002LL;
135 DECLARE_ASM_CONST(8, uint64_t, bm00001111)=0x00000000FFFFFFFFLL;
136 DECLARE_ASM_CONST(8, uint64_t, bm00000111)=0x0000000000FFFFFFLL;
137 DECLARE_ASM_CONST(8, uint64_t, bm11111000)=0xFFFFFFFFFF000000LL;
138 DECLARE_ASM_CONST(8, uint64_t, bm01010101)=0x00FF00FF00FF00FFLL;
139
140 const DECLARE_ALIGNED(8, uint64_t, ff_dither4)[2] = {
141         0x0103010301030103LL,
142         0x0200020002000200LL,};
143
144 const DECLARE_ALIGNED(8, uint64_t, ff_dither8)[2] = {
145         0x0602060206020602LL,
146         0x0004000400040004LL,};
147
148 DECLARE_ASM_CONST(8, uint64_t, b16Mask)=   0x001F001F001F001FLL;
149 DECLARE_ASM_CONST(8, uint64_t, g16Mask)=   0x07E007E007E007E0LL;
150 DECLARE_ASM_CONST(8, uint64_t, r16Mask)=   0xF800F800F800F800LL;
151 DECLARE_ASM_CONST(8, uint64_t, b15Mask)=   0x001F001F001F001FLL;
152 DECLARE_ASM_CONST(8, uint64_t, g15Mask)=   0x03E003E003E003E0LL;
153 DECLARE_ASM_CONST(8, uint64_t, r15Mask)=   0x7C007C007C007C00LL;
154
155 DECLARE_ALIGNED(8, const uint64_t, ff_M24A)         = 0x00FF0000FF0000FFLL;
156 DECLARE_ALIGNED(8, const uint64_t, ff_M24B)         = 0xFF0000FF0000FF00LL;
157 DECLARE_ALIGNED(8, const uint64_t, ff_M24C)         = 0x0000FF0000FF0000LL;
158
159 #ifdef FAST_BGR2YV12
160 DECLARE_ALIGNED(8, const uint64_t, ff_bgr2YCoeff)   = 0x000000210041000DULL;
161 DECLARE_ALIGNED(8, const uint64_t, ff_bgr2UCoeff)   = 0x0000FFEEFFDC0038ULL;
162 DECLARE_ALIGNED(8, const uint64_t, ff_bgr2VCoeff)   = 0x00000038FFD2FFF8ULL;
163 #else
164 DECLARE_ALIGNED(8, const uint64_t, ff_bgr2YCoeff)   = 0x000020E540830C8BULL;
165 DECLARE_ALIGNED(8, const uint64_t, ff_bgr2UCoeff)   = 0x0000ED0FDAC23831ULL;
166 DECLARE_ALIGNED(8, const uint64_t, ff_bgr2VCoeff)   = 0x00003831D0E6F6EAULL;
167 #endif /* FAST_BGR2YV12 */
168 DECLARE_ALIGNED(8, const uint64_t, ff_bgr2YOffset)  = 0x1010101010101010ULL;
169 DECLARE_ALIGNED(8, const uint64_t, ff_bgr2UVOffset) = 0x8080808080808080ULL;
170 DECLARE_ALIGNED(8, const uint64_t, ff_w1111)        = 0x0001000100010001ULL;
171
172 DECLARE_ASM_CONST(8, uint64_t, ff_bgr24toY1Coeff) = 0x0C88000040870C88ULL;
173 DECLARE_ASM_CONST(8, uint64_t, ff_bgr24toY2Coeff) = 0x20DE4087000020DEULL;
174 DECLARE_ASM_CONST(8, uint64_t, ff_rgb24toY1Coeff) = 0x20DE0000408720DEULL;
175 DECLARE_ASM_CONST(8, uint64_t, ff_rgb24toY2Coeff) = 0x0C88408700000C88ULL;
176 DECLARE_ASM_CONST(8, uint64_t, ff_bgr24toYOffset) = 0x0008400000084000ULL;
177
178 DECLARE_ASM_CONST(8, uint64_t, ff_bgr24toUV)[2][4] = {
179     {0x38380000DAC83838ULL, 0xECFFDAC80000ECFFULL, 0xF6E40000D0E3F6E4ULL, 0x3838D0E300003838ULL},
180     {0xECFF0000DAC8ECFFULL, 0x3838DAC800003838ULL, 0x38380000D0E33838ULL, 0xF6E4D0E30000F6E4ULL},
181 };
182
183 DECLARE_ASM_CONST(8, uint64_t, ff_bgr24toUVOffset)= 0x0040400000404000ULL;
184
185 #endif /* ARCH_X86 */
186
187 DECLARE_ALIGNED(8, static const uint8_t, dither_2x2_4)[2][8]={
188 {  1,   3,   1,   3,   1,   3,   1,   3, },
189 {  2,   0,   2,   0,   2,   0,   2,   0, },
190 };
191
192 DECLARE_ALIGNED(8, static const uint8_t, dither_2x2_8)[2][8]={
193 {  6,   2,   6,   2,   6,   2,   6,   2, },
194 {  0,   4,   0,   4,   0,   4,   0,   4, },
195 };
196
197 DECLARE_ALIGNED(8, const uint8_t, dither_4x4_16)[4][8]={
198 {  8,   4,  11,   7,   8,   4,  11,   7, },
199 {  2,  14,   1,  13,   2,  14,   1,  13, },
200 { 10,   6,   9,   5,  10,   6,   9,   5, },
201 {  0,  12,   3,  15,   0,  12,   3,  15, },
202 };
203
204 DECLARE_ALIGNED(8, const uint8_t, dither_8x8_32)[8][8]={
205 { 17,   9,  23,  15,  16,   8,  22,  14, },
206 {  5,  29,   3,  27,   4,  28,   2,  26, },
207 { 21,  13,  19,  11,  20,  12,  18,  10, },
208 {  0,  24,   6,  30,   1,  25,   7,  31, },
209 { 16,   8,  22,  14,  17,   9,  23,  15, },
210 {  4,  28,   2,  26,   5,  29,   3,  27, },
211 { 20,  12,  18,  10,  21,  13,  19,  11, },
212 {  1,  25,   7,  31,   0,  24,   6,  30, },
213 };
214
215 DECLARE_ALIGNED(8, const uint8_t, dither_8x8_73)[8][8]={
216 {  0,  55,  14,  68,   3,  58,  17,  72, },
217 { 37,  18,  50,  32,  40,  22,  54,  35, },
218 {  9,  64,   5,  59,  13,  67,   8,  63, },
219 { 46,  27,  41,  23,  49,  31,  44,  26, },
220 {  2,  57,  16,  71,   1,  56,  15,  70, },
221 { 39,  21,  52,  34,  38,  19,  51,  33, },
222 { 11,  66,   7,  62,  10,  65,   6,  60, },
223 { 48,  30,  43,  25,  47,  29,  42,  24, },
224 };
225
226 #if 1
227 DECLARE_ALIGNED(8, const uint8_t, dither_8x8_220)[8][8]={
228 {117,  62, 158, 103, 113,  58, 155, 100, },
229 { 34, 199,  21, 186,  31, 196,  17, 182, },
230 {144,  89, 131,  76, 141,  86, 127,  72, },
231 {  0, 165,  41, 206,  10, 175,  52, 217, },
232 {110,  55, 151,  96, 120,  65, 162, 107, },
233 { 28, 193,  14, 179,  38, 203,  24, 189, },
234 {138,  83, 124,  69, 148,  93, 134,  79, },
235 {  7, 172,  48, 213,   3, 168,  45, 210, },
236 };
237 #elif 1
238 // tries to correct a gamma of 1.5
239 DECLARE_ALIGNED(8, const uint8_t, dither_8x8_220)[8][8]={
240 {  0, 143,  18, 200,   2, 156,  25, 215, },
241 { 78,  28, 125,  64,  89,  36, 138,  74, },
242 { 10, 180,   3, 161,  16, 195,   8, 175, },
243 {109,  51,  93,  38, 121,  60, 105,  47, },
244 {  1, 152,  23, 210,   0, 147,  20, 205, },
245 { 85,  33, 134,  71,  81,  30, 130,  67, },
246 { 14, 190,   6, 171,  12, 185,   5, 166, },
247 {117,  57, 101,  44, 113,  54,  97,  41, },
248 };
249 #elif 1
250 // tries to correct a gamma of 2.0
251 DECLARE_ALIGNED(8, const uint8_t, dither_8x8_220)[8][8]={
252 {  0, 124,   8, 193,   0, 140,  12, 213, },
253 { 55,  14, 104,  42,  66,  19, 119,  52, },
254 {  3, 168,   1, 145,   6, 187,   3, 162, },
255 { 86,  31,  70,  21,  99,  39,  82,  28, },
256 {  0, 134,  11, 206,   0, 129,   9, 200, },
257 { 62,  17, 114,  48,  58,  16, 109,  45, },
258 {  5, 181,   2, 157,   4, 175,   1, 151, },
259 { 95,  36,  78,  26,  90,  34,  74,  24, },
260 };
261 #else
262 // tries to correct a gamma of 2.5
263 DECLARE_ALIGNED(8, const uint8_t, dither_8x8_220)[8][8]={
264 {  0, 107,   3, 187,   0, 125,   6, 212, },
265 { 39,   7,  86,  28,  49,  11, 102,  36, },
266 {  1, 158,   0, 131,   3, 180,   1, 151, },
267 { 68,  19,  52,  12,  81,  25,  64,  17, },
268 {  0, 119,   5, 203,   0, 113,   4, 195, },
269 { 45,   9,  96,  33,  42,   8,  91,  30, },
270 {  2, 172,   1, 144,   2, 165,   0, 137, },
271 { 77,  23,  60,  15,  72,  21,  56,  14, },
272 };
273 #endif
274
275 static av_always_inline void yuv2yuvX16inC_template(const int16_t *lumFilter, const int16_t **lumSrc, int lumFilterSize,
276                                                     const int16_t *chrFilter, const int16_t **chrSrc, int chrFilterSize,
277                                                     const int16_t **alpSrc, uint16_t *dest, uint16_t *uDest, uint16_t *vDest, uint16_t *aDest,
278                                                     int dstW, int chrDstW, int big_endian)
279 {
280     //FIXME Optimize (just quickly written not optimized..)
281     int i;
282
283     for (i = 0; i < dstW; i++) {
284         int val = 1 << 10;
285         int j;
286
287         for (j = 0; j < lumFilterSize; j++)
288             val += lumSrc[j][i] * lumFilter[j];
289
290         if (big_endian) {
291             AV_WB16(&dest[i], av_clip_uint16(val >> 11));
292         } else {
293             AV_WL16(&dest[i], av_clip_uint16(val >> 11));
294         }
295     }
296
297     if (uDest) {
298         for (i = 0; i < chrDstW; i++) {
299             int u = 1 << 10;
300             int v = 1 << 10;
301             int j;
302
303             for (j = 0; j < chrFilterSize; j++) {
304                 u += chrSrc[j][i       ] * chrFilter[j];
305                 v += chrSrc[j][i + VOFW] * chrFilter[j];
306             }
307
308             if (big_endian) {
309                 AV_WB16(&uDest[i], av_clip_uint16(u >> 11));
310                 AV_WB16(&vDest[i], av_clip_uint16(v >> 11));
311             } else {
312                 AV_WL16(&uDest[i], av_clip_uint16(u >> 11));
313                 AV_WL16(&vDest[i], av_clip_uint16(v >> 11));
314             }
315         }
316     }
317
318     if (CONFIG_SWSCALE_ALPHA && aDest) {
319         for (i = 0; i < dstW; i++) {
320             int val = 1 << 10;
321             int j;
322
323             for (j = 0; j < lumFilterSize; j++)
324                 val += alpSrc[j][i] * lumFilter[j];
325
326             if (big_endian) {
327                 AV_WB16(&aDest[i], av_clip_uint16(val >> 11));
328             } else {
329                 AV_WL16(&aDest[i], av_clip_uint16(val >> 11));
330             }
331         }
332     }
333 }
334
335 static inline void yuv2yuvX16inC(const int16_t *lumFilter, const int16_t **lumSrc, int lumFilterSize,
336                                  const int16_t *chrFilter, const int16_t **chrSrc, int chrFilterSize,
337                                  const int16_t **alpSrc, uint16_t *dest, uint16_t *uDest, uint16_t *vDest, uint16_t *aDest, int dstW, int chrDstW,
338                                  enum PixelFormat dstFormat)
339 {
340     if (isBE(dstFormat)) {
341         yuv2yuvX16inC_template(lumFilter, lumSrc, lumFilterSize,
342                                chrFilter, chrSrc, chrFilterSize,
343                                alpSrc,
344                                dest, uDest, vDest, aDest,
345                                dstW, chrDstW, 1);
346     } else {
347         yuv2yuvX16inC_template(lumFilter, lumSrc, lumFilterSize,
348                                chrFilter, chrSrc, chrFilterSize,
349                                alpSrc,
350                                dest, uDest, vDest, aDest,
351                                dstW, chrDstW, 0);
352     }
353 }
354
355 static inline void yuv2yuvXinC(const int16_t *lumFilter, const int16_t **lumSrc, int lumFilterSize,
356                                const int16_t *chrFilter, const int16_t **chrSrc, int chrFilterSize,
357                                const int16_t **alpSrc, uint8_t *dest, uint8_t *uDest, uint8_t *vDest, uint8_t *aDest, int dstW, int chrDstW)
358 {
359     //FIXME Optimize (just quickly written not optimized..)
360     int i;
361     for (i=0; i<dstW; i++) {
362         int val=1<<18;
363         int j;
364         for (j=0; j<lumFilterSize; j++)
365             val += lumSrc[j][i] * lumFilter[j];
366
367         dest[i]= av_clip_uint8(val>>19);
368     }
369
370     if (uDest)
371         for (i=0; i<chrDstW; i++) {
372             int u=1<<18;
373             int v=1<<18;
374             int j;
375             for (j=0; j<chrFilterSize; j++) {
376                 u += chrSrc[j][i] * chrFilter[j];
377                 v += chrSrc[j][i + VOFW] * chrFilter[j];
378             }
379
380             uDest[i]= av_clip_uint8(u>>19);
381             vDest[i]= av_clip_uint8(v>>19);
382         }
383
384     if (CONFIG_SWSCALE_ALPHA && aDest)
385         for (i=0; i<dstW; i++) {
386             int val=1<<18;
387             int j;
388             for (j=0; j<lumFilterSize; j++)
389                 val += alpSrc[j][i] * lumFilter[j];
390
391             aDest[i]= av_clip_uint8(val>>19);
392         }
393
394 }
395
396 static inline void yuv2nv12XinC(const int16_t *lumFilter, const int16_t **lumSrc, int lumFilterSize,
397                                 const int16_t *chrFilter, const int16_t **chrSrc, int chrFilterSize,
398                                 uint8_t *dest, uint8_t *uDest, int dstW, int chrDstW, int dstFormat)
399 {
400     //FIXME Optimize (just quickly written not optimized..)
401     int i;
402     for (i=0; i<dstW; i++) {
403         int val=1<<18;
404         int j;
405         for (j=0; j<lumFilterSize; j++)
406             val += lumSrc[j][i] * lumFilter[j];
407
408         dest[i]= av_clip_uint8(val>>19);
409     }
410
411     if (!uDest)
412         return;
413
414     if (dstFormat == PIX_FMT_NV12)
415         for (i=0; i<chrDstW; i++) {
416             int u=1<<18;
417             int v=1<<18;
418             int j;
419             for (j=0; j<chrFilterSize; j++) {
420                 u += chrSrc[j][i] * chrFilter[j];
421                 v += chrSrc[j][i + VOFW] * chrFilter[j];
422             }
423
424             uDest[2*i]= av_clip_uint8(u>>19);
425             uDest[2*i+1]= av_clip_uint8(v>>19);
426         }
427     else
428         for (i=0; i<chrDstW; i++) {
429             int u=1<<18;
430             int v=1<<18;
431             int j;
432             for (j=0; j<chrFilterSize; j++) {
433                 u += chrSrc[j][i] * chrFilter[j];
434                 v += chrSrc[j][i + VOFW] * chrFilter[j];
435             }
436
437             uDest[2*i]= av_clip_uint8(v>>19);
438             uDest[2*i+1]= av_clip_uint8(u>>19);
439         }
440 }
441
442 #define YSCALE_YUV_2_PACKEDX_NOCLIP_C(type,alpha) \
443     for (i=0; i<(dstW>>1); i++) {\
444         int j;\
445         int Y1 = 1<<18;\
446         int Y2 = 1<<18;\
447         int U  = 1<<18;\
448         int V  = 1<<18;\
449         int av_unused A1, A2;\
450         type av_unused *r, *b, *g;\
451         const int i2= 2*i;\
452         \
453         for (j=0; j<lumFilterSize; j++) {\
454             Y1 += lumSrc[j][i2] * lumFilter[j];\
455             Y2 += lumSrc[j][i2+1] * lumFilter[j];\
456         }\
457         for (j=0; j<chrFilterSize; j++) {\
458             U += chrSrc[j][i] * chrFilter[j];\
459             V += chrSrc[j][i+VOFW] * chrFilter[j];\
460         }\
461         Y1>>=19;\
462         Y2>>=19;\
463         U >>=19;\
464         V >>=19;\
465         if (alpha) {\
466             A1 = 1<<18;\
467             A2 = 1<<18;\
468             for (j=0; j<lumFilterSize; j++) {\
469                 A1 += alpSrc[j][i2  ] * lumFilter[j];\
470                 A2 += alpSrc[j][i2+1] * lumFilter[j];\
471             }\
472             A1>>=19;\
473             A2>>=19;\
474         }
475
476 #define YSCALE_YUV_2_PACKEDX_C(type,alpha) \
477         YSCALE_YUV_2_PACKEDX_NOCLIP_C(type,alpha)\
478         if ((Y1|Y2|U|V)&256) {\
479             if (Y1>255)   Y1=255; \
480             else if (Y1<0)Y1=0;   \
481             if (Y2>255)   Y2=255; \
482             else if (Y2<0)Y2=0;   \
483             if (U>255)    U=255;  \
484             else if (U<0) U=0;    \
485             if (V>255)    V=255;  \
486             else if (V<0) V=0;    \
487         }\
488         if (alpha && ((A1|A2)&256)) {\
489             A1=av_clip_uint8(A1);\
490             A2=av_clip_uint8(A2);\
491         }
492
493 #define YSCALE_YUV_2_PACKEDX_FULL_C(rnd,alpha) \
494     for (i=0; i<dstW; i++) {\
495         int j;\
496         int Y = 0;\
497         int U = -128<<19;\
498         int V = -128<<19;\
499         int av_unused A;\
500         int R,G,B;\
501         \
502         for (j=0; j<lumFilterSize; j++) {\
503             Y += lumSrc[j][i     ] * lumFilter[j];\
504         }\
505         for (j=0; j<chrFilterSize; j++) {\
506             U += chrSrc[j][i     ] * chrFilter[j];\
507             V += chrSrc[j][i+VOFW] * chrFilter[j];\
508         }\
509         Y >>=10;\
510         U >>=10;\
511         V >>=10;\
512         if (alpha) {\
513             A = rnd;\
514             for (j=0; j<lumFilterSize; j++)\
515                 A += alpSrc[j][i     ] * lumFilter[j];\
516             A >>=19;\
517             if (A&256)\
518                 A = av_clip_uint8(A);\
519         }
520
521 #define YSCALE_YUV_2_RGBX_FULL_C(rnd,alpha) \
522     YSCALE_YUV_2_PACKEDX_FULL_C(rnd>>3,alpha)\
523         Y-= c->yuv2rgb_y_offset;\
524         Y*= c->yuv2rgb_y_coeff;\
525         Y+= rnd;\
526         R= Y + V*c->yuv2rgb_v2r_coeff;\
527         G= Y + V*c->yuv2rgb_v2g_coeff + U*c->yuv2rgb_u2g_coeff;\
528         B= Y +                          U*c->yuv2rgb_u2b_coeff;\
529         if ((R|G|B)&(0xC0000000)) {\
530             if (R>=(256<<22))   R=(256<<22)-1; \
531             else if (R<0)R=0;   \
532             if (G>=(256<<22))   G=(256<<22)-1; \
533             else if (G<0)G=0;   \
534             if (B>=(256<<22))   B=(256<<22)-1; \
535             else if (B<0)B=0;   \
536         }
537
538 #define YSCALE_YUV_2_GRAY16_C \
539     for (i=0; i<(dstW>>1); i++) {\
540         int j;\
541         int Y1 = 1<<18;\
542         int Y2 = 1<<18;\
543         int U  = 1<<18;\
544         int V  = 1<<18;\
545         \
546         const int i2= 2*i;\
547         \
548         for (j=0; j<lumFilterSize; j++) {\
549             Y1 += lumSrc[j][i2] * lumFilter[j];\
550             Y2 += lumSrc[j][i2+1] * lumFilter[j];\
551         }\
552         Y1>>=11;\
553         Y2>>=11;\
554         if ((Y1|Y2|U|V)&65536) {\
555             if (Y1>65535)   Y1=65535; \
556             else if (Y1<0)Y1=0;   \
557             if (Y2>65535)   Y2=65535; \
558             else if (Y2<0)Y2=0;   \
559         }
560
561 #define YSCALE_YUV_2_RGBX_C(type,alpha) \
562     YSCALE_YUV_2_PACKEDX_C(type,alpha)  /* FIXME fix tables so that clipping is not needed and then use _NOCLIP*/\
563     r = (type *)c->table_rV[V];   \
564     g = (type *)(c->table_gU[U] + c->table_gV[V]); \
565     b = (type *)c->table_bU[U];
566
567 #define YSCALE_YUV_2_PACKED2_C(type,alpha)   \
568     for (i=0; i<(dstW>>1); i++) { \
569         const int i2= 2*i;       \
570         int Y1= (buf0[i2  ]*yalpha1+buf1[i2  ]*yalpha)>>19;           \
571         int Y2= (buf0[i2+1]*yalpha1+buf1[i2+1]*yalpha)>>19;           \
572         int U= (uvbuf0[i     ]*uvalpha1+uvbuf1[i     ]*uvalpha)>>19;  \
573         int V= (uvbuf0[i+VOFW]*uvalpha1+uvbuf1[i+VOFW]*uvalpha)>>19;  \
574         type av_unused *r, *b, *g;                                    \
575         int av_unused A1, A2;                                         \
576         if (alpha) {\
577             A1= (abuf0[i2  ]*yalpha1+abuf1[i2  ]*yalpha)>>19;         \
578             A2= (abuf0[i2+1]*yalpha1+abuf1[i2+1]*yalpha)>>19;         \
579         }
580
581 #define YSCALE_YUV_2_GRAY16_2_C   \
582     for (i=0; i<(dstW>>1); i++) { \
583         const int i2= 2*i;       \
584         int Y1= (buf0[i2  ]*yalpha1+buf1[i2  ]*yalpha)>>11;           \
585         int Y2= (buf0[i2+1]*yalpha1+buf1[i2+1]*yalpha)>>11;
586
587 #define YSCALE_YUV_2_RGB2_C(type,alpha) \
588     YSCALE_YUV_2_PACKED2_C(type,alpha)\
589     r = (type *)c->table_rV[V];\
590     g = (type *)(c->table_gU[U] + c->table_gV[V]);\
591     b = (type *)c->table_bU[U];
592
593 #define YSCALE_YUV_2_PACKED1_C(type,alpha) \
594     for (i=0; i<(dstW>>1); i++) {\
595         const int i2= 2*i;\
596         int Y1= buf0[i2  ]>>7;\
597         int Y2= buf0[i2+1]>>7;\
598         int U= (uvbuf1[i     ])>>7;\
599         int V= (uvbuf1[i+VOFW])>>7;\
600         type av_unused *r, *b, *g;\
601         int av_unused A1, A2;\
602         if (alpha) {\
603             A1= abuf0[i2  ]>>7;\
604             A2= abuf0[i2+1]>>7;\
605         }
606
607 #define YSCALE_YUV_2_GRAY16_1_C \
608     for (i=0; i<(dstW>>1); i++) {\
609         const int i2= 2*i;\
610         int Y1= buf0[i2  ]<<1;\
611         int Y2= buf0[i2+1]<<1;
612
613 #define YSCALE_YUV_2_RGB1_C(type,alpha) \
614     YSCALE_YUV_2_PACKED1_C(type,alpha)\
615     r = (type *)c->table_rV[V];\
616     g = (type *)(c->table_gU[U] + c->table_gV[V]);\
617     b = (type *)c->table_bU[U];
618
619 #define YSCALE_YUV_2_PACKED1B_C(type,alpha) \
620     for (i=0; i<(dstW>>1); i++) {\
621         const int i2= 2*i;\
622         int Y1= buf0[i2  ]>>7;\
623         int Y2= buf0[i2+1]>>7;\
624         int U= (uvbuf0[i     ] + uvbuf1[i     ])>>8;\
625         int V= (uvbuf0[i+VOFW] + uvbuf1[i+VOFW])>>8;\
626         type av_unused *r, *b, *g;\
627         int av_unused A1, A2;\
628         if (alpha) {\
629             A1= abuf0[i2  ]>>7;\
630             A2= abuf0[i2+1]>>7;\
631         }
632
633 #define YSCALE_YUV_2_RGB1B_C(type,alpha) \
634     YSCALE_YUV_2_PACKED1B_C(type,alpha)\
635     r = (type *)c->table_rV[V];\
636     g = (type *)(c->table_gU[U] + c->table_gV[V]);\
637     b = (type *)c->table_bU[U];
638
639 #define YSCALE_YUV_2_MONO2_C \
640     const uint8_t * const d128=dither_8x8_220[y&7];\
641     uint8_t *g= c->table_gU[128] + c->table_gV[128];\
642     for (i=0; i<dstW-7; i+=8) {\
643         int acc;\
644         acc =       g[((buf0[i  ]*yalpha1+buf1[i  ]*yalpha)>>19) + d128[0]];\
645         acc+= acc + g[((buf0[i+1]*yalpha1+buf1[i+1]*yalpha)>>19) + d128[1]];\
646         acc+= acc + g[((buf0[i+2]*yalpha1+buf1[i+2]*yalpha)>>19) + d128[2]];\
647         acc+= acc + g[((buf0[i+3]*yalpha1+buf1[i+3]*yalpha)>>19) + d128[3]];\
648         acc+= acc + g[((buf0[i+4]*yalpha1+buf1[i+4]*yalpha)>>19) + d128[4]];\
649         acc+= acc + g[((buf0[i+5]*yalpha1+buf1[i+5]*yalpha)>>19) + d128[5]];\
650         acc+= acc + g[((buf0[i+6]*yalpha1+buf1[i+6]*yalpha)>>19) + d128[6]];\
651         acc+= acc + g[((buf0[i+7]*yalpha1+buf1[i+7]*yalpha)>>19) + d128[7]];\
652         ((uint8_t*)dest)[0]= c->dstFormat == PIX_FMT_MONOBLACK ? acc : ~acc;\
653         dest++;\
654     }
655
656 #define YSCALE_YUV_2_MONOX_C \
657     const uint8_t * const d128=dither_8x8_220[y&7];\
658     uint8_t *g= c->table_gU[128] + c->table_gV[128];\
659     int acc=0;\
660     for (i=0; i<dstW-1; i+=2) {\
661         int j;\
662         int Y1=1<<18;\
663         int Y2=1<<18;\
664 \
665         for (j=0; j<lumFilterSize; j++) {\
666             Y1 += lumSrc[j][i] * lumFilter[j];\
667             Y2 += lumSrc[j][i+1] * lumFilter[j];\
668         }\
669         Y1>>=19;\
670         Y2>>=19;\
671         if ((Y1|Y2)&256) {\
672             if (Y1>255)   Y1=255;\
673             else if (Y1<0)Y1=0;\
674             if (Y2>255)   Y2=255;\
675             else if (Y2<0)Y2=0;\
676         }\
677         acc+= acc + g[Y1+d128[(i+0)&7]];\
678         acc+= acc + g[Y2+d128[(i+1)&7]];\
679         if ((i&7)==6) {\
680             ((uint8_t*)dest)[0]= c->dstFormat == PIX_FMT_MONOBLACK ? acc : ~acc;\
681             dest++;\
682         }\
683     }
684
685 #define YSCALE_YUV_2_ANYRGB_C(func, func2, func_g16, func_monoblack)\
686     switch(c->dstFormat) {\
687     case PIX_FMT_RGB48BE:\
688     case PIX_FMT_RGB48LE:\
689         func(uint8_t,0)\
690             ((uint8_t*)dest)[ 0]= r[Y1];\
691             ((uint8_t*)dest)[ 1]= r[Y1];\
692             ((uint8_t*)dest)[ 2]= g[Y1];\
693             ((uint8_t*)dest)[ 3]= g[Y1];\
694             ((uint8_t*)dest)[ 4]= b[Y1];\
695             ((uint8_t*)dest)[ 5]= b[Y1];\
696             ((uint8_t*)dest)[ 6]= r[Y2];\
697             ((uint8_t*)dest)[ 7]= r[Y2];\
698             ((uint8_t*)dest)[ 8]= g[Y2];\
699             ((uint8_t*)dest)[ 9]= g[Y2];\
700             ((uint8_t*)dest)[10]= b[Y2];\
701             ((uint8_t*)dest)[11]= b[Y2];\
702             dest+=12;\
703         }\
704         break;\
705     case PIX_FMT_RGBA:\
706     case PIX_FMT_BGRA:\
707         if (CONFIG_SMALL) {\
708             int needAlpha = CONFIG_SWSCALE_ALPHA && c->alpPixBuf;\
709             func(uint32_t,needAlpha)\
710                 ((uint32_t*)dest)[i2+0]= r[Y1] + g[Y1] + b[Y1] + (needAlpha ? (A1<<24) : 0);\
711                 ((uint32_t*)dest)[i2+1]= r[Y2] + g[Y2] + b[Y2] + (needAlpha ? (A2<<24) : 0);\
712             }\
713         } else {\
714             if (CONFIG_SWSCALE_ALPHA && c->alpPixBuf) {\
715                 func(uint32_t,1)\
716                     ((uint32_t*)dest)[i2+0]= r[Y1] + g[Y1] + b[Y1] + (A1<<24);\
717                     ((uint32_t*)dest)[i2+1]= r[Y2] + g[Y2] + b[Y2] + (A2<<24);\
718                 }\
719             } else {\
720                 func(uint32_t,0)\
721                     ((uint32_t*)dest)[i2+0]= r[Y1] + g[Y1] + b[Y1];\
722                     ((uint32_t*)dest)[i2+1]= r[Y2] + g[Y2] + b[Y2];\
723                 }\
724             }\
725         }\
726         break;\
727     case PIX_FMT_ARGB:\
728     case PIX_FMT_ABGR:\
729         if (CONFIG_SMALL) {\
730             int needAlpha = CONFIG_SWSCALE_ALPHA && c->alpPixBuf;\
731             func(uint32_t,needAlpha)\
732                 ((uint32_t*)dest)[i2+0]= r[Y1] + g[Y1] + b[Y1] + (needAlpha ? A1 : 0);\
733                 ((uint32_t*)dest)[i2+1]= r[Y2] + g[Y2] + b[Y2] + (needAlpha ? A2 : 0);\
734             }\
735         } else {\
736             if (CONFIG_SWSCALE_ALPHA && c->alpPixBuf) {\
737                 func(uint32_t,1)\
738                     ((uint32_t*)dest)[i2+0]= r[Y1] + g[Y1] + b[Y1] + A1;\
739                     ((uint32_t*)dest)[i2+1]= r[Y2] + g[Y2] + b[Y2] + A2;\
740                 }\
741             } else {\
742                 func(uint32_t,0)\
743                     ((uint32_t*)dest)[i2+0]= r[Y1] + g[Y1] + b[Y1];\
744                     ((uint32_t*)dest)[i2+1]= r[Y2] + g[Y2] + b[Y2];\
745                 }\
746             }\
747         }                \
748         break;\
749     case PIX_FMT_RGB24:\
750         func(uint8_t,0)\
751             ((uint8_t*)dest)[0]= r[Y1];\
752             ((uint8_t*)dest)[1]= g[Y1];\
753             ((uint8_t*)dest)[2]= b[Y1];\
754             ((uint8_t*)dest)[3]= r[Y2];\
755             ((uint8_t*)dest)[4]= g[Y2];\
756             ((uint8_t*)dest)[5]= b[Y2];\
757             dest+=6;\
758         }\
759         break;\
760     case PIX_FMT_BGR24:\
761         func(uint8_t,0)\
762             ((uint8_t*)dest)[0]= b[Y1];\
763             ((uint8_t*)dest)[1]= g[Y1];\
764             ((uint8_t*)dest)[2]= r[Y1];\
765             ((uint8_t*)dest)[3]= b[Y2];\
766             ((uint8_t*)dest)[4]= g[Y2];\
767             ((uint8_t*)dest)[5]= r[Y2];\
768             dest+=6;\
769         }\
770         break;\
771     case PIX_FMT_RGB565BE:\
772     case PIX_FMT_RGB565LE:\
773     case PIX_FMT_BGR565BE:\
774     case PIX_FMT_BGR565LE:\
775         {\
776             const int dr1= dither_2x2_8[y&1    ][0];\
777             const int dg1= dither_2x2_4[y&1    ][0];\
778             const int db1= dither_2x2_8[(y&1)^1][0];\
779             const int dr2= dither_2x2_8[y&1    ][1];\
780             const int dg2= dither_2x2_4[y&1    ][1];\
781             const int db2= dither_2x2_8[(y&1)^1][1];\
782             func(uint16_t,0)\
783                 ((uint16_t*)dest)[i2+0]= r[Y1+dr1] + g[Y1+dg1] + b[Y1+db1];\
784                 ((uint16_t*)dest)[i2+1]= r[Y2+dr2] + g[Y2+dg2] + b[Y2+db2];\
785             }\
786         }\
787         break;\
788     case PIX_FMT_RGB555BE:\
789     case PIX_FMT_RGB555LE:\
790     case PIX_FMT_BGR555BE:\
791     case PIX_FMT_BGR555LE:\
792         {\
793             const int dr1= dither_2x2_8[y&1    ][0];\
794             const int dg1= dither_2x2_8[y&1    ][1];\
795             const int db1= dither_2x2_8[(y&1)^1][0];\
796             const int dr2= dither_2x2_8[y&1    ][1];\
797             const int dg2= dither_2x2_8[y&1    ][0];\
798             const int db2= dither_2x2_8[(y&1)^1][1];\
799             func(uint16_t,0)\
800                 ((uint16_t*)dest)[i2+0]= r[Y1+dr1] + g[Y1+dg1] + b[Y1+db1];\
801                 ((uint16_t*)dest)[i2+1]= r[Y2+dr2] + g[Y2+dg2] + b[Y2+db2];\
802             }\
803         }\
804         break;\
805     case PIX_FMT_RGB444BE:\
806     case PIX_FMT_RGB444LE:\
807     case PIX_FMT_BGR444BE:\
808     case PIX_FMT_BGR444LE:\
809         {\
810             const int dr1= dither_4x4_16[y&3    ][0];\
811             const int dg1= dither_4x4_16[y&3    ][1];\
812             const int db1= dither_4x4_16[(y&3)^3][0];\
813             const int dr2= dither_4x4_16[y&3    ][1];\
814             const int dg2= dither_4x4_16[y&3    ][0];\
815             const int db2= dither_4x4_16[(y&3)^3][1];\
816             func(uint16_t,0)\
817                 ((uint16_t*)dest)[i2+0]= r[Y1+dr1] + g[Y1+dg1] + b[Y1+db1];\
818                 ((uint16_t*)dest)[i2+1]= r[Y2+dr2] + g[Y2+dg2] + b[Y2+db2];\
819             }\
820         }\
821         break;\
822     case PIX_FMT_RGB8:\
823     case PIX_FMT_BGR8:\
824         {\
825             const uint8_t * const d64= dither_8x8_73[y&7];\
826             const uint8_t * const d32= dither_8x8_32[y&7];\
827             func(uint8_t,0)\
828                 ((uint8_t*)dest)[i2+0]= r[Y1+d32[(i2+0)&7]] + g[Y1+d32[(i2+0)&7]] + b[Y1+d64[(i2+0)&7]];\
829                 ((uint8_t*)dest)[i2+1]= r[Y2+d32[(i2+1)&7]] + g[Y2+d32[(i2+1)&7]] + b[Y2+d64[(i2+1)&7]];\
830             }\
831         }\
832         break;\
833     case PIX_FMT_RGB4:\
834     case PIX_FMT_BGR4:\
835         {\
836             const uint8_t * const d64= dither_8x8_73 [y&7];\
837             const uint8_t * const d128=dither_8x8_220[y&7];\
838             func(uint8_t,0)\
839                 ((uint8_t*)dest)[i]= r[Y1+d128[(i2+0)&7]] + g[Y1+d64[(i2+0)&7]] + b[Y1+d128[(i2+0)&7]]\
840                                  + ((r[Y2+d128[(i2+1)&7]] + g[Y2+d64[(i2+1)&7]] + b[Y2+d128[(i2+1)&7]])<<4);\
841             }\
842         }\
843         break;\
844     case PIX_FMT_RGB4_BYTE:\
845     case PIX_FMT_BGR4_BYTE:\
846         {\
847             const uint8_t * const d64= dither_8x8_73 [y&7];\
848             const uint8_t * const d128=dither_8x8_220[y&7];\
849             func(uint8_t,0)\
850                 ((uint8_t*)dest)[i2+0]= r[Y1+d128[(i2+0)&7]] + g[Y1+d64[(i2+0)&7]] + b[Y1+d128[(i2+0)&7]];\
851                 ((uint8_t*)dest)[i2+1]= r[Y2+d128[(i2+1)&7]] + g[Y2+d64[(i2+1)&7]] + b[Y2+d128[(i2+1)&7]];\
852             }\
853         }\
854         break;\
855     case PIX_FMT_MONOBLACK:\
856     case PIX_FMT_MONOWHITE:\
857         {\
858             func_monoblack\
859         }\
860         break;\
861     case PIX_FMT_YUYV422:\
862         func2\
863             ((uint8_t*)dest)[2*i2+0]= Y1;\
864             ((uint8_t*)dest)[2*i2+1]= U;\
865             ((uint8_t*)dest)[2*i2+2]= Y2;\
866             ((uint8_t*)dest)[2*i2+3]= V;\
867         }                \
868         break;\
869     case PIX_FMT_UYVY422:\
870         func2\
871             ((uint8_t*)dest)[2*i2+0]= U;\
872             ((uint8_t*)dest)[2*i2+1]= Y1;\
873             ((uint8_t*)dest)[2*i2+2]= V;\
874             ((uint8_t*)dest)[2*i2+3]= Y2;\
875         }                \
876         break;\
877     case PIX_FMT_GRAY16BE:\
878         func_g16\
879             ((uint8_t*)dest)[2*i2+0]= Y1>>8;\
880             ((uint8_t*)dest)[2*i2+1]= Y1;\
881             ((uint8_t*)dest)[2*i2+2]= Y2>>8;\
882             ((uint8_t*)dest)[2*i2+3]= Y2;\
883         }                \
884         break;\
885     case PIX_FMT_GRAY16LE:\
886         func_g16\
887             ((uint8_t*)dest)[2*i2+0]= Y1;\
888             ((uint8_t*)dest)[2*i2+1]= Y1>>8;\
889             ((uint8_t*)dest)[2*i2+2]= Y2;\
890             ((uint8_t*)dest)[2*i2+3]= Y2>>8;\
891         }                \
892         break;\
893     }
894
895 static inline void yuv2packedXinC(SwsContext *c, const int16_t *lumFilter, const int16_t **lumSrc, int lumFilterSize,
896                                   const int16_t *chrFilter, const int16_t **chrSrc, int chrFilterSize,
897                                   const int16_t **alpSrc, uint8_t *dest, int dstW, int y)
898 {
899     int i;
900     YSCALE_YUV_2_ANYRGB_C(YSCALE_YUV_2_RGBX_C, YSCALE_YUV_2_PACKEDX_C(void,0), YSCALE_YUV_2_GRAY16_C, YSCALE_YUV_2_MONOX_C)
901 }
902
903 static inline void yuv2rgbXinC_full(SwsContext *c, const int16_t *lumFilter, const int16_t **lumSrc, int lumFilterSize,
904                                     const int16_t *chrFilter, const int16_t **chrSrc, int chrFilterSize,
905                                     const int16_t **alpSrc, uint8_t *dest, int dstW, int y)
906 {
907     int i;
908     int step= c->dstFormatBpp/8;
909     int aidx= 3;
910
911     switch(c->dstFormat) {
912     case PIX_FMT_ARGB:
913         dest++;
914         aidx= 0;
915     case PIX_FMT_RGB24:
916         aidx--;
917     case PIX_FMT_RGBA:
918         if (CONFIG_SMALL) {
919             int needAlpha = CONFIG_SWSCALE_ALPHA && c->alpPixBuf;
920             YSCALE_YUV_2_RGBX_FULL_C(1<<21, needAlpha)
921                 dest[aidx]= needAlpha ? A : 255;
922                 dest[0]= R>>22;
923                 dest[1]= G>>22;
924                 dest[2]= B>>22;
925                 dest+= step;
926             }
927         } else {
928             if (CONFIG_SWSCALE_ALPHA && c->alpPixBuf) {
929                 YSCALE_YUV_2_RGBX_FULL_C(1<<21, 1)
930                     dest[aidx]= A;
931                     dest[0]= R>>22;
932                     dest[1]= G>>22;
933                     dest[2]= B>>22;
934                     dest+= step;
935                 }
936             } else {
937                 YSCALE_YUV_2_RGBX_FULL_C(1<<21, 0)
938                     dest[aidx]= 255;
939                     dest[0]= R>>22;
940                     dest[1]= G>>22;
941                     dest[2]= B>>22;
942                     dest+= step;
943                 }
944             }
945         }
946         break;
947     case PIX_FMT_ABGR:
948         dest++;
949         aidx= 0;
950     case PIX_FMT_BGR24:
951         aidx--;
952     case PIX_FMT_BGRA:
953         if (CONFIG_SMALL) {
954             int needAlpha = CONFIG_SWSCALE_ALPHA && c->alpPixBuf;
955             YSCALE_YUV_2_RGBX_FULL_C(1<<21, needAlpha)
956                 dest[aidx]= needAlpha ? A : 255;
957                 dest[0]= B>>22;
958                 dest[1]= G>>22;
959                 dest[2]= R>>22;
960                 dest+= step;
961             }
962         } else {
963             if (CONFIG_SWSCALE_ALPHA && c->alpPixBuf) {
964                 YSCALE_YUV_2_RGBX_FULL_C(1<<21, 1)
965                     dest[aidx]= A;
966                     dest[0]= B>>22;
967                     dest[1]= G>>22;
968                     dest[2]= R>>22;
969                     dest+= step;
970                 }
971             } else {
972                 YSCALE_YUV_2_RGBX_FULL_C(1<<21, 0)
973                     dest[aidx]= 255;
974                     dest[0]= B>>22;
975                     dest[1]= G>>22;
976                     dest[2]= R>>22;
977                     dest+= step;
978                 }
979             }
980         }
981         break;
982     default:
983         assert(0);
984     }
985 }
986
987 static void fillPlane(uint8_t* plane, int stride, int width, int height, int y, uint8_t val)
988 {
989     int i;
990     uint8_t *ptr = plane + stride*y;
991     for (i=0; i<height; i++) {
992         memset(ptr, val, width);
993         ptr += stride;
994     }
995 }
996
997 static inline void rgb48ToY(uint8_t *dst, const uint8_t *src, int width,
998                             uint32_t *unused)
999 {
1000     int i;
1001     for (i = 0; i < width; i++) {
1002         int r = src[i*6+0];
1003         int g = src[i*6+2];
1004         int b = src[i*6+4];
1005
1006         dst[i] = (RY*r + GY*g + BY*b + (33<<(RGB2YUV_SHIFT-1))) >> RGB2YUV_SHIFT;
1007     }
1008 }
1009
1010 static inline void rgb48ToUV(uint8_t *dstU, uint8_t *dstV,
1011                              const uint8_t *src1, const uint8_t *src2,
1012                              int width, uint32_t *unused)
1013 {
1014     int i;
1015     assert(src1==src2);
1016     for (i = 0; i < width; i++) {
1017         int r = src1[6*i + 0];
1018         int g = src1[6*i + 2];
1019         int b = src1[6*i + 4];
1020
1021         dstU[i] = (RU*r + GU*g + BU*b + (257<<(RGB2YUV_SHIFT-1))) >> RGB2YUV_SHIFT;
1022         dstV[i] = (RV*r + GV*g + BV*b + (257<<(RGB2YUV_SHIFT-1))) >> RGB2YUV_SHIFT;
1023     }
1024 }
1025
1026 static inline void rgb48ToUV_half(uint8_t *dstU, uint8_t *dstV,
1027                                   const uint8_t *src1, const uint8_t *src2,
1028                                   int width, uint32_t *unused)
1029 {
1030     int i;
1031     assert(src1==src2);
1032     for (i = 0; i < width; i++) {
1033         int r= src1[12*i + 0] + src1[12*i + 6];
1034         int g= src1[12*i + 2] + src1[12*i + 8];
1035         int b= src1[12*i + 4] + src1[12*i + 10];
1036
1037         dstU[i]= (RU*r + GU*g + BU*b + (257<<RGB2YUV_SHIFT)) >> (RGB2YUV_SHIFT+1);
1038         dstV[i]= (RV*r + GV*g + BV*b + (257<<RGB2YUV_SHIFT)) >> (RGB2YUV_SHIFT+1);
1039     }
1040 }
1041
1042 #define BGR2Y(type, name, shr, shg, shb, maskr, maskg, maskb, RY, GY, BY, S)\
1043 static inline void name(uint8_t *dst, const uint8_t *src, long width, uint32_t *unused)\
1044 {\
1045     int i;\
1046     for (i=0; i<width; i++) {\
1047         int b= (((const type*)src)[i]>>shb)&maskb;\
1048         int g= (((const type*)src)[i]>>shg)&maskg;\
1049         int r= (((const type*)src)[i]>>shr)&maskr;\
1050 \
1051         dst[i]= (((RY)*r + (GY)*g + (BY)*b + (33<<((S)-1)))>>(S));\
1052     }\
1053 }
1054
1055 BGR2Y(uint32_t, bgr32ToY,16, 0, 0, 0x00FF, 0xFF00, 0x00FF, RY<< 8, GY   , BY<< 8, RGB2YUV_SHIFT+8)
1056 BGR2Y(uint32_t, rgb32ToY, 0, 0,16, 0x00FF, 0xFF00, 0x00FF, RY<< 8, GY   , BY<< 8, RGB2YUV_SHIFT+8)
1057 BGR2Y(uint16_t, bgr16ToY, 0, 0, 0, 0x001F, 0x07E0, 0xF800, RY<<11, GY<<5, BY    , RGB2YUV_SHIFT+8)
1058 BGR2Y(uint16_t, bgr15ToY, 0, 0, 0, 0x001F, 0x03E0, 0x7C00, RY<<10, GY<<5, BY    , RGB2YUV_SHIFT+7)
1059 BGR2Y(uint16_t, rgb16ToY, 0, 0, 0, 0xF800, 0x07E0, 0x001F, RY    , GY<<5, BY<<11, RGB2YUV_SHIFT+8)
1060 BGR2Y(uint16_t, rgb15ToY, 0, 0, 0, 0x7C00, 0x03E0, 0x001F, RY    , GY<<5, BY<<10, RGB2YUV_SHIFT+7)
1061
1062 static inline void abgrToA(uint8_t *dst, const uint8_t *src, long width, uint32_t *unused)
1063 {
1064     int i;
1065     for (i=0; i<width; i++) {
1066         dst[i]= src[4*i];
1067     }
1068 }
1069
1070 #define BGR2UV(type, name, shr, shg, shb, maska, maskr, maskg, maskb, RU, GU, BU, RV, GV, BV, S)\
1071 static inline void name(uint8_t *dstU, uint8_t *dstV, const uint8_t *src, const uint8_t *dummy, long width, uint32_t *unused)\
1072 {\
1073     int i;\
1074     for (i=0; i<width; i++) {\
1075         int b= (((const type*)src)[i]&maskb)>>shb;\
1076         int g= (((const type*)src)[i]&maskg)>>shg;\
1077         int r= (((const type*)src)[i]&maskr)>>shr;\
1078 \
1079         dstU[i]= ((RU)*r + (GU)*g + (BU)*b + (257<<((S)-1)))>>(S);\
1080         dstV[i]= ((RV)*r + (GV)*g + (BV)*b + (257<<((S)-1)))>>(S);\
1081     }\
1082 }\
1083 static inline void name ## _half(uint8_t *dstU, uint8_t *dstV, const uint8_t *src, const uint8_t *dummy, long width, uint32_t *unused)\
1084 {\
1085     int i;\
1086     for (i=0; i<width; i++) {\
1087         int pix0= ((const type*)src)[2*i+0];\
1088         int pix1= ((const type*)src)[2*i+1];\
1089         int g= (pix0&~(maskr|maskb))+(pix1&~(maskr|maskb));\
1090         int b= ((pix0+pix1-g)&(maskb|(2*maskb)))>>shb;\
1091         int r= ((pix0+pix1-g)&(maskr|(2*maskr)))>>shr;\
1092         g&= maskg|(2*maskg);\
1093 \
1094         g>>=shg;\
1095 \
1096         dstU[i]= ((RU)*r + (GU)*g + (BU)*b + (257<<(S)))>>((S)+1);\
1097         dstV[i]= ((RV)*r + (GV)*g + (BV)*b + (257<<(S)))>>((S)+1);\
1098     }\
1099 }
1100
1101 BGR2UV(uint32_t, bgr32ToUV,16, 0, 0, 0xFF000000, 0xFF0000, 0xFF00,   0x00FF, RU<< 8, GU   , BU<< 8, RV<< 8, GV   , BV<< 8, RGB2YUV_SHIFT+8)
1102 BGR2UV(uint32_t, rgb32ToUV, 0, 0,16, 0xFF000000,   0x00FF, 0xFF00, 0xFF0000, RU<< 8, GU   , BU<< 8, RV<< 8, GV   , BV<< 8, RGB2YUV_SHIFT+8)
1103 BGR2UV(uint16_t, bgr16ToUV, 0, 0, 0,          0,   0x001F, 0x07E0,   0xF800, RU<<11, GU<<5, BU    , RV<<11, GV<<5, BV    , RGB2YUV_SHIFT+8)
1104 BGR2UV(uint16_t, bgr15ToUV, 0, 0, 0,          0,   0x001F, 0x03E0,   0x7C00, RU<<10, GU<<5, BU    , RV<<10, GV<<5, BV    , RGB2YUV_SHIFT+7)
1105 BGR2UV(uint16_t, rgb16ToUV, 0, 0, 0,          0,   0xF800, 0x07E0,   0x001F, RU    , GU<<5, BU<<11, RV    , GV<<5, BV<<11, RGB2YUV_SHIFT+8)
1106 BGR2UV(uint16_t, rgb15ToUV, 0, 0, 0,          0,   0x7C00, 0x03E0,   0x001F, RU    , GU<<5, BU<<10, RV    , GV<<5, BV<<10, RGB2YUV_SHIFT+7)
1107
1108 static inline void palToY(uint8_t *dst, const uint8_t *src, long width, uint32_t *pal)
1109 {
1110     int i;
1111     for (i=0; i<width; i++) {
1112         int d= src[i];
1113
1114         dst[i]= pal[d] & 0xFF;
1115     }
1116 }
1117
1118 static inline void palToUV(uint8_t *dstU, uint8_t *dstV,
1119                            const uint8_t *src1, const uint8_t *src2,
1120                            long width, uint32_t *pal)
1121 {
1122     int i;
1123     assert(src1 == src2);
1124     for (i=0; i<width; i++) {
1125         int p= pal[src1[i]];
1126
1127         dstU[i]= p>>8;
1128         dstV[i]= p>>16;
1129     }
1130 }
1131
1132 static inline void monowhite2Y(uint8_t *dst, const uint8_t *src, long width, uint32_t *unused)
1133 {
1134     int i, j;
1135     for (i=0; i<width/8; i++) {
1136         int d= ~src[i];
1137         for(j=0; j<8; j++)
1138             dst[8*i+j]= ((d>>(7-j))&1)*255;
1139     }
1140 }
1141
1142 static inline void monoblack2Y(uint8_t *dst, const uint8_t *src, long width, uint32_t *unused)
1143 {
1144     int i, j;
1145     for (i=0; i<width/8; i++) {
1146         int d= src[i];
1147         for(j=0; j<8; j++)
1148             dst[8*i+j]= ((d>>(7-j))&1)*255;
1149     }
1150 }
1151
1152 //Note: we have C, MMX, MMX2, 3DNOW versions, there is no 3DNOW+MMX2 one
1153 //Plain C versions
1154 #if (!HAVE_MMX && !HAVE_ALTIVEC) || CONFIG_RUNTIME_CPUDETECT
1155 #define COMPILE_C
1156 #endif
1157
1158 #if ARCH_PPC
1159 #if HAVE_ALTIVEC
1160 #define COMPILE_ALTIVEC
1161 #endif
1162 #endif //ARCH_PPC
1163
1164 #if ARCH_X86
1165
1166 #if (HAVE_MMX && !HAVE_AMD3DNOW && !HAVE_MMX2) || CONFIG_RUNTIME_CPUDETECT
1167 #define COMPILE_MMX
1168 #endif
1169
1170 #if HAVE_MMX2 || CONFIG_RUNTIME_CPUDETECT
1171 #define COMPILE_MMX2
1172 #endif
1173
1174 #if (HAVE_AMD3DNOW && !HAVE_MMX2) || CONFIG_RUNTIME_CPUDETECT
1175 #define COMPILE_3DNOW
1176 #endif
1177 #endif //ARCH_X86
1178
1179 #define COMPILE_TEMPLATE_MMX 0
1180 #define COMPILE_TEMPLATE_MMX2 0
1181 #define COMPILE_TEMPLATE_AMD3DNOW 0
1182 #define COMPILE_TEMPLATE_ALTIVEC 0
1183
1184 #ifdef COMPILE_C
1185 #define RENAME(a) a ## _C
1186 #include "swscale_template.c"
1187 #endif
1188
1189 #ifdef COMPILE_ALTIVEC
1190 #undef RENAME
1191 #undef COMPILE_TEMPLATE_ALTIVEC
1192 #define COMPILE_TEMPLATE_ALTIVEC 1
1193 #define RENAME(a) a ## _altivec
1194 #include "swscale_template.c"
1195 #endif
1196
1197 #if ARCH_X86
1198
1199 //MMX versions
1200 #ifdef COMPILE_MMX
1201 #undef RENAME
1202 #undef COMPILE_TEMPLATE_MMX
1203 #undef COMPILE_TEMPLATE_MMX2
1204 #undef COMPILE_TEMPLATE_AMD3DNOW
1205 #define COMPILE_TEMPLATE_MMX 1
1206 #define COMPILE_TEMPLATE_MMX2 0
1207 #define COMPILE_TEMPLATE_AMD3DNOW 0
1208 #define RENAME(a) a ## _MMX
1209 #include "swscale_template.c"
1210 #endif
1211
1212 //MMX2 versions
1213 #ifdef COMPILE_MMX2
1214 #undef RENAME
1215 #undef COMPILE_TEMPLATE_MMX
1216 #undef COMPILE_TEMPLATE_MMX2
1217 #undef COMPILE_TEMPLATE_AMD3DNOW
1218 #define COMPILE_TEMPLATE_MMX 1
1219 #define COMPILE_TEMPLATE_MMX2 1
1220 #define COMPILE_TEMPLATE_AMD3DNOW 0
1221 #define RENAME(a) a ## _MMX2
1222 #include "swscale_template.c"
1223 #endif
1224
1225 //3DNOW versions
1226 #ifdef COMPILE_3DNOW
1227 #undef RENAME
1228 #undef COMPILE_TEMPLATE_MMX
1229 #undef COMPILE_TEMPLATE_MMX2
1230 #undef COMPILE_TEMPLATE_AMD3DNOW
1231 #define COMPILE_TEMPLATE_MMX 1
1232 #define COMPILE_TEMPLATE_MMX2 0
1233 #define COMPILE_TEMPLATE_AMD3DNOW 1
1234 #define RENAME(a) a ## _3DNow
1235 #include "swscale_template.c"
1236 #endif
1237
1238 #endif //ARCH_X86
1239
1240 SwsFunc ff_getSwsFunc(SwsContext *c)
1241 {
1242 #if CONFIG_RUNTIME_CPUDETECT
1243     int flags = c->flags;
1244
1245 #if ARCH_X86
1246     // ordered per speed fastest first
1247     if (flags & SWS_CPU_CAPS_MMX2) {
1248         sws_init_swScale_MMX2(c);
1249         return swScale_MMX2;
1250     } else if (flags & SWS_CPU_CAPS_3DNOW) {
1251         sws_init_swScale_3DNow(c);
1252         return swScale_3DNow;
1253     } else if (flags & SWS_CPU_CAPS_MMX) {
1254         sws_init_swScale_MMX(c);
1255         return swScale_MMX;
1256     } else {
1257         sws_init_swScale_C(c);
1258         return swScale_C;
1259     }
1260
1261 #else
1262 #ifdef COMPILE_ALTIVEC
1263     if (flags & SWS_CPU_CAPS_ALTIVEC) {
1264         sws_init_swScale_altivec(c);
1265         return swScale_altivec;
1266     } else {
1267         sws_init_swScale_C(c);
1268         return swScale_C;
1269     }
1270 #endif
1271     sws_init_swScale_C(c);
1272     return swScale_C;
1273 #endif /* ARCH_X86 */
1274 #else //CONFIG_RUNTIME_CPUDETECT
1275 #if   COMPILE_TEMPLATE_MMX2
1276     sws_init_swScale_MMX2(c);
1277     return swScale_MMX2;
1278 #elif COMPILE_TEMPLATE_AMD3DNOW
1279     sws_init_swScale_3DNow(c);
1280     return swScale_3DNow;
1281 #elif COMPILE_TEMPLATE_MMX
1282     sws_init_swScale_MMX(c);
1283     return swScale_MMX;
1284 #elif COMPILE_TEMPLATE_ALTIVEC
1285     sws_init_swScale_altivec(c);
1286     return swScale_altivec;
1287 #else
1288     sws_init_swScale_C(c);
1289     return swScale_C;
1290 #endif
1291 #endif //!CONFIG_RUNTIME_CPUDETECT
1292 }
1293
1294 static int planarToNv12Wrapper(SwsContext *c, const uint8_t* src[], int srcStride[], int srcSliceY,
1295                                int srcSliceH, uint8_t* dstParam[], int dstStride[])
1296 {
1297     uint8_t *dst=dstParam[0] + dstStride[0]*srcSliceY;
1298     /* Copy Y plane */
1299     if (dstStride[0]==srcStride[0] && srcStride[0] > 0)
1300         memcpy(dst, src[0], srcSliceH*dstStride[0]);
1301     else {
1302         int i;
1303         const uint8_t *srcPtr= src[0];
1304         uint8_t *dstPtr= dst;
1305         for (i=0; i<srcSliceH; i++) {
1306             memcpy(dstPtr, srcPtr, c->srcW);
1307             srcPtr+= srcStride[0];
1308             dstPtr+= dstStride[0];
1309         }
1310     }
1311     dst = dstParam[1] + dstStride[1]*srcSliceY/2;
1312     if (c->dstFormat == PIX_FMT_NV12)
1313         interleaveBytes(src[1], src[2], dst, c->srcW/2, srcSliceH/2, srcStride[1], srcStride[2], dstStride[0]);
1314     else
1315         interleaveBytes(src[2], src[1], dst, c->srcW/2, srcSliceH/2, srcStride[2], srcStride[1], dstStride[0]);
1316
1317     return srcSliceH;
1318 }
1319
1320 static int planarToYuy2Wrapper(SwsContext *c, const uint8_t* src[], int srcStride[], int srcSliceY,
1321                                int srcSliceH, uint8_t* dstParam[], int dstStride[])
1322 {
1323     uint8_t *dst=dstParam[0] + dstStride[0]*srcSliceY;
1324
1325     yv12toyuy2(src[0], src[1], src[2], dst, c->srcW, srcSliceH, srcStride[0], srcStride[1], dstStride[0]);
1326
1327     return srcSliceH;
1328 }
1329
1330 static int planarToUyvyWrapper(SwsContext *c, const uint8_t* src[], int srcStride[], int srcSliceY,
1331                                int srcSliceH, uint8_t* dstParam[], int dstStride[])
1332 {
1333     uint8_t *dst=dstParam[0] + dstStride[0]*srcSliceY;
1334
1335     yv12touyvy(src[0], src[1], src[2], dst, c->srcW, srcSliceH, srcStride[0], srcStride[1], dstStride[0]);
1336
1337     return srcSliceH;
1338 }
1339
1340 static int yuv422pToYuy2Wrapper(SwsContext *c, const uint8_t* src[], int srcStride[], int srcSliceY,
1341                                 int srcSliceH, uint8_t* dstParam[], int dstStride[])
1342 {
1343     uint8_t *dst=dstParam[0] + dstStride[0]*srcSliceY;
1344
1345     yuv422ptoyuy2(src[0],src[1],src[2],dst,c->srcW,srcSliceH,srcStride[0],srcStride[1],dstStride[0]);
1346
1347     return srcSliceH;
1348 }
1349
1350 static int yuv422pToUyvyWrapper(SwsContext *c, const uint8_t* src[], int srcStride[], int srcSliceY,
1351                                 int srcSliceH, uint8_t* dstParam[], int dstStride[])
1352 {
1353     uint8_t *dst=dstParam[0] + dstStride[0]*srcSliceY;
1354
1355     yuv422ptouyvy(src[0],src[1],src[2],dst,c->srcW,srcSliceH,srcStride[0],srcStride[1],dstStride[0]);
1356
1357     return srcSliceH;
1358 }
1359
1360 static int yuyvToYuv420Wrapper(SwsContext *c, const uint8_t* src[], int srcStride[], int srcSliceY,
1361                                int srcSliceH, uint8_t* dstParam[], int dstStride[])
1362 {
1363     uint8_t *ydst=dstParam[0] + dstStride[0]*srcSliceY;
1364     uint8_t *udst=dstParam[1] + dstStride[1]*srcSliceY/2;
1365     uint8_t *vdst=dstParam[2] + dstStride[2]*srcSliceY/2;
1366
1367     yuyvtoyuv420(ydst, udst, vdst, src[0], c->srcW, srcSliceH, dstStride[0], dstStride[1], srcStride[0]);
1368
1369     if (dstParam[3])
1370         fillPlane(dstParam[3], dstStride[3], c->srcW, srcSliceH, srcSliceY, 255);
1371
1372     return srcSliceH;
1373 }
1374
1375 static int yuyvToYuv422Wrapper(SwsContext *c, const uint8_t* src[], int srcStride[], int srcSliceY,
1376                                int srcSliceH, uint8_t* dstParam[], int dstStride[])
1377 {
1378     uint8_t *ydst=dstParam[0] + dstStride[0]*srcSliceY;
1379     uint8_t *udst=dstParam[1] + dstStride[1]*srcSliceY;
1380     uint8_t *vdst=dstParam[2] + dstStride[2]*srcSliceY;
1381
1382     yuyvtoyuv422(ydst, udst, vdst, src[0], c->srcW, srcSliceH, dstStride[0], dstStride[1], srcStride[0]);
1383
1384     return srcSliceH;
1385 }
1386
1387 static int uyvyToYuv420Wrapper(SwsContext *c, const uint8_t* src[], int srcStride[], int srcSliceY,
1388                                int srcSliceH, uint8_t* dstParam[], int dstStride[])
1389 {
1390     uint8_t *ydst=dstParam[0] + dstStride[0]*srcSliceY;
1391     uint8_t *udst=dstParam[1] + dstStride[1]*srcSliceY/2;
1392     uint8_t *vdst=dstParam[2] + dstStride[2]*srcSliceY/2;
1393
1394     uyvytoyuv420(ydst, udst, vdst, src[0], c->srcW, srcSliceH, dstStride[0], dstStride[1], srcStride[0]);
1395
1396     if (dstParam[3])
1397         fillPlane(dstParam[3], dstStride[3], c->srcW, srcSliceH, srcSliceY, 255);
1398
1399     return srcSliceH;
1400 }
1401
1402 static int uyvyToYuv422Wrapper(SwsContext *c, const uint8_t* src[], int srcStride[], int srcSliceY,
1403                                int srcSliceH, uint8_t* dstParam[], int dstStride[])
1404 {
1405     uint8_t *ydst=dstParam[0] + dstStride[0]*srcSliceY;
1406     uint8_t *udst=dstParam[1] + dstStride[1]*srcSliceY;
1407     uint8_t *vdst=dstParam[2] + dstStride[2]*srcSliceY;
1408
1409     uyvytoyuv422(ydst, udst, vdst, src[0], c->srcW, srcSliceH, dstStride[0], dstStride[1], srcStride[0]);
1410
1411     return srcSliceH;
1412 }
1413
1414 static int palToRgbWrapper(SwsContext *c, const uint8_t* src[], int srcStride[], int srcSliceY,
1415                            int srcSliceH, uint8_t* dst[], int dstStride[])
1416 {
1417     const enum PixelFormat srcFormat= c->srcFormat;
1418     const enum PixelFormat dstFormat= c->dstFormat;
1419     void (*conv)(const uint8_t *src, uint8_t *dst, long num_pixels,
1420                  const uint8_t *palette)=NULL;
1421     int i;
1422     uint8_t *dstPtr= dst[0] + dstStride[0]*srcSliceY;
1423     const uint8_t *srcPtr= src[0];
1424
1425     if (usePal(srcFormat)) {
1426         switch (dstFormat) {
1427         case PIX_FMT_RGB32  : conv = palette8topacked32; break;
1428         case PIX_FMT_BGR32  : conv = palette8topacked32; break;
1429         case PIX_FMT_BGR32_1: conv = palette8topacked32; break;
1430         case PIX_FMT_RGB32_1: conv = palette8topacked32; break;
1431         case PIX_FMT_RGB24  : conv = palette8topacked24; break;
1432         case PIX_FMT_BGR24  : conv = palette8topacked24; break;
1433         }
1434     }
1435
1436     if (!conv)
1437         av_log(c, AV_LOG_ERROR, "internal error %s -> %s converter\n",
1438                sws_format_name(srcFormat), sws_format_name(dstFormat));
1439     else {
1440         for (i=0; i<srcSliceH; i++) {
1441             conv(srcPtr, dstPtr, c->srcW, (uint8_t *) c->pal_rgb);
1442             srcPtr+= srcStride[0];
1443             dstPtr+= dstStride[0];
1444         }
1445     }
1446
1447     return srcSliceH;
1448 }
1449
1450 #define isRGBA32(x) (            \
1451            (x) == PIX_FMT_ARGB   \
1452         || (x) == PIX_FMT_RGBA   \
1453         || (x) == PIX_FMT_BGRA   \
1454         || (x) == PIX_FMT_ABGR   \
1455         )
1456
1457 /* {RGB,BGR}{15,16,24,32,32_1} -> {RGB,BGR}{15,16,24,32} */
1458 static int rgbToRgbWrapper(SwsContext *c, const uint8_t* src[], int srcStride[], int srcSliceY,
1459                            int srcSliceH, uint8_t* dst[], int dstStride[])
1460 {
1461     const enum PixelFormat srcFormat= c->srcFormat;
1462     const enum PixelFormat dstFormat= c->dstFormat;
1463     const int srcBpp= (c->srcFormatBpp + 7) >> 3;
1464     const int dstBpp= (c->dstFormatBpp + 7) >> 3;
1465     const int srcId= c->srcFormatBpp >> 2; /* 1:0, 4:1, 8:2, 15:3, 16:4, 24:6, 32:8 */
1466     const int dstId= c->dstFormatBpp >> 2;
1467     void (*conv)(const uint8_t *src, uint8_t *dst, long src_size)=NULL;
1468
1469 #define CONV_IS(src, dst) (srcFormat == PIX_FMT_##src && dstFormat == PIX_FMT_##dst)
1470
1471     if (isRGBA32(srcFormat) && isRGBA32(dstFormat)) {
1472         if (     CONV_IS(ABGR, RGBA)
1473               || CONV_IS(ARGB, BGRA)
1474               || CONV_IS(BGRA, ARGB)
1475               || CONV_IS(RGBA, ABGR)) conv = shuffle_bytes_3210;
1476         else if (CONV_IS(ABGR, ARGB)
1477               || CONV_IS(ARGB, ABGR)) conv = shuffle_bytes_0321;
1478         else if (CONV_IS(ABGR, BGRA)
1479               || CONV_IS(ARGB, RGBA)) conv = shuffle_bytes_1230;
1480         else if (CONV_IS(BGRA, RGBA)
1481               || CONV_IS(RGBA, BGRA)) conv = shuffle_bytes_2103;
1482         else if (CONV_IS(BGRA, ABGR)
1483               || CONV_IS(RGBA, ARGB)) conv = shuffle_bytes_3012;
1484     } else
1485     /* BGR -> BGR */
1486     if (  (isBGRinInt(srcFormat) && isBGRinInt(dstFormat))
1487        || (isRGBinInt(srcFormat) && isRGBinInt(dstFormat))) {
1488         switch(srcId | (dstId<<4)) {
1489         case 0x34: conv= rgb16to15; break;
1490         case 0x36: conv= rgb24to15; break;
1491         case 0x38: conv= rgb32to15; break;
1492         case 0x43: conv= rgb15to16; break;
1493         case 0x46: conv= rgb24to16; break;
1494         case 0x48: conv= rgb32to16; break;
1495         case 0x63: conv= rgb15to24; break;
1496         case 0x64: conv= rgb16to24; break;
1497         case 0x68: conv= rgb32to24; break;
1498         case 0x83: conv= rgb15to32; break;
1499         case 0x84: conv= rgb16to32; break;
1500         case 0x86: conv= rgb24to32; break;
1501         }
1502     } else if (  (isBGRinInt(srcFormat) && isRGBinInt(dstFormat))
1503              || (isRGBinInt(srcFormat) && isBGRinInt(dstFormat))) {
1504         switch(srcId | (dstId<<4)) {
1505         case 0x33: conv= rgb15tobgr15; break;
1506         case 0x34: conv= rgb16tobgr15; break;
1507         case 0x36: conv= rgb24tobgr15; break;
1508         case 0x38: conv= rgb32tobgr15; break;
1509         case 0x43: conv= rgb15tobgr16; break;
1510         case 0x44: conv= rgb16tobgr16; break;
1511         case 0x46: conv= rgb24tobgr16; break;
1512         case 0x48: conv= rgb32tobgr16; break;
1513         case 0x63: conv= rgb15tobgr24; break;
1514         case 0x64: conv= rgb16tobgr24; break;
1515         case 0x66: conv= rgb24tobgr24; break;
1516         case 0x68: conv= rgb32tobgr24; break;
1517         case 0x83: conv= rgb15tobgr32; break;
1518         case 0x84: conv= rgb16tobgr32; break;
1519         case 0x86: conv= rgb24tobgr32; break;
1520         }
1521     }
1522
1523     if (!conv) {
1524         av_log(c, AV_LOG_ERROR, "internal error %s -> %s converter\n",
1525                sws_format_name(srcFormat), sws_format_name(dstFormat));
1526     } else {
1527         const uint8_t *srcPtr= src[0];
1528               uint8_t *dstPtr= dst[0];
1529         if ((srcFormat == PIX_FMT_RGB32_1 || srcFormat == PIX_FMT_BGR32_1) && !isRGBA32(dstFormat))
1530             srcPtr += ALT32_CORR;
1531
1532         if ((dstFormat == PIX_FMT_RGB32_1 || dstFormat == PIX_FMT_BGR32_1) && !isRGBA32(srcFormat))
1533             dstPtr += ALT32_CORR;
1534
1535         if (dstStride[0]*srcBpp == srcStride[0]*dstBpp && srcStride[0] > 0)
1536             conv(srcPtr, dstPtr + dstStride[0]*srcSliceY, srcSliceH*srcStride[0]);
1537         else {
1538             int i;
1539             dstPtr += dstStride[0]*srcSliceY;
1540
1541             for (i=0; i<srcSliceH; i++) {
1542                 conv(srcPtr, dstPtr, c->srcW*srcBpp);
1543                 srcPtr+= srcStride[0];
1544                 dstPtr+= dstStride[0];
1545             }
1546         }
1547     }
1548     return srcSliceH;
1549 }
1550
1551 static int bgr24ToYv12Wrapper(SwsContext *c, const uint8_t* src[], int srcStride[], int srcSliceY,
1552                               int srcSliceH, uint8_t* dst[], int dstStride[])
1553 {
1554     rgb24toyv12(
1555         src[0],
1556         dst[0]+ srcSliceY    *dstStride[0],
1557         dst[1]+(srcSliceY>>1)*dstStride[1],
1558         dst[2]+(srcSliceY>>1)*dstStride[2],
1559         c->srcW, srcSliceH,
1560         dstStride[0], dstStride[1], srcStride[0]);
1561     if (dst[3])
1562         fillPlane(dst[3], dstStride[3], c->srcW, srcSliceH, srcSliceY, 255);
1563     return srcSliceH;
1564 }
1565
1566 static int yvu9ToYv12Wrapper(SwsContext *c, const uint8_t* src[], int srcStride[], int srcSliceY,
1567                              int srcSliceH, uint8_t* dst[], int dstStride[])
1568 {
1569     int i;
1570
1571     /* copy Y */
1572     if (srcStride[0]==dstStride[0] && srcStride[0] > 0)
1573         memcpy(dst[0]+ srcSliceY*dstStride[0], src[0], srcStride[0]*srcSliceH);
1574     else {
1575         const uint8_t *srcPtr= src[0];
1576         uint8_t *dstPtr= dst[0] + dstStride[0]*srcSliceY;
1577
1578         for (i=0; i<srcSliceH; i++) {
1579             memcpy(dstPtr, srcPtr, c->srcW);
1580             srcPtr+= srcStride[0];
1581             dstPtr+= dstStride[0];
1582         }
1583     }
1584
1585     if (c->dstFormat==PIX_FMT_YUV420P || c->dstFormat==PIX_FMT_YUVA420P) {
1586         planar2x(src[1], dst[1] + dstStride[1]*(srcSliceY >> 1), c->chrSrcW,
1587                  srcSliceH >> 2, srcStride[1], dstStride[1]);
1588         planar2x(src[2], dst[2] + dstStride[2]*(srcSliceY >> 1), c->chrSrcW,
1589                  srcSliceH >> 2, srcStride[2], dstStride[2]);
1590     } else {
1591         planar2x(src[1], dst[2] + dstStride[2]*(srcSliceY >> 1), c->chrSrcW,
1592                  srcSliceH >> 2, srcStride[1], dstStride[2]);
1593         planar2x(src[2], dst[1] + dstStride[1]*(srcSliceY >> 1), c->chrSrcW,
1594                  srcSliceH >> 2, srcStride[2], dstStride[1]);
1595     }
1596     if (dst[3])
1597         fillPlane(dst[3], dstStride[3], c->srcW, srcSliceH, srcSliceY, 255);
1598     return srcSliceH;
1599 }
1600
1601 /* unscaled copy like stuff (assumes nearly identical formats) */
1602 static int packedCopyWrapper(SwsContext *c, const uint8_t* src[], int srcStride[], int srcSliceY,
1603                              int srcSliceH, uint8_t* dst[], int dstStride[])
1604 {
1605     if (dstStride[0]==srcStride[0] && srcStride[0] > 0)
1606         memcpy(dst[0] + dstStride[0]*srcSliceY, src[0], srcSliceH*dstStride[0]);
1607     else {
1608         int i;
1609         const uint8_t *srcPtr= src[0];
1610         uint8_t *dstPtr= dst[0] + dstStride[0]*srcSliceY;
1611         int length=0;
1612
1613         /* universal length finder */
1614         while(length+c->srcW <= FFABS(dstStride[0])
1615            && length+c->srcW <= FFABS(srcStride[0])) length+= c->srcW;
1616         assert(length!=0);
1617
1618         for (i=0; i<srcSliceH; i++) {
1619             memcpy(dstPtr, srcPtr, length);
1620             srcPtr+= srcStride[0];
1621             dstPtr+= dstStride[0];
1622         }
1623     }
1624     return srcSliceH;
1625 }
1626
1627 static int planarCopyWrapper(SwsContext *c, const uint8_t* src[], int srcStride[], int srcSliceY,
1628                              int srcSliceH, uint8_t* dst[], int dstStride[])
1629 {
1630     int plane, i, j;
1631     for (plane=0; plane<4; plane++) {
1632         int length= (plane==0 || plane==3) ? c->srcW  : -((-c->srcW  )>>c->chrDstHSubSample);
1633         int y=      (plane==0 || plane==3) ? srcSliceY: -((-srcSliceY)>>c->chrDstVSubSample);
1634         int height= (plane==0 || plane==3) ? srcSliceH: -((-srcSliceH)>>c->chrDstVSubSample);
1635         const uint8_t *srcPtr= src[plane];
1636         uint8_t *dstPtr= dst[plane] + dstStride[plane]*y;
1637
1638         if (!dst[plane]) continue;
1639         // ignore palette for GRAY8
1640         if (plane == 1 && !dst[2]) continue;
1641         if (!src[plane] || (plane == 1 && !src[2])) {
1642             if(is16BPS(c->dstFormat))
1643                 length*=2;
1644             fillPlane(dst[plane], dstStride[plane], length, height, y, (plane==3) ? 255 : 128);
1645         } else {
1646             if(is16BPS(c->srcFormat) && !is16BPS(c->dstFormat)) {
1647                 if (!isBE(c->srcFormat)) srcPtr++;
1648                 for (i=0; i<height; i++) {
1649                     for (j=0; j<length; j++) dstPtr[j] = srcPtr[j<<1];
1650                     srcPtr+= srcStride[plane];
1651                     dstPtr+= dstStride[plane];
1652                 }
1653             } else if(!is16BPS(c->srcFormat) && is16BPS(c->dstFormat)) {
1654                 for (i=0; i<height; i++) {
1655                     for (j=0; j<length; j++) {
1656                         dstPtr[ j<<1   ] = srcPtr[j];
1657                         dstPtr[(j<<1)+1] = srcPtr[j];
1658                     }
1659                     srcPtr+= srcStride[plane];
1660                     dstPtr+= dstStride[plane];
1661                 }
1662             } else if(is16BPS(c->srcFormat) && is16BPS(c->dstFormat)
1663                   && isBE(c->srcFormat) != isBE(c->dstFormat)) {
1664
1665                 for (i=0; i<height; i++) {
1666                     for (j=0; j<length; j++)
1667                         ((uint16_t*)dstPtr)[j] = bswap_16(((const uint16_t*)srcPtr)[j]);
1668                     srcPtr+= srcStride[plane];
1669                     dstPtr+= dstStride[plane];
1670                 }
1671             } else if (dstStride[plane]==srcStride[plane] && srcStride[plane] > 0)
1672                 memcpy(dst[plane] + dstStride[plane]*y, src[plane], height*dstStride[plane]);
1673             else {
1674                 if(is16BPS(c->srcFormat) && is16BPS(c->dstFormat))
1675                     length*=2;
1676                 for (i=0; i<height; i++) {
1677                     memcpy(dstPtr, srcPtr, length);
1678                     srcPtr+= srcStride[plane];
1679                     dstPtr+= dstStride[plane];
1680                 }
1681             }
1682         }
1683     }
1684     return srcSliceH;
1685 }
1686
1687 int ff_hardcodedcpuflags(void)
1688 {
1689     int flags = 0;
1690 #if   COMPILE_TEMPLATE_MMX2
1691     flags |= SWS_CPU_CAPS_MMX|SWS_CPU_CAPS_MMX2;
1692 #elif COMPILE_TEMPLATE_AMD3DNOW
1693     flags |= SWS_CPU_CAPS_MMX|SWS_CPU_CAPS_3DNOW;
1694 #elif COMPILE_TEMPLATE_MMX
1695     flags |= SWS_CPU_CAPS_MMX;
1696 #elif COMPILE_TEMPLATE_ALTIVEC
1697     flags |= SWS_CPU_CAPS_ALTIVEC;
1698 #elif ARCH_BFIN
1699     flags |= SWS_CPU_CAPS_BFIN;
1700 #endif
1701     return flags;
1702 }
1703
1704 void ff_get_unscaled_swscale(SwsContext *c)
1705 {
1706     const enum PixelFormat srcFormat = c->srcFormat;
1707     const enum PixelFormat dstFormat = c->dstFormat;
1708     const int flags = c->flags;
1709     const int dstH = c->dstH;
1710     int needsDither;
1711
1712     needsDither= isAnyRGB(dstFormat)
1713         &&  c->dstFormatBpp < 24
1714         && (c->dstFormatBpp < c->srcFormatBpp || (!isAnyRGB(srcFormat)));
1715
1716     /* yv12_to_nv12 */
1717     if ((srcFormat == PIX_FMT_YUV420P || srcFormat == PIX_FMT_YUVA420P) && (dstFormat == PIX_FMT_NV12 || dstFormat == PIX_FMT_NV21)) {
1718         c->swScale= planarToNv12Wrapper;
1719     }
1720     /* yuv2bgr */
1721     if ((srcFormat==PIX_FMT_YUV420P || srcFormat==PIX_FMT_YUV422P || srcFormat==PIX_FMT_YUVA420P) && isAnyRGB(dstFormat)
1722         && !(flags & SWS_ACCURATE_RND) && !(dstH&1)) {
1723         c->swScale= ff_yuv2rgb_get_func_ptr(c);
1724     }
1725
1726     if (srcFormat==PIX_FMT_YUV410P && (dstFormat==PIX_FMT_YUV420P || dstFormat==PIX_FMT_YUVA420P) && !(flags & SWS_BITEXACT)) {
1727         c->swScale= yvu9ToYv12Wrapper;
1728     }
1729
1730     /* bgr24toYV12 */
1731     if (srcFormat==PIX_FMT_BGR24 && (dstFormat==PIX_FMT_YUV420P || dstFormat==PIX_FMT_YUVA420P) && !(flags & SWS_ACCURATE_RND))
1732         c->swScale= bgr24ToYv12Wrapper;
1733
1734     /* RGB/BGR -> RGB/BGR (no dither needed forms) */
1735     if (   isAnyRGB(srcFormat)
1736         && isAnyRGB(dstFormat)
1737         && srcFormat != PIX_FMT_BGR8      && dstFormat != PIX_FMT_BGR8
1738         && srcFormat != PIX_FMT_RGB8      && dstFormat != PIX_FMT_RGB8
1739         && srcFormat != PIX_FMT_BGR4      && dstFormat != PIX_FMT_BGR4
1740         && srcFormat != PIX_FMT_RGB4      && dstFormat != PIX_FMT_RGB4
1741         && srcFormat != PIX_FMT_BGR4_BYTE && dstFormat != PIX_FMT_BGR4_BYTE
1742         && srcFormat != PIX_FMT_RGB4_BYTE && dstFormat != PIX_FMT_RGB4_BYTE
1743         && srcFormat != PIX_FMT_MONOBLACK && dstFormat != PIX_FMT_MONOBLACK
1744         && srcFormat != PIX_FMT_MONOWHITE && dstFormat != PIX_FMT_MONOWHITE
1745         && srcFormat != PIX_FMT_RGB48LE   && dstFormat != PIX_FMT_RGB48LE
1746         && srcFormat != PIX_FMT_RGB48BE   && dstFormat != PIX_FMT_RGB48BE
1747         && (!needsDither || (c->flags&(SWS_FAST_BILINEAR|SWS_POINT))))
1748         c->swScale= rgbToRgbWrapper;
1749
1750     if ((usePal(srcFormat) && (
1751         dstFormat == PIX_FMT_RGB32   ||
1752         dstFormat == PIX_FMT_RGB32_1 ||
1753         dstFormat == PIX_FMT_RGB24   ||
1754         dstFormat == PIX_FMT_BGR32   ||
1755         dstFormat == PIX_FMT_BGR32_1 ||
1756         dstFormat == PIX_FMT_BGR24)))
1757         c->swScale= palToRgbWrapper;
1758
1759     if (srcFormat == PIX_FMT_YUV422P) {
1760         if (dstFormat == PIX_FMT_YUYV422)
1761             c->swScale= yuv422pToYuy2Wrapper;
1762         else if (dstFormat == PIX_FMT_UYVY422)
1763             c->swScale= yuv422pToUyvyWrapper;
1764     }
1765
1766     /* LQ converters if -sws 0 or -sws 4*/
1767     if (c->flags&(SWS_FAST_BILINEAR|SWS_POINT)) {
1768         /* yv12_to_yuy2 */
1769         if (srcFormat == PIX_FMT_YUV420P || srcFormat == PIX_FMT_YUVA420P) {
1770             if (dstFormat == PIX_FMT_YUYV422)
1771                 c->swScale= planarToYuy2Wrapper;
1772             else if (dstFormat == PIX_FMT_UYVY422)
1773                 c->swScale= planarToUyvyWrapper;
1774         }
1775     }
1776     if(srcFormat == PIX_FMT_YUYV422 && (dstFormat == PIX_FMT_YUV420P || dstFormat == PIX_FMT_YUVA420P))
1777         c->swScale= yuyvToYuv420Wrapper;
1778     if(srcFormat == PIX_FMT_UYVY422 && (dstFormat == PIX_FMT_YUV420P || dstFormat == PIX_FMT_YUVA420P))
1779         c->swScale= uyvyToYuv420Wrapper;
1780     if(srcFormat == PIX_FMT_YUYV422 && dstFormat == PIX_FMT_YUV422P)
1781         c->swScale= yuyvToYuv422Wrapper;
1782     if(srcFormat == PIX_FMT_UYVY422 && dstFormat == PIX_FMT_YUV422P)
1783         c->swScale= uyvyToYuv422Wrapper;
1784
1785 #ifdef COMPILE_ALTIVEC
1786     if ((c->flags & SWS_CPU_CAPS_ALTIVEC) &&
1787         !(c->flags & SWS_BITEXACT) &&
1788         srcFormat == PIX_FMT_YUV420P) {
1789         // unscaled YV12 -> packed YUV, we want speed
1790         if (dstFormat == PIX_FMT_YUYV422)
1791             c->swScale= yv12toyuy2_unscaled_altivec;
1792         else if (dstFormat == PIX_FMT_UYVY422)
1793             c->swScale= yv12touyvy_unscaled_altivec;
1794     }
1795 #endif
1796
1797     /* simple copy */
1798     if (  srcFormat == dstFormat
1799         || (srcFormat == PIX_FMT_YUVA420P && dstFormat == PIX_FMT_YUV420P)
1800         || (srcFormat == PIX_FMT_YUV420P && dstFormat == PIX_FMT_YUVA420P)
1801         || (isPlanarYUV(srcFormat) && isGray(dstFormat))
1802         || (isPlanarYUV(dstFormat) && isGray(srcFormat))
1803         || (isGray(dstFormat) && isGray(srcFormat))
1804         || (isPlanarYUV(srcFormat) && isPlanarYUV(dstFormat)
1805             && c->chrDstHSubSample == c->chrSrcHSubSample
1806             && c->chrDstVSubSample == c->chrSrcVSubSample
1807             && dstFormat != PIX_FMT_NV12 && dstFormat != PIX_FMT_NV21
1808             && srcFormat != PIX_FMT_NV12 && srcFormat != PIX_FMT_NV21))
1809     {
1810         if (isPacked(c->srcFormat))
1811             c->swScale= packedCopyWrapper;
1812         else /* Planar YUV or gray */
1813             c->swScale= planarCopyWrapper;
1814     }
1815 #if ARCH_BFIN
1816     if (flags & SWS_CPU_CAPS_BFIN)
1817         ff_bfin_get_unscaled_swscale (c);
1818 #endif
1819 }
1820
1821 static void reset_ptr(const uint8_t* src[], int format)
1822 {
1823     if(!isALPHA(format))
1824         src[3]=NULL;
1825     if(!isPlanarYUV(format)) {
1826         src[3]=src[2]=NULL;
1827
1828         if (!usePal(format))
1829             src[1]= NULL;
1830     }
1831 }
1832
1833 /**
1834  * swscale wrapper, so we don't need to export the SwsContext.
1835  * Assumes planar YUV to be in YUV order instead of YVU.
1836  */
1837 int sws_scale(SwsContext *c, const uint8_t* const src[], const int srcStride[], int srcSliceY,
1838               int srcSliceH, uint8_t* const dst[], const int dstStride[])
1839 {
1840     int i;
1841     const uint8_t* src2[4]= {src[0], src[1], src[2], src[3]};
1842     uint8_t* dst2[4]= {dst[0], dst[1], dst[2], dst[3]};
1843
1844     // do not mess up sliceDir if we have a "trailing" 0-size slice
1845     if (srcSliceH == 0)
1846         return 0;
1847
1848     if (c->sliceDir == 0 && srcSliceY != 0 && srcSliceY + srcSliceH != c->srcH) {
1849         av_log(c, AV_LOG_ERROR, "Slices start in the middle!\n");
1850         return 0;
1851     }
1852     if (c->sliceDir == 0) {
1853         if (srcSliceY == 0) c->sliceDir = 1; else c->sliceDir = -1;
1854     }
1855
1856     if (usePal(c->srcFormat)) {
1857         for (i=0; i<256; i++) {
1858             int p, r, g, b,y,u,v;
1859             if(c->srcFormat == PIX_FMT_PAL8) {
1860                 p=((const uint32_t*)(src[1]))[i];
1861                 r= (p>>16)&0xFF;
1862                 g= (p>> 8)&0xFF;
1863                 b=  p     &0xFF;
1864             } else if(c->srcFormat == PIX_FMT_RGB8) {
1865                 r= (i>>5    )*36;
1866                 g= ((i>>2)&7)*36;
1867                 b= (i&3     )*85;
1868             } else if(c->srcFormat == PIX_FMT_BGR8) {
1869                 b= (i>>6    )*85;
1870                 g= ((i>>3)&7)*36;
1871                 r= (i&7     )*36;
1872             } else if(c->srcFormat == PIX_FMT_RGB4_BYTE) {
1873                 r= (i>>3    )*255;
1874                 g= ((i>>1)&3)*85;
1875                 b= (i&1     )*255;
1876             } else if(c->srcFormat == PIX_FMT_GRAY8) {
1877                 r = g = b = i;
1878             } else {
1879                 assert(c->srcFormat == PIX_FMT_BGR4_BYTE);
1880                 b= (i>>3    )*255;
1881                 g= ((i>>1)&3)*85;
1882                 r= (i&1     )*255;
1883             }
1884             y= av_clip_uint8((RY*r + GY*g + BY*b + ( 33<<(RGB2YUV_SHIFT-1)))>>RGB2YUV_SHIFT);
1885             u= av_clip_uint8((RU*r + GU*g + BU*b + (257<<(RGB2YUV_SHIFT-1)))>>RGB2YUV_SHIFT);
1886             v= av_clip_uint8((RV*r + GV*g + BV*b + (257<<(RGB2YUV_SHIFT-1)))>>RGB2YUV_SHIFT);
1887             c->pal_yuv[i]= y + (u<<8) + (v<<16);
1888
1889             switch(c->dstFormat) {
1890             case PIX_FMT_BGR32:
1891 #if !HAVE_BIGENDIAN
1892             case PIX_FMT_RGB24:
1893 #endif
1894                 c->pal_rgb[i]=  r + (g<<8) + (b<<16);
1895                 break;
1896             case PIX_FMT_BGR32_1:
1897 #if HAVE_BIGENDIAN
1898             case PIX_FMT_BGR24:
1899 #endif
1900                 c->pal_rgb[i]= (r + (g<<8) + (b<<16)) << 8;
1901                 break;
1902             case PIX_FMT_RGB32_1:
1903 #if HAVE_BIGENDIAN
1904             case PIX_FMT_RGB24:
1905 #endif
1906                 c->pal_rgb[i]= (b + (g<<8) + (r<<16)) << 8;
1907                 break;
1908             case PIX_FMT_RGB32:
1909 #if !HAVE_BIGENDIAN
1910             case PIX_FMT_BGR24:
1911 #endif
1912             default:
1913                 c->pal_rgb[i]=  b + (g<<8) + (r<<16);
1914             }
1915         }
1916     }
1917
1918     // copy strides, so they can safely be modified
1919     if (c->sliceDir == 1) {
1920         // slices go from top to bottom
1921         int srcStride2[4]= {srcStride[0], srcStride[1], srcStride[2], srcStride[3]};
1922         int dstStride2[4]= {dstStride[0], dstStride[1], dstStride[2], dstStride[3]};
1923
1924         reset_ptr(src2, c->srcFormat);
1925         reset_ptr((const uint8_t**)dst2, c->dstFormat);
1926
1927         /* reset slice direction at end of frame */
1928         if (srcSliceY + srcSliceH == c->srcH)
1929             c->sliceDir = 0;
1930
1931         return c->swScale(c, src2, srcStride2, srcSliceY, srcSliceH, dst2, dstStride2);
1932     } else {
1933         // slices go from bottom to top => we flip the image internally
1934         int srcStride2[4]= {-srcStride[0], -srcStride[1], -srcStride[2], -srcStride[3]};
1935         int dstStride2[4]= {-dstStride[0], -dstStride[1], -dstStride[2], -dstStride[3]};
1936
1937         src2[0] += (srcSliceH-1)*srcStride[0];
1938         if (!usePal(c->srcFormat))
1939             src2[1] += ((srcSliceH>>c->chrSrcVSubSample)-1)*srcStride[1];
1940         src2[2] += ((srcSliceH>>c->chrSrcVSubSample)-1)*srcStride[2];
1941         src2[3] += (srcSliceH-1)*srcStride[3];
1942         dst2[0] += ( c->dstH                      -1)*dstStride[0];
1943         dst2[1] += ((c->dstH>>c->chrDstVSubSample)-1)*dstStride[1];
1944         dst2[2] += ((c->dstH>>c->chrDstVSubSample)-1)*dstStride[2];
1945         dst2[3] += ( c->dstH                      -1)*dstStride[3];
1946
1947         reset_ptr(src2, c->srcFormat);
1948         reset_ptr((const uint8_t**)dst2, c->dstFormat);
1949
1950         /* reset slice direction at end of frame */
1951         if (!srcSliceY)
1952             c->sliceDir = 0;
1953
1954         return c->swScale(c, src2, srcStride2, c->srcH-srcSliceY-srcSliceH, srcSliceH, dst2, dstStride2);
1955     }
1956 }
1957
1958 #if LIBSWSCALE_VERSION_MAJOR < 1
1959 int sws_scale_ordered(SwsContext *c, const uint8_t* const src[], int srcStride[], int srcSliceY,
1960                       int srcSliceH, uint8_t* dst[], int dstStride[])
1961 {
1962     return sws_scale(c, src, srcStride, srcSliceY, srcSliceH, dst, dstStride);
1963 }
1964 #endif