11176bdada62cabc6ec4b0308a930e83b679d5d36John Reck/* 21176bdada62cabc6ec4b0308a930e83b679d5d36John Reck * Copyright © 2012 Siarhei Siamashka <siarhei.siamashka@gmail.com> 31176bdada62cabc6ec4b0308a930e83b679d5d36John Reck * 41176bdada62cabc6ec4b0308a930e83b679d5d36John Reck * Permission is hereby granted, free of charge, to any person obtaining a 51176bdada62cabc6ec4b0308a930e83b679d5d36John Reck * copy of this software and associated documentation files (the "Software"), 61176bdada62cabc6ec4b0308a930e83b679d5d36John Reck * to deal in the Software without restriction, including without limitation 71176bdada62cabc6ec4b0308a930e83b679d5d36John Reck * the rights to use, copy, modify, merge, publish, distribute, sublicense, 81176bdada62cabc6ec4b0308a930e83b679d5d36John Reck * and/or sell copies of the Software, and to permit persons to whom the 91176bdada62cabc6ec4b0308a930e83b679d5d36John Reck * Software is furnished to do so, subject to the following conditions: 101176bdada62cabc6ec4b0308a930e83b679d5d36John Reck * 111176bdada62cabc6ec4b0308a930e83b679d5d36John Reck * The above copyright notice and this permission notice (including the next 121176bdada62cabc6ec4b0308a930e83b679d5d36John Reck * paragraph) shall be included in all copies or substantial portions of the 131176bdada62cabc6ec4b0308a930e83b679d5d36John Reck * Software. 141176bdada62cabc6ec4b0308a930e83b679d5d36John Reck * 151176bdada62cabc6ec4b0308a930e83b679d5d36John Reck * THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR 161176bdada62cabc6ec4b0308a930e83b679d5d36John Reck * IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY, 171176bdada62cabc6ec4b0308a930e83b679d5d36John Reck * FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL 181176bdada62cabc6ec4b0308a930e83b679d5d36John Reck * THE AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER 191176bdada62cabc6ec4b0308a930e83b679d5d36John Reck * LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING 201176bdada62cabc6ec4b0308a930e83b679d5d36John Reck * FROM, OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER 211176bdada62cabc6ec4b0308a930e83b679d5d36John Reck * DEALINGS IN THE SOFTWARE. 221176bdada62cabc6ec4b0308a930e83b679d5d36John Reck */ 231176bdada62cabc6ec4b0308a930e83b679d5d36John Reck 241176bdada62cabc6ec4b0308a930e83b679d5d36John Reck#include "utils.h" 251176bdada62cabc6ec4b0308a930e83b679d5d36John Reck#include <assert.h> 261176bdada62cabc6ec4b0308a930e83b679d5d36John Reck#include <stdlib.h> 271176bdada62cabc6ec4b0308a930e83b679d5d36John Reck#include <stdio.h> 281176bdada62cabc6ec4b0308a930e83b679d5d36John Reck#include <math.h> 291176bdada62cabc6ec4b0308a930e83b679d5d36John Reck 301176bdada62cabc6ec4b0308a930e83b679d5d36John Reck#ifdef HAVE_FLOAT128 311176bdada62cabc6ec4b0308a930e83b679d5d36John Reck 321176bdada62cabc6ec4b0308a930e83b679d5d36John Reck#define pixman_fixed_to_float128(x) (((__float128)(x)) / 65536.0Q) 331176bdada62cabc6ec4b0308a930e83b679d5d36John Reck 341176bdada62cabc6ec4b0308a930e83b679d5d36John Recktypedef struct { __float128 v[3]; } pixman_vector_f128_t; 351176bdada62cabc6ec4b0308a930e83b679d5d36John Recktypedef struct { __float128 m[3][3]; } pixman_transform_f128_t; 361176bdada62cabc6ec4b0308a930e83b679d5d36John Reck 371176bdada62cabc6ec4b0308a930e83b679d5d36John Reckpixman_bool_t 381176bdada62cabc6ec4b0308a930e83b679d5d36John Reckpixman_transform_point_f128 (const pixman_transform_f128_t *t, 391176bdada62cabc6ec4b0308a930e83b679d5d36John Reck const pixman_vector_f128_t *v, 401176bdada62cabc6ec4b0308a930e83b679d5d36John Reck pixman_vector_f128_t *result) 411176bdada62cabc6ec4b0308a930e83b679d5d36John Reck{ 421176bdada62cabc6ec4b0308a930e83b679d5d36John Reck int i; 431176bdada62cabc6ec4b0308a930e83b679d5d36John Reck for (i = 0; i < 3; i++) 441176bdada62cabc6ec4b0308a930e83b679d5d36John Reck { 451176bdada62cabc6ec4b0308a930e83b679d5d36John Reck result->v[i] = t->m[i][0] * v->v[0] + 461176bdada62cabc6ec4b0308a930e83b679d5d36John Reck t->m[i][1] * v->v[1] + 471176bdada62cabc6ec4b0308a930e83b679d5d36John Reck t->m[i][2] * v->v[2]; 481176bdada62cabc6ec4b0308a930e83b679d5d36John Reck } 491176bdada62cabc6ec4b0308a930e83b679d5d36John Reck if (result->v[2] != 0) 501176bdada62cabc6ec4b0308a930e83b679d5d36John Reck { 511176bdada62cabc6ec4b0308a930e83b679d5d36John Reck result->v[0] /= result->v[2]; 521176bdada62cabc6ec4b0308a930e83b679d5d36John Reck result->v[1] /= result->v[2]; 531176bdada62cabc6ec4b0308a930e83b679d5d36John Reck result->v[2] = 1; 541176bdada62cabc6ec4b0308a930e83b679d5d36John Reck return TRUE; 551176bdada62cabc6ec4b0308a930e83b679d5d36John Reck } 561176bdada62cabc6ec4b0308a930e83b679d5d36John Reck else 571176bdada62cabc6ec4b0308a930e83b679d5d36John Reck { 581176bdada62cabc6ec4b0308a930e83b679d5d36John Reck return FALSE; 591176bdada62cabc6ec4b0308a930e83b679d5d36John Reck } 601176bdada62cabc6ec4b0308a930e83b679d5d36John Reck} 611176bdada62cabc6ec4b0308a930e83b679d5d36John Reck 621176bdada62cabc6ec4b0308a930e83b679d5d36John Reckpixman_bool_t does_it_fit_fixed_48_16 (__float128 x) 631176bdada62cabc6ec4b0308a930e83b679d5d36John Reck{ 641176bdada62cabc6ec4b0308a930e83b679d5d36John Reck if (x >= 65536.0Q * 65536.0Q * 32768.0Q) 651176bdada62cabc6ec4b0308a930e83b679d5d36John Reck return FALSE; 661176bdada62cabc6ec4b0308a930e83b679d5d36John Reck if (x <= -65536.0Q * 65536.0Q * 32768.0Q) 671176bdada62cabc6ec4b0308a930e83b679d5d36John Reck return FALSE; 681176bdada62cabc6ec4b0308a930e83b679d5d36John Reck return TRUE; 691176bdada62cabc6ec4b0308a930e83b679d5d36John Reck} 701176bdada62cabc6ec4b0308a930e83b679d5d36John Reck 711176bdada62cabc6ec4b0308a930e83b679d5d36John Reck#endif 721176bdada62cabc6ec4b0308a930e83b679d5d36John Reck 731176bdada62cabc6ec4b0308a930e83b679d5d36John Reckuint32_t 741176bdada62cabc6ec4b0308a930e83b679d5d36John Recktest_matrix (int testnum, int verbose) 751176bdada62cabc6ec4b0308a930e83b679d5d36John Reck{ 761176bdada62cabc6ec4b0308a930e83b679d5d36John Reck uint32_t crc32 = 0; 771176bdada62cabc6ec4b0308a930e83b679d5d36John Reck int i, j, k; 781176bdada62cabc6ec4b0308a930e83b679d5d36John Reck pixman_bool_t is_affine; 791176bdada62cabc6ec4b0308a930e83b679d5d36John Reck 801176bdada62cabc6ec4b0308a930e83b679d5d36John Reck prng_srand (testnum); 811176bdada62cabc6ec4b0308a930e83b679d5d36John Reck 821176bdada62cabc6ec4b0308a930e83b679d5d36John Reck for (i = 0; i < 100; i++) 831176bdada62cabc6ec4b0308a930e83b679d5d36John Reck { 841176bdada62cabc6ec4b0308a930e83b679d5d36John Reck pixman_bool_t transform_ok; 851176bdada62cabc6ec4b0308a930e83b679d5d36John Reck pixman_transform_t ti; 861176bdada62cabc6ec4b0308a930e83b679d5d36John Reck pixman_vector_48_16_t vi, result_i; 871176bdada62cabc6ec4b0308a930e83b679d5d36John Reck#ifdef HAVE_FLOAT128 881176bdada62cabc6ec4b0308a930e83b679d5d36John Reck pixman_transform_f128_t tf; 891176bdada62cabc6ec4b0308a930e83b679d5d36John Reck pixman_vector_f128_t vf, result_f; 901176bdada62cabc6ec4b0308a930e83b679d5d36John Reck#endif 911176bdada62cabc6ec4b0308a930e83b679d5d36John Reck prng_randmemset (&ti, sizeof(ti), 0); 921176bdada62cabc6ec4b0308a930e83b679d5d36John Reck prng_randmemset (&vi, sizeof(vi), 0); 931176bdada62cabc6ec4b0308a930e83b679d5d36John Reck 941176bdada62cabc6ec4b0308a930e83b679d5d36John Reck for (j = 0; j < 3; j++) 951176bdada62cabc6ec4b0308a930e83b679d5d36John Reck { 961176bdada62cabc6ec4b0308a930e83b679d5d36John Reck /* make sure that "vi" contains 31.16 fixed point data */ 971176bdada62cabc6ec4b0308a930e83b679d5d36John Reck vi.v[j] >>= 17; 981176bdada62cabc6ec4b0308a930e83b679d5d36John Reck /* and apply random shift */ 991176bdada62cabc6ec4b0308a930e83b679d5d36John Reck if (prng_rand_n (3) == 0) 1001176bdada62cabc6ec4b0308a930e83b679d5d36John Reck vi.v[j] >>= prng_rand_n (46); 1011176bdada62cabc6ec4b0308a930e83b679d5d36John Reck } 1021176bdada62cabc6ec4b0308a930e83b679d5d36John Reck 1031176bdada62cabc6ec4b0308a930e83b679d5d36John Reck if (prng_rand_n (2)) 1041176bdada62cabc6ec4b0308a930e83b679d5d36John Reck { 1051176bdada62cabc6ec4b0308a930e83b679d5d36John Reck /* random shift for the matrix */ 1061176bdada62cabc6ec4b0308a930e83b679d5d36John Reck for (j = 0; j < 3; j++) 1071176bdada62cabc6ec4b0308a930e83b679d5d36John Reck for (k = 0; k < 3; k++) 1081176bdada62cabc6ec4b0308a930e83b679d5d36John Reck ti.matrix[j][k] >>= prng_rand_n (30); 1091176bdada62cabc6ec4b0308a930e83b679d5d36John Reck } 1101176bdada62cabc6ec4b0308a930e83b679d5d36John Reck 1111176bdada62cabc6ec4b0308a930e83b679d5d36John Reck if (prng_rand_n (2)) 1121176bdada62cabc6ec4b0308a930e83b679d5d36John Reck { 1131176bdada62cabc6ec4b0308a930e83b679d5d36John Reck /* affine matrix */ 1141176bdada62cabc6ec4b0308a930e83b679d5d36John Reck ti.matrix[2][0] = 0; 1151176bdada62cabc6ec4b0308a930e83b679d5d36John Reck ti.matrix[2][1] = 0; 1161176bdada62cabc6ec4b0308a930e83b679d5d36John Reck ti.matrix[2][2] = pixman_fixed_1; 1171176bdada62cabc6ec4b0308a930e83b679d5d36John Reck } 1181176bdada62cabc6ec4b0308a930e83b679d5d36John Reck 1191176bdada62cabc6ec4b0308a930e83b679d5d36John Reck if (prng_rand_n (2)) 1201176bdada62cabc6ec4b0308a930e83b679d5d36John Reck { 1211176bdada62cabc6ec4b0308a930e83b679d5d36John Reck /* cartesian coordinates */ 1221176bdada62cabc6ec4b0308a930e83b679d5d36John Reck vi.v[2] = pixman_fixed_1; 1231176bdada62cabc6ec4b0308a930e83b679d5d36John Reck } 1241176bdada62cabc6ec4b0308a930e83b679d5d36John Reck 1251176bdada62cabc6ec4b0308a930e83b679d5d36John Reck is_affine = (ti.matrix[2][0] == 0 && ti.matrix[2][1] == 0 && 1261176bdada62cabc6ec4b0308a930e83b679d5d36John Reck ti.matrix[2][2] == pixman_fixed_1 && 1271176bdada62cabc6ec4b0308a930e83b679d5d36John Reck vi.v[2] == pixman_fixed_1); 1281176bdada62cabc6ec4b0308a930e83b679d5d36John Reck 1291176bdada62cabc6ec4b0308a930e83b679d5d36John Reck transform_ok = TRUE; 1301176bdada62cabc6ec4b0308a930e83b679d5d36John Reck if (is_affine && prng_rand_n (2)) 1311176bdada62cabc6ec4b0308a930e83b679d5d36John Reck pixman_transform_point_31_16_affine (&ti, &vi, &result_i); 1321176bdada62cabc6ec4b0308a930e83b679d5d36John Reck else 1331176bdada62cabc6ec4b0308a930e83b679d5d36John Reck transform_ok = pixman_transform_point_31_16 (&ti, &vi, &result_i); 1341176bdada62cabc6ec4b0308a930e83b679d5d36John Reck 1351176bdada62cabc6ec4b0308a930e83b679d5d36John Reck crc32 = compute_crc32 (crc32, &result_i, sizeof(result_i)); 1361176bdada62cabc6ec4b0308a930e83b679d5d36John Reck 1371176bdada62cabc6ec4b0308a930e83b679d5d36John Reck#ifdef HAVE_FLOAT128 1381176bdada62cabc6ec4b0308a930e83b679d5d36John Reck /* compare with a reference 128-bit floating point implementation */ 1391176bdada62cabc6ec4b0308a930e83b679d5d36John Reck for (j = 0; j < 3; j++) 1401176bdada62cabc6ec4b0308a930e83b679d5d36John Reck { 1411176bdada62cabc6ec4b0308a930e83b679d5d36John Reck vf.v[j] = pixman_fixed_to_float128 (vi.v[j]); 1421176bdada62cabc6ec4b0308a930e83b679d5d36John Reck for (k = 0; k < 3; k++) 1431176bdada62cabc6ec4b0308a930e83b679d5d36John Reck { 1441176bdada62cabc6ec4b0308a930e83b679d5d36John Reck tf.m[j][k] = pixman_fixed_to_float128 (ti.matrix[j][k]); 1451176bdada62cabc6ec4b0308a930e83b679d5d36John Reck } 1461176bdada62cabc6ec4b0308a930e83b679d5d36John Reck } 1471176bdada62cabc6ec4b0308a930e83b679d5d36John Reck 1481176bdada62cabc6ec4b0308a930e83b679d5d36John Reck if (pixman_transform_point_f128 (&tf, &vf, &result_f)) 1491176bdada62cabc6ec4b0308a930e83b679d5d36John Reck { 1501176bdada62cabc6ec4b0308a930e83b679d5d36John Reck if (transform_ok || 1511176bdada62cabc6ec4b0308a930e83b679d5d36John Reck (does_it_fit_fixed_48_16 (result_f.v[0]) && 1521176bdada62cabc6ec4b0308a930e83b679d5d36John Reck does_it_fit_fixed_48_16 (result_f.v[1]) && 1531176bdada62cabc6ec4b0308a930e83b679d5d36John Reck does_it_fit_fixed_48_16 (result_f.v[2]))) 1541176bdada62cabc6ec4b0308a930e83b679d5d36John Reck { 1551176bdada62cabc6ec4b0308a930e83b679d5d36John Reck for (j = 0; j < 3; j++) 1561176bdada62cabc6ec4b0308a930e83b679d5d36John Reck { 1571176bdada62cabc6ec4b0308a930e83b679d5d36John Reck double diff = fabs (result_f.v[j] - 1581176bdada62cabc6ec4b0308a930e83b679d5d36John Reck pixman_fixed_to_float128 (result_i.v[j])); 1591176bdada62cabc6ec4b0308a930e83b679d5d36John Reck 1601176bdada62cabc6ec4b0308a930e83b679d5d36John Reck if (is_affine && diff > (0.51 / 65536.0)) 1611176bdada62cabc6ec4b0308a930e83b679d5d36John Reck { 1621176bdada62cabc6ec4b0308a930e83b679d5d36John Reck printf ("%d:%d: bad precision for affine (%.12f)\n", 1631176bdada62cabc6ec4b0308a930e83b679d5d36John Reck testnum, i, diff); 1641176bdada62cabc6ec4b0308a930e83b679d5d36John Reck abort (); 1651176bdada62cabc6ec4b0308a930e83b679d5d36John Reck } 1661176bdada62cabc6ec4b0308a930e83b679d5d36John Reck else if (diff > (0.71 / 65536.0)) 1671176bdada62cabc6ec4b0308a930e83b679d5d36John Reck { 1681176bdada62cabc6ec4b0308a930e83b679d5d36John Reck printf ("%d:%d: bad precision for projective (%.12f)\n", 1691176bdada62cabc6ec4b0308a930e83b679d5d36John Reck testnum, i, diff); 1701176bdada62cabc6ec4b0308a930e83b679d5d36John Reck abort (); 1711176bdada62cabc6ec4b0308a930e83b679d5d36John Reck } 1721176bdada62cabc6ec4b0308a930e83b679d5d36John Reck } 1731176bdada62cabc6ec4b0308a930e83b679d5d36John Reck } 1741176bdada62cabc6ec4b0308a930e83b679d5d36John Reck } 1751176bdada62cabc6ec4b0308a930e83b679d5d36John Reck#endif 1761176bdada62cabc6ec4b0308a930e83b679d5d36John Reck } 1771176bdada62cabc6ec4b0308a930e83b679d5d36John Reck return crc32; 1781176bdada62cabc6ec4b0308a930e83b679d5d36John Reck} 1791176bdada62cabc6ec4b0308a930e83b679d5d36John Reck 1801176bdada62cabc6ec4b0308a930e83b679d5d36John Reckint 1811176bdada62cabc6ec4b0308a930e83b679d5d36John Reckmain (int argc, const char *argv[]) 1821176bdada62cabc6ec4b0308a930e83b679d5d36John Reck{ 1831176bdada62cabc6ec4b0308a930e83b679d5d36John Reck return fuzzer_test_main ("matrix", 20000, 1841176bdada62cabc6ec4b0308a930e83b679d5d36John Reck 0xBEBF98C3, 1851176bdada62cabc6ec4b0308a930e83b679d5d36John Reck test_matrix, argc, argv); 1861176bdada62cabc6ec4b0308a930e83b679d5d36John Reck} 187