1 // Copyright (C) 2016-2017 Free Software Foundation, Inc.
3 // This file is part of the GNU ISO C++ Library. This library is free
4 // software; you can redistribute it and/or modify it under the
5 // terms of the GNU General Public License as published by the
6 // Free Software Foundation; either version 3, or (at your option)
9 // This library is distributed in the hope that it will be useful,
10 // but WITHOUT ANY WARRANTY; without even the implied warranty of
11 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
12 // GNU General Public License for more details.
14 // You should have received a copy of the GNU General Public License along
15 // with this library; see the file COPYING3. If not see
16 // <http://www.gnu.org/licenses/>.
18 // { dg-options "-std=gnu++17" }
19 // { dg-do run { target c++17 } }
20 // { dg-xfail-run-if "AIX long double" { powerpc-ibm-aix* } }
23 #include <type_traits>
24 #if defined(__TEST_DEBUG)
29 std::cout << "line " << __LINE__ \
30 << " max_abs_frac = " << max_abs_frac \
31 << " tolerance = " << toler \
35 #include <testsuite_hooks.h>
39 static_assert(is_same_v
<double, decltype(std::hypot(0.0, 0.0, 0.0))>);
40 static_assert(is_same_v
<double, decltype(std::hypot(0.0f
, 0.0, 0.0))>);
41 static_assert(is_same_v
<double, decltype(std::hypot(0.0, 0.0f
, 0.0))>);
42 static_assert(is_same_v
<double, decltype(std::hypot(0.0, 0.0, 0.0f
))>);
43 static_assert(is_same_v
<double, decltype(std::hypot(0.0f
, 0.0f
, 0.0))>);
44 static_assert(is_same_v
<double, decltype(std::hypot(0.0f
, 0.0, 0))>);
45 static_assert(is_same_v
<long double, decltype(std::hypot(0.0f
, 0.0, 0.0l))>);
46 static_assert(is_same_v
<long double, decltype(std::hypot(0, 0.0, 0.0l))>);
48 template<typename T
> struct testcase_hypot
{ T x
, y
, z
, f0
; };
50 template<typename Tp
, unsigned int Num
>
52 test(const testcase_hypot
<Tp
> (&data
)[Num
], Tp toler
)
54 const Tp eps
= std::numeric_limits
<Tp
>::epsilon();
55 Tp max_abs_diff
= -Tp(1);
56 Tp max_abs_frac
= -Tp(1);
57 unsigned int num_datum
= Num
;
58 for (unsigned int i
= 0; i
< num_datum
; ++i
)
60 const Tp f
= std::hypot(data
[i
].x
, data
[i
].y
, data
[i
].z
);
61 const Tp f0
= data
[i
].f0
;
62 const Tp diff
= f
- f0
;
63 if (std::abs(diff
) > max_abs_diff
)
64 max_abs_diff
= std::abs(diff
);
65 if (std::abs(f0
) > Tp(10) * eps
&& std::abs(f
) > Tp(10) * eps
)
67 const Tp frac
= diff
/ f0
;
68 if (std::abs(frac
) > max_abs_frac
)
69 max_abs_frac
= std::abs(frac
);
72 VERIFY(max_abs_frac
< toler
);
75 const testcase_hypot
<double> data1
[] = {
76 { 0.0, 0.0, 0.0, 0.0 },
77 { 0.0, 1.0, 1.0, std::sqrt(2.0) },
78 { 1.0, 1.0, 1.0, std::sqrt(3.0) },
79 { 1.0, 2.0, 2.0, 3.0 },
80 { 2.0, 3.0, 6.0, 7.0 },
81 { 1.0, 4.0, 8.0, 9.0 },
82 { 4.0, 4.0, 7.0, 9.0 },
83 { 12.0, 16.0, 21.0, 29.0 },
84 { 1e8
, 1., 1e-8, 1e8
},
85 { 1., 1e8
, 1e-8, 1e8
},
86 { 1e-8, 1., 1e8
, 1e8
},
87 { 1e-2, 1e-4, 1e-4, 0.01000099995 },
88 { 214748364., 214748364., 214748364., 371955077.2902952 }
90 const double toler1
= 1e-12;
92 const testcase_hypot
<float> data2
[] = {
93 { 0.0f
, 0.0f
, 0.0f
, 0.0f
},
94 { 0.0f
, 1.0f
, 1.0f
, std::sqrt(2.0f
) },
95 { 1.0f
, 1.0f
, 1.0f
, std::sqrt(3.0f
) },
96 { 1.0f
, 2.0f
, 2.0f
, 3.0f
},
97 { 2.0f
, 3.0f
, 6.0f
, 7.0f
},
98 { 1.0f
, 4.0f
, 8.0f
, 9.0f
},
99 { 4.0f
, 4.0f
, 7.0f
, 9.0f
},
100 { 12.0f
, 16.0f
, 21.0f
, 29.0f
},
101 { 1e8f
, 1.f
, 1e-8f
, 1e8f
},
102 { 1.f
, 1e8f
, 1e-8f
, 1e8f
},
103 { 1e-8f
, 1.f
, 1e8f
, 1e8f
},
104 { 1e-2f
, 1e-4f
, 1e-4f
, 0.010001f
},
105 { 214748364.f
, 214748364.f
, 214748364.f
, 371955072.f
}
107 const float toler2
= 1e-7f
;
109 const testcase_hypot
<long double> data3
[] = {
110 { 0.0l, 0.0l, 0.0l, 0.0l },
111 { 0.0l, 1.0l, 1.0l, std::sqrt(2.0l) },
112 { 1.0l, 1.0l, 1.0l, std::sqrt(3.0l) },
113 { 1.0l, 2.0l, 2.0l, 3.0l },
114 { 2.0l, 3.0l, 6.0l, 7.0l },
115 { 1.0l, 4.0l, 8.0l, 9.0l },
116 { 4.0l, 4.0l, 7.0l, 9.0l },
117 { 12.0l, 16.0l, 21.0l, 29.0l },
118 { 1e8l
, 1.l
, 1e-8l, 1e8l
},
119 { 1.l
, 1e8l
, 1e-8l, 1e8l
},
120 { 1e-8l, 1.l
, 1e8l
, 1e8l
},
121 { 1e-2l, 1e-4l, 1e-4l, 0.010000999950004999375l },
122 { 2147483647.l
, 2147483647.l
, 2147483647.l
, 3719550785.027307813987l }
124 const long double toler3
= 1e-16l;