swscale-test.c 13.4 KB
Newer Older
1
/*
2 3
 * Copyright (C) 2003 Michael Niedermayer <michaelni@gmx.at>
 *
4
 * This file is part of Libav.
5
 *
6
 * Libav is free software; you can redistribute it and/or
7 8 9
 * 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
 * Libav is distributed in the hope that it will be useful,
12
 * 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
 * You should have received a copy of the GNU Lesser General Public
17
 * License along with Libav; if not, write to the Free Software
18
 * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
19
 */
20 21 22 23 24 25 26

#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#include <inttypes.h>
#include <stdarg.h>

27
#undef HAVE_AV_CONFIG_H
28
#include "libavutil/imgutils.h"
29
#include "libavutil/mem.h"
30
#include "libavutil/avutil.h"
31
#include "libavutil/crc.h"
32
#include "libavutil/pixdesc.h"
33
#include "libavutil/lfg.h"
34 35
#include "swscale.h"

36 37
/* HACK Duplicated from swscale_internal.h.
 * Should be removed when a cleaner pixel format system exists. */
38
#define isGray(x)                      \
39
    ((x) == AV_PIX_FMT_GRAY8       ||     \
40
     (x) == AV_PIX_FMT_YA8         ||     \
41
     (x) == AV_PIX_FMT_GRAY16BE    ||     \
42 43 44
     (x) == AV_PIX_FMT_GRAY16LE    ||     \
     (x) == AV_PIX_FMT_YA16BE      ||     \
     (x) == AV_PIX_FMT_YA16LE)
45 46
#define hasChroma(x)                   \
    (!(isGray(x)                ||     \
47 48
       (x) == AV_PIX_FMT_MONOBLACK ||     \
       (x) == AV_PIX_FMT_MONOWHITE))
49
#define isALPHA(x)                     \
50 51 52 53 54
    ((x) == AV_PIX_FMT_BGR32   ||         \
     (x) == AV_PIX_FMT_BGR32_1 ||         \
     (x) == AV_PIX_FMT_RGB32   ||         \
     (x) == AV_PIX_FMT_RGB32_1 ||         \
     (x) == AV_PIX_FMT_YUVA420P)
55 56 57

static uint64_t getSSD(uint8_t *src1, uint8_t *src2, int stride1,
                       int stride2, int w, int h)
Ramiro Polla's avatar
Ramiro Polla committed
58
{
59 60
    int x, y;
    uint64_t ssd = 0;
61

62 63 64 65
    for (y = 0; y < h; y++) {
        for (x = 0; x < w; x++) {
            int d = src1[x + y * stride1] - src2[x + y * stride2];
            ssd += d * d;
66 67 68
        }
    }
    return ssd;
69 70
}

71 72 73 74 75 76 77 78
struct Results {
    uint64_t ssdY;
    uint64_t ssdU;
    uint64_t ssdV;
    uint64_t ssdA;
    uint32_t crc;
};

79 80
// test by ref -> src -> dst -> out & compare out against ref
// ref & out are YV12
Stefano Sabatini's avatar
Stefano Sabatini committed
81
static int doTest(uint8_t *ref[4], int refStride[4], int w, int h,
82
                  enum AVPixelFormat srcFormat, enum AVPixelFormat dstFormat,
83 84
                  int srcW, int srcH, int dstW, int dstH, int flags,
                  struct Results *r)
Ramiro Polla's avatar
Ramiro Polla committed
85
{
86 87 88
    const AVPixFmtDescriptor *desc_yuva420p = av_pix_fmt_desc_get(AV_PIX_FMT_YUVA420P);
    const AVPixFmtDescriptor *desc_src      = av_pix_fmt_desc_get(srcFormat);
    const AVPixFmtDescriptor *desc_dst      = av_pix_fmt_desc_get(dstFormat);
89
    static enum AVPixelFormat cur_srcFormat;
90
    static int cur_srcW, cur_srcH;
91 92
    static uint8_t *src[4];
    static int srcStride[4];
93 94
    uint8_t *dst[4] = { 0 };
    uint8_t *out[4] = { 0 };
95
    int dstStride[4];
96
    int i;
97
    uint64_t ssdY, ssdU = 0, ssdV = 0, ssdA = 0;
98
    struct SwsContext *dstContext = NULL, *outContext = NULL;
99
    uint32_t crc = 0;
100
    int res      = 0;
101

102
    if (cur_srcFormat != srcFormat || cur_srcW != srcW || cur_srcH != srcH) {
103 104 105 106
        struct SwsContext *srcContext = NULL;
        int p;

        for (p = 0; p < 4; p++)
107
            av_freep(&src[p]);
108 109 110 111

        av_image_fill_linesizes(srcStride, srcFormat, srcW);
        for (p = 0; p < 4; p++) {
            if (srcStride[p])
112
                src[p] = av_mallocz(srcStride[p] * srcH + 16);
113 114 115 116 117 118
            if (srcStride[p] && !src[p]) {
                perror("Malloc");
                res = -1;
                goto end;
            }
        }
119
        srcContext = sws_getContext(w, h, AV_PIX_FMT_YUVA420P, srcW, srcH,
120
                                    srcFormat, SWS_BILINEAR, NULL, NULL, NULL);
121 122
        if (!srcContext) {
            fprintf(stderr, "Failed to get %s ---> %s\n",
123 124
                    desc_yuva420p->name,
                    desc_src->name);
125 126 127 128 129 130 131
            res = -1;
            goto end;
        }
        sws_scale(srcContext, ref, refStride, 0, h, src, srcStride);
        sws_freeContext(srcContext);

        cur_srcFormat = srcFormat;
132 133
        cur_srcW      = srcW;
        cur_srcH      = srcH;
134 135
    }

136
    av_image_fill_linesizes(dstStride, dstFormat, dstW);
137
    for (i = 0; i < 4; i++) {
138 139 140 141
        /* Image buffers passed into libswscale can be allocated any way you
         * prefer, as long as they're aligned enough for the architecture, and
         * they're freed appropriately (such as using av_free for buffers
         * allocated with av_malloc). */
142 143
        /* An extra 16 bytes is being allocated because some scalers may write
         * out of bounds. */
144
        if (dstStride[i])
145
            dst[i] = av_mallocz(dstStride[i] * dstH + 16);
146
        if (dstStride[i] && !dst[i]) {
147 148 149 150 151 152 153
            perror("Malloc");
            res = -1;

            goto end;
        }
    }

154 155
    dstContext = sws_getContext(srcW, srcH, srcFormat, dstW, dstH, dstFormat,
                                flags, NULL, NULL, NULL);
156
    if (!dstContext) {
157
        fprintf(stderr, "Failed to get %s ---> %s\n",
158
                desc_src->name, desc_dst->name);
159 160 161 162
        res = -1;
        goto end;
    }

163
    printf(" %s %dx%d -> %s %3dx%3d flags=%2d",
164 165
           desc_src->name, srcW, srcH,
           desc_dst->name, dstW, dstH,
166 167 168
           flags);
    fflush(stdout);

169
    sws_scale(dstContext, src, srcStride, 0, srcH, dst, dstStride);
Alex Beregszaszi's avatar
Alex Beregszaszi committed
170

171 172 173
    for (i = 0; i < 4 && dstStride[i]; i++)
        crc = av_crc(av_crc_get_table(AV_CRC_32_IEEE), crc, dst[i],
                     dstStride[i] * dstH);
174

175 176 177 178 179 180
    if (r && crc == r->crc) {
        ssdY = r->ssdY;
        ssdU = r->ssdU;
        ssdV = r->ssdV;
        ssdA = r->ssdA;
    } else {
181
        for (i = 0; i < 4; i++) {
182
            if (refStride[i])
183
                out[i] = av_mallocz(refStride[i] * h);
184 185 186 187 188 189
            if (refStride[i] && !out[i]) {
                perror("Malloc");
                res = -1;
                goto end;
            }
        }
190
        outContext = sws_getContext(dstW, dstH, dstFormat, w, h,
191
                                    AV_PIX_FMT_YUVA420P, SWS_BILINEAR,
192
                                    NULL, NULL, NULL);
193 194
        if (!outContext) {
            fprintf(stderr, "Failed to get %s ---> %s\n",
195 196
                    desc_dst->name,
                    desc_yuva420p->name);
197 198 199 200 201
            res = -1;
            goto end;
        }
        sws_scale(outContext, dst, dstStride, 0, dstH, out, refStride);

202
        ssdY = getSSD(ref[0], out[0], refStride[0], refStride[0], w, h);
Ramiro Polla's avatar
Ramiro Polla committed
203 204
        if (hasChroma(srcFormat) && hasChroma(dstFormat)) {
            //FIXME check that output is really gray
205 206 207 208
            ssdU = getSSD(ref[1], out[1], refStride[1], refStride[1],
                          (w + 1) >> 1, (h + 1) >> 1);
            ssdV = getSSD(ref[2], out[2], refStride[2], refStride[2],
                          (w + 1) >> 1, (h + 1) >> 1);
Ramiro Polla's avatar
Ramiro Polla committed
209 210
        }
        if (isALPHA(srcFormat) && isALPHA(dstFormat))
211
            ssdA = getSSD(ref[3], out[3], refStride[3], refStride[3], w, h);
212

213 214 215 216
        ssdY /= w * h;
        ssdU /= w * h / 4;
        ssdV /= w * h / 4;
        ssdA /= w * h;
217

218 219
        sws_freeContext(outContext);

220
        for (i = 0; i < 4; i++)
221 222 223 224
            if (refStride[i])
                av_free(out[i]);
    }

225
    printf(" CRC=%08x SSD=%5"PRId64 ",%5"PRId64 ",%5"PRId64 ",%5"PRId64 "\n",
226
           crc, ssdY, ssdU, ssdV, ssdA);
227

228
end:
229
    sws_freeContext(dstContext);
230

231
    for (i = 0; i < 4; i++)
232
        if (dstStride[i])
Ramiro Polla's avatar
Ramiro Polla committed
233
            av_free(dst[i]);
234

235
    return res;
236 237
}

238
static void selfTest(uint8_t *ref[4], int refStride[4], int w, int h,
239 240
                     enum AVPixelFormat srcFormat_in,
                     enum AVPixelFormat dstFormat_in)
Ramiro Polla's avatar
Ramiro Polla committed
241
{
242 243 244 245 246 247
    const int flags[] = { SWS_FAST_BILINEAR, SWS_BILINEAR, SWS_BICUBIC,
                          SWS_X, SWS_POINT, SWS_AREA, 0 };
    const int srcW   = w;
    const int srcH   = h;
    const int dstW[] = { srcW - srcW / 3, srcW, srcW + srcW / 3, 0 };
    const int dstH[] = { srcH - srcH / 3, srcH, srcH + srcH / 3, 0 };
248
    enum AVPixelFormat srcFormat, dstFormat;
249
    const AVPixFmtDescriptor *desc_src, *desc_dst;
250

251 252
    for (srcFormat = srcFormat_in != AV_PIX_FMT_NONE ? srcFormat_in : 0;
         srcFormat < AV_PIX_FMT_NB; srcFormat++) {
253 254
        if (!sws_isSupportedInput(srcFormat) ||
            !sws_isSupportedOutput(srcFormat))
255 256
            continue;

257 258
        desc_src = av_pix_fmt_desc_get(srcFormat);

259 260
        for (dstFormat = dstFormat_in != AV_PIX_FMT_NONE ? dstFormat_in : 0;
             dstFormat < AV_PIX_FMT_NB; dstFormat++) {
261 262 263
            int i, j, k;
            int res = 0;

264 265
            if (!sws_isSupportedInput(dstFormat) ||
                !sws_isSupportedOutput(dstFormat))
266
                continue;
267

268 269 270
            desc_dst = av_pix_fmt_desc_get(dstFormat);

            printf("%s -> %s\n", desc_src->name, desc_dst->name);
271
            fflush(stdout);
272

273
            for (k = 0; flags[k] && !res; k++)
Ramiro Polla's avatar
Ramiro Polla committed
274 275
                for (i = 0; dstW[i] && !res; i++)
                    for (j = 0; dstH[j] && !res; j++)
276
                        res = doTest(ref, refStride, w, h,
277
                                     srcFormat, dstFormat,
278 279
                                     srcW, srcH, dstW[i], dstH[j], flags[k],
                                     NULL);
280
            if (dstFormat_in != AV_PIX_FMT_NONE)
281
                break;
282
        }
283
        if (srcFormat_in != AV_PIX_FMT_NONE)
284
            break;
285
    }
286 287
}

288
static int fileTest(uint8_t *ref[4], int refStride[4], int w, int h, FILE *fp,
289 290
                    enum AVPixelFormat srcFormat_in,
                    enum AVPixelFormat dstFormat_in)
291 292 293 294 295
{
    char buf[256];

    while (fgets(buf, sizeof(buf), fp)) {
        struct Results r;
296
        enum AVPixelFormat srcFormat;
297 298
        char srcStr[12];
        int srcW, srcH;
299
        enum AVPixelFormat dstFormat;
300 301 302 303 304
        char dstStr[12];
        int dstW, dstH;
        int flags;
        int ret;

305 306 307 308 309
        ret = sscanf(buf,
                     " %12s %dx%d -> %12s %dx%d flags=%d CRC=%x"
                     " SSD=%"PRId64 ", %"PRId64 ", %"PRId64 ", %"PRId64 "\n",
                     srcStr, &srcW, &srcH, dstStr, &dstW, &dstH,
                     &flags, &r.crc, &r.ssdY, &r.ssdU, &r.ssdV, &r.ssdA);
310 311
        if (ret != 12) {
            srcStr[0] = dstStr[0] = 0;
312
            ret       = sscanf(buf, "%12s -> %12s\n", srcStr, dstStr);
313 314 315 316 317
        }

        srcFormat = av_get_pix_fmt(srcStr);
        dstFormat = av_get_pix_fmt(dstStr);

318
        if (srcFormat == AV_PIX_FMT_NONE || dstFormat == AV_PIX_FMT_NONE) {
319 320 321
            fprintf(stderr, "malformed input file\n");
            return -1;
        }
322 323
        if ((srcFormat_in != AV_PIX_FMT_NONE && srcFormat_in != srcFormat) ||
            (dstFormat_in != AV_PIX_FMT_NONE && dstFormat_in != dstFormat))
324
            continue;
325 326 327 328 329 330 331 332 333 334 335 336 337 338
        if (ret != 12) {
            printf("%s", buf);
            continue;
        }

        doTest(ref, refStride, w, h,
               srcFormat, dstFormat,
               srcW, srcH, dstW, dstH, flags,
               &r);
    }

    return 0;
}

339 340 341
#define W 96
#define H 96

Ramiro Polla's avatar
Ramiro Polla committed
342 343
int main(int argc, char **argv)
{
344 345
    enum AVPixelFormat srcFormat = AV_PIX_FMT_NONE;
    enum AVPixelFormat dstFormat = AV_PIX_FMT_NONE;
346
    uint8_t *rgb_data   = av_malloc(W * H * 4);
347 348
    uint8_t *rgb_src[4] = { rgb_data, NULL, NULL, NULL };
    int rgb_stride[4]   = { 4 * W, 0, 0, 0 };
349 350 351
    uint8_t *data       = av_malloc(4 * W * H);
    uint8_t *src[4]     = { data, data + W * H, data + W * H * 2, data + W * H * 3 };
    int stride[4]       = { W, W, W, W };
352 353
    int x, y;
    struct SwsContext *sws;
354
    AVLFG rand;
355 356
    int res = -1;
    int i;
357

358 359 360
    if (!rgb_data || !data)
        return -1;

361 362
    sws = sws_getContext(W / 12, H / 12, AV_PIX_FMT_RGB32, W, H,
                         AV_PIX_FMT_YUVA420P, SWS_BILINEAR, NULL, NULL, NULL);
363

364 365
    av_lfg_init(&rand, 1);

366 367 368
    for (y = 0; y < H; y++)
        for (x = 0; x < W * 4; x++)
            rgb_data[ x + y * 4 * W] = av_lfg_get(&rand);
369
    sws_scale(sws, rgb_src, rgb_stride, 0, H, src, stride);
370
    sws_freeContext(sws);
371
    av_free(rgb_data);
Alex Beregszaszi's avatar
Alex Beregszaszi committed
372

373
    for (i = 1; i < argc; i += 2) {
374
        if (argv[i][0] != '-' || i + 1 == argc)
375 376
            goto bad_option;
        if (!strcmp(argv[i], "-ref")) {
377
            FILE *fp = fopen(argv[i + 1], "r");
378
            if (!fp) {
379
                fprintf(stderr, "could not open '%s'\n", argv[i + 1]);
380 381
                goto error;
            }
382
            res = fileTest(src, stride, W, H, fp, srcFormat, dstFormat);
383 384
            fclose(fp);
            goto end;
385
        } else if (!strcmp(argv[i], "-src")) {
386
            srcFormat = av_get_pix_fmt(argv[i + 1]);
387
            if (srcFormat == AV_PIX_FMT_NONE) {
388
                fprintf(stderr, "invalid pixel format %s\n", argv[i + 1]);
389 390 391
                return -1;
            }
        } else if (!strcmp(argv[i], "-dst")) {
392
            dstFormat = av_get_pix_fmt(argv[i + 1]);
393
            if (dstFormat == AV_PIX_FMT_NONE) {
394
                fprintf(stderr, "invalid pixel format %s\n", argv[i + 1]);
395 396
                return -1;
            }
397 398 399 400 401 402 403
        } else {
bad_option:
            fprintf(stderr, "bad option or argument missing (%s)\n", argv[i]);
            goto error;
        }
    }

404
    selfTest(src, stride, W, H, srcFormat, dstFormat);
405 406 407
end:
    res = 0;
error:
408
    av_free(data);
409

410
    return res;
411
}