yading@11: /*
yading@11:  * Copyright (C) 2003-2011 Michael Niedermayer <michaelni@gmx.at>
yading@11:  *
yading@11:  * This file is part of FFmpeg.
yading@11:  *
yading@11:  * FFmpeg is free software; you can redistribute it and/or
yading@11:  * modify it under the terms of the GNU Lesser General Public
yading@11:  * License as published by the Free Software Foundation; either
yading@11:  * version 2.1 of the License, or (at your option) any later version.
yading@11:  *
yading@11:  * FFmpeg is distributed in the hope that it will be useful,
yading@11:  * but WITHOUT ANY WARRANTY; without even the implied warranty of
yading@11:  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
yading@11:  * Lesser General Public License for more details.
yading@11:  *
yading@11:  * You should have received a copy of the GNU Lesser General Public
yading@11:  * License along with FFmpeg; if not, write to the Free Software
yading@11:  * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
yading@11:  */
yading@11: 
yading@11: #include <stdio.h>
yading@11: #include <stdlib.h>
yading@11: #include <string.h>
yading@11: #include <inttypes.h>
yading@11: #include <stdarg.h>
yading@11: 
yading@11: #undef HAVE_AV_CONFIG_H
yading@11: #include "libavutil/imgutils.h"
yading@11: #include "libavutil/mem.h"
yading@11: #include "libavutil/avutil.h"
yading@11: #include "libavutil/crc.h"
yading@11: #include "libavutil/pixdesc.h"
yading@11: #include "libavutil/lfg.h"
yading@11: #include "swscale.h"
yading@11: 
yading@11: /* HACK Duplicated from swscale_internal.h.
yading@11:  * Should be removed when a cleaner pixel format system exists. */
yading@11: #define isGray(x)                      \
yading@11:     ((x) == AV_PIX_FMT_GRAY8       ||     \
yading@11:      (x) == AV_PIX_FMT_Y400A       ||     \
yading@11:      (x) == AV_PIX_FMT_GRAY16BE    ||     \
yading@11:      (x) == AV_PIX_FMT_GRAY16LE)
yading@11: #define hasChroma(x)                   \
yading@11:     (!(isGray(x)                ||     \
yading@11:        (x) == AV_PIX_FMT_MONOBLACK ||     \
yading@11:        (x) == AV_PIX_FMT_MONOWHITE))
yading@11: #define isALPHA(x)                     \
yading@11:     ((x) == AV_PIX_FMT_BGR32   ||         \
yading@11:      (x) == AV_PIX_FMT_BGR32_1 ||         \
yading@11:      (x) == AV_PIX_FMT_RGB32   ||         \
yading@11:      (x) == AV_PIX_FMT_RGB32_1 ||         \
yading@11:      (x) == AV_PIX_FMT_YUVA420P)
yading@11: 
yading@11: static uint64_t getSSD(const uint8_t *src1, const uint8_t *src2, int stride1,
yading@11:                        int stride2, int w, int h)
yading@11: {
yading@11:     int x, y;
yading@11:     uint64_t ssd = 0;
yading@11: 
yading@11:     for (y = 0; y < h; y++) {
yading@11:         for (x = 0; x < w; x++) {
yading@11:             int d = src1[x + y * stride1] - src2[x + y * stride2];
yading@11:             ssd += d * d;
yading@11:         }
yading@11:     }
yading@11:     return ssd;
yading@11: }
yading@11: 
yading@11: struct Results {
yading@11:     uint64_t ssdY;
yading@11:     uint64_t ssdU;
yading@11:     uint64_t ssdV;
yading@11:     uint64_t ssdA;
yading@11:     uint32_t crc;
yading@11: };
yading@11: 
yading@11: // test by ref -> src -> dst -> out & compare out against ref
yading@11: // ref & out are YV12
yading@11: static int doTest(uint8_t *ref[4], int refStride[4], int w, int h,
yading@11:                   enum AVPixelFormat srcFormat, enum AVPixelFormat dstFormat,
yading@11:                   int srcW, int srcH, int dstW, int dstH, int flags,
yading@11:                   struct Results *r)
yading@11: {
yading@11:     const AVPixFmtDescriptor *desc_yuva420p = av_pix_fmt_desc_get(AV_PIX_FMT_YUVA420P);
yading@11:     const AVPixFmtDescriptor *desc_src      = av_pix_fmt_desc_get(srcFormat);
yading@11:     const AVPixFmtDescriptor *desc_dst      = av_pix_fmt_desc_get(dstFormat);
yading@11:     static enum AVPixelFormat cur_srcFormat;
yading@11:     static int cur_srcW, cur_srcH;
yading@11:     static uint8_t *src[4];
yading@11:     static int srcStride[4];
yading@11:     uint8_t *dst[4] = { 0 };
yading@11:     uint8_t *out[4] = { 0 };
yading@11:     int dstStride[4] = {0};
yading@11:     int i;
yading@11:     uint64_t ssdY, ssdU = 0, ssdV = 0, ssdA = 0;
yading@11:     struct SwsContext *dstContext = NULL, *outContext = NULL;
yading@11:     uint32_t crc = 0;
yading@11:     int res      = 0;
yading@11: 
yading@11:     if (cur_srcFormat != srcFormat || cur_srcW != srcW || cur_srcH != srcH) {
yading@11:         struct SwsContext *srcContext = NULL;
yading@11:         int p;
yading@11: 
yading@11:         for (p = 0; p < 4; p++)
yading@11:             av_freep(&src[p]);
yading@11: 
yading@11:         av_image_fill_linesizes(srcStride, srcFormat, srcW);
yading@11:         for (p = 0; p < 4; p++) {
yading@11:             srcStride[p] = FFALIGN(srcStride[p], 16);
yading@11:             if (srcStride[p])
yading@11:                 src[p] = av_mallocz(srcStride[p] * srcH + 16);
yading@11:             if (srcStride[p] && !src[p]) {
yading@11:                 perror("Malloc");
yading@11:                 res = -1;
yading@11:                 goto end;
yading@11:             }
yading@11:         }
yading@11:         srcContext = sws_getContext(w, h, AV_PIX_FMT_YUVA420P, srcW, srcH,
yading@11:                                     srcFormat, SWS_BILINEAR, NULL, NULL, NULL);
yading@11:         if (!srcContext) {
yading@11:             fprintf(stderr, "Failed to get %s ---> %s\n",
yading@11:                     desc_yuva420p->name,
yading@11:                     desc_src->name);
yading@11:             res = -1;
yading@11:             goto end;
yading@11:         }
yading@11:         sws_scale(srcContext, (const uint8_t * const*)ref, refStride, 0, h, src, srcStride);
yading@11:         sws_freeContext(srcContext);
yading@11: 
yading@11:         cur_srcFormat = srcFormat;
yading@11:         cur_srcW      = srcW;
yading@11:         cur_srcH      = srcH;
yading@11:     }
yading@11: 
yading@11:     av_image_fill_linesizes(dstStride, dstFormat, dstW);
yading@11:     for (i = 0; i < 4; i++) {
yading@11:         /* Image buffers passed into libswscale can be allocated any way you
yading@11:          * prefer, as long as they're aligned enough for the architecture, and
yading@11:          * they're freed appropriately (such as using av_free for buffers
yading@11:          * allocated with av_malloc). */
yading@11:         /* An extra 16 bytes is being allocated because some scalers may write
yading@11:          * out of bounds. */
yading@11:         dstStride[i] = FFALIGN(dstStride[i], 16);
yading@11:         if (dstStride[i])
yading@11:             dst[i] = av_mallocz(dstStride[i] * dstH + 16);
yading@11:         if (dstStride[i] && !dst[i]) {
yading@11:             perror("Malloc");
yading@11:             res = -1;
yading@11: 
yading@11:             goto end;
yading@11:         }
yading@11:     }
yading@11: 
yading@11:     dstContext = sws_getContext(srcW, srcH, srcFormat, dstW, dstH, dstFormat,
yading@11:                                 flags, NULL, NULL, NULL);
yading@11:     if (!dstContext) {
yading@11:         fprintf(stderr, "Failed to get %s ---> %s\n",
yading@11:                 desc_src->name, desc_dst->name);
yading@11:         res = -1;
yading@11:         goto end;
yading@11:     }
yading@11: 
yading@11:     printf(" %s %dx%d -> %s %3dx%3d flags=%2d",
yading@11:            desc_src->name, srcW, srcH,
yading@11:            desc_dst->name, dstW, dstH,
yading@11:            flags);
yading@11:     fflush(stdout);
yading@11: 
yading@11:     sws_scale(dstContext, (const uint8_t * const*)src, srcStride, 0, srcH, dst, dstStride);
yading@11: 
yading@11:     for (i = 0; i < 4 && dstStride[i]; i++)
yading@11:         crc = av_crc(av_crc_get_table(AV_CRC_32_IEEE), crc, dst[i],
yading@11:                      dstStride[i] * dstH);
yading@11: 
yading@11:     if (r && crc == r->crc) {
yading@11:         ssdY = r->ssdY;
yading@11:         ssdU = r->ssdU;
yading@11:         ssdV = r->ssdV;
yading@11:         ssdA = r->ssdA;
yading@11:     } else {
yading@11:         for (i = 0; i < 4; i++) {
yading@11:             refStride[i] = FFALIGN(refStride[i], 16);
yading@11:             if (refStride[i])
yading@11:                 out[i] = av_mallocz(refStride[i] * h);
yading@11:             if (refStride[i] && !out[i]) {
yading@11:                 perror("Malloc");
yading@11:                 res = -1;
yading@11:                 goto end;
yading@11:             }
yading@11:         }
yading@11:         outContext = sws_getContext(dstW, dstH, dstFormat, w, h,
yading@11:                                     AV_PIX_FMT_YUVA420P, SWS_BILINEAR,
yading@11:                                     NULL, NULL, NULL);
yading@11:         if (!outContext) {
yading@11:             fprintf(stderr, "Failed to get %s ---> %s\n",
yading@11:                     desc_dst->name,
yading@11:                     desc_yuva420p->name);
yading@11:             res = -1;
yading@11:             goto end;
yading@11:         }
yading@11:         sws_scale(outContext, (const uint8_t * const*)dst, dstStride, 0, dstH, out, refStride);
yading@11: 
yading@11:         ssdY = getSSD(ref[0], out[0], refStride[0], refStride[0], w, h);
yading@11:         if (hasChroma(srcFormat) && hasChroma(dstFormat)) {
yading@11:             //FIXME check that output is really gray
yading@11:             ssdU = getSSD(ref[1], out[1], refStride[1], refStride[1],
yading@11:                           (w + 1) >> 1, (h + 1) >> 1);
yading@11:             ssdV = getSSD(ref[2], out[2], refStride[2], refStride[2],
yading@11:                           (w + 1) >> 1, (h + 1) >> 1);
yading@11:         }
yading@11:         if (isALPHA(srcFormat) && isALPHA(dstFormat))
yading@11:             ssdA = getSSD(ref[3], out[3], refStride[3], refStride[3], w, h);
yading@11: 
yading@11:         ssdY /= w * h;
yading@11:         ssdU /= w * h / 4;
yading@11:         ssdV /= w * h / 4;
yading@11:         ssdA /= w * h;
yading@11: 
yading@11:         sws_freeContext(outContext);
yading@11: 
yading@11:         for (i = 0; i < 4; i++)
yading@11:             if (refStride[i])
yading@11:                 av_free(out[i]);
yading@11:     }
yading@11: 
yading@11:     printf(" CRC=%08x SSD=%5"PRId64 ",%5"PRId64 ",%5"PRId64 ",%5"PRId64 "\n",
yading@11:            crc, ssdY, ssdU, ssdV, ssdA);
yading@11: 
yading@11: end:
yading@11:     sws_freeContext(dstContext);
yading@11: 
yading@11:     for (i = 0; i < 4; i++)
yading@11:         if (dstStride[i])
yading@11:             av_free(dst[i]);
yading@11: 
yading@11:     return res;
yading@11: }
yading@11: 
yading@11: static void selfTest(uint8_t *ref[4], int refStride[4], int w, int h,
yading@11:                      enum AVPixelFormat srcFormat_in,
yading@11:                      enum AVPixelFormat dstFormat_in)
yading@11: {
yading@11:     const int flags[] = { SWS_FAST_BILINEAR, SWS_BILINEAR, SWS_BICUBIC,
yading@11:                           SWS_X, SWS_POINT, SWS_AREA, 0 };
yading@11:     const int srcW   = w;
yading@11:     const int srcH   = h;
yading@11:     const int dstW[] = { srcW - srcW / 3, srcW, srcW + srcW / 3, 0 };
yading@11:     const int dstH[] = { srcH - srcH / 3, srcH, srcH + srcH / 3, 0 };
yading@11:     enum AVPixelFormat srcFormat, dstFormat;
yading@11:     const AVPixFmtDescriptor *desc_src, *desc_dst;
yading@11: 
yading@11:     for (srcFormat = srcFormat_in != AV_PIX_FMT_NONE ? srcFormat_in : 0;
yading@11:          srcFormat < AV_PIX_FMT_NB; srcFormat++) {
yading@11:         if (!sws_isSupportedInput(srcFormat) ||
yading@11:             !sws_isSupportedOutput(srcFormat))
yading@11:             continue;
yading@11: 
yading@11:         desc_src = av_pix_fmt_desc_get(srcFormat);
yading@11: 
yading@11:         for (dstFormat = dstFormat_in != AV_PIX_FMT_NONE ? dstFormat_in : 0;
yading@11:              dstFormat < AV_PIX_FMT_NB; dstFormat++) {
yading@11:             int i, j, k;
yading@11:             int res = 0;
yading@11: 
yading@11:             if (!sws_isSupportedInput(dstFormat) ||
yading@11:                 !sws_isSupportedOutput(dstFormat))
yading@11:                 continue;
yading@11: 
yading@11:             desc_dst = av_pix_fmt_desc_get(dstFormat);
yading@11: 
yading@11:             printf("%s -> %s\n", desc_src->name, desc_dst->name);
yading@11:             fflush(stdout);
yading@11: 
yading@11:             for (k = 0; flags[k] && !res; k++)
yading@11:                 for (i = 0; dstW[i] && !res; i++)
yading@11:                     for (j = 0; dstH[j] && !res; j++)
yading@11:                         res = doTest(ref, refStride, w, h,
yading@11:                                      srcFormat, dstFormat,
yading@11:                                      srcW, srcH, dstW[i], dstH[j], flags[k],
yading@11:                                      NULL);
yading@11:             if (dstFormat_in != AV_PIX_FMT_NONE)
yading@11:                 break;
yading@11:         }
yading@11:         if (srcFormat_in != AV_PIX_FMT_NONE)
yading@11:             break;
yading@11:     }
yading@11: }
yading@11: 
yading@11: static int fileTest(uint8_t *ref[4], int refStride[4], int w, int h, FILE *fp,
yading@11:                     enum AVPixelFormat srcFormat_in,
yading@11:                     enum AVPixelFormat dstFormat_in)
yading@11: {
yading@11:     char buf[256];
yading@11: 
yading@11:     while (fgets(buf, sizeof(buf), fp)) {
yading@11:         struct Results r;
yading@11:         enum AVPixelFormat srcFormat;
yading@11:         char srcStr[12];
yading@11:         int srcW, srcH;
yading@11:         enum AVPixelFormat dstFormat;
yading@11:         char dstStr[12];
yading@11:         int dstW, dstH;
yading@11:         int flags;
yading@11:         int ret;
yading@11: 
yading@11:         ret = sscanf(buf,
yading@11:                      " %12s %dx%d -> %12s %dx%d flags=%d CRC=%x"
yading@11:                      " SSD=%"SCNd64 ", %"SCNd64 ", %"SCNd64 ", %"SCNd64 "\n",
yading@11:                      srcStr, &srcW, &srcH, dstStr, &dstW, &dstH,
yading@11:                      &flags, &r.crc, &r.ssdY, &r.ssdU, &r.ssdV, &r.ssdA);
yading@11:         if (ret != 12) {
yading@11:             srcStr[0] = dstStr[0] = 0;
yading@11:             ret       = sscanf(buf, "%12s -> %12s\n", srcStr, dstStr);
yading@11:         }
yading@11: 
yading@11:         srcFormat = av_get_pix_fmt(srcStr);
yading@11:         dstFormat = av_get_pix_fmt(dstStr);
yading@11: 
yading@11:         if (srcFormat == AV_PIX_FMT_NONE || dstFormat == AV_PIX_FMT_NONE ||
yading@11:             srcW > 8192U || srcH > 8192U || dstW > 8192U || dstH > 8192U) {
yading@11:             fprintf(stderr, "malformed input file\n");
yading@11:             return -1;
yading@11:         }
yading@11:         if ((srcFormat_in != AV_PIX_FMT_NONE && srcFormat_in != srcFormat) ||
yading@11:             (dstFormat_in != AV_PIX_FMT_NONE && dstFormat_in != dstFormat))
yading@11:             continue;
yading@11:         if (ret != 12) {
yading@11:             printf("%s", buf);
yading@11:             continue;
yading@11:         }
yading@11: 
yading@11:         doTest(ref, refStride, w, h,
yading@11:                srcFormat, dstFormat,
yading@11:                srcW, srcH, dstW, dstH, flags,
yading@11:                &r);
yading@11:     }
yading@11: 
yading@11:     return 0;
yading@11: }
yading@11: 
yading@11: #define W 96
yading@11: #define H 96
yading@11: 
yading@11: int main(int argc, char **argv)
yading@11: {
yading@11:     enum AVPixelFormat srcFormat = AV_PIX_FMT_NONE;
yading@11:     enum AVPixelFormat dstFormat = AV_PIX_FMT_NONE;
yading@11:     uint8_t *rgb_data   = av_malloc(W * H * 4);
yading@11:     const uint8_t * const rgb_src[4] = { rgb_data, NULL, NULL, NULL };
yading@11:     int rgb_stride[4]   = { 4 * W, 0, 0, 0 };
yading@11:     uint8_t *data       = av_malloc(4 * W * H);
yading@11:     uint8_t *src[4]     = { data, data + W * H, data + W * H * 2, data + W * H * 3 };
yading@11:     int stride[4]       = { W, W, W, W };
yading@11:     int x, y;
yading@11:     struct SwsContext *sws;
yading@11:     AVLFG rand;
yading@11:     int res = -1;
yading@11:     int i;
yading@11:     FILE *fp = NULL;
yading@11: 
yading@11:     if (!rgb_data || !data)
yading@11:         return -1;
yading@11: 
yading@11:     for (i = 1; i < argc; i += 2) {
yading@11:         if (argv[i][0] != '-' || i + 1 == argc)
yading@11:             goto bad_option;
yading@11:         if (!strcmp(argv[i], "-ref")) {
yading@11:             fp = fopen(argv[i + 1], "r");
yading@11:             if (!fp) {
yading@11:                 fprintf(stderr, "could not open '%s'\n", argv[i + 1]);
yading@11:                 goto error;
yading@11:             }
yading@11:         } else if (!strcmp(argv[i], "-src")) {
yading@11:             srcFormat = av_get_pix_fmt(argv[i + 1]);
yading@11:             if (srcFormat == AV_PIX_FMT_NONE) {
yading@11:                 fprintf(stderr, "invalid pixel format %s\n", argv[i + 1]);
yading@11:                 return -1;
yading@11:             }
yading@11:         } else if (!strcmp(argv[i], "-dst")) {
yading@11:             dstFormat = av_get_pix_fmt(argv[i + 1]);
yading@11:             if (dstFormat == AV_PIX_FMT_NONE) {
yading@11:                 fprintf(stderr, "invalid pixel format %s\n", argv[i + 1]);
yading@11:                 return -1;
yading@11:             }
yading@11:         } else {
yading@11: bad_option:
yading@11:             fprintf(stderr, "bad option or argument missing (%s)\n", argv[i]);
yading@11:             goto error;
yading@11:         }
yading@11:     }
yading@11: 
yading@11:     sws = sws_getContext(W / 12, H / 12, AV_PIX_FMT_RGB32, W, H,
yading@11:                          AV_PIX_FMT_YUVA420P, SWS_BILINEAR, NULL, NULL, NULL);
yading@11: 
yading@11:     av_lfg_init(&rand, 1);
yading@11: 
yading@11:     for (y = 0; y < H; y++)
yading@11:         for (x = 0; x < W * 4; x++)
yading@11:             rgb_data[ x + y * 4 * W] = av_lfg_get(&rand);
yading@11:     sws_scale(sws, rgb_src, rgb_stride, 0, H, src, stride);
yading@11:     sws_freeContext(sws);
yading@11:     av_free(rgb_data);
yading@11: 
yading@11:     if(fp) {
yading@11:         res = fileTest(src, stride, W, H, fp, srcFormat, dstFormat);
yading@11:         fclose(fp);
yading@11:     } else {
yading@11:         selfTest(src, stride, W, H, srcFormat, dstFormat);
yading@11:         res = 0;
yading@11:     }
yading@11: error:
yading@11:     av_free(data);
yading@11: 
yading@11:     return res;
yading@11: }