More accurate value for log2(10)
[libav.git] / libswscale / swscale-test.c
CommitLineData
c9b99ea6 1/*
d026b45e
DB
2 * Copyright (C) 2003 Michael Niedermayer <michaelni@gmx.at>
3 *
4 * This file is part of FFmpeg.
5 *
54e3ae2e
DB
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.
d026b45e
DB
10 *
11 * FFmpeg is distributed in the hope that it will be useful,
12 * but WITHOUT ANY WARRANTY; without even the implied warranty of
54e3ae2e
DB
13 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
14 * Lesser General Public License for more details.
d026b45e 15 *
54e3ae2e
DB
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
b19bcbaa 18 * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
d026b45e 19 */
c9b99ea6
MN
20
21#include <stdio.h>
22#include <stdlib.h>
23#include <string.h>
24#include <inttypes.h>
25#include <stdarg.h>
26
e9e12f0e 27#undef HAVE_AV_CONFIG_H
51d83986 28#include "libavutil/mem.h"
83da2c6f 29#include "libavutil/avutil.h"
9a034dde 30#include "libavutil/lfg.h"
c9b99ea6 31#include "swscale.h"
c9b99ea6 32
71ed6c8f
RP
33/* HACK Duplicated from swscale_internal.h.
34 * Should be removed when a cleaner pixel format system exists. */
35const char *sws_format_name(enum PixelFormat format);
f5a9c4ee
RP
36#define isGray(x) ( \
37 (x)==PIX_FMT_GRAY8 \
38 || (x)==PIX_FMT_GRAY16BE \
39 || (x)==PIX_FMT_GRAY16LE \
40 )
a48ce2c3 41#define hasChroma(x) (!( \
d426ea20
RP
42 isGray(x) \
43 || (x)==PIX_FMT_MONOBLACK \
44 || (x)==PIX_FMT_MONOWHITE \
45 ))
71ed6c8f
RP
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 )
dc822901 53
dd68318c
RP
54static uint64_t getSSD(uint8_t *src1, uint8_t *src2, int stride1, int stride2, int w, int h)
55{
221b804f
DB
56 int x,y;
57 uint64_t ssd=0;
c9b99ea6
MN
58
59//printf("%d %d\n", w, h);
6a4970ab 60
dd68318c
RP
61 for (y=0; y<h; y++) {
62 for (x=0; x<w; x++) {
221b804f
DB
63 int d= src1[x + y*stride1] - src2[x + y*stride2];
64 ssd+= d*d;
c9b99ea6 65//printf("%d", abs(src1[x + y*stride1] - src2[x + y*stride2])/26 );
221b804f 66 }
c9b99ea6 67//printf("\n");
221b804f
DB
68 }
69 return ssd;
c9b99ea6
MN
70}
71
72// test by ref -> src -> dst -> out & compare out against ref
73// ref & out are YV12
2d8d885b
SS
74static int doTest(uint8_t *ref[4], int refStride[4], int w, int h,
75 enum PixelFormat srcFormat, enum PixelFormat dstFormat,
dd68318c
RP
76 int srcW, int srcH, int dstW, int dstH, int flags)
77{
1cd98da8
BF
78 uint8_t *src[4] = {0};
79 uint8_t *dst[4] = {0};
80 uint8_t *out[4] = {0};
110f147a 81 int srcStride[4], dstStride[4];
221b804f 82 int i;
f5a9c4ee 83 uint64_t ssdY, ssdU=0, ssdV=0, ssdA=0;
54ea7736
DB
84 struct SwsContext *srcContext = NULL, *dstContext = NULL,
85 *outContext = NULL;
221b804f
DB
86 int res;
87
88 res = 0;
dd68318c 89 for (i=0; i<4; i++) {
221b804f
DB
90 // avoid stride % bpp != 0
91 if (srcFormat==PIX_FMT_RGB24 || srcFormat==PIX_FMT_BGR24)
92 srcStride[i]= srcW*3;
5d5a8813
KS
93 else if (srcFormat==PIX_FMT_RGB48BE || srcFormat==PIX_FMT_RGB48LE)
94 srcStride[i]= srcW*6;
221b804f
DB
95 else
96 srcStride[i]= srcW*4;
97
98 if (dstFormat==PIX_FMT_RGB24 || dstFormat==PIX_FMT_BGR24)
99 dstStride[i]= dstW*3;
5d5a8813
KS
100 else if (dstFormat==PIX_FMT_RGB48BE || dstFormat==PIX_FMT_RGB48LE)
101 dstStride[i]= dstW*6;
221b804f
DB
102 else
103 dstStride[i]= dstW*4;
104
9e5f164d
RP
105 /* Image buffers passed into libswscale can be allocated any way you
106 * prefer, as long as they're aligned enough for the architecture, and
107 * they're freed appropriately (such as using av_free for buffers
108 * allocated with av_malloc). */
7dac0c43
SS
109 src[i]= av_mallocz(srcStride[i]*srcH);
110 dst[i]= av_mallocz(dstStride[i]*dstH);
111 out[i]= av_mallocz(refStride[i]*h);
1b0a4572 112 if (!src[i] || !dst[i] || !out[i]) {
221b804f
DB
113 perror("Malloc");
114 res = -1;
115
116 goto end;
117 }
118 }
119
110f147a 120 srcContext= sws_getContext(w, h, PIX_FMT_YUVA420P, srcW, srcH, srcFormat, flags, NULL, NULL, NULL);
1b0a4572 121 if (!srcContext) {
221b804f 122 fprintf(stderr, "Failed to get %s ---> %s\n",
110f147a 123 sws_format_name(PIX_FMT_YUVA420P),
221b804f
DB
124 sws_format_name(srcFormat));
125 res = -1;
126
127 goto end;
128 }
129 dstContext= sws_getContext(srcW, srcH, srcFormat, dstW, dstH, dstFormat, flags, NULL, NULL, NULL);
1b0a4572 130 if (!dstContext) {
221b804f
DB
131 fprintf(stderr, "Failed to get %s ---> %s\n",
132 sws_format_name(srcFormat),
133 sws_format_name(dstFormat));
134 res = -1;
135
136 goto end;
137 }
110f147a 138 outContext= sws_getContext(dstW, dstH, dstFormat, w, h, PIX_FMT_YUVA420P, flags, NULL, NULL, NULL);
1b0a4572 139 if (!outContext) {
221b804f
DB
140 fprintf(stderr, "Failed to get %s ---> %s\n",
141 sws_format_name(dstFormat),
110f147a 142 sws_format_name(PIX_FMT_YUVA420P));
221b804f
DB
143 res = -1;
144
145 goto end;
146 }
147// printf("test %X %X %X -> %X %X %X\n", (int)ref[0], (int)ref[1], (int)ref[2],
148// (int)src[0], (int)src[1], (int)src[2]);
149
150 sws_scale(srcContext, ref, refStride, 0, h , src, srcStride);
151 sws_scale(dstContext, src, srcStride, 0, srcH, dst, dstStride);
152 sws_scale(outContext, dst, dstStride, 0, dstH, out, refStride);
27161c06 153
221b804f 154 ssdY= getSSD(ref[0], out[0], refStride[0], refStride[0], w, h);
d426ea20 155 if (hasChroma(srcFormat) && hasChroma(dstFormat)) {
f5a9c4ee 156 //FIXME check that output is really gray
475f3613
RP
157 ssdU= getSSD(ref[1], out[1], refStride[1], refStride[1], (w+1)>>1, (h+1)>>1);
158 ssdV= getSSD(ref[2], out[2], refStride[2], refStride[2], (w+1)>>1, (h+1)>>1);
f5a9c4ee 159 }
110f147a
CS
160 if (isALPHA(srcFormat) && isALPHA(dstFormat))
161 ssdA= getSSD(ref[3], out[3], refStride[3], refStride[3], w, h);
6a4970ab 162
221b804f
DB
163 ssdY/= w*h;
164 ssdU/= w*h/4;
165 ssdV/= w*h/4;
110f147a 166 ssdA/= w*h;
6a4970ab 167
4ed88643 168 printf(" %s %dx%d -> %s %4dx%4d flags=%2d SSD=%5"PRId64",%5"PRId64",%5"PRId64",%5"PRId64"\n",
ed2164f3
DB
169 sws_format_name(srcFormat), srcW, srcH,
170 sws_format_name(dstFormat), dstW, dstH,
110f147a 171 flags, ssdY, ssdU, ssdV, ssdA);
ed2164f3 172 fflush(stdout);
c9b99ea6 173
a48ce2c3 174end:
6a4970ab 175
221b804f
DB
176 sws_freeContext(srcContext);
177 sws_freeContext(dstContext);
178 sws_freeContext(outContext);
c9b99ea6 179
dd68318c 180 for (i=0; i<4; i++) {
51d83986
MR
181 av_free(src[i]);
182 av_free(dst[i]);
183 av_free(out[i]);
221b804f 184 }
c9b99ea6 185
221b804f 186 return res;
c9b99ea6
MN
187}
188
89fb9516 189static void selfTest(uint8_t *ref[4], int refStride[4], int w, int h)
dd68318c 190{
6b052ba8 191 const int flags[] = { SWS_FAST_BILINEAR,
dd653355
RP
192 SWS_BILINEAR, SWS_BICUBIC,
193 SWS_X , SWS_POINT , SWS_AREA, 0 };
6b052ba8
RP
194 const int srcW = w;
195 const int srcH = h;
196 const int dstW[] = { srcW - srcW/3, srcW, srcW + srcW/3, 0 };
197 const int dstH[] = { srcH - srcH/3, srcH, srcH + srcH/3, 0 };
221b804f 198 enum PixelFormat srcFormat, dstFormat;
221b804f
DB
199
200 for (srcFormat = 0; srcFormat < PIX_FMT_NB; srcFormat++) {
613f2dc4 201 if (!sws_isSupportedInput(srcFormat) || !sws_isSupportedOutput(srcFormat))
731c04ad
SS
202 continue;
203
221b804f 204 for (dstFormat = 0; dstFormat < PIX_FMT_NB; dstFormat++) {
6b052ba8 205 int i, j, k;
e55ed689
RP
206 int res = 0;
207
613f2dc4 208 if (!sws_isSupportedInput(dstFormat) || !sws_isSupportedOutput(dstFormat))
731c04ad
SS
209 continue;
210
221b804f
DB
211 printf("%s -> %s\n",
212 sws_format_name(srcFormat),
213 sws_format_name(dstFormat));
ae9e0e83 214 fflush(stdout);
221b804f 215
6b052ba8
RP
216 for (i = 0; dstW[i] && !res; i++)
217 for (j = 0; dstH[j] && !res; j++)
218 for (k = 0; flags[k] && !res; k++)
89fb9516 219 res = doTest(ref, refStride, w, h, srcFormat, dstFormat,
6b052ba8 220 srcW, srcH, dstW[i], dstH[j], flags[k]);
221b804f
DB
221 }
222 }
c9b99ea6
MN
223}
224
225#define W 96
226#define H 96
227
dd68318c
RP
228int main(int argc, char **argv)
229{
51d83986 230 uint8_t *rgb_data = av_malloc (W*H*4);
221b804f
DB
231 uint8_t *rgb_src[3]= {rgb_data, NULL, NULL};
232 int rgb_stride[3]={4*W, 0, 0};
51d83986 233 uint8_t *data = av_malloc (4*W*H);
110f147a
CS
234 uint8_t *src[4]= {data, data+W*H, data+W*H*2, data+W*H*3};
235 int stride[4]={W, W, W, W};
221b804f
DB
236 int x, y;
237 struct SwsContext *sws;
9a034dde 238 AVLFG rand;
221b804f 239
7248797c
RP
240 if (!rgb_data || !data)
241 return -1;
242
7da35acd 243 sws= sws_getContext(W/12, H/12, PIX_FMT_RGB32, W, H, PIX_FMT_YUVA420P, SWS_BILINEAR, NULL, NULL, NULL);
221b804f 244
9a034dde
RP
245 av_lfg_init(&rand, 1);
246
dd68318c
RP
247 for (y=0; y<H; y++) {
248 for (x=0; x<W*4; x++) {
9a034dde 249 rgb_data[ x + y*4*W]= av_lfg_get(&rand);
221b804f
DB
250 }
251 }
e5091488 252 sws_scale(sws, rgb_src, rgb_stride, 0, H, src, stride);
64d39b7c 253 sws_freeContext(sws);
51d83986 254 av_free(rgb_data);
27161c06 255
e5091488 256 selfTest(src, stride, W, H);
51d83986 257 av_free(data);
516b1f82 258
d7cf0a6f 259 return 0;
c9b99ea6 260}