utils.c 86.1 KB
Newer Older
1 2 3 4 5
/*
 * Copyright (C) 2001-2003 Michael Niedermayer <michaelni@gmx.at>
 *
 * This file is part of FFmpeg.
 *
6 7 8 9
 * FFmpeg is free software; you can redistribute it and/or
 * modify it under the terms of the GNU Lesser General Public
 * License as published by the Free Software Foundation; either
 * version 2.1 of the License, or (at your option) any later version.
10 11 12
 *
 * FFmpeg is distributed in the hope that it will be useful,
 * but WITHOUT ANY WARRANTY; without even the implied warranty of
13 14
 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
 * Lesser General Public License for more details.
15
 *
16 17
 * You should have received a copy of the GNU Lesser General Public
 * License along with FFmpeg; if not, write to the Free Software
18 19 20
 * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
 */

21 22
#include "config.h"

23
#define _DEFAULT_SOURCE
24
#define _SVID_SOURCE // needed for MAP_ANONYMOUS
25
#define _DARWIN_C_SOURCE // needed for MAP_ANON
26 27 28
#include <inttypes.h>
#include <math.h>
#include <stdio.h>
29
#include <string.h>
30 31 32 33 34 35 36 37 38 39
#if HAVE_SYS_MMAN_H
#include <sys/mman.h>
#if defined(MAP_ANON) && !defined(MAP_ANONYMOUS)
#define MAP_ANONYMOUS MAP_ANON
#endif
#endif
#if HAVE_VIRTUALALLOC
#define WIN32_LEAN_AND_MEAN
#include <windows.h>
#endif
40

41
#include "libavutil/attributes.h"
42
#include "libavutil/avassert.h"
43 44
#include "libavutil/avutil.h"
#include "libavutil/bswap.h"
45
#include "libavutil/cpu.h"
46
#include "libavutil/imgutils.h"
47
#include "libavutil/intreadwrite.h"
48
#include "libavutil/libm.h"
49
#include "libavutil/mathematics.h"
50
#include "libavutil/opt.h"
51
#include "libavutil/pixdesc.h"
52
#include "libavutil/ppc/cpu.h"
53
#include "libavutil/x86/asm.h"
54
#include "libavutil/x86/cpu.h"
55 56 57
#include "rgb2rgb.h"
#include "swscale.h"
#include "swscale_internal.h"
58

59 60
static void handle_formats(SwsContext *c);

61 62
unsigned swscale_version(void)
{
63
    av_assert0(LIBSWSCALE_VERSION_MICRO >= 100);
64 65 66 67 68
    return LIBSWSCALE_VERSION_INT;
}

const char *swscale_configuration(void)
{
69
    return FFMPEG_CONFIGURATION;
70 71 72 73 74
}

const char *swscale_license(void)
{
#define LICENSE_PREFIX "libswscale license: "
75
    return LICENSE_PREFIX FFMPEG_LICENSE + sizeof(LICENSE_PREFIX) - 1;
76 77
}

78
typedef struct FormatEntry {
79 80 81
    uint8_t is_supported_in         :1;
    uint8_t is_supported_out        :1;
    uint8_t is_supported_endianness :1;
82 83
} FormatEntry;

84 85 86 87 88 89 90 91 92 93 94 95 96 97
static const FormatEntry format_entries[AV_PIX_FMT_NB] = {
    [AV_PIX_FMT_YUV420P]     = { 1, 1 },
    [AV_PIX_FMT_YUYV422]     = { 1, 1 },
    [AV_PIX_FMT_RGB24]       = { 1, 1 },
    [AV_PIX_FMT_BGR24]       = { 1, 1 },
    [AV_PIX_FMT_YUV422P]     = { 1, 1 },
    [AV_PIX_FMT_YUV444P]     = { 1, 1 },
    [AV_PIX_FMT_YUV410P]     = { 1, 1 },
    [AV_PIX_FMT_YUV411P]     = { 1, 1 },
    [AV_PIX_FMT_GRAY8]       = { 1, 1 },
    [AV_PIX_FMT_MONOWHITE]   = { 1, 1 },
    [AV_PIX_FMT_MONOBLACK]   = { 1, 1 },
    [AV_PIX_FMT_PAL8]        = { 1, 0 },
    [AV_PIX_FMT_YUVJ420P]    = { 1, 1 },
Michael Niedermayer's avatar
Michael Niedermayer committed
98
    [AV_PIX_FMT_YUVJ411P]    = { 1, 1 },
99 100
    [AV_PIX_FMT_YUVJ422P]    = { 1, 1 },
    [AV_PIX_FMT_YUVJ444P]    = { 1, 1 },
101
    [AV_PIX_FMT_YVYU422]     = { 1, 1 },
102 103 104 105 106 107 108 109 110 111 112 113 114 115
    [AV_PIX_FMT_UYVY422]     = { 1, 1 },
    [AV_PIX_FMT_UYYVYY411]   = { 0, 0 },
    [AV_PIX_FMT_BGR8]        = { 1, 1 },
    [AV_PIX_FMT_BGR4]        = { 0, 1 },
    [AV_PIX_FMT_BGR4_BYTE]   = { 1, 1 },
    [AV_PIX_FMT_RGB8]        = { 1, 1 },
    [AV_PIX_FMT_RGB4]        = { 0, 1 },
    [AV_PIX_FMT_RGB4_BYTE]   = { 1, 1 },
    [AV_PIX_FMT_NV12]        = { 1, 1 },
    [AV_PIX_FMT_NV21]        = { 1, 1 },
    [AV_PIX_FMT_ARGB]        = { 1, 1 },
    [AV_PIX_FMT_RGBA]        = { 1, 1 },
    [AV_PIX_FMT_ABGR]        = { 1, 1 },
    [AV_PIX_FMT_BGRA]        = { 1, 1 },
116 117 118 119
    [AV_PIX_FMT_0RGB]        = { 1, 1 },
    [AV_PIX_FMT_RGB0]        = { 1, 1 },
    [AV_PIX_FMT_0BGR]        = { 1, 1 },
    [AV_PIX_FMT_BGR0]        = { 1, 1 },
120 121 122 123
    [AV_PIX_FMT_GRAY16BE]    = { 1, 1 },
    [AV_PIX_FMT_GRAY16LE]    = { 1, 1 },
    [AV_PIX_FMT_YUV440P]     = { 1, 1 },
    [AV_PIX_FMT_YUVJ440P]    = { 1, 1 },
124 125 126 127
    [AV_PIX_FMT_YUV440P10LE] = { 1, 1 },
    [AV_PIX_FMT_YUV440P10BE] = { 1, 1 },
    [AV_PIX_FMT_YUV440P12LE] = { 1, 1 },
    [AV_PIX_FMT_YUV440P12BE] = { 1, 1 },
128
    [AV_PIX_FMT_YUVA420P]    = { 1, 1 },
129 130
    [AV_PIX_FMT_YUVA422P]    = { 1, 1 },
    [AV_PIX_FMT_YUVA444P]    = { 1, 1 },
131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148
    [AV_PIX_FMT_YUVA420P9BE] = { 1, 1 },
    [AV_PIX_FMT_YUVA420P9LE] = { 1, 1 },
    [AV_PIX_FMT_YUVA422P9BE] = { 1, 1 },
    [AV_PIX_FMT_YUVA422P9LE] = { 1, 1 },
    [AV_PIX_FMT_YUVA444P9BE] = { 1, 1 },
    [AV_PIX_FMT_YUVA444P9LE] = { 1, 1 },
    [AV_PIX_FMT_YUVA420P10BE]= { 1, 1 },
    [AV_PIX_FMT_YUVA420P10LE]= { 1, 1 },
    [AV_PIX_FMT_YUVA422P10BE]= { 1, 1 },
    [AV_PIX_FMT_YUVA422P10LE]= { 1, 1 },
    [AV_PIX_FMT_YUVA444P10BE]= { 1, 1 },
    [AV_PIX_FMT_YUVA444P10LE]= { 1, 1 },
    [AV_PIX_FMT_YUVA420P16BE]= { 1, 1 },
    [AV_PIX_FMT_YUVA420P16LE]= { 1, 1 },
    [AV_PIX_FMT_YUVA422P16BE]= { 1, 1 },
    [AV_PIX_FMT_YUVA422P16LE]= { 1, 1 },
    [AV_PIX_FMT_YUVA444P16BE]= { 1, 1 },
    [AV_PIX_FMT_YUVA444P16LE]= { 1, 1 },
149 150
    [AV_PIX_FMT_RGB48BE]     = { 1, 1 },
    [AV_PIX_FMT_RGB48LE]     = { 1, 1 },
151 152
    [AV_PIX_FMT_RGBA64BE]    = { 1, 1, 1 },
    [AV_PIX_FMT_RGBA64LE]    = { 1, 1, 1 },
153 154 155 156 157 158 159 160 161 162 163 164 165 166 167 168 169 170
    [AV_PIX_FMT_RGB565BE]    = { 1, 1 },
    [AV_PIX_FMT_RGB565LE]    = { 1, 1 },
    [AV_PIX_FMT_RGB555BE]    = { 1, 1 },
    [AV_PIX_FMT_RGB555LE]    = { 1, 1 },
    [AV_PIX_FMT_BGR565BE]    = { 1, 1 },
    [AV_PIX_FMT_BGR565LE]    = { 1, 1 },
    [AV_PIX_FMT_BGR555BE]    = { 1, 1 },
    [AV_PIX_FMT_BGR555LE]    = { 1, 1 },
    [AV_PIX_FMT_YUV420P16LE] = { 1, 1 },
    [AV_PIX_FMT_YUV420P16BE] = { 1, 1 },
    [AV_PIX_FMT_YUV422P16LE] = { 1, 1 },
    [AV_PIX_FMT_YUV422P16BE] = { 1, 1 },
    [AV_PIX_FMT_YUV444P16LE] = { 1, 1 },
    [AV_PIX_FMT_YUV444P16BE] = { 1, 1 },
    [AV_PIX_FMT_RGB444LE]    = { 1, 1 },
    [AV_PIX_FMT_RGB444BE]    = { 1, 1 },
    [AV_PIX_FMT_BGR444LE]    = { 1, 1 },
    [AV_PIX_FMT_BGR444BE]    = { 1, 1 },
171
    [AV_PIX_FMT_YA8]         = { 1, 1 },
172 173
    [AV_PIX_FMT_YA16BE]      = { 1, 0 },
    [AV_PIX_FMT_YA16LE]      = { 1, 0 },
174 175
    [AV_PIX_FMT_BGR48BE]     = { 1, 1 },
    [AV_PIX_FMT_BGR48LE]     = { 1, 1 },
176 177
    [AV_PIX_FMT_BGRA64BE]    = { 1, 1, 1 },
    [AV_PIX_FMT_BGRA64LE]    = { 1, 1, 1 },
178 179 180 181
    [AV_PIX_FMT_YUV420P9BE]  = { 1, 1 },
    [AV_PIX_FMT_YUV420P9LE]  = { 1, 1 },
    [AV_PIX_FMT_YUV420P10BE] = { 1, 1 },
    [AV_PIX_FMT_YUV420P10LE] = { 1, 1 },
182 183 184 185
    [AV_PIX_FMT_YUV420P12BE] = { 1, 1 },
    [AV_PIX_FMT_YUV420P12LE] = { 1, 1 },
    [AV_PIX_FMT_YUV420P14BE] = { 1, 1 },
    [AV_PIX_FMT_YUV420P14LE] = { 1, 1 },
186 187 188 189
    [AV_PIX_FMT_YUV422P9BE]  = { 1, 1 },
    [AV_PIX_FMT_YUV422P9LE]  = { 1, 1 },
    [AV_PIX_FMT_YUV422P10BE] = { 1, 1 },
    [AV_PIX_FMT_YUV422P10LE] = { 1, 1 },
190 191 192 193
    [AV_PIX_FMT_YUV422P12BE] = { 1, 1 },
    [AV_PIX_FMT_YUV422P12LE] = { 1, 1 },
    [AV_PIX_FMT_YUV422P14BE] = { 1, 1 },
    [AV_PIX_FMT_YUV422P14LE] = { 1, 1 },
194 195 196 197
    [AV_PIX_FMT_YUV444P9BE]  = { 1, 1 },
    [AV_PIX_FMT_YUV444P9LE]  = { 1, 1 },
    [AV_PIX_FMT_YUV444P10BE] = { 1, 1 },
    [AV_PIX_FMT_YUV444P10LE] = { 1, 1 },
198 199 200 201
    [AV_PIX_FMT_YUV444P12BE] = { 1, 1 },
    [AV_PIX_FMT_YUV444P12LE] = { 1, 1 },
    [AV_PIX_FMT_YUV444P14BE] = { 1, 1 },
    [AV_PIX_FMT_YUV444P14LE] = { 1, 1 },
202
    [AV_PIX_FMT_GBRP]        = { 1, 1 },
203 204 205 206 207 208 209 210
    [AV_PIX_FMT_GBRP9LE]     = { 1, 1 },
    [AV_PIX_FMT_GBRP9BE]     = { 1, 1 },
    [AV_PIX_FMT_GBRP10LE]    = { 1, 1 },
    [AV_PIX_FMT_GBRP10BE]    = { 1, 1 },
    [AV_PIX_FMT_GBRP12LE]    = { 1, 1 },
    [AV_PIX_FMT_GBRP12BE]    = { 1, 1 },
    [AV_PIX_FMT_GBRP14LE]    = { 1, 1 },
    [AV_PIX_FMT_GBRP14BE]    = { 1, 1 },
211 212
    [AV_PIX_FMT_GBRP16LE]    = { 1, 0 },
    [AV_PIX_FMT_GBRP16BE]    = { 1, 0 },
213 214 215
    [AV_PIX_FMT_GBRAP]       = { 1, 1 },
    [AV_PIX_FMT_GBRAP16LE]   = { 1, 0 },
    [AV_PIX_FMT_GBRAP16BE]   = { 1, 0 },
216 217 218 219 220 221 222 223 224 225 226 227
    [AV_PIX_FMT_BAYER_BGGR8] = { 1, 0 },
    [AV_PIX_FMT_BAYER_RGGB8] = { 1, 0 },
    [AV_PIX_FMT_BAYER_GBRG8] = { 1, 0 },
    [AV_PIX_FMT_BAYER_GRBG8] = { 1, 0 },
    [AV_PIX_FMT_BAYER_BGGR16LE] = { 1, 0 },
    [AV_PIX_FMT_BAYER_BGGR16BE] = { 1, 0 },
    [AV_PIX_FMT_BAYER_RGGB16LE] = { 1, 0 },
    [AV_PIX_FMT_BAYER_RGGB16BE] = { 1, 0 },
    [AV_PIX_FMT_BAYER_GBRG16LE] = { 1, 0 },
    [AV_PIX_FMT_BAYER_GBRG16BE] = { 1, 0 },
    [AV_PIX_FMT_BAYER_GRBG16LE] = { 1, 0 },
    [AV_PIX_FMT_BAYER_GRBG16BE] = { 1, 0 },
228 229
    [AV_PIX_FMT_XYZ12BE]     = { 1, 1, 1 },
    [AV_PIX_FMT_XYZ12LE]     = { 1, 1, 1 },
230
    [AV_PIX_FMT_AYUV64LE]    = { 1, 1},
231 232
};

233
int sws_isSupportedInput(enum AVPixelFormat pix_fmt)
234
{
235
    return (unsigned)pix_fmt < AV_PIX_FMT_NB ?
236
           format_entries[pix_fmt].is_supported_in : 0;
237 238
}

239
int sws_isSupportedOutput(enum AVPixelFormat pix_fmt)
240
{
241
    return (unsigned)pix_fmt < AV_PIX_FMT_NB ?
242
           format_entries[pix_fmt].is_supported_out : 0;
243 244
}

245 246 247 248 249 250
int sws_isSupportedEndiannessConversion(enum AVPixelFormat pix_fmt)
{
    return (unsigned)pix_fmt < AV_PIX_FMT_NB ?
           format_entries[pix_fmt].is_supported_endianness : 0;
}

251 252
static double getSplineCoeff(double a, double b, double c, double d,
                             double dist)
253
{
254 255 256 257 258 259 260 261
    if (dist <= 1.0)
        return ((d * dist + c) * dist + b) * dist + a;
    else
        return getSplineCoeff(0.0,
                               b + 2.0 * c + 3.0 * d,
                               c + 3.0 * d,
                              -b - 3.0 * c - 6.0 * d,
                              dist - 1.0);
262 263
}

264 265
static av_cold int get_local_pos(SwsContext *s, int chr_subsample, int pos, int dir)
{
266
    if (pos == -1 || pos <= -513) {
267 268 269 270 271 272
        pos = (128 << chr_subsample) - 128;
    }
    pos += 128; // relative to ideal left edge
    return pos >> chr_subsample;
}

273 274 275 276 277 278 279 280 281 282 283 284 285 286 287 288 289 290 291 292
typedef struct {
    int flag;                   ///< flag associated to the algorithm
    const char *description;    ///< human-readable description
    int size_factor;            ///< size factor used when initing the filters
} ScaleAlgorithm;

static const ScaleAlgorithm scale_algorithms[] = {
    { SWS_AREA,          "area averaging",                  1 /* downscale only, for upscale it is bilinear */ },
    { SWS_BICUBIC,       "bicubic",                         4 },
    { SWS_BICUBLIN,      "luma bicubic / chroma bilinear", -1 },
    { SWS_BILINEAR,      "bilinear",                        2 },
    { SWS_FAST_BILINEAR, "fast bilinear",                  -1 },
    { SWS_GAUSS,         "Gaussian",                        8 /* infinite ;) */ },
    { SWS_LANCZOS,       "Lanczos",                        -1 /* custom */ },
    { SWS_POINT,         "nearest neighbor / point",       -1 },
    { SWS_SINC,          "sinc",                           20 /* infinite ;) */ },
    { SWS_SPLINE,        "bicubic spline",                 20 /* infinite :)*/ },
    { SWS_X,             "experimental",                    8 },
};

293 294 295 296 297
static av_cold int initFilter(int16_t **outFilter, int32_t **filterPos,
                              int *outFilterSize, int xInc, int srcW,
                              int dstW, int filterAlign, int one,
                              int flags, int cpu_flags,
                              SwsVector *srcFilter, SwsVector *dstFilter,
298
                              double param[2], int srcPos, int dstPos)
299 300 301 302 303
{
    int i;
    int filterSize;
    int filter2Size;
    int minFilterSize;
304 305
    int64_t *filter    = NULL;
    int64_t *filter2   = NULL;
306
    const int64_t fone = 1LL << (54 - FFMIN(av_log2(srcW/dstW), 8));
307
    int ret            = -1;
308

309
    emms_c(); // FIXME should not be required but IS (even for non-MMX versions)
310

311
    // NOTE: the +3 is for the MMX(+1) / SSE(+3) scaler which reads over the end
312
    FF_ALLOC_ARRAY_OR_GOTO(NULL, *filterPos, (dstW + 3), sizeof(**filterPos), fail);
313

314
    if (FFABS(xInc - 0x10000) < 10 && srcPos == dstPos) { // unscaled
315
        int i;
316
        filterSize = 1;
317 318
        FF_ALLOCZ_ARRAY_OR_GOTO(NULL, filter,
                                dstW, sizeof(*filter) * filterSize, fail);
319

320 321 322
        for (i = 0; i < dstW; i++) {
            filter[i * filterSize] = fone;
            (*filterPos)[i]        = i;
323
        }
324
    } else if (flags & SWS_POINT) { // lame looking point sampling mode
325
        int i;
326
        int64_t xDstInSrc;
327
        filterSize = 1;
328 329
        FF_ALLOC_ARRAY_OR_GOTO(NULL, filter,
                               dstW, sizeof(*filter) * filterSize, fail);
330

331
        xDstInSrc = ((dstPos*(int64_t)xInc)>>8) - ((srcPos*0x8000LL)>>7);
332 333
        for (i = 0; i < dstW; i++) {
            int xx = (xDstInSrc - ((filterSize - 1) << 15) + (1 << 15)) >> 16;
334

335 336 337
            (*filterPos)[i] = xx;
            filter[i]       = fone;
            xDstInSrc      += xInc;
338
        }
339 340
    } else if ((xInc <= (1 << 16) && (flags & SWS_AREA)) ||
               (flags & SWS_FAST_BILINEAR)) { // bilinear upscale
341
        int i;
342
        int64_t xDstInSrc;
343
        filterSize = 2;
344 345
        FF_ALLOC_ARRAY_OR_GOTO(NULL, filter,
                               dstW, sizeof(*filter) * filterSize, fail);
346

347
        xDstInSrc = ((dstPos*(int64_t)xInc)>>8) - ((srcPos*0x8000LL)>>7);
348 349
        for (i = 0; i < dstW; i++) {
            int xx = (xDstInSrc - ((filterSize - 1) << 15) + (1 << 15)) >> 16;
350 351
            int j;

352 353 354
            (*filterPos)[i] = xx;
            // bilinear upscale / linear interpolate / area averaging
            for (j = 0; j < filterSize; j++) {
355
                int64_t coeff= fone - FFABS(((int64_t)xx<<16) - xDstInSrc)*(fone>>16);
356 357 358
                if (coeff < 0)
                    coeff = 0;
                filter[i * filterSize + j] = coeff;
359 360
                xx++;
            }
361
            xDstInSrc += xInc;
362 363
        }
    } else {
364
        int64_t xDstInSrc;
365 366 367
        int sizeFactor = -1;

        for (i = 0; i < FF_ARRAY_ELEMS(scale_algorithms); i++) {
368
            if (flags & scale_algorithms[i].flag && scale_algorithms[i].size_factor > 0) {
369 370 371
                sizeFactor = scale_algorithms[i].size_factor;
                break;
            }
372
        }
373 374 375
        if (flags & SWS_LANCZOS)
            sizeFactor = param[0] != SWS_PARAM_DEFAULT ? ceil(2 * param[0]) : 6;
        av_assert0(sizeFactor > 0);
376

377 378 379 380
        if (xInc <= 1 << 16)
            filterSize = 1 + sizeFactor;    // upscale
        else
            filterSize = 1 + (sizeFactor * srcW + dstW - 1) / dstW;
381

382 383
        filterSize = FFMIN(filterSize, srcW - 2);
        filterSize = FFMAX(filterSize, 1);
384

385 386
        FF_ALLOC_ARRAY_OR_GOTO(NULL, filter,
                               dstW, sizeof(*filter) * filterSize, fail);
387

388
        xDstInSrc = ((dstPos*(int64_t)xInc)>>7) - ((srcPos*0x10000LL)>>7);
389
        for (i = 0; i < dstW; i++) {
390
            int xx = (xDstInSrc - (filterSize - 2) * (1LL<<16)) / (1 << 17);
391
            int j;
392 393
            (*filterPos)[i] = xx;
            for (j = 0; j < filterSize; j++) {
394
                int64_t d = (FFABS(((int64_t)xx * (1 << 17)) - xDstInSrc)) << 13;
395 396 397
                double floatd;
                int64_t coeff;

398 399 400
                if (xInc > 1 << 16)
                    d = d * dstW / srcW;
                floatd = d * (1.0 / (1 << 30));
401 402

                if (flags & SWS_BICUBIC) {
403 404
                    int64_t B = (param[0] != SWS_PARAM_DEFAULT ? param[0] :   0) * (1 << 24);
                    int64_t C = (param[1] != SWS_PARAM_DEFAULT ? param[1] : 0.6) * (1 << 24);
405

406
                    if (d >= 1LL << 31) {
407 408 409 410 411
                        coeff = 0.0;
                    } else {
                        int64_t dd  = (d  * d) >> 30;
                        int64_t ddd = (dd * d) >> 30;

412 413 414 415
                        if (d < 1LL << 30)
                            coeff =  (12 * (1 << 24) -  9 * B - 6 * C) * ddd +
                                    (-18 * (1 << 24) + 12 * B + 6 * C) *  dd +
                                      (6 * (1 << 24) -  2 * B)         * (1 << 30);
416
                        else
417 418 419 420
                            coeff =      (-B -  6 * C) * ddd +
                                      (6 * B + 30 * C) * dd  +
                                    (-12 * B - 48 * C) * d   +
                                      (8 * B + 24 * C) * (1 << 30);
421
                    }
422
                    coeff /= (1LL<<54)/fone;
423
                }
424 425 426 427 428
#if 0
                else if (flags & SWS_X) {
                    double p  = param ? param * 0.01 : 0.3;
                    coeff     = d ? sin(d * M_PI) / (d * M_PI) : 1.0;
                    coeff    *= pow(2.0, -p * d * d);
429
                }
430
#endif
431
                else if (flags & SWS_X) {
432
                    double A = param[0] != SWS_PARAM_DEFAULT ? param[0] : 1.0;
433 434
                    double c;

435 436
                    if (floatd < 1.0)
                        c = cos(floatd * M_PI);
437
                    else
438 439 440
                        c = -1.0;
                    if (c < 0.0)
                        c = -pow(-c, A);
441
                    else
442 443
                        c = pow(c, A);
                    coeff = (c * 0.5 + 0.5) * fone;
444
                } else if (flags & SWS_AREA) {
445 446 447 448 449 450 451 452
                    int64_t d2 = d - (1 << 29);
                    if (d2 * xInc < -(1LL << (29 + 16)))
                        coeff = 1.0 * (1LL << (30 + 16));
                    else if (d2 * xInc < (1LL << (29 + 16)))
                        coeff = -d2 * xInc + (1LL << (29 + 16));
                    else
                        coeff = 0.0;
                    coeff *= fone >> (30 + 16);
453
                } else if (flags & SWS_GAUSS) {
454
                    double p = param[0] != SWS_PARAM_DEFAULT ? param[0] : 3.0;
455
                    coeff = exp2(-p * floatd * floatd) * fone;
456
                } else if (flags & SWS_SINC) {
457
                    coeff = (d ? sin(floatd * M_PI) / (floatd * M_PI) : 1.0) * fone;
458
                } else if (flags & SWS_LANCZOS) {
459 460 461 462 463
                    double p = param[0] != SWS_PARAM_DEFAULT ? param[0] : 3.0;
                    coeff = (d ? sin(floatd * M_PI) * sin(floatd * M_PI / p) /
                             (floatd * floatd * M_PI * M_PI / p) : 1.0) * fone;
                    if (floatd > p)
                        coeff = 0;
464
                } else if (flags & SWS_BILINEAR) {
465 466 467
                    coeff = (1 << 30) - d;
                    if (coeff < 0)
                        coeff = 0;
468 469
                    coeff *= fone >> 30;
                } else if (flags & SWS_SPLINE) {
470 471
                    double p = -2.196152422706632;
                    coeff = getSplineCoeff(1.0, 0.0, p, -p - 1.0, floatd) * fone;
472
                } else {
473
                    av_assert0(0);
474 475
                }

476
                filter[i * filterSize + j] = coeff;
477 478
                xx++;
            }
479
            xDstInSrc += 2 * xInc;
480 481 482 483
        }
    }

    /* apply src & dst Filter to filter -> filter2
484 485
     * av_free(filter);
     */
486
    av_assert0(filterSize > 0);
487 488 489 490 491
    filter2Size = filterSize;
    if (srcFilter)
        filter2Size += srcFilter->length - 1;
    if (dstFilter)
        filter2Size += dstFilter->length - 1;
492
    av_assert0(filter2Size > 0);
493
    FF_ALLOCZ_ARRAY_OR_GOTO(NULL, filter2, dstW, filter2Size * sizeof(*filter2), fail);
494 495

    for (i = 0; i < dstW; i++) {
496 497
        int j, k;

498 499 500 501 502
        if (srcFilter) {
            for (k = 0; k < srcFilter->length; k++) {
                for (j = 0; j < filterSize; j++)
                    filter2[i * filter2Size + k + j] +=
                        srcFilter->coeff[k] * filter[i * filterSize + j];
503 504
            }
        } else {
505 506
            for (j = 0; j < filterSize; j++)
                filter2[i * filter2Size + j] = filter[i * filterSize + j];
507
        }
508
        // FIXME dstFilter
509

510
        (*filterPos)[i] += (filterSize - 1) / 2 - (filter2Size - 1) / 2;
511 512 513 514 515
    }
    av_freep(&filter);

    /* try to reduce the filter-size (step1 find size and shift left) */
    // Assume it is near normalized (*0.5 or *2.0 is OK but * 0.001 is not).
516 517 518
    minFilterSize = 0;
    for (i = dstW - 1; i >= 0; i--) {
        int min = filter2Size;
519
        int j;
520
        int64_t cutOff = 0.0;
521 522

        /* get rid of near zero elements on the left by shifting left */
523
        for (j = 0; j < filter2Size; j++) {
524
            int k;
525
            cutOff += FFABS(filter2[i * filter2Size]);
526

527 528
            if (cutOff > SWS_MAX_REDUCE_CUTOFF * fone)
                break;
529

530 531 532 533
            /* preserve monotonicity because the core can't handle the
             * filter otherwise */
            if (i < dstW - 1 && (*filterPos)[i] >= (*filterPos)[i + 1])
                break;
534 535

            // move filter coefficients left
536 537 538
            for (k = 1; k < filter2Size; k++)
                filter2[i * filter2Size + k - 1] = filter2[i * filter2Size + k];
            filter2[i * filter2Size + k - 1] = 0;
539 540 541
            (*filterPos)[i]++;
        }

542
        cutOff = 0;
543
        /* count near zeros on the right */
544 545
        for (j = filter2Size - 1; j > 0; j--) {
            cutOff += FFABS(filter2[i * filter2Size + j]);
546

547 548
            if (cutOff > SWS_MAX_REDUCE_CUTOFF * fone)
                break;
549 550 551
            min--;
        }

552 553
        if (min > minFilterSize)
            minFilterSize = min;
554 555
    }

556
    if (PPC_ALTIVEC(cpu_flags)) {
557
        // we can handle the special case 4, so we don't want to go the full 8
558 559 560
        if (minFilterSize < 5)
            filterAlign = 4;

561 562 563
        /* We really don't want to waste our time doing useless computation, so
         * fall back on the scalar C code for very small filters.
         * Vectorizing is worth it only if you have a decent-sized vector. */
564 565 566 567
        if (minFilterSize < 3)
            filterAlign = 1;
    }

568
    if (HAVE_MMX && cpu_flags & AV_CPU_FLAG_MMX) {
569 570
        // special case for unscaled vertical filtering
        if (minFilterSize == 1 && filterAlign == 2)
571
            filterAlign = 1;
572 573
    }

574
    av_assert0(minFilterSize > 0);
575
    filterSize = (minFilterSize + (filterAlign - 1)) & (~(filterAlign - 1));
576
    av_assert0(filterSize > 0);
577
    filter = av_malloc_array(dstW, filterSize * sizeof(*filter));
578 579
    if (!filter)
        goto fail;
580
    if (filterSize >= MAX_FILTER_SIZE * 16 /
581
                      ((flags & SWS_ACCURATE_RND) ? APCK_SIZE : 16)) {
582
        ret = RETCODE_USE_CASCADE;
583
        goto fail;
584
    }
585
    *outFilterSize = filterSize;
586

587 588 589 590
    if (flags & SWS_PRINT_INFO)
        av_log(NULL, AV_LOG_VERBOSE,
               "SwScaler: reducing / aligning filtersize %d -> %d\n",
               filter2Size, filterSize);
591
    /* try to reduce the filter-size (step2 reduce it) */
592
    for (i = 0; i < dstW; i++) {
593 594
        int j;

595 596 597 598 599 600 601
        for (j = 0; j < filterSize; j++) {
            if (j >= filter2Size)
                filter[i * filterSize + j] = 0;
            else
                filter[i * filterSize + j] = filter2[i * filter2Size + j];
            if ((flags & SWS_BITEXACT) && j >= minFilterSize)
                filter[i * filterSize + j] = 0;
602 603 604
        }
    }

605
    // FIXME try to align filterPos if possible
606

607
    // fix borders
608
    for (i = 0; i < dstW; i++) {
609 610 611
        int j;
        if ((*filterPos)[i] < 0) {
            // move filter coefficients left to compensate for filterPos
612 613 614 615
            for (j = 1; j < filterSize; j++) {
                int left = FFMAX(j + (*filterPos)[i], 0);
                filter[i * filterSize + left] += filter[i * filterSize + j];
                filter[i * filterSize + j]     = 0;
616 617 618 619 620
            }
            (*filterPos)[i]= 0;
        }

        if ((*filterPos)[i] + filterSize > srcW) {
621
            int shift = (*filterPos)[i] + FFMIN(filterSize - srcW, 0);
622
            int64_t acc = 0;
623

624 625 626 627 628
            for (j = filterSize - 1; j >= 0; j--) {
                if ((*filterPos)[i] + j >= srcW) {
                    acc += filter[i * filterSize + j];
                    filter[i * filterSize + j] = 0;
                }
629
            }
630 631 632 633 634 635 636 637
            for (j = filterSize - 1; j >= 0; j--) {
                if (j < shift) {
                    filter[i * filterSize + j] = 0;
                } else {
                    filter[i * filterSize + j] = filter[i * filterSize + j - shift];
                }
            }

638
            (*filterPos)[i]-= shift;
639
            filter[i * filterSize + srcW - 1 - (*filterPos)[i]] += acc;
640
        }
641 642 643 644 645 646 647
        av_assert0((*filterPos)[i] >= 0);
        av_assert0((*filterPos)[i] < srcW);
        if ((*filterPos)[i] + filterSize > srcW) {
            for (j = 0; j < filterSize; j++) {
                av_assert0((*filterPos)[i] + j < srcW || !filter[i * filterSize + j]);
            }
        }
648 649 650 651
    }

    // Note the +1 is for the MMX scaler which reads over the end
    /* align at 16 for AltiVec (needed by hScale_altivec_real) */
652 653
    FF_ALLOCZ_ARRAY_OR_GOTO(NULL, *outFilter,
                            (dstW + 3), *outFilterSize * sizeof(int16_t), fail);
654 655

    /* normalize & store in outFilter */
656
    for (i = 0; i < dstW; i++) {
657
        int j;
658 659
        int64_t error = 0;
        int64_t sum   = 0;
660

661 662
        for (j = 0; j < filterSize; j++) {
            sum += filter[i * filterSize + j];
663
        }
664
        sum = (sum + one / 2) / one;
665 666 667 668
        if (!sum) {
            av_log(NULL, AV_LOG_WARNING, "SwScaler: zero vector in scaling\n");
            sum = 1;
        }
669 670 671 672 673
        for (j = 0; j < *outFilterSize; j++) {
            int64_t v = filter[i * filterSize + j] + error;
            int intV  = ROUNDED_DIV(v, sum);
            (*outFilter)[i * (*outFilterSize) + j] = intV;
            error                                  = v - intV * sum;
674 675 676
        }
    }

677 678 679 680 681 682
    (*filterPos)[dstW + 0] =
    (*filterPos)[dstW + 1] =
    (*filterPos)[dstW + 2] = (*filterPos)[dstW - 1]; /* the MMX/SSE scaler will
                                                      * read over the end */
    for (i = 0; i < *outFilterSize; i++) {
        int k = (dstW - 1) * (*outFilterSize) + i;
683 684 685
        (*outFilter)[k + 1 * (*outFilterSize)] =
        (*outFilter)[k + 2 * (*outFilterSize)] =
        (*outFilter)[k + 3 * (*outFilterSize)] = (*outFilter)[k];
686 687
    }

688 689
    ret = 0;

690
fail:
691
    if(ret < 0)
692
        av_log(NULL, ret == RETCODE_USE_CASCADE ? AV_LOG_DEBUG : AV_LOG_ERROR, "sws: initFilter failed\n");
693 694 695 696 697
    av_free(filter);
    av_free(filter2);
    return ret;
}

698 699
static void fill_rgb2yuv_table(SwsContext *c, const int table[4], int dstRange)
{
700
    int64_t W, V, Z, Cy, Cu, Cv;
701 702 703 704 705 706
    int64_t vr =  table[0];
    int64_t ub =  table[1];
    int64_t ug = -table[2];
    int64_t vg = -table[3];
    int64_t ONE = 65536;
    int64_t cy = ONE;
707 708 709 710 711 712 713 714 715 716 717 718 719 720 721
    uint8_t *p = (uint8_t*)c->input_rgb2yuv_table;
    int i;
    static const int8_t map[] = {
    BY_IDX, GY_IDX, -1    , BY_IDX, BY_IDX, GY_IDX, -1    , BY_IDX,
    RY_IDX, -1    , GY_IDX, RY_IDX, RY_IDX, -1    , GY_IDX, RY_IDX,
    RY_IDX, GY_IDX, -1    , RY_IDX, RY_IDX, GY_IDX, -1    , RY_IDX,
    BY_IDX, -1    , GY_IDX, BY_IDX, BY_IDX, -1    , GY_IDX, BY_IDX,
    BU_IDX, GU_IDX, -1    , BU_IDX, BU_IDX, GU_IDX, -1    , BU_IDX,
    RU_IDX, -1    , GU_IDX, RU_IDX, RU_IDX, -1    , GU_IDX, RU_IDX,
    RU_IDX, GU_IDX, -1    , RU_IDX, RU_IDX, GU_IDX, -1    , RU_IDX,
    BU_IDX, -1    , GU_IDX, BU_IDX, BU_IDX, -1    , GU_IDX, BU_IDX,
    BV_IDX, GV_IDX, -1    , BV_IDX, BV_IDX, GV_IDX, -1    , BV_IDX,
    RV_IDX, -1    , GV_IDX, RV_IDX, RV_IDX, -1    , GV_IDX, RV_IDX,
    RV_IDX, GV_IDX, -1    , RV_IDX, RV_IDX, GV_IDX, -1    , RV_IDX,
    BV_IDX, -1    , GV_IDX, BV_IDX, BV_IDX, -1    , GV_IDX, BV_IDX,
722 723 724 725 726 727 728 729 730 731 732
    RY_IDX, BY_IDX, RY_IDX, BY_IDX, RY_IDX, BY_IDX, RY_IDX, BY_IDX,
    BY_IDX, RY_IDX, BY_IDX, RY_IDX, BY_IDX, RY_IDX, BY_IDX, RY_IDX,
    GY_IDX, -1    , GY_IDX, -1    , GY_IDX, -1    , GY_IDX, -1    ,
    -1    , GY_IDX, -1    , GY_IDX, -1    , GY_IDX, -1    , GY_IDX,
    RU_IDX, BU_IDX, RU_IDX, BU_IDX, RU_IDX, BU_IDX, RU_IDX, BU_IDX,
    BU_IDX, RU_IDX, BU_IDX, RU_IDX, BU_IDX, RU_IDX, BU_IDX, RU_IDX,
    GU_IDX, -1    , GU_IDX, -1    , GU_IDX, -1    , GU_IDX, -1    ,
    -1    , GU_IDX, -1    , GU_IDX, -1    , GU_IDX, -1    , GU_IDX,
    RV_IDX, BV_IDX, RV_IDX, BV_IDX, RV_IDX, BV_IDX, RV_IDX, BV_IDX,
    BV_IDX, RV_IDX, BV_IDX, RV_IDX, BV_IDX, RV_IDX, BV_IDX, RV_IDX,
    GV_IDX, -1    , GV_IDX, -1    , GV_IDX, -1    , GV_IDX, -1    ,
733 734 735 736 737 738 739 740 741 742 743 744
    -1    , GV_IDX, -1    , GV_IDX, -1    , GV_IDX, -1    , GV_IDX, //23
    -1    , -1    , -1    , -1    , -1    , -1    , -1    , -1    , //24
    -1    , -1    , -1    , -1    , -1    , -1    , -1    , -1    , //25
    -1    , -1    , -1    , -1    , -1    , -1    , -1    , -1    , //26
    -1    , -1    , -1    , -1    , -1    , -1    , -1    , -1    , //27
    -1    , -1    , -1    , -1    , -1    , -1    , -1    , -1    , //28
    -1    , -1    , -1    , -1    , -1    , -1    , -1    , -1    , //29
    -1    , -1    , -1    , -1    , -1    , -1    , -1    , -1    , //30
    -1    , -1    , -1    , -1    , -1    , -1    , -1    , -1    , //31
    BY_IDX, GY_IDX, RY_IDX, -1    , -1    , -1    , -1    , -1    , //32
    BU_IDX, GU_IDX, RU_IDX, -1    , -1    , -1    , -1    , -1    , //33
    BV_IDX, GV_IDX, RV_IDX, -1    , -1    , -1    , -1    , -1    , //34
745
    };
746 747 748 749 750 751 752 753 754 755 756

    dstRange = 0; //FIXME range = 1 is handled elsewhere

    if (!dstRange) {
        cy = cy * 255 / 219;
    } else {
        vr = vr * 224 / 255;
        ub = ub * 224 / 255;
        ug = ug * 224 / 255;
        vg = vg * 224 / 255;
    }
757 758
    W = ROUNDED_DIV(ONE*ONE*ug, ub);
    V = ROUNDED_DIV(ONE*ONE*vg, vr);
759
    Z = ONE*ONE-W-V;
760

761 762 763
    Cy = ROUNDED_DIV(cy*Z, ONE);
    Cu = ROUNDED_DIV(ub*Z, ONE);
    Cv = ROUNDED_DIV(vr*Z, ONE);
764

765 766 767
    c->input_rgb2yuv_table[RY_IDX] = -ROUNDED_DIV((1 << RGB2YUV_SHIFT)*V        , Cy);
    c->input_rgb2yuv_table[GY_IDX] =  ROUNDED_DIV((1 << RGB2YUV_SHIFT)*ONE*ONE  , Cy);
    c->input_rgb2yuv_table[BY_IDX] = -ROUNDED_DIV((1 << RGB2YUV_SHIFT)*W        , Cy);
768

769 770 771
    c->input_rgb2yuv_table[RU_IDX] =  ROUNDED_DIV((1 << RGB2YUV_SHIFT)*V        , Cu);
    c->input_rgb2yuv_table[GU_IDX] = -ROUNDED_DIV((1 << RGB2YUV_SHIFT)*ONE*ONE  , Cu);
    c->input_rgb2yuv_table[BU_IDX] =  ROUNDED_DIV((1 << RGB2YUV_SHIFT)*(Z+W)    , Cu);
772

773 774 775
    c->input_rgb2yuv_table[RV_IDX] =  ROUNDED_DIV((1 << RGB2YUV_SHIFT)*(V+Z)    , Cv);
    c->input_rgb2yuv_table[GV_IDX] = -ROUNDED_DIV((1 << RGB2YUV_SHIFT)*ONE*ONE  , Cv);
    c->input_rgb2yuv_table[BV_IDX] =  ROUNDED_DIV((1 << RGB2YUV_SHIFT)*W        , Cv);
776

777
    if(/*!dstRange && */!memcmp(table, ff_yuv2rgb_coeffs[SWS_CS_DEFAULT], sizeof(ff_yuv2rgb_coeffs[SWS_CS_DEFAULT]))) {
778 779 780 781 782 783 784 785 786 787
        c->input_rgb2yuv_table[BY_IDX] =  ((int)(0.114 * 219 / 255 * (1 << RGB2YUV_SHIFT) + 0.5));
        c->input_rgb2yuv_table[BV_IDX] = (-(int)(0.081 * 224 / 255 * (1 << RGB2YUV_SHIFT) + 0.5));
        c->input_rgb2yuv_table[BU_IDX] =  ((int)(0.500 * 224 / 255 * (1 << RGB2YUV_SHIFT) + 0.5));
        c->input_rgb2yuv_table[GY_IDX] =  ((int)(0.587 * 219 / 255 * (1 << RGB2YUV_SHIFT) + 0.5));
        c->input_rgb2yuv_table[GV_IDX] = (-(int)(0.419 * 224 / 255 * (1 << RGB2YUV_SHIFT) + 0.5));
        c->input_rgb2yuv_table[GU_IDX] = (-(int)(0.331 * 224 / 255 * (1 << RGB2YUV_SHIFT) + 0.5));
        c->input_rgb2yuv_table[RY_IDX] =  ((int)(0.299 * 219 / 255 * (1 << RGB2YUV_SHIFT) + 0.5));
        c->input_rgb2yuv_table[RV_IDX] =  ((int)(0.500 * 224 / 255 * (1 << RGB2YUV_SHIFT) + 0.5));
        c->input_rgb2yuv_table[RU_IDX] = (-(int)(0.169 * 224 / 255 * (1 << RGB2YUV_SHIFT) + 0.5));
    }
788 789
    for(i=0; i<FF_ARRAY_ELEMS(map); i++)
        AV_WL16(p + 16*4 + 2*i, map[i] >= 0 ? c->input_rgb2yuv_table[map[i]] : 0);
790 791
}

792 793 794 795 796
static void fill_xyztables(struct SwsContext *c)
{
    int i;
    double xyzgamma = XYZ_GAMMA;
    double rgbgamma = 1.0 / RGB_GAMMA;
797 798
    double xyzgammainv = 1.0 / XYZ_GAMMA;
    double rgbgammainv = RGB_GAMMA;
799 800 801 802
    static const int16_t xyz2rgb_matrix[3][4] = {
        {13270, -6295, -2041},
        {-3969,  7682,   170},
        {  228,  -835,  4329} };
803 804 805 806 807
    static const int16_t rgb2xyz_matrix[3][4] = {
        {1689, 1464,  739},
        { 871, 2929,  296},
        {  79,  488, 3891} };
    static int16_t xyzgamma_tab[4096], rgbgamma_tab[4096], xyzgammainv_tab[4096], rgbgammainv_tab[4096];
808 809

    memcpy(c->xyz2rgb_matrix, xyz2rgb_matrix, sizeof(c->xyz2rgb_matrix));
810
    memcpy(c->rgb2xyz_matrix, rgb2xyz_matrix, sizeof(c->rgb2xyz_matrix));
811 812
    c->xyzgamma = xyzgamma_tab;
    c->rgbgamma = rgbgamma_tab;
813 814
    c->xyzgammainv = xyzgammainv_tab;
    c->rgbgammainv = rgbgammainv_tab;
815 816 817

    if (rgbgamma_tab[4095])
        return;
818 819 820

    /* set gamma vectors */
    for (i = 0; i < 4096; i++) {
821 822
        xyzgamma_tab[i] = lrint(pow(i / 4095.0, xyzgamma) * 4095.0);
        rgbgamma_tab[i] = lrint(pow(i / 4095.0, rgbgamma) * 4095.0);
823 824
        xyzgammainv_tab[i] = lrint(pow(i / 4095.0, xyzgammainv) * 4095.0);
        rgbgammainv_tab[i] = lrint(pow(i / 4095.0, rgbgammainv) * 4095.0);
825 826 827
    }
}

828 829 830
int sws_setColorspaceDetails(struct SwsContext *c, const int inv_table[4],
                             int srcRange, const int table[4], int dstRange,
                             int brightness, int contrast, int saturation)
831
{
832 833
    const AVPixFmtDescriptor *desc_dst;
    const AVPixFmtDescriptor *desc_src;
834
    int need_reinit = 0;
835

836 837 838 839
    handle_formats(c);
    desc_dst = av_pix_fmt_desc_get(c->dstFormat);
    desc_src = av_pix_fmt_desc_get(c->srcFormat);

840 841 842 843 844
    if(!isYUV(c->dstFormat) && !isGray(c->dstFormat))
        dstRange = 0;
    if(!isYUV(c->srcFormat) && !isGray(c->srcFormat))
        srcRange = 0;

845 846 847 848 849 850 851 852 853 854 855 856 857 858 859
    if (c->srcRange != srcRange ||
        c->dstRange != dstRange ||
        c->brightness != brightness ||
        c->contrast   != contrast ||
        c->saturation != saturation ||
        memcmp(c->srcColorspaceTable, inv_table, sizeof(int) * 4) ||
        memcmp(c->dstColorspaceTable,     table, sizeof(int) * 4)
    )
        need_reinit = 1;

    memmove(c->srcColorspaceTable, inv_table, sizeof(int) * 4);
    memmove(c->dstColorspaceTable, table, sizeof(int) * 4);



860 861 862 863 864
    c->brightness = brightness;
    c->contrast   = contrast;
    c->saturation = saturation;
    c->srcRange   = srcRange;
    c->dstRange   = dstRange;
865

866 867
    //The srcBpc check is possibly wrong but we seem to lack a definitive reference to test this
    //and what we have in ticket 2939 looks better with this check
868
    if (need_reinit && (c->srcBpc == 8 || !isYUV(c->srcFormat)))
869 870
        ff_sws_init_range_convert(c);

871 872 873
    c->dstFormatBpp = av_get_bits_per_pixel(desc_dst);
    c->srcFormatBpp = av_get_bits_per_pixel(desc_src);

874 875
    if (c->cascaded_context[c->cascaded_mainindex])
        return sws_setColorspaceDetails(c->cascaded_context[c->cascaded_mainindex],inv_table, srcRange,table, dstRange, brightness,  contrast, saturation);
876

877 878 879
    if (!need_reinit)
        return 0;

880 881 882 883 884 885 886 887 888 889 890 891 892 893 894 895 896 897 898 899 900 901 902 903 904 905 906 907 908 909 910 911 912 913 914 915 916 917 918 919 920 921 922 923 924 925 926 927 928 929 930 931 932 933 934 935 936 937 938 939 940 941 942 943 944
    if ((isYUV(c->dstFormat) || isGray(c->dstFormat)) && (isYUV(c->srcFormat) || isGray(c->srcFormat))) {
        if (!c->cascaded_context[0] &&
            memcmp(c->dstColorspaceTable, c->srcColorspaceTable, sizeof(int) * 4) &&
            c->srcW && c->srcH && c->dstW && c->dstH) {
            enum AVPixelFormat tmp_format;
            int tmp_width, tmp_height;
            int srcW = c->srcW;
            int srcH = c->srcH;
            int dstW = c->dstW;
            int dstH = c->dstH;
            int ret;
            av_log(c, AV_LOG_VERBOSE, "YUV color matrix differs for YUV->YUV, using intermediate RGB to convert\n");

            if (isNBPS(c->dstFormat) || is16BPS(c->dstFormat)) {
                if (isALPHA(c->srcFormat) && isALPHA(c->dstFormat)) {
                    tmp_format = AV_PIX_FMT_BGRA64;
                } else {
                    tmp_format = AV_PIX_FMT_BGR48;
                }
            } else {
                if (isALPHA(c->srcFormat) && isALPHA(c->dstFormat)) {
                    tmp_format = AV_PIX_FMT_BGRA;
                } else {
                    tmp_format = AV_PIX_FMT_BGR24;
                }
            }

            if (srcW*srcH > dstW*dstH) {
                tmp_width  = dstW;
                tmp_height = dstH;
            } else {
                tmp_width  = srcW;
                tmp_height = srcH;
            }

            ret = av_image_alloc(c->cascaded_tmp, c->cascaded_tmpStride,
                                tmp_width, tmp_height, tmp_format, 64);
            if (ret < 0)
                return ret;

            c->cascaded_context[0] = sws_alloc_set_opts(srcW, srcH, c->srcFormat,
                                                        tmp_width, tmp_height, tmp_format,
                                                        c->flags, c->param);
            if (!c->cascaded_context[0])
                return -1;

            c->cascaded_context[0]->alphablend = c->alphablend;
            ret = sws_init_context(c->cascaded_context[0], NULL , NULL);
            if (ret < 0)
                return ret;
            //we set both src and dst depending on that the RGB side will be ignored
            sws_setColorspaceDetails(c->cascaded_context[0], inv_table,
                                     srcRange, table, dstRange,
                                     brightness, contrast, saturation);

            c->cascaded_context[1] = sws_getContext(tmp_width, tmp_height, tmp_format,
                                                    dstW, dstH, c->dstFormat,
                                                    c->flags, NULL, NULL, c->param);
            if (!c->cascaded_context[1])
                return -1;
            sws_setColorspaceDetails(c->cascaded_context[1], inv_table,
                                     srcRange, table, dstRange,
                                     0, 1 << 16, 1 << 16);
            return 0;
        }
945
        return -1;
946
    }
947

948
    if (!isYUV(c->dstFormat) && !isGray(c->dstFormat)) {
949 950 951
        ff_yuv2rgb_c_init_tables(c, inv_table, srcRange, brightness,
                                 contrast, saturation);
        // FIXME factorize
952

953 954 955
        if (ARCH_PPC)
            ff_yuv2rgb_init_tables_ppc(c, inv_table, brightness,
                                       contrast, saturation);
956
    }
957 958 959

    fill_rgb2yuv_table(c, table, dstRange);

960 961 962
    return 0;
}

963 964 965
int sws_getColorspaceDetails(struct SwsContext *c, int **inv_table,
                             int *srcRange, int **table, int *dstRange,
                             int *brightness, int *contrast, int *saturation)
966
{
967
    if (!c )
968
        return -1;
969

970 971 972 973 974 975 976
    *inv_table  = c->srcColorspaceTable;
    *table      = c->dstColorspaceTable;
    *srcRange   = c->srcRange;
    *dstRange   = c->dstRange;
    *brightness = c->brightness;
    *contrast   = c->contrast;
    *saturation = c->saturation;
977 978 979 980

    return 0;
}

981
static int handle_jpeg(enum AVPixelFormat *format)
982 983
{
    switch (*format) {
984 985
    case AV_PIX_FMT_YUVJ420P:
        *format = AV_PIX_FMT_YUV420P;
986
        return 1;
Michael Niedermayer's avatar
Michael Niedermayer committed
987 988 989
    case AV_PIX_FMT_YUVJ411P:
        *format = AV_PIX_FMT_YUV411P;
        return 1;
990 991
    case AV_PIX_FMT_YUVJ422P:
        *format = AV_PIX_FMT_YUV422P;
992
        return 1;
993 994
    case AV_PIX_FMT_YUVJ444P:
        *format = AV_PIX_FMT_YUV444P;
995
        return 1;
996 997
    case AV_PIX_FMT_YUVJ440P:
        *format = AV_PIX_FMT_YUV440P;
998
        return 1;
999
    case AV_PIX_FMT_GRAY8:
1000
    case AV_PIX_FMT_YA8:
1001 1002
    case AV_PIX_FMT_GRAY16LE:
    case AV_PIX_FMT_GRAY16BE:
1003 1004
    case AV_PIX_FMT_YA16BE:
    case AV_PIX_FMT_YA16LE:
1005
        return 1;
1006 1007
    default:
        return 0;
1008 1009 1010
    }
}

1011
static int handle_0alpha(enum AVPixelFormat *format)
1012 1013
{
    switch (*format) {
1014 1015 1016 1017
    case AV_PIX_FMT_0BGR    : *format = AV_PIX_FMT_ABGR   ; return 1;
    case AV_PIX_FMT_BGR0    : *format = AV_PIX_FMT_BGRA   ; return 4;
    case AV_PIX_FMT_0RGB    : *format = AV_PIX_FMT_ARGB   ; return 1;
    case AV_PIX_FMT_RGB0    : *format = AV_PIX_FMT_RGBA   ; return 4;
1018
    default:                                          return 0;
1019 1020 1021
    }
}

1022 1023 1024 1025 1026 1027 1028 1029 1030
static int handle_xyz(enum AVPixelFormat *format)
{
    switch (*format) {
    case AV_PIX_FMT_XYZ12BE : *format = AV_PIX_FMT_RGB48BE; return 1;
    case AV_PIX_FMT_XYZ12LE : *format = AV_PIX_FMT_RGB48LE; return 1;
    default:                                                return 0;
    }
}

1031 1032
static void handle_formats(SwsContext *c)
{
1033 1034 1035 1036
    c->src0Alpha |= handle_0alpha(&c->srcFormat);
    c->dst0Alpha |= handle_0alpha(&c->dstFormat);
    c->srcXYZ    |= handle_xyz(&c->srcFormat);
    c->dstXYZ    |= handle_xyz(&c->dstFormat);
1037 1038
    if (c->srcXYZ || c->dstXYZ)
        fill_xyztables(c);
1039 1040
}

1041 1042
SwsContext *sws_alloc_context(void)
{
1043
    SwsContext *c = av_mallocz(sizeof(SwsContext));
1044

1045 1046
    av_assert0(offsetof(SwsContext, redDither) + DITHER32_INT == offsetof(SwsContext, dither32));

1047
    if (c) {
1048
        c->av_class = &ff_sws_context_class;
1049 1050
        av_opt_set_defaults(c);
    }
1051 1052 1053 1054

    return c;
}

1055 1056 1057 1058 1059 1060 1061 1062 1063 1064 1065 1066 1067 1068
static uint16_t * alloc_gamma_tbl(double e)
{
    int i = 0;
    uint16_t * tbl;
    tbl = (uint16_t*)av_malloc(sizeof(uint16_t) * 1 << 16);
    if (!tbl)
        return NULL;

    for (i = 0; i < 65536; ++i) {
        tbl[i] = pow(i / 65535.0, e) * 65535.0;
    }
    return tbl;
}

1069 1070 1071
static enum AVPixelFormat alphaless_fmt(enum AVPixelFormat fmt)
{
    switch(fmt) {
1072 1073 1074 1075 1076
    case AV_PIX_FMT_ARGB:       return AV_PIX_FMT_RGB24;
    case AV_PIX_FMT_RGBA:       return AV_PIX_FMT_RGB24;
    case AV_PIX_FMT_ABGR:       return AV_PIX_FMT_BGR24;
    case AV_PIX_FMT_BGRA:       return AV_PIX_FMT_BGR24;
    case AV_PIX_FMT_YA8:        return AV_PIX_FMT_GRAY8;
1077 1078 1079

    case AV_PIX_FMT_YUVA420P:   return AV_PIX_FMT_YUV420P;
    case AV_PIX_FMT_YUVA422P:   return AV_PIX_FMT_YUV422P;
1080 1081 1082 1083 1084 1085 1086
    case AV_PIX_FMT_YUVA444P:           return AV_PIX_FMT_YUV444P;

    case AV_PIX_FMT_GBRAP:              return AV_PIX_FMT_GBRP;

    case AV_PIX_FMT_GBRAP16LE:          return AV_PIX_FMT_GBRP16;
    case AV_PIX_FMT_GBRAP16BE:          return AV_PIX_FMT_GBRP16;

1087 1088 1089 1090
    case AV_PIX_FMT_RGBA64LE:   return AV_PIX_FMT_RGB48;
    case AV_PIX_FMT_RGBA64BE:   return AV_PIX_FMT_RGB48;
    case AV_PIX_FMT_BGRA64LE:   return AV_PIX_FMT_BGR48;
    case AV_PIX_FMT_BGRA64BE:   return AV_PIX_FMT_BGR48;
1091

1092 1093
    case AV_PIX_FMT_YA16BE:             return AV_PIX_FMT_GRAY16;
    case AV_PIX_FMT_YA16LE:             return AV_PIX_FMT_GRAY16;
1094

1095 1096
    case AV_PIX_FMT_YUVA420P9BE:        return AV_PIX_FMT_YUV420P9;
    case AV_PIX_FMT_YUVA422P9BE:        return AV_PIX_FMT_YUV422P9;
1097
    case AV_PIX_FMT_YUVA444P9BE:        return AV_PIX_FMT_YUV444P9;
1098 1099
    case AV_PIX_FMT_YUVA420P9LE:        return AV_PIX_FMT_YUV420P9;
    case AV_PIX_FMT_YUVA422P9LE:        return AV_PIX_FMT_YUV422P9;
1100
    case AV_PIX_FMT_YUVA444P9LE:        return AV_PIX_FMT_YUV444P9;
1101 1102
    case AV_PIX_FMT_YUVA420P10BE:       return AV_PIX_FMT_YUV420P10;
    case AV_PIX_FMT_YUVA422P10BE:       return AV_PIX_FMT_YUV422P10;
1103
    case AV_PIX_FMT_YUVA444P10BE:       return AV_PIX_FMT_YUV444P10;
1104 1105
    case AV_PIX_FMT_YUVA420P10LE:       return AV_PIX_FMT_YUV420P10;
    case AV_PIX_FMT_YUVA422P10LE:       return AV_PIX_FMT_YUV422P10;
1106
    case AV_PIX_FMT_YUVA444P10LE:       return AV_PIX_FMT_YUV444P10;
1107 1108
    case AV_PIX_FMT_YUVA420P16BE:       return AV_PIX_FMT_YUV420P16;
    case AV_PIX_FMT_YUVA422P16BE:       return AV_PIX_FMT_YUV422P16;
1109
    case AV_PIX_FMT_YUVA444P16BE:       return AV_PIX_FMT_YUV444P16;
1110 1111
    case AV_PIX_FMT_YUVA420P16LE:       return AV_PIX_FMT_YUV420P16;
    case AV_PIX_FMT_YUVA422P16LE:       return AV_PIX_FMT_YUV422P16;
1112 1113 1114 1115 1116 1117 1118 1119 1120
    case AV_PIX_FMT_YUVA444P16LE:       return AV_PIX_FMT_YUV444P16;

//     case AV_PIX_FMT_AYUV64LE:
//     case AV_PIX_FMT_AYUV64BE:
//     case AV_PIX_FMT_PAL8:
    default: return AV_PIX_FMT_NONE;
    }
}

1121 1122
av_cold int sws_init_context(SwsContext *c, SwsFilter *srcFilter,
                             SwsFilter *dstFilter)
1123
{
1124
    int i, j;
1125 1126
    int usesVFilter, usesHFilter;
    int unscaled;
1127 1128 1129 1130 1131
    SwsFilter dummyFilter = { NULL, NULL, NULL, NULL };
    int srcW              = c->srcW;
    int srcH              = c->srcH;
    int dstW              = c->dstW;
    int dstH              = c->dstH;
1132
    int dst_stride        = FFALIGN(dstW * sizeof(int16_t) + 66, 16);
1133
    int flags, cpu_flags;
1134 1135
    enum AVPixelFormat srcFormat = c->srcFormat;
    enum AVPixelFormat dstFormat = c->dstFormat;
1136 1137
    const AVPixFmtDescriptor *desc_src;
    const AVPixFmtDescriptor *desc_dst;
1138
    int ret = 0;
1139
    enum AVPixelFormat tmpFmt;
1140

1141 1142
    cpu_flags = av_get_cpu_flags();
    flags     = c->flags;
1143
    emms_c();
1144
    if (!rgb15to16)
1145
        ff_sws_rgb2rgb_init();
1146 1147 1148

    unscaled = (srcW == dstW && srcH == dstH);

1149 1150 1151
    c->srcRange |= handle_jpeg(&c->srcFormat);
    c->dstRange |= handle_jpeg(&c->dstFormat);

1152 1153 1154
    if(srcFormat!=c->srcFormat || dstFormat!=c->dstFormat)
        av_log(c, AV_LOG_WARNING, "deprecated pixel format used, make sure you did set range correctly\n");

1155 1156 1157 1158 1159
    if (!c->contrast && !c->saturation && !c->dstFormatBpp)
        sws_setColorspaceDetails(c, ff_yuv2rgb_coeffs[SWS_CS_DEFAULT], c->srcRange,
                                 ff_yuv2rgb_coeffs[SWS_CS_DEFAULT],
                                 c->dstRange, 0, 1 << 16, 1 << 16);

1160 1161 1162 1163 1164
    handle_formats(c);
    srcFormat = c->srcFormat;
    dstFormat = c->dstFormat;
    desc_src = av_pix_fmt_desc_get(srcFormat);
    desc_dst = av_pix_fmt_desc_get(dstFormat);
1165

1166 1167 1168 1169
    // If the source has no alpha then disable alpha blendaway
    if (c->src0Alpha)
        c->alphablend = SWS_ALPHA_BLEND_NONE;

1170 1171
    if (!(unscaled && sws_isSupportedEndiannessConversion(srcFormat) &&
          av_pix_fmt_swap_endianness(srcFormat) == dstFormat)) {
1172
    if (!sws_isSupportedInput(srcFormat)) {
1173
        av_log(c, AV_LOG_ERROR, "%s is not supported as input pixel format\n",
1174
               av_get_pix_fmt_name(srcFormat));
1175
        return AVERROR(EINVAL);
1176
    }
1177
    if (!sws_isSupportedOutput(dstFormat)) {
1178
        av_log(c, AV_LOG_ERROR, "%s is not supported as output pixel format\n",
1179
               av_get_pix_fmt_name(dstFormat));
1180
        return AVERROR(EINVAL);
1181
    }
1182
    }
1183
    av_assert2(desc_src && desc_dst);
1184

1185 1186 1187 1188 1189 1190 1191 1192 1193 1194 1195
    i = flags & (SWS_POINT         |
                 SWS_AREA          |
                 SWS_BILINEAR      |
                 SWS_FAST_BILINEAR |
                 SWS_BICUBIC       |
                 SWS_X             |
                 SWS_GAUSS         |
                 SWS_LANCZOS       |
                 SWS_SINC          |
                 SWS_SPLINE        |
                 SWS_BICUBLIN);
1196 1197 1198 1199

    /* provide a default scaler if not set by caller */
    if (!i) {
        if (dstW < srcW && dstH < srcH)
1200
            flags |= SWS_BICUBIC;
1201
        else if (dstW > srcW && dstH > srcH)
1202
            flags |= SWS_BICUBIC;
1203
        else
1204
            flags |= SWS_BICUBIC;
1205 1206
        c->flags = flags;
    } else if (i & (i - 1)) {
1207
        av_log(c, AV_LOG_ERROR,
1208
               "Exactly one scaler algorithm must be chosen, got %X\n", i);
1209
        return AVERROR(EINVAL);
1210 1211
    }
    /* sanity check */
1212
    if (srcW < 1 || srcH < 1 || dstW < 1 || dstH < 1) {
1213 1214
        /* FIXME check if these are enough and try to lower them after
         * fixing the relevant parts of the code */
1215
        av_log(c, AV_LOG_ERROR, "%dx%d -> %dx%d is invalid scaling dimension\n",
1216
               srcW, srcH, dstW, dstH);
1217
        return AVERROR(EINVAL);
1218
    }
1219 1220 1221 1222 1223 1224
    if (flags & SWS_FAST_BILINEAR) {
        if (srcW < 8 || dstW < 8) {
            flags ^= SWS_FAST_BILINEAR | SWS_BILINEAR;
            c->flags = flags;
        }
    }
1225

1226 1227 1228 1229
    if (!dstFilter)
        dstFilter = &dummyFilter;
    if (!srcFilter)
        srcFilter = &dummyFilter;
1230

1231 1232
    c->lumXInc      = (((int64_t)srcW << 16) + (dstW >> 1)) / dstW;
    c->lumYInc      = (((int64_t)srcH << 16) + (dstH >> 1)) / dstH;
1233 1234
    c->dstFormatBpp = av_get_bits_per_pixel(desc_dst);
    c->srcFormatBpp = av_get_bits_per_pixel(desc_src);
1235
    c->vRounder     = 4 * 0x0001000100010001ULL;
1236

1237 1238 1239 1240 1241 1242 1243 1244
    usesVFilter = (srcFilter->lumV && srcFilter->lumV->length > 1) ||
                  (srcFilter->chrV && srcFilter->chrV->length > 1) ||
                  (dstFilter->lumV && dstFilter->lumV->length > 1) ||
                  (dstFilter->chrV && dstFilter->chrV->length > 1);
    usesHFilter = (srcFilter->lumH && srcFilter->lumH->length > 1) ||
                  (srcFilter->chrH && srcFilter->chrH->length > 1) ||
                  (dstFilter->lumH && dstFilter->lumH->length > 1) ||
                  (dstFilter->chrH && dstFilter->chrH->length > 1);
1245

1246 1247
    av_pix_fmt_get_chroma_sub_sample(srcFormat, &c->chrSrcHSubSample, &c->chrSrcVSubSample);
    av_pix_fmt_get_chroma_sub_sample(dstFormat, &c->chrDstHSubSample, &c->chrDstVSubSample);
1248

1249 1250 1251 1252 1253
    if (isAnyRGB(dstFormat) && !(flags&SWS_FULL_CHR_H_INT)) {
        if (dstW&1) {
            av_log(c, AV_LOG_DEBUG, "Forcing full internal H chroma due to odd output size\n");
            flags |= SWS_FULL_CHR_H_INT;
            c->flags = flags;
1254
        }
1255 1256 1257 1258 1259 1260 1261 1262 1263 1264

        if (   c->chrSrcHSubSample == 0
            && c->chrSrcVSubSample == 0
            && c->dither != SWS_DITHER_BAYER //SWS_FULL_CHR_H_INT is currently not supported with SWS_DITHER_BAYER
            && !(c->flags & SWS_FAST_BILINEAR)
        ) {
            av_log(c, AV_LOG_DEBUG, "Forcing full internal H chroma due to input having non subsampled chroma\n");
            flags |= SWS_FULL_CHR_H_INT;
            c->flags = flags;
        }
1265
    }
1266

1267 1268 1269 1270 1271
    if (c->dither == SWS_DITHER_AUTO) {
        if (flags & SWS_ERROR_DIFFUSION)
            c->dither = SWS_DITHER_ED;
    }

1272 1273 1274 1275
    if(dstFormat == AV_PIX_FMT_BGR4_BYTE ||
       dstFormat == AV_PIX_FMT_RGB4_BYTE ||
       dstFormat == AV_PIX_FMT_BGR8 ||
       dstFormat == AV_PIX_FMT_RGB8) {
1276 1277
        if (c->dither == SWS_DITHER_AUTO)
            c->dither = (flags & SWS_FULL_CHR_H_INT) ? SWS_DITHER_ED : SWS_DITHER_BAYER;
1278
        if (!(flags & SWS_FULL_CHR_H_INT)) {
1279
            if (c->dither == SWS_DITHER_ED || c->dither == SWS_DITHER_A_DITHER || c->dither == SWS_DITHER_X_DITHER) {
1280 1281 1282 1283 1284 1285
                av_log(c, AV_LOG_DEBUG,
                    "Desired dithering only supported in full chroma interpolation for destination format '%s'\n",
                    av_get_pix_fmt_name(dstFormat));
                flags   |= SWS_FULL_CHR_H_INT;
                c->flags = flags;
            }
1286
        }
1287 1288 1289 1290 1291 1292 1293
        if (flags & SWS_FULL_CHR_H_INT) {
            if (c->dither == SWS_DITHER_BAYER) {
                av_log(c, AV_LOG_DEBUG,
                    "Ordered dither is not supported in full chroma interpolation for destination format '%s'\n",
                    av_get_pix_fmt_name(dstFormat));
                c->dither = SWS_DITHER_ED;
            }
1294
        }
1295
    }
1296
    if (isPlanarRGB(dstFormat)) {
1297 1298
        if (!(flags & SWS_FULL_CHR_H_INT)) {
            av_log(c, AV_LOG_DEBUG,
1299 1300
                   "%s output is not supported with half chroma resolution, switching to full\n",
                   av_get_pix_fmt_name(dstFormat));
1301 1302 1303 1304
            flags   |= SWS_FULL_CHR_H_INT;
            c->flags = flags;
        }
    }
1305

1306 1307
    /* reuse chroma for 2 pixels RGB/BGR unless user wants full
     * chroma interpolation */
1308
    if (flags & SWS_FULL_CHR_H_INT &&
1309
        isAnyRGB(dstFormat)        &&
1310
        !isPlanarRGB(dstFormat)    &&
1311 1312 1313 1314 1315 1316 1317 1318
        dstFormat != AV_PIX_FMT_RGBA64LE &&
        dstFormat != AV_PIX_FMT_RGBA64BE &&
        dstFormat != AV_PIX_FMT_BGRA64LE &&
        dstFormat != AV_PIX_FMT_BGRA64BE &&
        dstFormat != AV_PIX_FMT_RGB48LE &&
        dstFormat != AV_PIX_FMT_RGB48BE &&
        dstFormat != AV_PIX_FMT_BGR48LE &&
        dstFormat != AV_PIX_FMT_BGR48BE &&
1319 1320 1321 1322 1323
        dstFormat != AV_PIX_FMT_RGBA  &&
        dstFormat != AV_PIX_FMT_ARGB  &&
        dstFormat != AV_PIX_FMT_BGRA  &&
        dstFormat != AV_PIX_FMT_ABGR  &&
        dstFormat != AV_PIX_FMT_RGB24 &&
1324 1325 1326 1327 1328 1329
        dstFormat != AV_PIX_FMT_BGR24 &&
        dstFormat != AV_PIX_FMT_BGR4_BYTE &&
        dstFormat != AV_PIX_FMT_RGB4_BYTE &&
        dstFormat != AV_PIX_FMT_BGR8 &&
        dstFormat != AV_PIX_FMT_RGB8
    ) {
1330 1331
        av_log(c, AV_LOG_WARNING,
               "full chroma interpolation for destination format '%s' not yet implemented\n",
1332
               av_get_pix_fmt_name(dstFormat));
1333
        flags   &= ~SWS_FULL_CHR_H_INT;
1334
        c->flags = flags;
1335
    }
1336 1337
    if (isAnyRGB(dstFormat) && !(flags & SWS_FULL_CHR_H_INT))
        c->chrDstHSubSample = 1;
1338 1339

    // drop some chroma lines if the user wants it
1340 1341 1342 1343 1344 1345 1346
    c->vChrDrop          = (flags & SWS_SRC_V_CHR_DROP_MASK) >>
                           SWS_SRC_V_CHR_DROP_SHIFT;
    c->chrSrcVSubSample += c->vChrDrop;

    /* drop every other pixel for chroma calculation unless user
     * wants full chroma */
    if (isAnyRGB(srcFormat) && !(flags & SWS_FULL_CHR_H_INP)   &&
1347 1348 1349
        srcFormat != AV_PIX_FMT_RGB8 && srcFormat != AV_PIX_FMT_BGR8 &&
        srcFormat != AV_PIX_FMT_RGB4 && srcFormat != AV_PIX_FMT_BGR4 &&
        srcFormat != AV_PIX_FMT_RGB4_BYTE && srcFormat != AV_PIX_FMT_BGR4_BYTE &&
1350 1351 1352 1353
        srcFormat != AV_PIX_FMT_GBRP9BE   && srcFormat != AV_PIX_FMT_GBRP9LE  &&
        srcFormat != AV_PIX_FMT_GBRP10BE  && srcFormat != AV_PIX_FMT_GBRP10LE &&
        srcFormat != AV_PIX_FMT_GBRP12BE  && srcFormat != AV_PIX_FMT_GBRP12LE &&
        srcFormat != AV_PIX_FMT_GBRP14BE  && srcFormat != AV_PIX_FMT_GBRP14LE &&
1354
        srcFormat != AV_PIX_FMT_GBRP16BE  && srcFormat != AV_PIX_FMT_GBRP16LE &&
1355 1356 1357
        ((dstW >> c->chrDstHSubSample) <= (srcW >> 1) ||
         (flags & SWS_FAST_BILINEAR)))
        c->chrSrcHSubSample = 1;
1358

1359 1360 1361 1362 1363
    // Note the FF_CEIL_RSHIFT is so that we always round toward +inf.
    c->chrSrcW = FF_CEIL_RSHIFT(srcW, c->chrSrcHSubSample);
    c->chrSrcH = FF_CEIL_RSHIFT(srcH, c->chrSrcVSubSample);
    c->chrDstW = FF_CEIL_RSHIFT(dstW, c->chrDstHSubSample);
    c->chrDstH = FF_CEIL_RSHIFT(dstH, c->chrDstVSubSample);
1364

1365
    FF_ALLOCZ_OR_GOTO(c, c->formatConvBuffer, FFALIGN(srcW*2+78, 16) * 2, fail);
1366

1367
    c->srcBpc = desc_src->comp[0].depth;
1368 1369
    if (c->srcBpc < 8)
        c->srcBpc = 8;
1370
    c->dstBpc = desc_dst->comp[0].depth;
1371 1372
    if (c->dstBpc < 8)
        c->dstBpc = 8;
1373
    if (isAnyRGB(srcFormat) || srcFormat == AV_PIX_FMT_PAL8)
1374
        c->srcBpc = 16;
1375
    if (c->dstBpc == 16)
1376
        dst_stride <<= 1;
1377

1378
    if (INLINE_MMXEXT(cpu_flags) && c->srcBpc == 8 && c->dstBpc <= 14) {
1379
        c->canMMXEXTBeUsed = dstW >= srcW && (dstW & 31) == 0 &&
1380
                             c->chrDstW >= c->chrSrcW &&
1381
                             (srcW & 15) == 0;
1382
        if (!c->canMMXEXTBeUsed && dstW >= srcW && c->chrDstW >= c->chrSrcW && (srcW & 15) == 0
1383

1384 1385 1386
            && (flags & SWS_FAST_BILINEAR)) {
            if (flags & SWS_PRINT_INFO)
                av_log(c, AV_LOG_INFO,
1387
                       "output width is not a multiple of 32 -> no MMXEXT scaler\n");
1388
        }
1389
        if (usesHFilter || isNBPS(c->srcFormat) || is16BPS(c->srcFormat) || isAnyRGB(c->srcFormat))
1390
            c->canMMXEXTBeUsed = 0;
1391
    } else
1392
        c->canMMXEXTBeUsed = 0;
1393 1394 1395 1396 1397 1398 1399 1400 1401 1402 1403 1404

    c->chrXInc = (((int64_t)c->chrSrcW << 16) + (c->chrDstW >> 1)) / c->chrDstW;
    c->chrYInc = (((int64_t)c->chrSrcH << 16) + (c->chrDstH >> 1)) / c->chrDstH;

    /* Match pixel 0 of the src to pixel 0 of dst and match pixel n-2 of src
     * to pixel n-2 of dst, but only for the FAST_BILINEAR mode otherwise do
     * correct scaling.
     * n-2 is the last chrominance sample available.
     * This is not perfect, but no one should notice the difference, the more
     * correct variant would be like the vertical one, but that would require
     * some special code for the first and last pixel */
    if (flags & SWS_FAST_BILINEAR) {
1405
        if (c->canMMXEXTBeUsed) {
1406 1407
            c->lumXInc += 20;
            c->chrXInc += 20;
1408
        }
1409
        // we don't use the x86 asm scaler if MMX is available
1410
        else if (INLINE_MMX(cpu_flags) && c->dstBpc <= 14) {
1411 1412
            c->lumXInc = ((int64_t)(srcW       - 2) << 16) / (dstW       - 2) - 20;
            c->chrXInc = ((int64_t)(c->chrSrcW - 2) << 16) / (c->chrDstW - 2) - 20;
1413 1414 1415
        }
    }

1416 1417 1418 1419 1420 1421
    // hardcoded for now
    c->gamma_value = 2.2;
    tmpFmt = AV_PIX_FMT_RGBA64LE;


    if (!unscaled && c->gamma_flag && (srcFormat != tmpFmt || dstFormat != tmpFmt)) {
1422
        SwsContext *c2;
1423 1424 1425 1426 1427 1428 1429 1430 1431 1432 1433 1434 1435 1436 1437 1438
        c->cascaded_context[0] = NULL;

        ret = av_image_alloc(c->cascaded_tmp, c->cascaded_tmpStride,
                            srcW, srcH, tmpFmt, 64);
        if (ret < 0)
            return ret;

        c->cascaded_context[0] = sws_getContext(srcW, srcH, srcFormat,
                                                srcW, srcH, tmpFmt,
                                                flags, NULL, NULL, c->param);
        if (!c->cascaded_context[0]) {
            return -1;
        }

        c->cascaded_context[1] = sws_getContext(srcW, srcH, tmpFmt,
                                                dstW, dstH, tmpFmt,
1439
                                                flags, srcFilter, dstFilter, c->param);
1440 1441 1442 1443

        if (!c->cascaded_context[1])
            return -1;

1444 1445 1446 1447 1448 1449 1450
        c2 = c->cascaded_context[1];
        c2->is_internal_gamma = 1;
        c2->gamma     = alloc_gamma_tbl(    c->gamma_value);
        c2->inv_gamma = alloc_gamma_tbl(1.f/c->gamma_value);
        if (!c2->gamma || !c2->inv_gamma)
            return AVERROR(ENOMEM);

1451 1452 1453 1454 1455 1456 1457 1458 1459
        // is_internal_flag is set after creating the context
        // to properly create the gamma convert FilterDescriptor
        // we have to re-initialize it
        ff_free_filters(c2);
        if (ff_init_filters(c2) < 0) {
            sws_freeContext(c2);
            return -1;
        }

1460 1461 1462 1463 1464 1465 1466 1467 1468 1469 1470 1471 1472 1473 1474 1475
        c->cascaded_context[2] = NULL;
        if (dstFormat != tmpFmt) {
            ret = av_image_alloc(c->cascaded1_tmp, c->cascaded1_tmpStride,
                                dstW, dstH, tmpFmt, 64);
            if (ret < 0)
                return ret;

            c->cascaded_context[2] = sws_getContext(dstW, dstH, tmpFmt,
                                                dstW, dstH, dstFormat,
                                                flags, NULL, NULL, c->param);
            if (!c->cascaded_context[2])
                return -1;
        }
        return 0;
    }

1476 1477 1478 1479 1480 1481 1482 1483 1484 1485 1486 1487 1488 1489 1490 1491 1492 1493 1494 1495 1496 1497 1498 1499 1500
    if (isBayer(srcFormat)) {
        if (!unscaled ||
            (dstFormat != AV_PIX_FMT_RGB24 && dstFormat != AV_PIX_FMT_YUV420P)) {
            enum AVPixelFormat tmpFormat = AV_PIX_FMT_RGB24;

            ret = av_image_alloc(c->cascaded_tmp, c->cascaded_tmpStride,
                                srcW, srcH, tmpFormat, 64);
            if (ret < 0)
                return ret;

            c->cascaded_context[0] = sws_getContext(srcW, srcH, srcFormat,
                                                    srcW, srcH, tmpFormat,
                                                    flags, srcFilter, NULL, c->param);
            if (!c->cascaded_context[0])
                return -1;

            c->cascaded_context[1] = sws_getContext(srcW, srcH, tmpFormat,
                                                    dstW, dstH, dstFormat,
                                                    flags, NULL, dstFilter, c->param);
            if (!c->cascaded_context[1])
                return -1;
            return 0;
        }
    }

1501 1502 1503 1504 1505 1506 1507 1508 1509
    if (CONFIG_SWSCALE_ALPHA && isALPHA(srcFormat) && !isALPHA(dstFormat)) {
        enum AVPixelFormat tmpFormat = alphaless_fmt(srcFormat);

        if (tmpFormat != AV_PIX_FMT_NONE && c->alphablend != SWS_ALPHA_BLEND_NONE)
        if (!unscaled ||
            dstFormat != tmpFormat ||
            usesHFilter || usesVFilter ||
            c->srcRange != c->dstRange
        ) {
1510
            c->cascaded_mainindex = 1;
1511 1512 1513 1514 1515 1516 1517 1518 1519 1520 1521 1522 1523 1524 1525
            ret = av_image_alloc(c->cascaded_tmp, c->cascaded_tmpStride,
                                srcW, srcH, tmpFormat, 64);
            if (ret < 0)
                return ret;

            c->cascaded_context[0] = sws_alloc_set_opts(srcW, srcH, srcFormat,
                                                        srcW, srcH, tmpFormat,
                                                        flags, c->param);
            if (!c->cascaded_context[0])
                return -1;
            c->cascaded_context[0]->alphablend = c->alphablend;
            ret = sws_init_context(c->cascaded_context[0], NULL , NULL);
            if (ret < 0)
                return ret;

1526 1527 1528
            c->cascaded_context[1] = sws_alloc_set_opts(srcW, srcH, tmpFormat,
                                                        dstW, dstH, dstFormat,
                                                        flags, c->param);
1529 1530
            if (!c->cascaded_context[1])
                return -1;
1531 1532 1533 1534 1535 1536 1537

            c->cascaded_context[1]->srcRange = c->srcRange;
            c->cascaded_context[1]->dstRange = c->dstRange;
            ret = sws_init_context(c->cascaded_context[1], srcFilter , dstFilter);
            if (ret < 0)
                return ret;

1538 1539 1540 1541
            return 0;
        }
    }

1542 1543
#define USE_MMAP (HAVE_MMAP && HAVE_MPROTECT && defined MAP_ANONYMOUS)

1544 1545
    /* precalculate horizontal scaler filter coefficients */
    {
1546
#if HAVE_MMXEXT_INLINE
1547
// can't downscale !!!
1548
        if (c->canMMXEXTBeUsed && (flags & SWS_FAST_BILINEAR)) {
1549
            c->lumMmxextFilterCodeSize = ff_init_hscaler_mmxext(dstW, c->lumXInc, NULL,
1550
                                                             NULL, NULL, 8);
1551
            c->chrMmxextFilterCodeSize = ff_init_hscaler_mmxext(c->chrDstW, c->chrXInc,
1552
                                                             NULL, NULL, NULL, 4);
1553

1554
#if USE_MMAP
1555 1556 1557 1558 1559 1560 1561 1562
            c->lumMmxextFilterCode = mmap(NULL, c->lumMmxextFilterCodeSize,
                                          PROT_READ | PROT_WRITE,
                                          MAP_PRIVATE | MAP_ANONYMOUS,
                                          -1, 0);
            c->chrMmxextFilterCode = mmap(NULL, c->chrMmxextFilterCodeSize,
                                          PROT_READ | PROT_WRITE,
                                          MAP_PRIVATE | MAP_ANONYMOUS,
                                          -1, 0);
1563
#elif HAVE_VIRTUALALLOC
1564 1565 1566 1567 1568 1569 1570 1571
            c->lumMmxextFilterCode = VirtualAlloc(NULL,
                                                  c->lumMmxextFilterCodeSize,
                                                  MEM_COMMIT,
                                                  PAGE_EXECUTE_READWRITE);
            c->chrMmxextFilterCode = VirtualAlloc(NULL,
                                                  c->chrMmxextFilterCodeSize,
                                                  MEM_COMMIT,
                                                  PAGE_EXECUTE_READWRITE);
1572
#else
1573 1574
            c->lumMmxextFilterCode = av_malloc(c->lumMmxextFilterCodeSize);
            c->chrMmxextFilterCode = av_malloc(c->chrMmxextFilterCodeSize);
1575 1576
#endif

1577
#ifdef MAP_ANONYMOUS
1578
            if (c->lumMmxextFilterCode == MAP_FAILED || c->chrMmxextFilterCode == MAP_FAILED)
1579
#else
1580
            if (!c->lumMmxextFilterCode || !c->chrMmxextFilterCode)
1581
#endif
1582 1583
            {
                av_log(c, AV_LOG_ERROR, "Failed to allocate MMX2FilterCode\n");
1584
                return AVERROR(ENOMEM);
1585
            }
1586

1587 1588 1589 1590
            FF_ALLOCZ_OR_GOTO(c, c->hLumFilter,    (dstW           / 8 + 8) * sizeof(int16_t), fail);
            FF_ALLOCZ_OR_GOTO(c, c->hChrFilter,    (c->chrDstW     / 4 + 8) * sizeof(int16_t), fail);
            FF_ALLOCZ_OR_GOTO(c, c->hLumFilterPos, (dstW       / 2 / 8 + 8) * sizeof(int32_t), fail);
            FF_ALLOCZ_OR_GOTO(c, c->hChrFilterPos, (c->chrDstW / 2 / 4 + 8) * sizeof(int32_t), fail);
1591

1592
            ff_init_hscaler_mmxext(      dstW, c->lumXInc, c->lumMmxextFilterCode,
1593
                                c->hLumFilter, (uint32_t*)c->hLumFilterPos, 8);
1594
            ff_init_hscaler_mmxext(c->chrDstW, c->chrXInc, c->chrMmxextFilterCode,
1595
                                c->hChrFilter, (uint32_t*)c->hChrFilterPos, 4);
1596

1597
#if USE_MMAP
1598 1599 1600 1601 1602
            if (   mprotect(c->lumMmxextFilterCode, c->lumMmxextFilterCodeSize, PROT_EXEC | PROT_READ) == -1
                || mprotect(c->chrMmxextFilterCode, c->chrMmxextFilterCodeSize, PROT_EXEC | PROT_READ) == -1) {
                av_log(c, AV_LOG_ERROR, "mprotect failed, cannot use fast bilinear scaler\n");
                goto fail;
            }
1603 1604
#endif
        } else
1605
#endif /* HAVE_MMXEXT_INLINE */
1606
        {
1607 1608
            const int filterAlign = X86_MMX(cpu_flags)     ? 4 :
                                    PPC_ALTIVEC(cpu_flags) ? 8 : 1;
1609

1610
            if ((ret = initFilter(&c->hLumFilter, &c->hLumFilterPos,
1611 1612 1613 1614
                           &c->hLumFilterSize, c->lumXInc,
                           srcW, dstW, filterAlign, 1 << 14,
                           (flags & SWS_BICUBLIN) ? (flags | SWS_BICUBIC) : flags,
                           cpu_flags, srcFilter->lumH, dstFilter->lumH,
1615 1616
                           c->param,
                           get_local_pos(c, 0, 0, 0),
1617
                           get_local_pos(c, 0, 0, 0))) < 0)
1618
                goto fail;
1619
            if ((ret = initFilter(&c->hChrFilter, &c->hChrFilterPos,
1620 1621 1622 1623
                           &c->hChrFilterSize, c->chrXInc,
                           c->chrSrcW, c->chrDstW, filterAlign, 1 << 14,
                           (flags & SWS_BICUBLIN) ? (flags | SWS_BILINEAR) : flags,
                           cpu_flags, srcFilter->chrH, dstFilter->chrH,
1624 1625
                           c->param,
                           get_local_pos(c, c->chrSrcHSubSample, c->src_h_chr_pos, 0),
1626
                           get_local_pos(c, c->chrDstHSubSample, c->dst_h_chr_pos, 0))) < 0)
1627 1628 1629 1630 1631 1632
                goto fail;
        }
    } // initialize horizontal stuff

    /* precalculate vertical scaler filter coefficients */
    {
1633 1634
        const int filterAlign = X86_MMX(cpu_flags)     ? 2 :
                                PPC_ALTIVEC(cpu_flags) ? 8 : 1;
1635

1636
        if ((ret = initFilter(&c->vLumFilter, &c->vLumFilterPos, &c->vLumFilterSize,
1637 1638 1639
                       c->lumYInc, srcH, dstH, filterAlign, (1 << 12),
                       (flags & SWS_BICUBLIN) ? (flags | SWS_BICUBIC) : flags,
                       cpu_flags, srcFilter->lumV, dstFilter->lumV,
1640 1641
                       c->param,
                       get_local_pos(c, 0, 0, 1),
1642
                       get_local_pos(c, 0, 0, 1))) < 0)
1643
            goto fail;
1644
        if ((ret = initFilter(&c->vChrFilter, &c->vChrFilterPos, &c->vChrFilterSize,
1645 1646 1647 1648
                       c->chrYInc, c->chrSrcH, c->chrDstH,
                       filterAlign, (1 << 12),
                       (flags & SWS_BICUBLIN) ? (flags | SWS_BILINEAR) : flags,
                       cpu_flags, srcFilter->chrV, dstFilter->chrV,
1649 1650
                       c->param,
                       get_local_pos(c, c->chrSrcVSubSample, c->src_v_chr_pos, 1),
1651
                       get_local_pos(c, c->chrDstVSubSample, c->dst_v_chr_pos, 1))) < 0)
1652

1653 1654
            goto fail;

1655
#if HAVE_ALTIVEC
1656 1657
        FF_ALLOC_OR_GOTO(c, c->vYCoeffsBank, sizeof(vector signed short) * c->vLumFilterSize * c->dstH,    fail);
        FF_ALLOC_OR_GOTO(c, c->vCCoeffsBank, sizeof(vector signed short) * c->vChrFilterSize * c->chrDstH, fail);
1658

1659
        for (i = 0; i < c->vLumFilterSize * c->dstH; i++) {
1660 1661
            int j;
            short *p = (short *)&c->vYCoeffsBank[i];
1662
            for (j = 0; j < 8; j++)
1663 1664 1665
                p[j] = c->vLumFilter[i];
        }

1666
        for (i = 0; i < c->vChrFilterSize * c->chrDstH; i++) {
1667 1668
            int j;
            short *p = (short *)&c->vCCoeffsBank[i];
1669
            for (j = 0; j < 8; j++)
1670 1671 1672 1673 1674 1675
                p[j] = c->vChrFilter[i];
        }
#endif
    }

    // calculate buffer sizes so that they won't run out while handling these damn slices
1676 1677 1678 1679 1680 1681 1682 1683 1684 1685 1686 1687 1688 1689 1690 1691
    c->vLumBufSize = c->vLumFilterSize;
    c->vChrBufSize = c->vChrFilterSize;
    for (i = 0; i < dstH; i++) {
        int chrI      = (int64_t)i * c->chrDstH / dstH;
        int nextSlice = FFMAX(c->vLumFilterPos[i] + c->vLumFilterSize - 1,
                              ((c->vChrFilterPos[chrI] + c->vChrFilterSize - 1)
                               << c->chrSrcVSubSample));

        nextSlice >>= c->chrSrcVSubSample;
        nextSlice <<= c->chrSrcVSubSample;
        if (c->vLumFilterPos[i] + c->vLumBufSize < nextSlice)
            c->vLumBufSize = nextSlice - c->vLumFilterPos[i];
        if (c->vChrFilterPos[chrI] + c->vChrBufSize <
            (nextSlice >> c->chrSrcVSubSample))
            c->vChrBufSize = (nextSlice >> c->chrSrcVSubSample) -
                             c->vChrFilterPos[chrI];
1692 1693
    }

1694 1695 1696
    for (i = 0; i < 4; i++)
        FF_ALLOCZ_OR_GOTO(c, c->dither_error[i], (c->dstW+2) * sizeof(int), fail);

1697 1698
    /* Allocate pixbufs (we use dynamic allocation because otherwise we would
     * need to allocate several megabytes to handle all possible cases) */
1699 1700 1701
    FF_ALLOCZ_OR_GOTO(c, c->lumPixBuf,  c->vLumBufSize * 3 * sizeof(int16_t *), fail);
    FF_ALLOCZ_OR_GOTO(c, c->chrUPixBuf, c->vChrBufSize * 3 * sizeof(int16_t *), fail);
    FF_ALLOCZ_OR_GOTO(c, c->chrVPixBuf, c->vChrBufSize * 3 * sizeof(int16_t *), fail);
1702
    if (CONFIG_SWSCALE_ALPHA && isALPHA(c->srcFormat) && isALPHA(c->dstFormat))
1703 1704 1705
        FF_ALLOCZ_OR_GOTO(c, c->alpPixBuf, c->vLumBufSize * 3 * sizeof(int16_t *), fail);
    /* Note we need at least one pixel more at the end because of the MMX code
     * (just in case someone wants to replace the 4000/8000). */
1706
    /* align at 16 bytes for AltiVec */
1707 1708 1709 1710
    for (i = 0; i < c->vLumBufSize; i++) {
        FF_ALLOCZ_OR_GOTO(c, c->lumPixBuf[i + c->vLumBufSize],
                          dst_stride + 16, fail);
        c->lumPixBuf[i] = c->lumPixBuf[i + c->vLumBufSize];
1711
    }
1712
    // 64 / c->scalingBpp is the same as 16 / sizeof(scaling_intermediate)
1713
    c->uv_off   = (dst_stride>>1) + 64 / (c->dstBpc &~ 7);
1714
    c->uv_offx2 = dst_stride + 16;
1715 1716 1717 1718 1719 1720
    for (i = 0; i < c->vChrBufSize; i++) {
        FF_ALLOC_OR_GOTO(c, c->chrUPixBuf[i + c->vChrBufSize],
                         dst_stride * 2 + 32, fail);
        c->chrUPixBuf[i] = c->chrUPixBuf[i + c->vChrBufSize];
        c->chrVPixBuf[i] = c->chrVPixBuf[i + c->vChrBufSize]
                         = c->chrUPixBuf[i] + (dst_stride >> 1) + 8;
1721 1722
    }
    if (CONFIG_SWSCALE_ALPHA && c->alpPixBuf)
1723 1724 1725 1726
        for (i = 0; i < c->vLumBufSize; i++) {
            FF_ALLOCZ_OR_GOTO(c, c->alpPixBuf[i + c->vLumBufSize],
                              dst_stride + 16, fail);
            c->alpPixBuf[i] = c->alpPixBuf[i + c->vLumBufSize];
1727 1728
        }

1729 1730
    // try to avoid drawing green stuff between the right end and the stride end
    for (i = 0; i < c->vChrBufSize; i++)
1731
        if(desc_dst->comp[0].depth == 16){
1732
            av_assert0(c->dstBpc > 14);
1733 1734 1735 1736 1737
            for(j=0; j<dst_stride/2+1; j++)
                ((int32_t*)(c->chrUPixBuf[i]))[j] = 1<<18;
        } else
            for(j=0; j<dst_stride+1; j++)
                ((int16_t*)(c->chrUPixBuf[i]))[j] = 1<<14;
1738

1739
    av_assert0(c->chrDstH <= dstH);
1740

1741
    if (flags & SWS_PRINT_INFO) {
1742
        const char *scaler = NULL, *cpucaps;
1743

1744 1745 1746 1747 1748 1749 1750 1751 1752
        for (i = 0; i < FF_ARRAY_ELEMS(scale_algorithms); i++) {
            if (flags & scale_algorithms[i].flag) {
                scaler = scale_algorithms[i].description;
                break;
            }
        }
        if (!scaler)
            scaler =  "ehh flags invalid?!";
        av_log(c, AV_LOG_INFO, "%s scaler, from %s to %s%s ",
1753
               scaler,
1754
               av_get_pix_fmt_name(srcFormat),
1755
#ifdef DITHER1XBPP
1756 1757 1758
               dstFormat == AV_PIX_FMT_BGR555   || dstFormat == AV_PIX_FMT_BGR565   ||
               dstFormat == AV_PIX_FMT_RGB444BE || dstFormat == AV_PIX_FMT_RGB444LE ||
               dstFormat == AV_PIX_FMT_BGR444BE || dstFormat == AV_PIX_FMT_BGR444LE ?
1759
                                                             "dithered " : "",
1760 1761 1762
#else
               "",
#endif
1763
               av_get_pix_fmt_name(dstFormat));
1764

1765
        if (INLINE_MMXEXT(cpu_flags))
1766
            cpucaps = "MMXEXT";
1767
        else if (INLINE_AMD3DNOW(cpu_flags))
1768
            cpucaps = "3DNOW";
1769
        else if (INLINE_MMX(cpu_flags))
1770
            cpucaps = "MMX";
1771
        else if (PPC_ALTIVEC(cpu_flags))
1772
            cpucaps = "AltiVec";
1773
        else
1774 1775 1776
            cpucaps = "C";

        av_log(c, AV_LOG_INFO, "using %s\n", cpucaps);
1777 1778

        av_log(c, AV_LOG_VERBOSE, "%dx%d -> %dx%d\n", srcW, srcH, dstW, dstH);
1779 1780
        av_log(c, AV_LOG_DEBUG,
               "lum srcW=%d srcH=%d dstW=%d dstH=%d xInc=%d yInc=%d\n",
1781
               c->srcW, c->srcH, c->dstW, c->dstH, c->lumXInc, c->lumYInc);
1782 1783 1784 1785
        av_log(c, AV_LOG_DEBUG,
               "chr srcW=%d srcH=%d dstW=%d dstH=%d xInc=%d yInc=%d\n",
               c->chrSrcW, c->chrSrcH, c->chrDstW, c->chrDstH,
               c->chrXInc, c->chrYInc);
1786 1787
    }

1788 1789 1790 1791 1792 1793 1794 1795 1796 1797 1798 1799 1800 1801 1802 1803
    /* alpha blend special case, note this has been split via cascaded contexts if its scaled */
    if (unscaled && !usesHFilter && !usesVFilter &&
        c->alphablend != SWS_ALPHA_BLEND_NONE &&
        isALPHA(srcFormat) &&
        (c->srcRange == c->dstRange || isAnyRGB(dstFormat)) &&
        alphaless_fmt(srcFormat) == dstFormat
    ) {
        c->swscale = ff_sws_alphablendaway;

        if (flags & SWS_PRINT_INFO)
            av_log(c, AV_LOG_INFO,
                    "using alpha blendaway %s -> %s special converter\n",
                    av_get_pix_fmt_name(srcFormat), av_get_pix_fmt_name(dstFormat));
        return 0;
    }

1804 1805 1806 1807 1808 1809 1810 1811 1812 1813 1814 1815 1816 1817
    /* unscaled special cases */
    if (unscaled && !usesHFilter && !usesVFilter &&
        (c->srcRange == c->dstRange || isAnyRGB(dstFormat))) {
        ff_get_unscaled_swscale(c);

        if (c->swscale) {
            if (flags & SWS_PRINT_INFO)
                av_log(c, AV_LOG_INFO,
                       "using unscaled %s -> %s special converter\n",
                       av_get_pix_fmt_name(srcFormat), av_get_pix_fmt_name(dstFormat));
            return 0;
        }
    }

1818
    c->swscale = ff_getSwsFunc(c);
1819
    return ff_init_filters(c);
1820
fail: // FIXME replace things by appropriate error codes
1821 1822 1823 1824 1825
    if (ret == RETCODE_USE_CASCADE)  {
        int tmpW = sqrt(srcW * (int64_t)dstW);
        int tmpH = sqrt(srcH * (int64_t)dstH);
        enum AVPixelFormat tmpFormat = AV_PIX_FMT_YUV420P;

1826 1827 1828
        if (isALPHA(srcFormat))
            tmpFormat = AV_PIX_FMT_YUVA420P;

1829 1830 1831 1832 1833 1834 1835 1836 1837 1838 1839 1840 1841 1842 1843 1844 1845 1846 1847 1848 1849
        if (srcW*(int64_t)srcH <= 4LL*dstW*dstH)
            return AVERROR(EINVAL);

        ret = av_image_alloc(c->cascaded_tmp, c->cascaded_tmpStride,
                             tmpW, tmpH, tmpFormat, 64);
        if (ret < 0)
            return ret;

        c->cascaded_context[0] = sws_getContext(srcW, srcH, srcFormat,
                                                tmpW, tmpH, tmpFormat,
                                                flags, srcFilter, NULL, c->param);
        if (!c->cascaded_context[0])
            return -1;

        c->cascaded_context[1] = sws_getContext(tmpW, tmpH, tmpFormat,
                                                dstW, dstH, dstFormat,
                                                flags, NULL, dstFilter, c->param);
        if (!c->cascaded_context[1])
            return -1;
        return 0;
    }
1850 1851
    return -1;
}
1852

1853 1854 1855
SwsContext *sws_alloc_set_opts(int srcW, int srcH, enum AVPixelFormat srcFormat,
                               int dstW, int dstH, enum AVPixelFormat dstFormat,
                               int flags, const double *param)
1856 1857 1858
{
    SwsContext *c;

1859
    if (!(c = sws_alloc_context()))
1860 1861
        return NULL;

1862 1863 1864 1865 1866 1867 1868
    c->flags     = flags;
    c->srcW      = srcW;
    c->srcH      = srcH;
    c->dstW      = dstW;
    c->dstH      = dstH;
    c->srcFormat = srcFormat;
    c->dstFormat = dstFormat;
1869 1870 1871 1872 1873 1874

    if (param) {
        c->param[0] = param[0];
        c->param[1] = param[1];
    }

1875 1876 1877 1878 1879 1880 1881 1882 1883 1884 1885 1886 1887 1888 1889 1890
    return c;
}

SwsContext *sws_getContext(int srcW, int srcH, enum AVPixelFormat srcFormat,
                           int dstW, int dstH, enum AVPixelFormat dstFormat,
                           int flags, SwsFilter *srcFilter,
                           SwsFilter *dstFilter, const double *param)
{
    SwsContext *c;

    c = sws_alloc_set_opts(srcW, srcH, srcFormat,
                           dstW, dstH, dstFormat,
                           flags, param);
    if (!c)
        return NULL;

1891
    if (sws_init_context(c, srcFilter, dstFilter) < 0) {
1892 1893 1894 1895 1896
        sws_freeContext(c);
        return NULL;
    }

    return c;
1897 1898
}

1899 1900 1901 1902 1903 1904 1905 1906 1907 1908 1909 1910 1911 1912 1913 1914
static int isnan_vec(SwsVector *a)
{
    int i;
    for (i=0; i<a->length; i++)
        if (isnan(a->coeff[i]))
            return 1;
    return 0;
}

static void makenan_vec(SwsVector *a)
{
    int i;
    for (i=0; i<a->length; i++)
        a->coeff[i] = NAN;
}

1915 1916 1917 1918 1919
SwsFilter *sws_getDefaultFilter(float lumaGBlur, float chromaGBlur,
                                float lumaSharpen, float chromaSharpen,
                                float chromaHShift, float chromaVShift,
                                int verbose)
{
1920
    SwsFilter *filter = av_malloc(sizeof(SwsFilter));
1921 1922 1923
    if (!filter)
        return NULL;

1924 1925 1926
    if (lumaGBlur != 0.0) {
        filter->lumH = sws_getGaussianVec(lumaGBlur, 3.0);
        filter->lumV = sws_getGaussianVec(lumaGBlur, 3.0);
1927
    } else {
1928 1929
        filter->lumH = sws_getIdentityVec();
        filter->lumV = sws_getIdentityVec();
1930 1931
    }

1932 1933 1934
    if (chromaGBlur != 0.0) {
        filter->chrH = sws_getGaussianVec(chromaGBlur, 3.0);
        filter->chrV = sws_getGaussianVec(chromaGBlur, 3.0);
1935
    } else {
1936 1937
        filter->chrH = sws_getIdentityVec();
        filter->chrV = sws_getIdentityVec();
1938 1939
    }

1940 1941
    if (!filter->lumH || !filter->lumV || !filter->chrH || !filter->chrV)
        goto fail;
1942

1943 1944
    if (chromaSharpen != 0.0) {
        SwsVector *id = sws_getIdentityVec();
1945 1946
        if (!id)
            goto fail;
1947 1948 1949 1950 1951 1952 1953
        sws_scaleVec(filter->chrH, -chromaSharpen);
        sws_scaleVec(filter->chrV, -chromaSharpen);
        sws_addVec(filter->chrH, id);
        sws_addVec(filter->chrV, id);
        sws_freeVec(id);
    }

1954 1955
    if (lumaSharpen != 0.0) {
        SwsVector *id = sws_getIdentityVec();
1956 1957
        if (!id)
            goto fail;
1958 1959 1960 1961 1962 1963 1964 1965
        sws_scaleVec(filter->lumH, -lumaSharpen);
        sws_scaleVec(filter->lumV, -lumaSharpen);
        sws_addVec(filter->lumH, id);
        sws_addVec(filter->lumV, id);
        sws_freeVec(id);
    }

    if (chromaHShift != 0.0)
1966
        sws_shiftVec(filter->chrH, (int)(chromaHShift + 0.5));
1967 1968

    if (chromaVShift != 0.0)
1969
        sws_shiftVec(filter->chrV, (int)(chromaVShift + 0.5));
1970 1971 1972 1973 1974 1975

    sws_normalizeVec(filter->chrH, 1.0);
    sws_normalizeVec(filter->chrV, 1.0);
    sws_normalizeVec(filter->lumH, 1.0);
    sws_normalizeVec(filter->lumV, 1.0);

1976 1977 1978 1979 1980 1981
    if (isnan_vec(filter->chrH) ||
        isnan_vec(filter->chrV) ||
        isnan_vec(filter->lumH) ||
        isnan_vec(filter->lumV))
        goto fail;

1982 1983 1984 1985
    if (verbose)
        sws_printVec2(filter->chrH, NULL, AV_LOG_DEBUG);
    if (verbose)
        sws_printVec2(filter->lumH, NULL, AV_LOG_DEBUG);
1986 1987

    return filter;
1988 1989 1990 1991 1992 1993 1994 1995

fail:
    sws_freeVec(filter->lumH);
    sws_freeVec(filter->lumV);
    sws_freeVec(filter->chrH);
    sws_freeVec(filter->chrV);
    av_freep(&filter);
    return NULL;
1996 1997 1998 1999
}

SwsVector *sws_allocVec(int length)
{
2000 2001 2002 2003 2004 2005
    SwsVector *vec;

    if(length <= 0 || length > INT_MAX/ sizeof(double))
        return NULL;

    vec = av_malloc(sizeof(SwsVector));
2006 2007 2008 2009 2010 2011 2012 2013 2014 2015 2016
    if (!vec)
        return NULL;
    vec->length = length;
    vec->coeff  = av_malloc(sizeof(double) * length);
    if (!vec->coeff)
        av_freep(&vec);
    return vec;
}

SwsVector *sws_getGaussianVec(double variance, double quality)
{
2017
    const int length = (int)(variance * quality + 0.5) | 1;
2018
    int i;
2019
    double middle  = (length - 1) * 0.5;
2020 2021 2022 2023 2024 2025
    SwsVector *vec;

    if(variance < 0 || quality < 0)
        return NULL;

    vec = sws_allocVec(length);
2026 2027 2028 2029

    if (!vec)
        return NULL;

2030 2031 2032 2033
    for (i = 0; i < length; i++) {
        double dist = i - middle;
        vec->coeff[i] = exp(-dist * dist / (2 * variance * variance)) /
                        sqrt(2 * variance * M_PI);
2034 2035 2036 2037 2038 2039 2040 2041 2042 2043
    }

    sws_normalizeVec(vec, 1.0);

    return vec;
}

SwsVector *sws_getConstVec(double c, int length)
{
    int i;
2044
    SwsVector *vec = sws_allocVec(length);
2045 2046 2047 2048

    if (!vec)
        return NULL;

2049 2050
    for (i = 0; i < length; i++)
        vec->coeff[i] = c;
2051 2052 2053 2054 2055 2056 2057 2058 2059

    return vec;
}

SwsVector *sws_getIdentityVec(void)
{
    return sws_getConstVec(1.0, 1);
}

2060
static double sws_dcVec(SwsVector *a)
2061 2062
{
    int i;
2063
    double sum = 0;
2064

2065 2066
    for (i = 0; i < a->length; i++)
        sum += a->coeff[i];
2067 2068 2069 2070 2071 2072 2073 2074

    return sum;
}

void sws_scaleVec(SwsVector *a, double scalar)
{
    int i;

2075 2076
    for (i = 0; i < a->length; i++)
        a->coeff[i] *= scalar;
2077 2078 2079 2080
}

void sws_normalizeVec(SwsVector *a, double height)
{
2081
    sws_scaleVec(a, height / sws_dcVec(a));
2082 2083 2084 2085
}

static SwsVector *sws_getConvVec(SwsVector *a, SwsVector *b)
{
2086
    int length = a->length + b->length - 1;
2087
    int i, j;
2088
    SwsVector *vec = sws_getConstVec(0.0, length);
2089 2090 2091 2092

    if (!vec)
        return NULL;

2093 2094 2095
    for (i = 0; i < a->length; i++) {
        for (j = 0; j < b->length; j++) {
            vec->coeff[i + j] += a->coeff[i] * b->coeff[j];
2096 2097 2098 2099 2100 2101 2102 2103
        }
    }

    return vec;
}

static SwsVector *sws_sumVec(SwsVector *a, SwsVector *b)
{
2104
    int length = FFMAX(a->length, b->length);
2105
    int i;
2106
    SwsVector *vec = sws_getConstVec(0.0, length);
2107 2108 2109 2110

    if (!vec)
        return NULL;

2111 2112 2113 2114
    for (i = 0; i < a->length; i++)
        vec->coeff[i + (length - 1) / 2 - (a->length - 1) / 2] += a->coeff[i];
    for (i = 0; i < b->length; i++)
        vec->coeff[i + (length - 1) / 2 - (b->length - 1) / 2] += b->coeff[i];
2115 2116 2117 2118 2119 2120

    return vec;
}

static SwsVector *sws_diffVec(SwsVector *a, SwsVector *b)
{
2121
    int length = FFMAX(a->length, b->length);
2122
    int i;
2123
    SwsVector *vec = sws_getConstVec(0.0, length);
2124 2125 2126 2127

    if (!vec)
        return NULL;

2128 2129 2130 2131
    for (i = 0; i < a->length; i++)
        vec->coeff[i + (length - 1) / 2 - (a->length - 1) / 2] += a->coeff[i];
    for (i = 0; i < b->length; i++)
        vec->coeff[i + (length - 1) / 2 - (b->length - 1) / 2] -= b->coeff[i];
2132 2133 2134 2135 2136 2137 2138

    return vec;
}

/* shift left / or right if "shift" is negative */
static SwsVector *sws_getShiftedVec(SwsVector *a, int shift)
{
2139
    int length = a->length + FFABS(shift) * 2;
2140
    int i;
2141
    SwsVector *vec = sws_getConstVec(0.0, length);
2142 2143 2144 2145

    if (!vec)
        return NULL;

2146 2147 2148
    for (i = 0; i < a->length; i++) {
        vec->coeff[i + (length    - 1) / 2 -
                       (a->length - 1) / 2 - shift] = a->coeff[i];
2149 2150 2151 2152 2153 2154 2155
    }

    return vec;
}

void sws_shiftVec(SwsVector *a, int shift)
{
2156
    SwsVector *shifted = sws_getShiftedVec(a, shift);
2157 2158 2159 2160
    if (!shifted) {
        makenan_vec(a);
        return;
    }
2161
    av_free(a->coeff);
2162 2163
    a->coeff  = shifted->coeff;
    a->length = shifted->length;
2164 2165 2166 2167 2168
    av_free(shifted);
}

void sws_addVec(SwsVector *a, SwsVector *b)
{
2169
    SwsVector *sum = sws_sumVec(a, b);
2170 2171 2172 2173
    if (!sum) {
        makenan_vec(a);
        return;
    }
2174
    av_free(a->coeff);
2175 2176
    a->coeff  = sum->coeff;
    a->length = sum->length;
2177 2178 2179 2180 2181
    av_free(sum);
}

void sws_subVec(SwsVector *a, SwsVector *b)
{
2182
    SwsVector *diff = sws_diffVec(a, b);
2183 2184 2185 2186
    if (!diff) {
        makenan_vec(a);
        return;
    }
2187
    av_free(a->coeff);
2188 2189
    a->coeff  = diff->coeff;
    a->length = diff->length;
2190 2191 2192 2193 2194
    av_free(diff);
}

void sws_convVec(SwsVector *a, SwsVector *b)
{
2195
    SwsVector *conv = sws_getConvVec(a, b);
2196 2197 2198 2199
    if (!conv) {
        makenan_vec(a);
        return;
    }
2200
    av_free(a->coeff);
2201 2202
    a->coeff  = conv->coeff;
    a->length = conv->length;
2203 2204 2205 2206 2207
    av_free(conv);
}

SwsVector *sws_cloneVec(SwsVector *a)
{
2208
    SwsVector *vec = sws_allocVec(a->length);
2209 2210 2211 2212

    if (!vec)
        return NULL;

2213
    memcpy(vec->coeff, a->coeff, a->length * sizeof(*a->coeff));
2214 2215 2216 2217 2218 2219 2220

    return vec;
}

void sws_printVec2(SwsVector *a, AVClass *log_ctx, int log_level)
{
    int i;
2221 2222
    double max = 0;
    double min = 0;
2223 2224
    double range;

2225 2226 2227
    for (i = 0; i < a->length; i++)
        if (a->coeff[i] > max)
            max = a->coeff[i];
2228

2229 2230 2231
    for (i = 0; i < a->length; i++)
        if (a->coeff[i] < min)
            min = a->coeff[i];
2232

2233
    range = max - min;
2234

2235 2236
    for (i = 0; i < a->length; i++) {
        int x = (int)((a->coeff[i] - min) * 60.0 / range + 0.5);
2237
        av_log(log_ctx, log_level, "%1.3f ", a->coeff[i]);
2238 2239
        for (; x > 0; x--)
            av_log(log_ctx, log_level, " ");
2240 2241 2242 2243 2244 2245
        av_log(log_ctx, log_level, "|\n");
    }
}

void sws_freeVec(SwsVector *a)
{
2246 2247
    if (!a)
        return;
2248
    av_freep(&a->coeff);
2249
    a->length = 0;
2250 2251 2252 2253 2254
    av_free(a);
}

void sws_freeFilter(SwsFilter *filter)
{
2255 2256 2257
    if (!filter)
        return;

2258 2259 2260 2261
    sws_freeVec(filter->lumH);
    sws_freeVec(filter->lumV);
    sws_freeVec(filter->chrH);
    sws_freeVec(filter->chrV);
2262 2263 2264 2265 2266 2267
    av_free(filter);
}

void sws_freeContext(SwsContext *c)
{
    int i;
2268 2269
    if (!c)
        return;
2270 2271

    if (c->lumPixBuf) {
2272
        for (i = 0; i < c->vLumBufSize; i++)
2273 2274 2275 2276
            av_freep(&c->lumPixBuf[i]);
        av_freep(&c->lumPixBuf);
    }

2277
    if (c->chrUPixBuf) {
2278
        for (i = 0; i < c->vChrBufSize; i++)
2279 2280 2281
            av_freep(&c->chrUPixBuf[i]);
        av_freep(&c->chrUPixBuf);
        av_freep(&c->chrVPixBuf);
2282 2283 2284
    }

    if (CONFIG_SWSCALE_ALPHA && c->alpPixBuf) {
2285
        for (i = 0; i < c->vLumBufSize; i++)
2286 2287 2288 2289
            av_freep(&c->alpPixBuf[i]);
        av_freep(&c->alpPixBuf);
    }

2290 2291 2292
    for (i = 0; i < 4; i++)
        av_freep(&c->dither_error[i]);

2293 2294 2295 2296
    av_freep(&c->vLumFilter);
    av_freep(&c->vChrFilter);
    av_freep(&c->hLumFilter);
    av_freep(&c->hChrFilter);
2297
#if HAVE_ALTIVEC
2298 2299 2300 2301 2302 2303 2304 2305 2306
    av_freep(&c->vYCoeffsBank);
    av_freep(&c->vCCoeffsBank);
#endif

    av_freep(&c->vLumFilterPos);
    av_freep(&c->vChrFilterPos);
    av_freep(&c->hLumFilterPos);
    av_freep(&c->hChrFilterPos);

2307
#if HAVE_MMX_INLINE
2308
#if USE_MMAP
2309 2310 2311 2312
    if (c->lumMmxextFilterCode)
        munmap(c->lumMmxextFilterCode, c->lumMmxextFilterCodeSize);
    if (c->chrMmxextFilterCode)
        munmap(c->chrMmxextFilterCode, c->chrMmxextFilterCodeSize);
2313
#elif HAVE_VIRTUALALLOC
2314 2315 2316 2317
    if (c->lumMmxextFilterCode)
        VirtualFree(c->lumMmxextFilterCode, 0, MEM_RELEASE);
    if (c->chrMmxextFilterCode)
        VirtualFree(c->chrMmxextFilterCode, 0, MEM_RELEASE);
2318
#else
2319 2320
    av_free(c->lumMmxextFilterCode);
    av_free(c->chrMmxextFilterCode);
2321
#endif
2322 2323
    c->lumMmxextFilterCode = NULL;
    c->chrMmxextFilterCode = NULL;
2324
#endif /* HAVE_MMX_INLINE */
2325 2326

    av_freep(&c->yuvTable);
2327
    av_freep(&c->formatConvBuffer);
2328

2329 2330
    sws_freeContext(c->cascaded_context[0]);
    sws_freeContext(c->cascaded_context[1]);
2331
    sws_freeContext(c->cascaded_context[2]);
2332 2333
    memset(c->cascaded_context, 0, sizeof(c->cascaded_context));
    av_freep(&c->cascaded_tmp[0]);
2334 2335 2336 2337 2338
    av_freep(&c->cascaded1_tmp[0]);

    av_freep(&c->gamma);
    av_freep(&c->inv_gamma);

2339
    ff_free_filters(c);
2340

2341 2342 2343
    av_free(c);
}

2344
struct SwsContext *sws_getCachedContext(struct SwsContext *context, int srcW,
2345
                                        int srcH, enum AVPixelFormat srcFormat,
2346
                                        int dstW, int dstH,
2347
                                        enum AVPixelFormat dstFormat, int flags,
2348 2349 2350
                                        SwsFilter *srcFilter,
                                        SwsFilter *dstFilter,
                                        const double *param)
2351
{
2352 2353
    static const double default_param[2] = { SWS_PARAM_DEFAULT,
                                             SWS_PARAM_DEFAULT };
2354 2355
    int64_t src_h_chr_pos = -513, dst_h_chr_pos = -513,
            src_v_chr_pos = -513, dst_v_chr_pos = -513;
2356 2357 2358 2359

    if (!param)
        param = default_param;

2360
    if (context &&
2361 2362 2363 2364 2365 2366 2367 2368
        (context->srcW      != srcW      ||
         context->srcH      != srcH      ||
         context->srcFormat != srcFormat ||
         context->dstW      != dstW      ||
         context->dstH      != dstH      ||
         context->dstFormat != dstFormat ||
         context->flags     != flags     ||
         context->param[0]  != param[0]  ||
2369
         context->param[1]  != param[1])) {
2370 2371 2372 2373 2374

        av_opt_get_int(context, "src_h_chr_pos", 0, &src_h_chr_pos);
        av_opt_get_int(context, "src_v_chr_pos", 0, &src_v_chr_pos);
        av_opt_get_int(context, "dst_h_chr_pos", 0, &dst_h_chr_pos);
        av_opt_get_int(context, "dst_v_chr_pos", 0, &dst_v_chr_pos);
2375 2376 2377
        sws_freeContext(context);
        context = NULL;
    }
2378

2379
    if (!context) {
2380 2381 2382 2383 2384
        if (!(context = sws_alloc_context()))
            return NULL;
        context->srcW      = srcW;
        context->srcH      = srcH;
        context->srcFormat = srcFormat;
2385 2386
        context->dstW      = dstW;
        context->dstH      = dstH;
2387 2388 2389 2390
        context->dstFormat = dstFormat;
        context->flags     = flags;
        context->param[0]  = param[0];
        context->param[1]  = param[1];
2391 2392 2393 2394 2395 2396

        av_opt_set_int(context, "src_h_chr_pos", src_h_chr_pos, 0);
        av_opt_set_int(context, "src_v_chr_pos", src_v_chr_pos, 0);
        av_opt_set_int(context, "dst_h_chr_pos", dst_h_chr_pos, 0);
        av_opt_set_int(context, "dst_v_chr_pos", dst_v_chr_pos, 0);

2397 2398 2399 2400
        if (sws_init_context(context, srcFilter, dstFilter) < 0) {
            sws_freeContext(context);
            return NULL;
        }
2401 2402 2403
    }
    return context;
}