1 /*
2  * Copyright © 2015 Intel Corporation
3  *
4  * Permission is hereby granted, free of charge, to any person obtaining a
5  * copy of this software and associated documentation files (the "Software"),
6  * to deal in the Software without restriction, including without limitation
7  * the rights to use, copy, modify, merge, publish, distribute, sublicense,
8  * and/or sell copies of the Software, and to permit persons to whom the
9  * Software is furnished to do so, subject to the following conditions:
10  *
11  * The above copyright notice and this permission notice (including the next
12  * paragraph) shall be included in all copies or substantial portions of the
13  * Software.
14  *
15  * THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
16  * IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
17  * FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT.  IN NO EVENT SHALL
18  * THE AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
19  * LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING
20  * FROM, OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS
21  * IN THE SOFTWARE.
22  */
23 
24 #include <stdio.h>
25 #include <stdbool.h>
26 #include <string.h>
27 #include <math.h>
28 
29 #include "macros.h"
30 #include "rounding.h"
31 
32 int main(int argc, char *argv[])
33 {
34    const struct {
35       float input, expected;
36    } float_data[] = {
37       { 0.0,                  0.0 },
38       { nextafterf(0.5, 0.0), 0.0 },
39       { 0.5,                  0.0 },
40       { nextafterf(0.5, 1.0), 1.0 },
41       { 1.0,                  1.0 },
42       { nextafterf(1.5, 1.0), 1.0 },
43       { 1.5,                  2.0 },
44       { nextafterf(1.5, 2.0), 2.0 },
45       { 2.0,                  2.0 },
46       { nextafterf(2.5, 2.0), 2.0 },
47       { 2.5,                  2.0 },
48       { nextafterf(2.5, 3.0), 3.0 },
49    };
50 
51    const struct {
52       double input, expected;
53    } double_data[] = {
54       { 0.0,                 0.0 },
55       { nextafter(0.5, 0.0), 0.0 },
56       { 0.5,                 0.0 },
57       { nextafter(0.5, 1.0), 1.0 },
58       { 1.0,                 1.0 },
59       { nextafter(1.5, 1.0), 1.0 },
60       { 1.5,                 2.0 },
61       { nextafter(1.5, 2.0), 2.0 },
62       { 2.0,                 2.0 },
63       { nextafter(2.5, 2.0), 2.0 },
64       { 2.5,                 2.0 },
65       { nextafter(2.5, 3.0), 3.0 },
66    };
67 
68    bool failed = false;
69    int i;
70 
71    for (i = 0; i < ARRAY_SIZE(float_data); i++) {
72       float output = _mesa_roundevenf(float_data[i].input);
73       if (memcmp(&float_data[i].expected, &output, sizeof(float))) {
74          fprintf(stderr, "%d float: expected %f (%a) from "
75                          "_mesa_roundevenf(%f (%a)) but got %f (%a)\n",
76                  i,
77                  float_data[i].expected,
78                  float_data[i].expected,
79                  float_data[i].input,
80                  float_data[i].input,
81                  output,
82                  output);
83          failed = true;
84       }
85    }
86 
87    /* Test negated values */
88    for (i = 0; i < ARRAY_SIZE(float_data); i++) {
89       float output = _mesa_roundevenf(-float_data[i].input);
90       float negated_expected = -float_data[i].expected;
91       if (memcmp(&negated_expected, &output, sizeof(float))) {
92          fprintf(stderr, "%d float: expected %f (%a) from "
93                          "_mesa_roundevenf(%f (%a)) but got %f (%a)\n",
94                  i,
95                  negated_expected,
96                  negated_expected,
97                  -float_data[i].input,
98                  -float_data[i].input,
99                  output,
100                  output);
101          failed = true;
102       }
103    }
104 
105    for (i = 0; i < ARRAY_SIZE(double_data); i++) {
106       double output = _mesa_roundeven(double_data[i].input);
107       if (memcmp(&double_data[i].expected, &output, sizeof(double))) {
108          fprintf(stderr, "%d double: expected %f (%a) from "
109                          "_mesa_roundeven(%f (%a)) but got %f (%a)\n",
110                  i,
111                  double_data[i].expected,
112                  double_data[i].expected,
113                  double_data[i].input,
114                  double_data[i].input,
115                  output,
116                  output);
117          failed = true;
118       }
119    }
120 
121    /* Test negated values */
122    for (i = 0; i < ARRAY_SIZE(double_data); i++) {
123       double output = _mesa_roundeven(-double_data[i].input);
124       double negated_expected = -double_data[i].expected;
125       if (memcmp(&negated_expected, &output, sizeof(double))) {
126          fprintf(stderr, "%d double: expected %f (%a) from "
127                          "_mesa_roundeven(%f (%a)) but got %f (%a)\n",
128                  i,
129                  negated_expected,
130                  negated_expected,
131                  -double_data[i].input,
132                  -double_data[i].input,
133                  output,
134                  output);
135          failed = true;
136       }
137    }
138 
139    return failed;
140 }
141