scaling-test.c revision 317c648b 
1/*
2 * Test program, which can detect problems with nearest neighbout scaling
3 * implementation. Also SRC and OVER opetations tested for 16bpp and 32bpp
4 * images.
5 *
6 * Just run it without any command line arguments, and it will report either
7 *   "scaling test passed" - everything is ok
8 *   "scaling test failed!" - there is some problem
9 *
10 * In the case of failure, finding the problem involves the following steps:
11 * 1. Get the reference 'scaling-test' binary. It makes sense to disable all
12 *    the cpu specific optimizations in pixman and also configure it with
13 *    '--disable-shared' option. Those who are paranoid can also tweak the
14 *    sources to disable all fastpath functions. The resulting binary
15 *    can be renamed to something like 'scaling-test.ref'.
16 * 2. Compile the buggy binary (also with the '--disable-shared' option).
17 * 3. Run 'ruby scaling-test-bisect.rb ./scaling-test.ref ./scaling-test'
18 * 4. Look at the information about failed case (destination buffer content
19 *    will be shown) and try to figure out what is wrong. It is possible
20 *    to use debugging print to stderr in pixman to get more information,
21 *    this does not interfere with the testing script.
22 */
23#include <assert.h>
24#include <stdlib.h>
25#include <stdio.h>
26#include "pixman.h"
27
28/* A primitive pseudorandom number generator, taken from POSIX.1-2001 example */
29
30static uint32_t lcg_seed;
31
32uint32_t lcg_rand(void)
33{
34    lcg_seed = lcg_seed * 1103515245 + 12345;
35    return ((uint32_t)(lcg_seed / 65536) % 32768);
36}
37
38void lcg_srand(uint32_t seed)
39{
40    lcg_seed = seed;
41}
42
43uint32_t lcg_rand_n(int max)
44{
45    return lcg_rand() % max;
46}
47
48/*----------------------------------------------------------------------------*\
49 *  CRC-32 version 2.0.0 by Craig Bruce, 2006-04-29.
50 *
51 *  This program generates the CRC-32 values for the files named in the
52 *  command-line arguments.  These are the same CRC-32 values used by GZIP,
53 *  PKZIP, and ZMODEM.  The Crc32_ComputeBuf() can also be detached and
54 *  used independently.
55 *
56 *  THIS PROGRAM IS PUBLIC-DOMAIN SOFTWARE.
57 *
58 *  Based on the byte-oriented implementation "File Verification Using CRC"
59 *  by Mark R. Nelson in Dr. Dobb's Journal, May 1992, pp. 64-67.
60 *
61 *  v1.0.0: original release.
62 *  v1.0.1: fixed printf formats.
63 *  v1.0.2: fixed something else.
64 *  v1.0.3: replaced CRC constant table by generator function.
65 *  v1.0.4: reformatted code, made ANSI C.  1994-12-05.
66 *  v2.0.0: rewrote to use memory buffer & static table, 2006-04-29.
67\*----------------------------------------------------------------------------*/
68
69/*----------------------------------------------------------------------------*\
70 *  NAME:
71 *     Crc32_ComputeBuf() - computes the CRC-32 value of a memory buffer
72 *  DESCRIPTION:
73 *     Computes or accumulates the CRC-32 value for a memory buffer.
74 *     The 'inCrc32' gives a previously accumulated CRC-32 value to allow
75 *     a CRC to be generated for multiple sequential buffer-fuls of data.
76 *     The 'inCrc32' for the first buffer must be zero.
77 *  ARGUMENTS:
78 *     inCrc32 - accumulated CRC-32 value, must be 0 on first call
79 *     buf     - buffer to compute CRC-32 value for
80 *     bufLen  - number of bytes in buffer
81 *  RETURNS:
82 *     crc32 - computed CRC-32 value
83 *  ERRORS:
84 *     (no errors are possible)
85\*----------------------------------------------------------------------------*/
86
87static uint32_t Crc32_ComputeBuf( uint32_t inCrc32, const void *buf,
88                                       size_t bufLen )
89{
90    static const uint32_t crcTable[256] = {
91   0x00000000,0x77073096,0xEE0E612C,0x990951BA,0x076DC419,0x706AF48F,0xE963A535,
92   0x9E6495A3,0x0EDB8832,0x79DCB8A4,0xE0D5E91E,0x97D2D988,0x09B64C2B,0x7EB17CBD,
93   0xE7B82D07,0x90BF1D91,0x1DB71064,0x6AB020F2,0xF3B97148,0x84BE41DE,0x1ADAD47D,
94   0x6DDDE4EB,0xF4D4B551,0x83D385C7,0x136C9856,0x646BA8C0,0xFD62F97A,0x8A65C9EC,
95   0x14015C4F,0x63066CD9,0xFA0F3D63,0x8D080DF5,0x3B6E20C8,0x4C69105E,0xD56041E4,
96   0xA2677172,0x3C03E4D1,0x4B04D447,0xD20D85FD,0xA50AB56B,0x35B5A8FA,0x42B2986C,
97   0xDBBBC9D6,0xACBCF940,0x32D86CE3,0x45DF5C75,0xDCD60DCF,0xABD13D59,0x26D930AC,
98   0x51DE003A,0xC8D75180,0xBFD06116,0x21B4F4B5,0x56B3C423,0xCFBA9599,0xB8BDA50F,
99   0x2802B89E,0x5F058808,0xC60CD9B2,0xB10BE924,0x2F6F7C87,0x58684C11,0xC1611DAB,
100   0xB6662D3D,0x76DC4190,0x01DB7106,0x98D220BC,0xEFD5102A,0x71B18589,0x06B6B51F,
101   0x9FBFE4A5,0xE8B8D433,0x7807C9A2,0x0F00F934,0x9609A88E,0xE10E9818,0x7F6A0DBB,
102   0x086D3D2D,0x91646C97,0xE6635C01,0x6B6B51F4,0x1C6C6162,0x856530D8,0xF262004E,
103   0x6C0695ED,0x1B01A57B,0x8208F4C1,0xF50FC457,0x65B0D9C6,0x12B7E950,0x8BBEB8EA,
104   0xFCB9887C,0x62DD1DDF,0x15DA2D49,0x8CD37CF3,0xFBD44C65,0x4DB26158,0x3AB551CE,
105   0xA3BC0074,0xD4BB30E2,0x4ADFA541,0x3DD895D7,0xA4D1C46D,0xD3D6F4FB,0x4369E96A,
106   0x346ED9FC,0xAD678846,0xDA60B8D0,0x44042D73,0x33031DE5,0xAA0A4C5F,0xDD0D7CC9,
107   0x5005713C,0x270241AA,0xBE0B1010,0xC90C2086,0x5768B525,0x206F85B3,0xB966D409,
108   0xCE61E49F,0x5EDEF90E,0x29D9C998,0xB0D09822,0xC7D7A8B4,0x59B33D17,0x2EB40D81,
109   0xB7BD5C3B,0xC0BA6CAD,0xEDB88320,0x9ABFB3B6,0x03B6E20C,0x74B1D29A,0xEAD54739,
110   0x9DD277AF,0x04DB2615,0x73DC1683,0xE3630B12,0x94643B84,0x0D6D6A3E,0x7A6A5AA8,
111   0xE40ECF0B,0x9309FF9D,0x0A00AE27,0x7D079EB1,0xF00F9344,0x8708A3D2,0x1E01F268,
112   0x6906C2FE,0xF762575D,0x806567CB,0x196C3671,0x6E6B06E7,0xFED41B76,0x89D32BE0,
113   0x10DA7A5A,0x67DD4ACC,0xF9B9DF6F,0x8EBEEFF9,0x17B7BE43,0x60B08ED5,0xD6D6A3E8,
114   0xA1D1937E,0x38D8C2C4,0x4FDFF252,0xD1BB67F1,0xA6BC5767,0x3FB506DD,0x48B2364B,
115   0xD80D2BDA,0xAF0A1B4C,0x36034AF6,0x41047A60,0xDF60EFC3,0xA867DF55,0x316E8EEF,
116   0x4669BE79,0xCB61B38C,0xBC66831A,0x256FD2A0,0x5268E236,0xCC0C7795,0xBB0B4703,
117   0x220216B9,0x5505262F,0xC5BA3BBE,0xB2BD0B28,0x2BB45A92,0x5CB36A04,0xC2D7FFA7,
118   0xB5D0CF31,0x2CD99E8B,0x5BDEAE1D,0x9B64C2B0,0xEC63F226,0x756AA39C,0x026D930A,
119   0x9C0906A9,0xEB0E363F,0x72076785,0x05005713,0x95BF4A82,0xE2B87A14,0x7BB12BAE,
120   0x0CB61B38,0x92D28E9B,0xE5D5BE0D,0x7CDCEFB7,0x0BDBDF21,0x86D3D2D4,0xF1D4E242,
121   0x68DDB3F8,0x1FDA836E,0x81BE16CD,0xF6B9265B,0x6FB077E1,0x18B74777,0x88085AE6,
122   0xFF0F6A70,0x66063BCA,0x11010B5C,0x8F659EFF,0xF862AE69,0x616BFFD3,0x166CCF45,
123   0xA00AE278,0xD70DD2EE,0x4E048354,0x3903B3C2,0xA7672661,0xD06016F7,0x4969474D,
124   0x3E6E77DB,0xAED16A4A,0xD9D65ADC,0x40DF0B66,0x37D83BF0,0xA9BCAE53,0xDEBB9EC5,
125   0x47B2CF7F,0x30B5FFE9,0xBDBDF21C,0xCABAC28A,0x53B39330,0x24B4A3A6,0xBAD03605,
126   0xCDD70693,0x54DE5729,0x23D967BF,0xB3667A2E,0xC4614AB8,0x5D681B02,0x2A6F2B94,
127   0xB40BBE37,0xC30C8EA1,0x5A05DF1B,0x2D02EF8D };
128    uint32_t crc32;
129    unsigned char *byteBuf;
130    size_t i;
131
132    /** accumulate crc32 for buffer **/
133    crc32 = inCrc32 ^ 0xFFFFFFFF;
134    byteBuf = (unsigned char*) buf;
135    for (i=0; i < bufLen; i++) {
136        crc32 = (crc32 >> 8) ^ crcTable[ (crc32 ^ byteBuf[i]) & 0xFF ];
137    }
138    return( crc32 ^ 0xFFFFFFFF );
139}
140
141
142#define MAX_SRC_WIDTH  10
143#define MAX_SRC_HEIGHT 10
144#define MAX_DST_WIDTH  10
145#define MAX_DST_HEIGHT 10
146#define MAX_STRIDE     4
147
148/*
149 * Composite operation with pseudorandom images
150 */
151uint32_t test_composite(uint32_t initcrc, int testnum, int verbose)
152{
153    int i;
154    pixman_image_t *src_img;
155    pixman_image_t *dst_img;
156    pixman_transform_t transform;
157    pixman_region16_t clip;
158    int src_width, src_height;
159    int dst_width, dst_height;
160    int src_stride, dst_stride;
161    int src_x, src_y;
162    int dst_x, dst_y;
163    int src_bpp;
164    int dst_bpp;
165    int w, h;
166    int scale_x = 32768, scale_y = 32768;
167    int op;
168    int repeat = 0;
169    int src_fmt, dst_fmt;
170    uint32_t *srcbuf;
171    uint32_t *dstbuf;
172    uint32_t crc32;
173
174    lcg_srand(testnum);
175
176    src_bpp = (lcg_rand_n(2) == 0) ? 2 : 4;
177    dst_bpp = (lcg_rand_n(2) == 0) ? 2 : 4;
178    op = (lcg_rand_n(2) == 0) ? PIXMAN_OP_SRC : PIXMAN_OP_OVER;
179
180    src_width = lcg_rand_n(MAX_SRC_WIDTH) + 1;
181    src_height = lcg_rand_n(MAX_SRC_HEIGHT) + 1;
182    dst_width = lcg_rand_n(MAX_DST_WIDTH) + 1;
183    dst_height = lcg_rand_n(MAX_DST_HEIGHT) + 1;
184    src_stride = src_width * src_bpp + lcg_rand_n(MAX_STRIDE) * src_bpp;
185    dst_stride = dst_width * dst_bpp + lcg_rand_n(MAX_STRIDE) * dst_bpp;
186    if (src_stride & 3) src_stride += 2;
187    if (dst_stride & 3) dst_stride += 2;
188
189    src_x = -(src_width / 4) + lcg_rand_n(src_width * 3 / 2);
190    src_y = -(src_height / 4) + lcg_rand_n(src_height * 3 / 2);
191    dst_x = -(dst_width / 4) + lcg_rand_n(dst_width * 3 / 2);
192    dst_y = -(dst_height / 4) + lcg_rand_n(dst_height * 3 / 2);
193    w = lcg_rand_n(dst_width * 3 / 2 - dst_x);
194    h = lcg_rand_n(dst_height * 3 / 2 - dst_y);
195
196    srcbuf = (uint32_t *)malloc(src_stride * src_height);
197    dstbuf = (uint32_t *)malloc(dst_stride * dst_height);
198    for (i = 0; i < src_stride * src_height; i++)
199        *((uint8_t *)srcbuf + i) = lcg_rand_n(256);
200    for (i = 0; i < dst_stride * dst_height; i++)
201        *((uint8_t *)dstbuf + i) = lcg_rand_n(256);
202
203    src_fmt = src_bpp == 4 ? (lcg_rand_n(2) == 0 ?
204        PIXMAN_a8r8g8b8 : PIXMAN_x8r8g8b8) : PIXMAN_r5g6b5;
205
206    dst_fmt = dst_bpp == 4 ? (lcg_rand_n(2) == 0 ?
207        PIXMAN_a8r8g8b8 : PIXMAN_x8r8g8b8) : PIXMAN_r5g6b5;
208
209    src_img = pixman_image_create_bits(
210        src_fmt, src_width, src_height, srcbuf, src_stride);
211
212    dst_img = pixman_image_create_bits(
213        dst_fmt, dst_width, dst_height, dstbuf, dst_stride);
214
215    if (lcg_rand_n(8) > 0) {
216        scale_x = 32768 + lcg_rand_n(65536);
217        scale_y = 32768 + lcg_rand_n(65536);
218        pixman_transform_init_scale(&transform, scale_x, scale_y);
219        pixman_image_set_transform(src_img, &transform);
220    }
221
222    switch (lcg_rand_n(4)) {
223        case 0: repeat = PIXMAN_REPEAT_NONE; break;
224        case 1: repeat = PIXMAN_REPEAT_NORMAL; break;
225        case 2: repeat = PIXMAN_REPEAT_PAD; break;
226        case 3: repeat = PIXMAN_REPEAT_REFLECT; break;
227    }
228    pixman_image_set_repeat(src_img, repeat);
229
230    if (verbose) {
231        printf("src_fmt=%08X, dst_fmt=%08X\n", src_fmt, dst_fmt);
232        printf("op=%d, scale_x=%d, scale_y=%d, repeat=%d\n",
233            op, scale_x, scale_y, repeat);
234        printf("src_width=%d, src_height=%d, dst_width=%d, dst_height=%d\n",
235            src_width, src_height, dst_width, dst_height);
236        printf("src_x=%d, src_y=%d, dst_x=%d, dst_y=%d\n",
237            src_x, src_y, dst_x, dst_y);
238        printf("w=%d, h=%d\n", w, h);
239    }
240
241    if (lcg_rand_n(8) == 0) {
242        pixman_box16_t clip_boxes[2];
243        int n = lcg_rand_n(2) + 1;
244        for (i = 0; i < n; i++) {
245            clip_boxes[i].x1 = lcg_rand_n(src_width);
246            clip_boxes[i].y1 = lcg_rand_n(src_height);
247            clip_boxes[i].x2 = clip_boxes[i].x1 + lcg_rand_n(src_width - clip_boxes[i].x1);
248            clip_boxes[i].y2 = clip_boxes[i].y1 + lcg_rand_n(src_height - clip_boxes[i].y1);
249            if (verbose) {
250                printf("source clip box: [%d,%d-%d,%d]\n",
251                    clip_boxes[i].x1, clip_boxes[i].y1,
252                    clip_boxes[i].x2, clip_boxes[i].y2);
253            }
254        }
255        pixman_region_init_rects(&clip, clip_boxes, n);
256        pixman_image_set_clip_region(src_img, &clip);
257        pixman_image_set_source_clipping(src_img, 1);
258        pixman_region_fini(&clip);
259    }
260
261    if (lcg_rand_n(8) == 0) {
262        pixman_box16_t clip_boxes[2];
263        int n = lcg_rand_n(2) + 1;
264        for (i = 0; i < n; i++) {
265            clip_boxes[i].x1 = lcg_rand_n(dst_width);
266            clip_boxes[i].y1 = lcg_rand_n(dst_height);
267            clip_boxes[i].x2 = clip_boxes[i].x1 + lcg_rand_n(dst_width - clip_boxes[i].x1);
268            clip_boxes[i].y2 = clip_boxes[i].y1 + lcg_rand_n(dst_height - clip_boxes[i].y1);
269            if (verbose) {
270                printf("destination clip box: [%d,%d-%d,%d]\n",
271                    clip_boxes[i].x1, clip_boxes[i].y1,
272                    clip_boxes[i].x2, clip_boxes[i].y2);
273            }
274        }
275        pixman_region_init_rects(&clip, clip_boxes, n);
276        pixman_image_set_clip_region(dst_img, &clip);
277        pixman_region_fini(&clip);
278    }
279
280    pixman_image_composite (op, src_img, NULL, dst_img,
281                            src_x, src_y, 0, 0, dst_x, dst_y, w, h);
282
283    if (dst_fmt == PIXMAN_x8r8g8b8) {
284        /* ignore unused part */
285        for (i = 0; i < dst_stride * dst_height / 4; i++)
286            dstbuf[i] &= 0xFFFFFF;
287    }
288
289    if (verbose) {
290        int j;
291        for (i = 0; i < dst_height; i++) {
292            for (j = 0; j < dst_stride; j++) {
293                printf("%02X ", *((uint8_t *)dstbuf + i * dst_stride + j));
294            }
295            printf("\n");
296        }
297    }
298
299    pixman_image_unref (src_img);
300    pixman_image_unref (dst_img);
301
302    crc32 = Crc32_ComputeBuf(initcrc, dstbuf, dst_stride * dst_height);
303    free(srcbuf);
304    free(dstbuf);
305    return crc32;
306}
307
308int main(int argc, char *argv[])
309{
310    int i, n = 0;
311    uint32_t crc = 0;
312
313    if (argc >= 2)
314        n = atoi(argv[1]);
315
316    if (n == 0) n = 3000000;
317
318    if (n < 0) {
319        crc = test_composite(0, -n, 1);
320        printf("crc32=%08X\n", crc);
321    }
322    else {
323        for (i = 1; i <= n; i++)
324        {
325            crc = test_composite(crc, i, 0);
326        }
327        printf("crc32=%08X\n", crc);
328#ifdef LITTLE_ENDIAN
329        if (n == 3000000) {
330            /* predefined value for running with all the fastpath functions disabled  */
331            /* it needs to be updated every time changes are introduced to this program! */
332            if (crc == 0xC950E5BB) {
333                printf("scaling test passed\n");
334            } else {
335                printf("scaling test failed!\n");
336            }
337        }
338#endif
339    }
340    return 0;
341}
342