]> git.sesse.net Git - ffmpeg/blob - libswscale/swscale-test.c
swscale-test: scale from reference to source only once for each algorithm
[ffmpeg] / libswscale / swscale-test.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/mem.h"
29 #include "libavutil/avutil.h"
30 #include "libavutil/pixdesc.h"
31 #include "libavutil/lfg.h"
32 #include "swscale.h"
33
34 /* HACK Duplicated from swscale_internal.h.
35  * Should be removed when a cleaner pixel format system exists. */
36 #define isGray(x)       (           \
37            (x)==PIX_FMT_GRAY8       \
38         || (x)==PIX_FMT_GRAY16BE    \
39         || (x)==PIX_FMT_GRAY16LE    \
40     )
41 #define hasChroma(x)   (!(          \
42             isGray(x)               \
43         || (x)==PIX_FMT_MONOBLACK   \
44         || (x)==PIX_FMT_MONOWHITE   \
45     ))
46 #define isALPHA(x)      (           \
47            (x)==PIX_FMT_BGR32       \
48         || (x)==PIX_FMT_BGR32_1     \
49         || (x)==PIX_FMT_RGB32       \
50         || (x)==PIX_FMT_RGB32_1     \
51         || (x)==PIX_FMT_YUVA420P    \
52     )
53
54 static uint64_t getSSD(uint8_t *src1, uint8_t *src2, int stride1, int stride2, int w, int h)
55 {
56     int x,y;
57     uint64_t ssd=0;
58
59 //printf("%d %d\n", w, h);
60
61     for (y=0; y<h; y++) {
62         for (x=0; x<w; x++) {
63             int d= src1[x + y*stride1] - src2[x + y*stride2];
64             ssd+= d*d;
65 //printf("%d", abs(src1[x + y*stride1] - src2[x + y*stride2])/26 );
66         }
67 //printf("\n");
68     }
69     return ssd;
70 }
71
72 // test by ref -> src -> dst -> out & compare out against ref
73 // ref & out are YV12
74 static int doTest(uint8_t *ref[4], int refStride[4], int w, int h,
75                   uint8_t *src[4], int srcStride[4],
76                   enum PixelFormat srcFormat, enum PixelFormat dstFormat,
77                   int srcW, int srcH, int dstW, int dstH, int flags)
78 {
79     uint8_t *dst[4] = {0};
80     uint8_t *out[4] = {0};
81     int dstStride[4];
82     int i;
83     uint64_t ssdY, ssdU=0, ssdV=0, ssdA=0;
84     struct SwsContext *dstContext = NULL, *outContext = NULL;
85     int res;
86
87     res = 0;
88     for (i=0; i<4; i++) {
89         // avoid stride % bpp != 0
90         if (dstFormat==PIX_FMT_RGB24 || dstFormat==PIX_FMT_BGR24)
91             dstStride[i]= dstW*3;
92         else if (dstFormat==PIX_FMT_RGB48BE || dstFormat==PIX_FMT_RGB48LE)
93             dstStride[i]= dstW*6;
94         else
95             dstStride[i]= dstW*4;
96
97         /* Image buffers passed into libswscale can be allocated any way you
98          * prefer, as long as they're aligned enough for the architecture, and
99          * they're freed appropriately (such as using av_free for buffers
100          * allocated with av_malloc). */
101         /* An extra 16 bytes is being allocated because some scalers may write
102          * out of bounds. */
103         dst[i]= av_mallocz(dstStride[i]*dstH+16);
104         out[i]= av_mallocz(refStride[i]*h);
105         if (!dst[i] || !out[i]) {
106             perror("Malloc");
107             res = -1;
108
109             goto end;
110         }
111     }
112
113     dstContext= sws_getContext(srcW, srcH, srcFormat, dstW, dstH, dstFormat, flags, NULL, NULL, NULL);
114     if (!dstContext) {
115         fprintf(stderr, "Failed to get %s ---> %s\n",
116                 av_pix_fmt_descriptors[srcFormat].name,
117                 av_pix_fmt_descriptors[dstFormat].name);
118         res = -1;
119
120         goto end;
121     }
122     outContext= sws_getContext(dstW, dstH, dstFormat, w, h, PIX_FMT_YUVA420P, flags, NULL, NULL, NULL);
123     if (!outContext) {
124         fprintf(stderr, "Failed to get %s ---> %s\n",
125                 av_pix_fmt_descriptors[dstFormat].name,
126                 av_pix_fmt_descriptors[PIX_FMT_YUVA420P].name);
127         res = -1;
128
129         goto end;
130     }
131 //    printf("test %X %X %X -> %X %X %X\n", (int)ref[0], (int)ref[1], (int)ref[2],
132 //        (int)src[0], (int)src[1], (int)src[2]);
133
134     printf(" %s %dx%d -> %s %3dx%3d flags=%2d",
135            av_pix_fmt_descriptors[srcFormat].name, srcW, srcH,
136            av_pix_fmt_descriptors[dstFormat].name, dstW, dstH,
137            flags);
138     fflush(stdout);
139
140     sws_scale(dstContext, src, srcStride, 0, srcH, dst, dstStride);
141     sws_scale(outContext, dst, dstStride, 0, dstH, out, refStride);
142
143     ssdY= getSSD(ref[0], out[0], refStride[0], refStride[0], w, h);
144     if (hasChroma(srcFormat) && hasChroma(dstFormat)) {
145         //FIXME check that output is really gray
146         ssdU= getSSD(ref[1], out[1], refStride[1], refStride[1], (w+1)>>1, (h+1)>>1);
147         ssdV= getSSD(ref[2], out[2], refStride[2], refStride[2], (w+1)>>1, (h+1)>>1);
148     }
149     if (isALPHA(srcFormat) && isALPHA(dstFormat))
150         ssdA= getSSD(ref[3], out[3], refStride[3], refStride[3], w, h);
151
152     ssdY/= w*h;
153     ssdU/= w*h/4;
154     ssdV/= w*h/4;
155     ssdA/= w*h;
156
157     printf(" SSD=%5"PRId64",%5"PRId64",%5"PRId64",%5"PRId64"\n",
158            ssdY, ssdU, ssdV, ssdA);
159
160 end:
161
162     sws_freeContext(dstContext);
163     sws_freeContext(outContext);
164
165     for (i=0; i<4; i++) {
166         av_free(dst[i]);
167         av_free(out[i]);
168     }
169
170     return res;
171 }
172
173 static void selfTest(uint8_t *ref[4], int refStride[4], int w, int h)
174 {
175     const int flags[] = { SWS_FAST_BILINEAR,
176                           SWS_BILINEAR, SWS_BICUBIC,
177                           SWS_X       , SWS_POINT  , SWS_AREA, 0 };
178     const int srcW = w;
179     const int srcH = h;
180     const int dstW[] = { srcW - srcW/3, srcW, srcW + srcW/3, 0 };
181     const int dstH[] = { srcH - srcH/3, srcH, srcH + srcH/3, 0 };
182     enum PixelFormat srcFormat, dstFormat;
183
184     for (srcFormat = 0; srcFormat < PIX_FMT_NB; srcFormat++) {
185         if (!sws_isSupportedInput(srcFormat) || !sws_isSupportedOutput(srcFormat))
186             continue;
187
188         for (dstFormat = 0; dstFormat < PIX_FMT_NB; dstFormat++) {
189             int i, j, k;
190             int res = 0;
191
192             if (!sws_isSupportedInput(dstFormat) || !sws_isSupportedOutput(dstFormat))
193                 continue;
194
195             printf("%s -> %s\n",
196                    av_pix_fmt_descriptors[srcFormat].name,
197                    av_pix_fmt_descriptors[dstFormat].name);
198             fflush(stdout);
199
200             for (k = 0; flags[k] && !res; k++) {
201                 struct SwsContext *srcContext = NULL;
202                 uint8_t *src[4] = {0};
203                 int srcStride[4];
204                 int p;
205                 for (p = 0; p < 4; p++) {
206                     if (srcFormat == PIX_FMT_RGB24 ||
207                         srcFormat == PIX_FMT_BGR24)
208                         srcStride[p] = srcW*3;
209                     else if (srcFormat==PIX_FMT_RGB48BE ||
210                              srcFormat==PIX_FMT_RGB48LE)
211                         srcStride[p] = srcW*6;
212                     else
213                         srcStride[p] = srcW*4;
214                     src[p] = av_mallocz(srcStride[p]*srcH+16);
215                     if (!src[p]) {
216                         perror("Malloc");
217                         return;
218                     }
219                 }
220                 srcContext = sws_getContext(w, h, PIX_FMT_YUVA420P, srcW, srcH,
221                                             srcFormat, flags[k], NULL, NULL, NULL);
222                 if (!srcContext) {
223                    fprintf(stderr, "Failed to get %s ---> %s\n",
224                             av_pix_fmt_descriptors[PIX_FMT_YUVA420P].name,
225                             av_pix_fmt_descriptors[srcFormat].name);
226                    return;
227                 }
228                 sws_scale(srcContext, ref, refStride, 0, h, src, srcStride);
229
230                 for (i = 0; dstW[i] && !res; i++)
231                     for (j = 0; dstH[j] && !res; j++)
232                         res = doTest(ref, refStride, w, h, src, srcStride,
233                                      srcFormat, dstFormat,
234                                      srcW, srcH, dstW[i], dstH[j], flags[k]);
235
236                 sws_freeContext(srcContext);
237                 for (p = 0; p < 4; p++)
238                     av_free(src[p]);
239             }
240         }
241     }
242 }
243
244 #define W 96
245 #define H 96
246
247 int main(int argc, char **argv)
248 {
249     uint8_t *rgb_data = av_malloc (W*H*4);
250     uint8_t *rgb_src[3]= {rgb_data, NULL, NULL};
251     int rgb_stride[3]={4*W, 0, 0};
252     uint8_t *data = av_malloc (4*W*H);
253     uint8_t *src[4]= {data, data+W*H, data+W*H*2, data+W*H*3};
254     int stride[4]={W, W, W, W};
255     int x, y;
256     struct SwsContext *sws;
257     AVLFG rand;
258
259     if (!rgb_data || !data)
260         return -1;
261
262     sws= sws_getContext(W/12, H/12, PIX_FMT_RGB32, W, H, PIX_FMT_YUVA420P, SWS_BILINEAR, NULL, NULL, NULL);
263
264     av_lfg_init(&rand, 1);
265
266     for (y=0; y<H; y++) {
267         for (x=0; x<W*4; x++) {
268             rgb_data[ x + y*4*W]= av_lfg_get(&rand);
269         }
270     }
271     sws_scale(sws, rgb_src, rgb_stride, 0, H, src, stride);
272     sws_freeContext(sws);
273     av_free(rgb_data);
274
275     selfTest(src, stride, W, H);
276     av_free(data);
277
278     return 0;
279 }