libswscale/swscale-example.c

   1 /*
   2  * Copyright (C) 2003 Michael Niedermayer <michaelni@gmx.at>
   3  *
   4  * This file is part of FFmpeg.
   5  *
   6  * FFmpeg is free software; you can redistribute it and/or
   7  * modify it under the terms of the GNU Lesser General Public
   8  * License as published by the Free Software Foundation; either
   9  * version 2.1 of the License, or (at your option) any later version.
  10  *
  11  * FFmpeg is distributed in the hope that it will be useful,
  12  * but WITHOUT ANY WARRANTY; without even the implied warranty of
  13  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
  14  * Lesser General Public License for more details.
  15  *
  16  * You should have received a copy of the GNU Lesser General Public
  17  * License along with FFmpeg; if not, write to the Free Software
  18  * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
  19  */
  20
  21 #include <stdio.h>
  22 #include <stdlib.h>
  23 #include <string.h>
  24 #include <inttypes.h>
  25 #include <stdarg.h>
  26
  27 #undef HAVE_AV_CONFIG_H
  28 #include "libavutil/avutil.h"
  29 #include "swscale.h"
  30 #include "swscale_internal.h"
  31 #include "rgb2rgb.h"
  32
  33 static uint64_t getSSD(uint8_t *src1, uint8_t *src2, int stride1, int stride2, int w, int h){
  34     int x,y;
  35     uint64_t ssd=0;
  36
  37 //printf("%d %d\n", w, h);
  38
  39     for (y=0; y<h; y++){
  40         for (x=0; x<w; x++){
  41             int d= src1[x + y*stride1] - src2[x + y*stride2];
  42             ssd+= d*d;
  43 //printf("%d", abs(src1[x + y*stride1] - src2[x + y*stride2])/26 );
  44         }
  45 //printf("\n");
  46     }
  47     return ssd;
  48 }
  49
  50 // test by ref -> src -> dst -> out & compare out against ref
  51 // ref & out are YV12
  52 static int doTest(uint8_t *ref[3], int refStride[3], int w, int h, int srcFormat, int dstFormat,
  53                   int srcW, int srcH, int dstW, int dstH, int flags){
  54     uint8_t *src[3];
  55     uint8_t *dst[3];
  56     uint8_t *out[3];
  57     int srcStride[3], dstStride[3];
  58     int i;
  59     uint64_t ssdY, ssdU, ssdV;
  60     struct SwsContext *srcContext, *dstContext, *outContext;
  61     int res;
  62
  63     res = 0;
  64     for (i=0; i<3; i++){
  65         // avoid stride % bpp != 0
  66         if (srcFormat==PIX_FMT_RGB24 || srcFormat==PIX_FMT_BGR24)
  67             srcStride[i]= srcW*3;
  68         else
  69             srcStride[i]= srcW*4;
  70
  71         if (dstFormat==PIX_FMT_RGB24 || dstFormat==PIX_FMT_BGR24)
  72             dstStride[i]= dstW*3;
  73         else
  74             dstStride[i]= dstW*4;
  75
  76         src[i]= (uint8_t*) malloc(srcStride[i]*srcH);
  77         dst[i]= (uint8_t*) malloc(dstStride[i]*dstH);
  78         out[i]= (uint8_t*) malloc(refStride[i]*h);
  79         if (!src[i] || !dst[i] || !out[i]) {
  80             perror("Malloc");
  81             res = -1;
  82
  83             goto end;
  84         }
  85     }
  86
  87     dstContext = outContext = NULL;
  88     srcContext= sws_getContext(w, h, PIX_FMT_YUV420P, srcW, srcH, srcFormat, flags, NULL, NULL, NULL);
  89     if (!srcContext) {
  90         fprintf(stderr, "Failed to get %s ---> %s\n",
  91                 sws_format_name(PIX_FMT_YUV420P),
  92                 sws_format_name(srcFormat));
  93         res = -1;
  94
  95         goto end;
  96     }
  97     dstContext= sws_getContext(srcW, srcH, srcFormat, dstW, dstH, dstFormat, flags, NULL, NULL, NULL);
  98     if (!dstContext) {
  99         fprintf(stderr, "Failed to get %s ---> %s\n",
 100                 sws_format_name(srcFormat),
 101                 sws_format_name(dstFormat));
 102         res = -1;
 103
 104         goto end;
 105     }
 106     outContext= sws_getContext(dstW, dstH, dstFormat, w, h, PIX_FMT_YUV420P, flags, NULL, NULL, NULL);
 107     if (!outContext) {
 108         fprintf(stderr, "Failed to get %s ---> %s\n",
 109                 sws_format_name(dstFormat),
 110                 sws_format_name(PIX_FMT_YUV420P));
 111         res = -1;
 112
 113         goto end;
 114     }
 115 //    printf("test %X %X %X -> %X %X %X\n", (int)ref[0], (int)ref[1], (int)ref[2],
 116 //        (int)src[0], (int)src[1], (int)src[2]);
 117
 118     sws_scale(srcContext, ref, refStride, 0, h   , src, srcStride);
 119     sws_scale(dstContext, src, srcStride, 0, srcH, dst, dstStride);
 120     sws_scale(outContext, dst, dstStride, 0, dstH, out, refStride);
 121
 122 #if defined(ARCH_X86)
 123     asm volatile ("emms\n\t");
 124 #endif
 125
 126     ssdY= getSSD(ref[0], out[0], refStride[0], refStride[0], w, h);
 127     ssdU= getSSD(ref[1], out[1], refStride[1], refStride[1], (w+1)>>1, (h+1)>>1);
 128     ssdV= getSSD(ref[2], out[2], refStride[2], refStride[2], (w+1)>>1, (h+1)>>1);
 129
 130     if (srcFormat == PIX_FMT_GRAY8 || dstFormat==PIX_FMT_GRAY8) ssdU=ssdV=0; //FIXME check that output is really gray
 131
 132     ssdY/= w*h;
 133     ssdU/= w*h/4;
 134     ssdV/= w*h/4;
 135
 136         printf(" %s %dx%d -> %s %4dx%4d flags=%2d SSD=%5lld,%5lld,%5lld\n",
 137                sws_format_name(srcFormat), srcW, srcH,
 138                sws_format_name(dstFormat), dstW, dstH,
 139                flags,
 140                ssdY, ssdU, ssdV);
 141         fflush(stdout);
 142
 143     end:
 144
 145     sws_freeContext(srcContext);
 146     sws_freeContext(dstContext);
 147     sws_freeContext(outContext);
 148
 149     for (i=0; i<3; i++){
 150         free(src[i]);
 151         free(dst[i]);
 152         free(out[i]);
 153     }
 154
 155     return res;
 156 }
 157
 158 void fast_memcpy(void *a, void *b, int s){ //FIXME
 159     memcpy(a, b, s);
 160 }
 161
 162 static void selfTest(uint8_t *src[3], int stride[3], int w, int h){
 163     enum PixelFormat srcFormat, dstFormat;
 164     int srcW, srcH, dstW, dstH;
 165     int flags;
 166
 167     for (srcFormat = 0; srcFormat < PIX_FMT_NB; srcFormat++) {
 168         for (dstFormat = 0; dstFormat < PIX_FMT_NB; dstFormat++) {
 169             printf("%s -> %s\n",
 170                    sws_format_name(srcFormat),
 171                    sws_format_name(dstFormat));
 172             fflush(stdout);
 173
 174             srcW= w;
 175             srcH= h;
 176             for (dstW=w - w/3; dstW<= 4*w/3; dstW+= w/3){
 177                 for (dstH=h - h/3; dstH<= 4*h/3; dstH+= h/3){
 178                     for (flags=1; flags<33; flags*=2) {
 179                         int res;
 180
 181                         res = doTest(src, stride, w, h, srcFormat, dstFormat,
 182                                      srcW, srcH, dstW, dstH, flags);
 183                         if (res < 0) {
 184                             dstW = 4 * w / 3;
 185                             dstH = 4 * h / 3;
 186                             flags = 33;
 187                         }
 188                     }
 189                 }
 190             }
 191         }
 192     }
 193 }
 194
 195 #define W 96
 196 #define H 96
 197
 198 int main(int argc, char **argv){
 199     uint8_t *rgb_data = malloc (W*H*4);
 200     uint8_t *rgb_src[3]= {rgb_data, NULL, NULL};
 201     int rgb_stride[3]={4*W, 0, 0};
 202     uint8_t *data = malloc (3*W*H);
 203     uint8_t *src[3]= {data, data+W*H, data+W*H*2};
 204     int stride[3]={W, W, W};
 205     int x, y;
 206     struct SwsContext *sws;
 207
 208     sws= sws_getContext(W/12, H/12, PIX_FMT_RGB32, W, H, PIX_FMT_YUV420P, 2, NULL, NULL, NULL);
 209
 210     for (y=0; y<H; y++){
 211         for (x=0; x<W*4; x++){
 212             rgb_data[ x + y*4*W]= random();
 213         }
 214     }
 215 #if defined(ARCH_X86)
 216     sws_rgb2rgb_init(SWS_CPU_CAPS_MMX*0);
 217 #else
 218     sws_rgb2rgb_init(0);
 219 #endif
 220     sws_scale(sws, rgb_src, rgb_stride, 0, H, src, stride);
 221
 222 #if defined(ARCH_X86)
 223     asm volatile ("emms\n\t");
 224 #endif
 225
 226     selfTest(src, stride, W, H);
 227
 228     return 123;
 229 }