1*4882a593Smuzhiyun /*
2*4882a593Smuzhiyun * Copyright (C) 2007
3*4882a593Smuzhiyun * Wolfgang Denk, DENX Software Engineering, wd@denx.de.
4*4882a593Smuzhiyun *
5*4882a593Smuzhiyun * SPDX-License-Identifier: GPL-2.0+
6*4882a593Smuzhiyun */
7*4882a593Smuzhiyun /*
8*4882a593Smuzhiyun * This file is originally a part of the GCC testsuite.
9*4882a593Smuzhiyun * Check that certain subnormal numbers (formerly known as denormalized
10*4882a593Smuzhiyun * numbers) are rounded to within 0.5 ulp. PR other/14354.
11*4882a593Smuzhiyun */
12*4882a593Smuzhiyun
13*4882a593Smuzhiyun #include <common.h>
14*4882a593Smuzhiyun
15*4882a593Smuzhiyun #include <post.h>
16*4882a593Smuzhiyun
17*4882a593Smuzhiyun GNU_FPOST_ATTR
18*4882a593Smuzhiyun
19*4882a593Smuzhiyun #if CONFIG_POST & CONFIG_SYS_POST_FPU
20*4882a593Smuzhiyun
21*4882a593Smuzhiyun union uf
22*4882a593Smuzhiyun {
23*4882a593Smuzhiyun unsigned int u;
24*4882a593Smuzhiyun float f;
25*4882a593Smuzhiyun };
26*4882a593Smuzhiyun
27*4882a593Smuzhiyun static float
u2f(unsigned int v)28*4882a593Smuzhiyun u2f (unsigned int v)
29*4882a593Smuzhiyun {
30*4882a593Smuzhiyun union uf u;
31*4882a593Smuzhiyun u.u = v;
32*4882a593Smuzhiyun return u.f;
33*4882a593Smuzhiyun }
34*4882a593Smuzhiyun
35*4882a593Smuzhiyun static unsigned int
f2u(float v)36*4882a593Smuzhiyun f2u (float v)
37*4882a593Smuzhiyun {
38*4882a593Smuzhiyun union uf u;
39*4882a593Smuzhiyun u.f = v;
40*4882a593Smuzhiyun return u.u;
41*4882a593Smuzhiyun }
42*4882a593Smuzhiyun
43*4882a593Smuzhiyun static int ok = 1;
44*4882a593Smuzhiyun
45*4882a593Smuzhiyun static void
tstmul(unsigned int ux,unsigned int uy,unsigned int ur)46*4882a593Smuzhiyun tstmul (unsigned int ux, unsigned int uy, unsigned int ur)
47*4882a593Smuzhiyun {
48*4882a593Smuzhiyun float x = u2f (ux);
49*4882a593Smuzhiyun float y = u2f (uy);
50*4882a593Smuzhiyun
51*4882a593Smuzhiyun if (f2u (x * y) != ur)
52*4882a593Smuzhiyun /* Set a variable rather than aborting here, to simplify tracing when
53*4882a593Smuzhiyun several computations are wrong. */
54*4882a593Smuzhiyun ok = 0;
55*4882a593Smuzhiyun }
56*4882a593Smuzhiyun
57*4882a593Smuzhiyun /* We don't want to make this const and static, or else we risk inlining
58*4882a593Smuzhiyun causing the test to fold as constants at compile-time. */
59*4882a593Smuzhiyun struct
60*4882a593Smuzhiyun {
61*4882a593Smuzhiyun unsigned int p1, p2, res;
62*4882a593Smuzhiyun } static volatile expected[] =
63*4882a593Smuzhiyun {
64*4882a593Smuzhiyun {0xfff, 0x3f800400, 0xfff},
65*4882a593Smuzhiyun {0xf, 0x3fc88888, 0x17},
66*4882a593Smuzhiyun {0xf, 0x3f844444, 0xf}
67*4882a593Smuzhiyun };
68*4882a593Smuzhiyun
fpu_post_test_math7(void)69*4882a593Smuzhiyun int fpu_post_test_math7 (void)
70*4882a593Smuzhiyun {
71*4882a593Smuzhiyun unsigned int i;
72*4882a593Smuzhiyun
73*4882a593Smuzhiyun for (i = 0; i < ARRAY_SIZE(expected); i++)
74*4882a593Smuzhiyun {
75*4882a593Smuzhiyun tstmul (expected[i].p1, expected[i].p2, expected[i].res);
76*4882a593Smuzhiyun tstmul (expected[i].p2, expected[i].p1, expected[i].res);
77*4882a593Smuzhiyun }
78*4882a593Smuzhiyun
79*4882a593Smuzhiyun if (!ok) {
80*4882a593Smuzhiyun post_log ("Error in FPU math7 test\n");
81*4882a593Smuzhiyun return -1;
82*4882a593Smuzhiyun }
83*4882a593Smuzhiyun return 0;
84*4882a593Smuzhiyun }
85*4882a593Smuzhiyun
86*4882a593Smuzhiyun #endif /* CONFIG_POST & CONFIG_SYS_POST_FPU */
87